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.     

Structural analysis of zinc-substituted cytochrome c

Zinc-substituted cytochrome c has been widely used in studies of protein-protein interactions and photo-induced electron transfer reactions between proteins. However, the coordination geometry of zinc in zinc-substituted cyt c has not yet been determined; two different opinions about the coordination have been reached. Here the solution structures of zinc-substituted cytochrome c that might be five-coordinated and six-coordinated have been refined separately by using (1)H NMR spectroscopy, and the zinc coordination geometry was determined just by NOE distance constraints. Structural analysis of the energy-minimized average solution structures of both the pentacoordinated and hexacoordinated geometries indicate that that zinc in zinc-substituted cyt c should be bound to both His18 and Met80, which means that the zinc is six-coordinated. RMSD values of the family of 25 six-coordinated structures from the average structure are 0.66+/-0.13 A and 1.09+/-0.16 A for the backbone and all heavy atoms, respectively. A statistical analysis of the structure indicates its satisfactory quality. Comparison of the solution structure of the six-coordinated energy-minimized average structure of zinc-substituted cytochrome c with the solution structure of reduced cytochrome c reveals that for the overall folding the secondary structure elements are very close. The availability of the structure provides for a better understanding of the protein-protein complex and for electron transfer processes between Zn cyt c and other metalloproteins.     

MAD structure of Pseudomonas nautica dimeric cytochrome c552 mimicks the c4 Dihemic cytochrome domain association.

The monohemic cytochrome c552from Pseudomonas nautica (c552-Pn) is thought to be the electron donor to cytochrome cd1, the so-called nitrite reductase (NiR). It shows as high levels of activity and affinity for the P. nautica NiR (NiR-Pn), as the Pseudomonas aeruginosa enzyme (NiR-Pa). Since cytochrome c552is by far the most abundant electron carrier in the periplasm, it is probably involved in numerous other reactions. Its sequence is related to that of the c type cytochromes, but resembles that of the dihemic c4cytochromes even more closely.The three-dimensional structure of P. nautica cytochrome c552has been solved to 2.2 A resolution using the multiple wavelength anomalous dispersion (MAD) technique, taking advantage of the presence of the eight Fe heme ions in the asymmetric unit. Density modification procedures involving 4-fold non-crystallographic averaging yielded a model with an R -factor value of 17.8 % (Rfree=20.8 %). Cytochrome c552forms a tight dimer in the crystal, and the dimer interface area amounts to 19% of the total cytochrome surface area. Four tighly packed dimers form the eight molecules of the asymmetric unit.The c552dimer is superimposable on each domain of the monomeric cytochrome c4from Pseudomomas stutzeri (c4-Ps), a dihemic cytochrome, and on the dihemic c domain of flavocytochrome c of Chromatium vinosum (Fcd-Cv). The interacting residues which form the dimer are both similar in character and position, which is also true for the propionates. The dimer observed in the crystal also exists in solution. It has been hypothesised that the dihemic c4-Ps may have evolved via monohemic cytochrome c gene duplication followed by evolutionary divergence and the adjunction of a connecting linker. In this process, our dimeric c552structure might be said to constitute a "living fossile" occurring in the course of evolution between the formation of the dimer and the gene duplication and fusion. The availability of the structure of the cytochrome c552-Pn and that of NiR from P. aeruginosa made it possible to identify putative surface patches at which the docking of c552to NiR-Pn may occur.     

The structure of Paracoccus denitrificans cytochrome c550.

The crystal structure of Paracoccus (formerly Micrococcus) denitrificans cytochrome c550 has been solved by x-ray diffraction to a resolution of 2.45 A. In both amino acid sequence and molecular structure it is evolutionarily homologous with mitochondrial cytochrome c from eukaryotes and photosynthetic cytochrome c2 from purple non-sulfur bacteria. All of these cytochromes c have the same basic folding pattern, with surface insertions of extra amino acids in c550. Various strains of c2 have all, some, or none of the extra insertions observed in c550. The hydrophobic heme environment, position of aromatic rings, and structure and environment of the heme crevice, are virtually identical in cytochromes c55o, c, and c2. Radical changes observed at all regions on the molecular surface except the heme crevice argue for the importance of the crevice and the exposed edge of the heme in the transfer of electrons to and from the cytochrome molecule.     

Solution structure of mitochondrial cytochrome c. I. 1H nuclear magnetic resonance of ferricytochrome c

The 1H nuclear magnetic resonance spectra of tuna and horse ferricytochromes c have been investigated and the resonances of all amino acid methyl groups have been assigned to specific absorption lines. The assignment procedure involves principally the comparison of one-dimensional nuclear magnetic resonance spectra from a range of homologous ferricytochromes c and does not require a prior knowledge of the secondary or tertiary protein structure. Of the 49 methyl groups of tuna cytochrome c, the assignment of 33 is made without reference to the X-ray crystal structure. The method should therefore be applicable to other proteins of similar size where X-ray structures are unavailable. The assignments will be used to investigate the structure of cytochrome c in solution.     

Differences in the solution structures of oxidized and reduced cytochrome c measured by small-angle X-ray scattering

While X-ray crystallographic data on cytochrome c show the reduced and oxidized forms to have very similar structures, there is a considerable body of data, mostly from solution studies, that indicates the reduced form is more stable and that the interior of the protein is less accessible to solvent in this state. These observations have led to the hypothesis that while the time-averaged structure is preserved between the two forms, the dynamics of the two forms are different. The oxidized form has been proposed to undergo more large-amplitude, low-frequency motions than the reduced form. The crystal structure data were derived from crystals grown in high salt concentrations, but the solution studies were done at relatively low ionic strength. Small-angle X-ray scattering has been used to examine the effects of the ionic strength and oxidation state on the solution structure of cytochrome c. We find that the radius of gyration and the maximum linear dimension of oxidized cytochrome c are significantly larger than those for reduced cytochrome c, in 5 mM phosphate buffer at pH 7.3, and further that this difference is suppressed by addition of 200 mM sodium chloride. We conclude that there is a real structural difference between the two forms at low ionic strength in solution and that this difference is likely to contribute to the observed differences in accessibility and compressibility.     

Resonance Raman investigations of cytochrome c conformational change upon interaction with the membranes of intact and Ca2+-exposed mitochondria

The conformational states of cytochrome c inside intact and Ca(2+)-exposed mitochondria have been investigated using resonance Raman spectroscopy. Intact and swelling bovine heart and rat liver mitochondria were examined with an excitation wavelength (413.1 nm) in resonance with the Soret transition of ferrous cytochrome c. The different b- to c-type cytochrome concentration ratio in mitochondria from two different tissues was used to help assign the Raman spectral components. Resonance Raman spectra were also recorded for mitochondria fractions (supernatants and pellets) obtained from swollen (Ca(2+)-exposed) mitochondria after differential centrifugation. The results illustrate that cytochrome c has an altered vibrational spectrum in solution, in intact, and in swollen mitochondria. When cytochrome c is released from mitochondria, its Raman spectrum becomes identical to that of ferrous cytochrome c in solution. The spectra of mitochondrial pellets indicate that a small amount of structurally modified cytochrome c remains associated with the heavy membrane fraction. Indeed, spectroscopic shifts in the low-frequency fingerprint and the high-frequency marker-band regions suggest that membrane binding leads to a partial opening of the heme pocket and an alteration of the heme thioether bonds. The results support the conclusion that most cytochrome c molecules in mitochondria are membrane-bound and that the cytochrome c structure changes upon binding. Furthermore, changes in the resonance Raman active mode located at 675 cm(-)(1) in the spectra of intact, swollen, and fractionated mitochondria indicate that b-type cytochromes may also undergo structural alterations during mitochondrial swelling and disruption.     

Comparison of the solution and crystal structures of mitochondrial cytochrome c. Analysis of paramagnetic shifts in the nuclear magnetic resonance spectrum of ferricytochrome c

The two accompanying papers describe the assignment of methyl-containing spin-systems in the 1H nuclear magnetic resonance spectra of tuna ferricytochrome c and tuna ferrocytochrome c. At present, 104 resonances from 208 C-H protons are assigned in both oxidation states. In this paper, the difference in chemical shift of a resonance between the two oxidation states is used together with a dipolar model of the unpaired electron spin of ferricytochrome c to compare the structure of cytochrome c in solution with three high-resolution structures of cytochrome c obtained by X-ray diffraction in single crystals. The overall protein fold and the positions of most of the haem-packing residues are shown to be invariant between the crystal and solution. However, three regions of the protein, at the C terminus, around the haem propionic acid groups and at the haem crevice near thioether-2, are found to undergo conformational changes on the removal of crystal packing constraints.     

Family of cytochrome c7-type proteins from Geobacter sulfurreducens: structure of one cytochrome c7 at 1.45 A resolution

The structure of a cytochrome c(7) (PpcA) from Geobacter sulfurreducens was determined by X-ray diffraction at 1.45 A resolution; the R factor is 18.2%. The protein contains a three-heme core that is surrounded by 71 amino acid residues. An unusual feature of this cytochrome is that it has 17 lysine residues, but only nine hydrophobic residues that are larger than alanine. The details of the structure are described and compared with those of cytochrome c(7) from Desulfuromonas acetoxidans and with cytochromes c(3). The two cytochrome c(7) molecules have sequences that are 46% identical, but the arrangements of the hemes in the two structures differ; the rms deviation of all alpha-carbons is 2.5 A. These cytochromes can reduce various metal ions. The reduction site of the chromate ion in D. acetoxidans is occupied by a sulfate ion in the crystal structure of PpcA. We identified four additional homologues of cytochrome c(7) in the G. sulfurreducens genome and three polymers of c(7)-type domains. Of the polymers, two have four repeats and one has nine repeats. On the basis of sequence alignments, one of the hemes in each of the cytochrome c(7)-type domains does not have the bis-histidine coordination. The packing of the molecules in the crystal structure of PpcA suggests that the polymers have an elongated conformation and might form a "nanowire".     

Time-resolved fluorescence and circular dichroism of porphyrin cytochrome c and Zn-porphyrin cytochrome c incorporated in reversed micelles

Interactions between fluorescent horse heart cytochrome c derivatives (e. g. porphyrin cytochrome c and Zn-porphyrin cytochrome c) with surfactant interfaces in reversed micellar solutions have been studied, using different spectroscopic techniques. Anionic [sodium bis(2-ethylhexyl)sulfosuccinate, AOT] and cationic (cetyltrime-thylammonium bromide, CTAB) surfactant solutions have been used in order to investigate the effects of charge interactions between proteins and interfaces. Circular dichroism reveals that much of the protein secondary structure is lost in AOT-reversed micelles, especially when the molar water/surfactant ratio, wo, is high (wo = 40), whereas in CTAB-reversed micelles secondary structure seems to be preserved. Time-resolved fluorescence measurements of the porphyrin in the cytochrome c molecule yields information about the changes in structure and the dynamics of the protein upon interaction with surfactant assemblies both in aqueous and in hydrocarbon solutions. With AOT as surfactant a strong interaction between protein and interface can be observed. The effects found in aqueous AOT solution are of the same kind as in hydrocarbon solution. In the CTAB systems the interactions between protein and surfactant are much less pronounced. The measured effects on the fluorescence properties of the proteins are different in aqueous and hydrocarbon solutions. In general, the observations can be explained by an electrostatic attraction between the overall positively charged protein molecules and the anionic AOT interface. Electrostatic attraction can also occur between the cytochrome c derivatives and CTAB because there is a negatively charged zone on the surface of the proteins. From the fluorescence anisotropy decays it can be concluded that in the CTAB-reversed micellar system these interactions are not important, whereas in an aqueous CTAB solution the proteins interact with surfactant molecules.     

A study of the electron transfer properties of the heme undecapeptide from cytochrome c by 1H nmr spectroscopy

Nuclear magnetic resonance (nmr) spectroscopy has been used to investigate the heme undecapeptide from cytochrome c. Assignments of resonances to specific residues have been made based on spin decoupling, redox titration, and the pH and temperature dependence of resonance lines. An outline structure is presented based on the assignments, secondary shift data, and the x-ray crystal structure of cytochrome c. An equation is derived to relate the width of an nmr line during a redox titration to the percentage of each oxidation state. Using this equation the self-exchange rate constant for electron transfer for the heme peptide is 1.3 x 10(7) M-1 sec-1 at 330 degrees K. Discussion of the self-exchange rate constants of cytochrome c, cytochrome c3, and cytochrome c551 is related to this constant for the heme undecapeptide.     

Interaction of yeast iso-1-cytochrome c with cytochrome c peroxidase investigated by [15N, 1H] heteronuclear NMR spectroscopy

The interaction of yeast iso-1-cytochrome c with its physiological redox partner cytochrome c peroxidase has been investigated using heteronuclear NMR techniques. Chemical shift perturbations for both 15N and 1H nuclei arising from the interaction of isotopically enriched 15N cytochrome c with cytochrome c peroxidase have been observed. For the diamagnetic, ferrous cytochrome c, 34 amides are affected by binding, corresponding to residues at the front face of the protein and in agreement with the interface observed in the 1:1 crystal structure of the complex. In contrast, for the paramagnetic, ferric protein, 56 amides are affected, corresponding to residues both at the front and toward the rear of the protein. In addition, the chemical shift perturbations were larger for the ferric protein. Using experimentally observed pseudocontact shifts the magnetic susceptibility tensor of yeast iso-1-cytochrome c in both the free and bound forms has been calculated with HN nuclei as inputs. In contrast to an earlier study, the results indicate that there is no change in the geometry of the magnetic axes for cytochrome c upon binding to cytochrome c peroxidase. This leads us to conclude that the additional effects observed for the ferric protein arise either from a difference in binding mode or from the more flexible overall structure causing a transmittance effect upon binding.     

Effect of single-site charge-reversal mutations on the catalytic properties of yeast cytochrome c peroxidase: mutations near the high-affinity cytochrome c binding site

Fifteen single-site charge-reversal mutations of yeast cytochrome c peroxidase (CcP) have been constructed to determine the effect of localized charge on the catalytic properties of the enzyme. The mutations are located on the front face of CcP, near the cytochrome c binding site identified in the crystallographic structure of the yeast cytochrome c-CcP complex [Pelletier, H., and Kraut, J. (1992) Science 258, 1748-1755]. The mutants are characterized by absorption spectroscopy and hydrogen peroxide reactivity at both pH 6.0 and 7.5 and by steady-state kinetic studies using recombinant yeast iso-1-ferrocytochrome c(C102T) as a substrate at pH 7.5. Some of the charge-reversal mutations cause detectable changes in the absorption spectrum, especially at pH 7.5, reflecting changes in the equilibrium between penta- and hexacoordinate heme species in the enzyme. An increase in the amount of hexacoordinate heme in the mutant enzymes correlates with an increase in the fraction of enzyme that does not react with hydrogen peroxide. Steady-state velocity measurements indicate that five of the 15 mutations cause large increases in the Michaelis constant (R31E, D34K, D37K, E118K, and E290K). These data support the hypothesis that the cytochrome c-CcP complex observed in the crystal is the dominant catalytically active complex in solution.     

Change of cytochrome c structure during development of the mouse.

The structure of cytochrome c during mouse development is investigated. For this purpose the amino acid sequence of cytochrome c of the adult mouse had to be determined. The structure of cytochrome c of adult differentiated mouse cells differs in two amino acid residues from the known amino acid sequence of rabbit cytochrome c. No indication of different forms of cytochrome c in the adult differentiated cells was obtained. The structure of cytochrome c from 11.5-day-old mouse embryos is identical with that of adult mouse tissues. Since germ cells after meiotic division are the immediate precursors of a new individual, the structure of cytochrome c from sperm-containing mice testes was investigated. By means of chromatography of the cytochrome c and of peptide maps and amino acid analyses of its tryptic peptides, it is shown that mouse testis contains two isocytochromes c in about equal amount. The structure of one of these two isocytochromes c is identical with the structure of the adult-type cytochrome c of mouse. The testis-specific cytochrome c, which is assumed to be located in the sperm cells, differs in 13 of its 104 amino acid residues from the adult-type cytochrome c. From comparison of the primary and the spatial structures of the adult-type and the sperm-type isocytochromes c with the known structures of cytochrome c of more than 65 different species it is concluded that the duplication of the cytochrome c structural gene, causing the existence of the two ontogenetic-specific isocytochromes c in mouse, has occurred early in the evolution of eucaryotes.     

Crystal structure of oxidized cytochrome c(6A) from Arabidopsis thaliana

Compared with algal and cyanobacterial cytochrome c(6), cytochrome c(6A) from higher plants contains an additional loop of 12 amino acid residues. We have determined the first crystal structure of cytochrome c(6A) from Arabidopsis thaliana at 1.5 Angstrom resolution in order to help elucidate its function. The overall structure of cytochrome c(6A) follows the topology of class I c-type cytochromes in which the heme prosthetic group covalently binds to Cys16 and Cys19, and the iron has octahedral coordination with His20 and Met60 as the axial ligands. Two cysteine residues (Cys67 and Cys73) within the characteristic 12 amino acids loop form a disulfide bond, contributing to the structural stability of cytochrome c(6A). Our model provides a chemical basis for the known low redox potential of cytochrome c(6A) which makes it an unsuitable electron carrier between cytochrome b(6)f and PSI.     

Solution structure of cytochrome c6 from the thermophilic cyanobacterium Synechococcus elongatus.

Cytochrome c6 is a small, soluble electron carrier between the two membrane-bound complexes cytochrome b6f and photosystem I (PSI) in oxygenic photosynthesis. We determined the solution structure of cytochrome c6 from the thermophilic cyanobacterium Synechococcus elongatus by NMR spectroscopy and molecular dynamics calculations based on 1586 interresidual distance and 28 dihedral angle restraints. The overall fold exhibits four alpha-helices and a small antiparallel beta-sheet in the vicinity of Met58, one of the axial heme ligands. The flat hydrophobic area in this cytochrome c6 is conserved in other c6 cytochromes and even in plastocyanin of higher plants. This docking region includes the site of electron transfer to PSI and possibly to the cytochrome b6f complex. The binding of cytochrome c6 to PSI in green algae involves interaction of a negative patch with a positive domain of PSI. This positive domain has not been inserted at the evolutionary level of cyanobacteria, but the negatively charged surface region is already present in S. elongatus cytochrome c6 and may thus have been optimized during evolution to improve the interaction with the positively charged cytochrome f. As the structure of PSI is known in S.elongatus, the reported cytochrome c6 structure can provide a basis for mutagenesis studies to delineate the mechanism of electron transfer between both.     

Preliminary X-ray studies of the tetra-heme cytochrome c3 and the octa-heme cytochrome c3 from Desulfovibrio gigas

A tetra-heme and an octa-heme cytochrome c3 from the sulfate bacterium Desulfovibrio gigas have been crystallized. Diffraction quality crystals of the tetra-heme cytochrome are obtained from solution by the addition of polyethylene glycol at pH 6.5. The crystals are orthorhombic, space group P2(1)2(1)2 with unit cell parameters a = 42.27 A, b = 52.54 A and c = 52.83 A. The octa-heme cytochrome crystals develop from low ionic strength solutions of phosphate or Tris-Cl in the pH range 6.2-7.6. The crystals belong to the trigonal system, space group P3(1) or the enantiomorph P3(2), with unit cell parameters a = b = 57.4 A, c = 97.3 A, gamma = 120 degrees. Single crystal diffraction studies of the structures of these two low-potential cytochromes are in progress.     

Crystal structure of Azotobacter cytochrome c5 at 2.5 A resolution

The crystal structure of cytochrome c5 from Azotobacter vinelandii has been solved and refined to an R value of 0.29 at 2.5 A resolution. The structure of the oxidized protein was solved using a monoclinic crystal form. The structure was solved by multiple isomorphous replacements, re-fit to a solvent-leveled multiple isomorphous replacement map, and refined by restrained least squares. The structure reveals monomers associated about the crystallographic 2-fold axis by hydrophobic contacts at the "exposed heme edge". The overall conformation for the monomer is similar to that of Pseudomonas aeruginosa cytochrome c551. However, relative to a common heme conformation, c5 and c551 differ by an average of 6.8 A over 82 alpha-carbon positions and the propionates of c5 are much more exposed to solvent. The shortest heme--heme contact at the "dimer" interface is 6.3 A (Fe to Fe 16.4 A). Alignment of c5 and c551 shows that the two cytochromes, in spite of sequence differences, have remarkably similar charge distributions. A disulfide stacks on a tyrosine between the N- and C-terminal helices.     

Thermodynamic volume cycles for electron transfer in the cytochrome c oxidase and for the binding of cytochrome c to cytochrome c oxidase.

Dilatometry is a sensitive technique for measuring volume changes occurring during a chemical reaction. We applied it to the reduction-oxidation cycle of cytochrome c oxidase, and to the binding of cytochrome c to the oxidase. We measured the volume changes that occur during the interconversion of oxidase intermediates. The numerical values of these volume changes have allowed the construction of a thermodynamic cycle that includes many of the redox intermediates. The system volume for each of the intermediates is different. We suggest that these differences arise by two mechanisms that are not mutually exclusive: intermediates in the catalytic cycle could be hydrated to different extents, and/or small voids in the protein could open and close. Based on our experience with osmotic stress, we believe that at least a portion of the volume changes represent the obligatory movement of solvent into and out of the oxidase during the combined electron and proton transfer process. The volume changes associated with the binding of cytochrome c to cytochrome c oxidase have been studied as a function of the redox state of the two proteins. The volume changes determined by dilatometry are large and negative. The data indicate quite clearly that there are structural alterations in the two proteins that occur on complex formation.     

Effects of extrinsic imidazole ligation on the molecular and electronic structure of cytochrome c

Although imidazole ligand binding to cytochrome c is not directly related to its physiological function, it has the potential to provide valuable information on the molecular and electronic structure of the protein. The solution structure of the imidazole adduct of oxidized horse heart cytochrome c (Im-cyt c) has been determined through 2D NMR spectroscopy. The Im-cyt c, 8 mM in 1.2 M imidazole solution at pH 5.7 and 313 K, provided altogether 2,542 NOEs (1,901 meaningful NOEs) and 194 pseudocontact shifts. The 35 conformers of the family show the RMSD values to the average structure of 0.063+/-0.007 nm for the backbone and 0.107+/-0.007 nm for all heavy atoms, respectively. The characterization of Im-cyt c is discussed in detail both in terms of structure and electronic properties. The replacement of the axial ligand Met80 with the exogenous imidazole ligand induces significant conformation changes in both backbone and side chains of the residues located in the distal axial ligand regions. The imidazole ligand binds essentially parallel to the imidazole of the proximal histidine, the two planes forming an angle of 8+/-7 degrees. The electron delocalization on the heme moiety and the magnetic susceptibility tensor are consistent with these structural features.