Electronic and vibrational spectroscopy of the cytochrome c:cytochrome c oxidase complexes from bovine and Paracoccus denitrificans.

The 1:1 complex between horse heart cytochrome c and bovine cytochrome c oxidase, and between yeast cytochrome c and Paracoccus denitrificans cytochrome c oxidase have been studied by a combination of second derivative absorption, circular dichroism (CD), and resonance Raman spectroscopy. The second derivative absorption and CD spectra reveal changes in the electronic transitions of cytochrome a upon complex formation. These results could reflect changes in ground state heme structure or changes in the protein environment surrounding the chromophore that affect either the ground or excited electronic states. The resonance Raman spectrum, on the other hand, reflects the heme structure in the ground electronic state only and shows no significant difference between cytochrome a vibrations in the complex or free enzyme. The only major difference between the Raman spectra of the free enzyme and complex is a broadening of the cytochrome a3 formyl band of the complex that is relieved upon complex dissociation at high ionic strength. These data suggest that the differences observed in the second derivative and CD spectra are the result of changes in the protein environment around cytochrome a that affect the electronic excited state. By analogy to other protein-chromophore systems, we suggest that the energy of the Soret pi* state of cytochrome a may be affected by (1) changes in the local dielectric, possibly brought about by movement of a charged amino acid side chain in proximity to the heme group, or (2) pi-pi interactions between the heme and aromatic amino acid residues.     

Sequence context strongly modulates association of polar residues in transmembrane helices

Polar residues are capable of mediating the association of membrane-embedded helices through the formation of side-chain/side-chain inter-helical hydrogen bonds. However, the extent to which native van der Waals packing of the residues surrounding the polar locus can enhance, or interfere with, the interaction of polar residues has not yet been studied. We examined the propensities of four polar residues (aspartic acid, asparagine, glutamic acid, and glutamine) to promote self-association of transmembrane (TM) domains in several biologically derived sequence environments, including (i). four naturally occurring TM domains that contain a Glu or Gln residue (Tnf5/CD40 ligand, C79a/Ig-alpha, C79b/Ig-beta, and Fut3/alpha-fucosyltransferase); and (ii). variants of bacteriophage M13 major coat protein TM segment with Asp and Asn at interfacial and non-interfacial positions. Self-association was quantified by the TOXCAT assay, which measures TM helix self-oligomerization in the Escherichia coli inner membrane. While an appropriately placed polar residue was found in several cases to significantly stabilize TM helix-helix interactions through the formation of an interhelical hydrogen bond, in other cases the strongly polar residues did not enhance the association of the two helices. Overall, these results suggest that an innate structural mechanism may operate to control non-specific association of membrane-embedded polar residues.     

The role of weakly polar and H-bonding interactions in the stabilization of the conformers of FGG, WGG, and YGG: an aqueous phase computational study

The energetics of intramolecular interactions on the conformational potential energy surface of the terminally protected N-Ac-Phe-Gly-Gly-NHMe (FGG), N-Ac-Trp-Gly-Gly-NHMe (WGG), and N-Ac-Tyr-Gly-Gly-NHMe (YGG) tripeptides was investigated. To identify the representative conformations, simulated annealing molecular dynamics (MD) and density functional theory (DFT) methods were used. The interaction energies were calculated at the BHandHLYP/aug-cc-pVTZ level of theory. In the global minima, 10%, 31%, and 10% of the stabilization energy come from weakly polar interactions, respectively, in FGG, WGG, and YGG. In the prominent cases 46%, 62%, and 46% of the stabilization energy is from the weakly polar interactions, respectively, in FGG, WGG, and YGG. On average, weakly polar interactions account for 15%, 34%, and 9% of the stabilization energies of the FGG, WGG, and YGG conformers, respectively. Thus, weakly polar interactions can make an important energetic contribution to protein structure and function. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 1002-1011, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

Environmental influences on epistatic interactions: viabilities of cytochrome c genotypes in interpopulation crosses

Abstract .The genetic incompatibilities that underlie F₂ hybrid breakdown and reproductive isolation between al-lopatric populations may be susceptible to environmental interactions. Here we show that epistatic interactions between cytochrome c ( CYC ) alleles and mitochondrial DNA (mtDNA) variation are dramatically influenced by environmental temperature in interpopulation hybrids of the copepod Tigriopus californicus . CYC is a nuclear-encoded gene that functionally interacts with electron transport system (ETS) complexes composed in part of mtDNA-encoded proteins. Previous studies have provided evidence for functional coadaptation between CYC and ETS complex IV (cytochrome c oxidase) and for cytoplasmic effects on the fitness of CYC genotype in copepod hybrids. In this study, selection on CYC genotype is shown to continue into advanced generation hybrids (F₂-F₈) increasing the likelihood that CYC itself is involved in the interaction (and not a linked factor). Relative viabilities varied markedly between copepods raised in two different temperature/light regimes. These results suggest that both intrinsic coadaptation and extrinsic selection will influence the outcome of natural hybridizations between populations. Furthermore, the results indicate that the fitness of particular hybrid genotypes depends on additional non-mtDNA encoded genes that interact with CYC.     

Stopped-flow NMR measurement of hydrogen exchange rates in reduced horse cytochrome c under strongly destabilizing conditions

A procedure to measure exchange rates of fast exchanging protein amide hydrogens by time-resolved NMR spectroscopy following in situ initiation of the reaction by diluting a native protein solution into an exchanging deuterated buffer is described. The method has been used to measure exchange rates of a small set of amide hydrogens of reduced cytochrome c, maintained in a strictly anaerobic atmosphere, in the presence of an otherwise inaccessible range of guanidinium deuterochloride concentrations. The results for the measured protons indicate that hydrogen exchange in the unfolding transition region of cytochrome c reach the EX2 limit, but emphasize the difficulty in interpretation of the exchange mechanism in protein hydrogen exchange studies. Comparison of free energies of structure opening for the measured hydrogens with the global unfolding free energy monitored by far-UV CD measurements has indicated the presence of at least one partially unfolded equilibrium species of reduced cytochrome c. The results provide the first report of measurement of free energy of opening of structure to exchange in the 0–2-kcal/mol range. Proteins 32:241–247, 1998. © 1998 Wiley-Liss, Inc.     

Replacements in a conserved leucine cluster in the hydrophobic heme pocket of cytochrome c.

A cluster of highly conserved leucine side chains from residues 9, 68, 85, 94, and 98 is located in the hydrophobic heme pocket of cytochrome c. The contributions of two of these, Leu 85 and Leu 94, have been studied using a protein structure-function-mutagenesis approach to probe their roles in the maintenance of overall structural integrity and electron transfer activity. Structural studies of the L85C, L85F, L85M, and L94S mutant proteins show that, in each case, the overall fold of cytochrome c is retained, but that localized conformational shifts are required to accommodate the introduced side chains. In particular, the side chains of Cys 85 and Phe 85 form energetically favorable interactions with Phe 82, whereas Met 85 takes on a more remote conformation to prevent an unfavorable interaction with the phenyl ring of Phe 82. In the case of the L94S mutant protein, the new polar group introduced is found to form hydrogen bonds to nearby carbonyl groups. In all proteins with substitutions at Leu 85, the hydrophobic nature of the heme pocket is preserved and no significant decrease in heme reduction potential is observed. Despite earlier predictions that Leu 85 is an important determinant in cytochrome c electron transfer partner complexation, our studies show this is unlikely to be the case because the considerable surface contour perturbations made by substitutions at this residue do not correspondingly translate into significant changes in electron transfer rates. For the L94S mutant protein, the substitution of a polar hydroxyl group directly into the hydrophobic heme pocket has a larger effect on heme reduction potential, but this is mitigated by two factors. First, the side chain of Ser 94 is rotated away from the heme group and, second, the side chain of Leu 98 shifts into a portion of the new space available, partially shielding the heme group. The Leu 94 Ser substitution does not perturb the highly conserved interface formed by the nearly perpendicular packing of the N- and C-terminal helices of cytochrome c, ruling this out as the cause of this mutant protein becoming thermally labile and having a lower functional activity. Our results show these effects are most likely attributable to disruption of the heme pocket region. Much of the ability of cytochrome c to absorb the introduction of mutations at Leu 85 and Leu 94 appears to be a consequence of the conformational flexibility afforded by the leucine cluster in this region as well as the presence of a nearby internal cavity.(ABSTRACT TRUNCATED AT 400 WORDS)     

Fast folding of cytochrome c.

Native iso-2 cytochrome c contains two residues (His 18, Met 80) coordinated to the covalently attached heme. On unfolding of iso-2, the His 18 ligand remains coordinated to the heme iron, whereas Met 80 is displaced by a non-native heme ligand, His 33 or His 39. To test whether non-native His-heme ligation slows folding, we have constructed a double mutant protein in which the non-native ligands are replaced by asparagine and lysine, respectively (H33N,H39K iso-2). The double mutant protein, which cannot form non-native histidine-heme coordinate bonds, folds significantly faster than normal iso-2 cytochrome c: gamma = 14-26 ms for H33N,H39K iso-2 versus gamma = 200-1,100 ms for iso-2. These results with iso-2 cytochrome c strongly support the hypothesis that non-native His-heme ligation results in a kinetic barrier to fast folding of cytochrome c. Assuming that the maximum rate of a conformational search is about 10(11) s-1, the results imply that the direct folding pathway of iso-2 involves passage through on the order of 10(9) or fewer partially folded conformers.     

Axial ligand replacement in horse heart cytochrome c by semisynthesis

Semisynthesis has been employed to replace the axial methionine in horse heart cytochrome c with histidine. The reduction potential of the His-80 protein (cyt c-His-80) is 41 mV vs NHE (0.1 M phosphate; pH 7.0; 25 degrees C). The absorption spectra of oxidized and reduced cyt c-His-80 are very similar to those of the native protein in the porphyrin region, but the 695 nm band is absent in the oxidized His-80 protein.     

Determinants of protein hydrogen exchange studied in equine cytochrome c.

The exchange of a large number of amide hydrogens in oxidized equine cytochrome c was measured by NMR and compared with structural parameters. Hydrogens known to exchange through local structural fluctuations and through larger unfolding reactions were separately considered. All hydrogens protected from exchange by factors greater than 10(3) are in defined H-bonds, and almost all H-bonded hydrogens including those at the protein surface were measured to exchange slowly. H-exchange rates do not correlate with H-bond strength (length) or crystallographic B factors. It appears that the transient structural fluctuation necessary to bring an exchangeable hydrogen into H-bonding contact with the H-exchange catalyst (OH(-)-ion) involves a fairly large separation of the H-bond donor and acceptor, several angstroms at least, and therefore depends on the relative resistance to distortion of immediately neighboring structure. Accordingly, H-exchange by way of local fluctuational pathways tends to be very slow for hydrogens that are neighbored by tightly anchored structure and for hydrogens that are well buried. The slowing of buried hydrogens may also reflect the need for additional motions that allow solvent access once the protecting H-bond is separated, although it is noteworthy that burial in a protein like cytochrome c does not exceed 4 angstroms. When local fluctuational pathways are very slow, exchange can become dominated by a different category of larger, cooperative, segmental unfolding reactions reaching up to global unfolding.     

Oxygen and proton pathways in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump, which couples the reduction of oxygen to water to the translocation of protons across the membrane. The recently solved x-ray structures of cytochrome c oxidase permit molecular dynamics simulations of the underlying transport processes. To eventually establish the proton pump mechanism, we investigate the transport of the substrates, oxygen and protons, through the enzyme. Molecular dynamics simulations of oxygen diffusion through the protein reveal a well-defined pathway to the oxygen-binding site starting at a hydrophobic cavity near the membrane-exposed surface of subunit I, close to the interface to subunit III. A large number of water sites are predicted within the protein, which could play an essential role for the transfer of protons in cytochrome c oxidase. The water molecules form two channels along which protons can enter from the cytoplasmic (matrix) side of the protein and reach the binuclear center. A possible pumping mechanism is proposed that involves a shuttling motion of a glutamic acid side chain, which could then transfer a proton to a propionate group of heme α₃. Proteins 30:100–107, 1998. © 1998 Wiley-Liss, Inc.     

Solution structure of reduced horse heart cytochrome c

 In the frame of a broad study on the structural differences between the two redox forms of cytochromes to be related to the electron transfer process, the NMR solution structure of horse heart cytochrome c in the reduced form has been determined. The structural data obtained in the present work are compared to those already available in the literature on the same protein and the presence of conformational differences is discussed in the light of the experimental method employed for the structure determination. Redox-state dependent changes are analyzed and in particular they are related to the role of propionate-7 of the heme. Also some hydrogen bonds are changed upon reduction of the heme iron. A substantial similarity is observed for the backbone fold, independently of the oxidation state. At variance, some meaningful differences are observed in the orientation of a few side chains. These changes are related to those found in the case of the highly homologous cytochrome c from Saccharomyces cerevisiae. The exchangeability of the NH protons has been investigated and found to be smaller than in the case of the oxidized protein. We think that this is a characteristic of reduced cytochromes and that mobility is a medium for molecular recognition in vivo. Received: 8 June 1998 / Accepted: 13 October 1998     

Probing electrostatic interactions in cytochrome c using site-directed chemical modification.

This communication reports the generation of an electrostatic probe using chemical modification of methionine side chains. The alkylation of methionine by iodoacetamide was achieved in a set of Saccharomyces cerevisiae iso-1-cytochrome c mutants, introducing the nontitratable, nondelocalized positive charge of a carboxyamidomethylmethionine sulfonium (CAMMS) ion at five surface and one buried site in the protein. Changes in redox potential and its variation with temperature were used to calculate microscopic effective dielectric constants operating between the probe and the heme iron. Dielectric constants (epsilon) derived from deltadeltaG values were not useful due to entropic effects, but epsilon(deltadeltaH) gave results that supported the theory. The effect on biological activity of surface derivatization was interpreted in terms of protein-protein interactions. The introduction of an electrostatic probe in cytochrome c often resulted in marked effects on activity with one of two physiological partners: cytochrome c reductase, especially if introduced at position 65, and cytochrome c oxidase, if at position 28.     

The molecular structure of an unusual cytochrome c2 determined at 2.0 A; the cytochrome cH from Methylobacterium extorquens.

Cytochrome cH is the electron donor to the oxidase in methylotrophic bacteria. Its amino acid sequence suggests that it is a typical Class 1 cytochrome c, but some features of the sequence indicated that its structure might be of special interest. The structure of oxidized cytochrome cH has been solved to 2.0 A resolution by X-ray diffraction. It has the classical tertiary structure of the Class 1 cytochromes c but bears a closer gross resemblance to mitochondrial cytochrome c than to the bacterial cytochrome c2. The left-hand side of the haem cleft is unique; in particular, it is highly hydrophobic, the usual water is absent, and the "conserved" Tyr67 is replaced by tryptophan. A number of features of the structure demonstrate that the usual hydrogen bonding network involving water in the haem channel is not essential and that other mechanisms may exist for modulation of redox potentials in this cytochrome.     

Oxidative lipidomics of apoptosis: redox catalytic interactions of cytochrome c with cardiolipin and phosphatidylserine

The primary life-supporting function of cytochrome c (cyt c) is control of cellular energetic metabolism as a mobile shuttle in the electron transport chain of mitochondria. Recently, cyt c's equally important life-terminating function as a trigger and regulator of apoptosis was identified. This dreadful role is realized through the relocalization of mitochondrial cyt c to the cytoplasm where it interacts with Apaf-1 in forming apoptosomes and mediating caspase-9 activation. Although the presence of heme moiety of cyt c is essential for the latter function, cyt c's redox catalytic features are not required. Lately, two other essential functions of cyt c in apoptosis, that may rely heavily on its redox activity have been suggested. Both functions are directed toward oxidation of two negatively charged phospholipids, cardiolipin (CL) in the mitochondria and phosphatidylserine (PS) in the plasma membrane. In both cases, oxidized phospholipids seem to be essential for the transduction of two distinctive apoptotic signals: one is participation of oxidized CL in the formation of the mitochondrial permeability transition pore that facilitates release of cyt c into the cytosol and the other is the contribution of oxidized PS to the externalization and recognition of PS (and possibly oxidized PS) on the cell surface by specialized receptors of phagocytes. In this review, we present a new concept that cyt c actuates both of these oxidative roles through a uniform mechanism: its specific interactions with each of these phospholipids result in the conversion and activation of cyt c, transforming it from an innocuous electron transporter into a calamitous peroxidase capable of oxidizing the activating phospholipids. We also show that this new concept is compatible with a leading role for reactive oxygen species in the execution of the apoptotic program, with cyt c as the main executioner.     

Modulation of Bacillus pasteurii cytochrome c 553 reduction potential by structural and solution parameters

 Direct cyclic voltammetry and ¹H NMR spectroscopy have been combined to investigate the electrochemical and spectroscopic properties of cytochrome c ₅₅₃ isolated from the alkaliphilic soil bacterium Bacillus pasteurii. A quasi-reversible diffusion-controlled redox process is exhibited by cytochrome c ₅₅₃ at a pyrolitic graphite edge microelectrode. The temperature dependence of the reduction potential, measured using a non-isothermal electrochemical cell, revealed a discontinuity at 308 K. The thermodynamic parameters determined in the low-temperature range (275–308 K;ΔS°′=–162.7±1.2 J mol^–1 K^–1, ΔH°′=–53.0±0.5 kJ mol^–1, ΔG°′=–4.5±0.1 kJ mol^–1, E°′=+47.0±0.6 mV) indicate the presence of large enthalpic and entropic effects, leading, respectively, to stabilization and destabilization of the reduced form of cytochrome c ₅₅₃. Both effects are more accentuated in the high-temperature range (308–323 K;ΔS°′=–294.1±8.4 J mol^–1 K^–1, ΔH°′=–93.4±3.1 kJ mol^–1, ΔG°′=–5.8±0.6 kJ mol^–1, E°′=+60.3±5.8 mV), with the net result being a slight increase of the standard reduction potential. These thermodynamic parameters are interpreted using the compensation theory of hydration of biopolymers as indicating the extrusion, upon reduction, of water molecules from the hydration sphere of the cytochrome. The low-T and high-T conformers differ by the number of water molecules in the solvation sphere: in the high-T conformer, the number of water molecules extruded upon reduction increases, as compared to the low-T conformer. The ionic strength dependence of the reduction potential at 298 K, treated within the frame of extended Debye-Hückel theory, yields values of E ^°′ _(I=0) =–25.4±1.4 mV, z _red=–11.3, and z _ox=–10.3. The pH dependence of the reduction potential at 298 K shows a plateau in the pH range 7–10 and an increase at more acidic pH, allowing the calculation of pK _O=5.5 and pK _R=5.7, together with the estimate of the reduction potentials of completely protonated (+71 mV) and deprotonated (+58 mV) forms of cytochrome c ₅₅₃. ¹H NMR spectra of the oxidized paramagnetic cytochrome c ₅₅₃ indicate the presence of a His-Met axial coordination of the low-spin (S=1/2) heme iron, which is maintained in the temperature interval 288–340 K at pH 7 and in the pH range 4.8–10.0 at 298 K. The temperature dependence of the hyperfine-shifted signals shows both Curie-type and anti-Curie-type behavior, with marked deviations from linearity, interpreted as indicating the presence of a fast equilibrium between the low-T and high-T conformers, having slightly different heme electronic structures resulting from the T-induced conformational change. Increasing the NaCl concentration in the range 0–0.2 M causes a slight change of the ¹H NMR chemical shifts of the hyperfine-shifted signals, with no influence on their linewidth. The calculated lower limit value of the apparent affinity constant for specific ion binding is estimated as 5.2±1.1 M^–1. The pH dependence of the isotropically shifted ¹H NMR signals of the oxidized cytochrome displays at least one ionization step with pK _O=5.7. The thermodynamic and spectroscopic data indicate a large solvent-derived entropic effect as the main cause for the observed low reduction potential of B. pasteurii cytochrome c ₅₅₃. Received: 9 January 1998 / Accepted: 8 April 1998     

Secondary and tertiary structure of the A-state of cytochrome c from resonance Raman spectroscopy.

Ferricytochrome c can be converted to the partially folded A-state at pH 2.2 in the presence of 1.5 M NaCl. The structure of the A-state has been studied in comparison with the native and unfolded states, using resonance Raman spectroscopy with visible and ultraviolet excitation wavelengths. Spectra obtained with 200 nm excitation show a decrease in amide II intensity consistent with loss of structure for the 50s and 70s helices. The 230-nm spectra contain information on vibrational modes of the single Trp 59 side chain and the four tyrosine side chains (Tyr 48, 67, 74, and 97). The Trp 59 modes indicate that the side chain remains in a hydrophobic environment but loses its tertiary hydrogen bond and is rotationally disordered. The tyrosine modes Y8b and Y9a show disruption of tertiary hydrogen bonding for the Tyr 48, 67, and 74 side chains. The high-wavenumber region of the 406.7-nm resonance Raman spectrum reveals a mixed spin heme iron atom, which arises from axial coordination to His 18 and a water molecule. The low-frequency spectral region reports on heme distortions and indicates a reduced degree of interaction between the heme and the polypeptide chain. A structural model for the A-state is proposed in which a folded protein subdomain, consisting of the heme and the N-terminal, C-terminal, and 60s helices, is stabilized through nonbonding interactions between helices and with the heme.     

Equilibrium unfolding of a small low-potential cytochrome, cytochrome c553 from Desulfovibrio vulgaris.

To understand general aspects of stability and folding of c-type cytochromes, we have studied the folding characteristics of cytochrome c553 from Desulfovibrio vulgaris (Hildenborough). This cytochrome is structurally similar but lacks sequence homology to other heme proteins; moreover, it has an abnormally low reduction potential. Unfolding of oxidized and reduced cytochrome c553 by guanidine hydrochloride (GuHCl) was monitored by circular dichroism (CD) and Soret absorption; the same unfolding curves were obtained with both methods supporting that cytochrome c553 unfolds by an apparent two-state process. Reduced cytochrome c553 is 7(3) kJ/mol more stable than the oxidized form; accordingly, the reduction potential of unfolded cytochrome c553 is 100(20) mV more negative than that of the folded protein. In contrast to many other unfolded cytochrome c proteins, upon unfolding at pH 7.0 both oxidized and reduced heme in cytochrome c553 become high-spin. The lack of heme misligation in unfolded cytochrome c553 implies that its unfolded structure is less constrained than those of cytochromes c with low-spin, misligated hemes.     

Constitutive presence of cytochrome c in the cytosol of a chemoresistant leukemic cell line

The release of holocytochrome c (cyt c) from mitochondria into the cytosol is reportedly a landmark of the execution phase of apoptosis. As shown here, the P-glycoprotein- (P-gp) expressing K562/ADR cell line (but not the parental K562 cell line) exhibits both cytosolic and mitochondrial cyt c in the absence of any signs of apoptosis. K562/ADR cells were found to be relatively resistant to a variety of different inducers of apoptosis, and blocking the P-gp did not reverse this resistance. The release of cyt c in non-apoptotic K562/ADR cells was not accompanied by that of any other mitochondrial apoptogenic protein, such as AIF or Smac/DIABLO, and was inhibited by Bcl-2 over expression. In addition, using a cell-free system, we show that mitochondria isolated from K562/ADR cells spontaneously released cyt c. These data suggest that cyt c release may be compatible with the preservation of mitochondrial integrity and function, as well as cell proliferation.     

Protein--protein docking of electron transfer complexes: cytochrome c oxidase and cytochrome c

Electron transferring protein complexes form only transiently and the crystal structures of electron transfer protein--protein complexes involving cytochrome c could so far be determined only for the pairs of yeast cytochrome c peroxidase (CcP) with iso-1-cytochrome c (iso-1-cyt c) and with horse heart cytochrome c (cyt c). This article presents models from computational docking for complexes of cytochrome c oxidase (COX) from Paracoccus denitrificans with horse heart cytochrome c, and with its physiological counterpart cytochrome c552 (c552). Initial docking is performed with the FTDOCK program, which permits an exhaustive search of translational and rotational space. A filtering procedure is then applied to reduce the number of complexes to a manageable number. In a final step of structural and energetic refinement, the complexes are optimized by rigid-body energy minimization with the molecular mechanics package CHARMM. This methodology was first tested on the CcP:iso-1-cyt c complex, in which the complex with the lowest CHARMM energy has an RMSD from the crystal structure of only 1.8 A (C(alpha) carbon atoms). Notably, the crystal conformation has an even lower energy. The same procedure was then applied to COX:cyt c and COX:c552. The lowest-energy COX:cyt c complex is very similar to a docking model previously described for the complex of bovine cytochrome c oxidase with horse heart cytochrome c. For the COX:c552 complex, cytochrome c552 is found in two different orientations, depending on whether it is docked against COX from a two-subunit or from a four-subunit crystal structure, respectively. Both conformations are discussed critically in the light of the available experimental data.