ATP specifically drives refolding of non-native conformations of cytochrome c

An increasing body of evidence ascribes to misfolded forms of cytochrome c (cyt c) a role in pathophysiological events such as apoptosis and disease. Here, we examine the conformational changes induced by lipid binding to horse heart cyt c at pH 7 and study the ability of ATP (and other nucleotides) to refold several forms of unfolded cyt c such as oleic acid-bound cyt c, nicked cyt c, and acid denatured cyt c. The CD and fluorescence spectra demonstrate that cyt c unfolded by oleic acid has an intact secondary structure, and a disrupted tertiary structure and heme environment. Furthermore, evidence from the Soret CD, electronic absorption, and resonance Raman spectra indicates the presence of an equilibrium of at least two low-spin species having distinct heme-iron(III) coordination. As a whole, the data indicate that binding of cyt c to oleic acid leads to a partially unfolded conformation of the protein, resembling that typical of the molten globule state. Interestingly, the native conformation is almost fully recovered in the presence of ATP or dATP, while other nucleotides, such as GTP, are ineffective. Molecular modeling of ATP binding to cyt c and mutagenesis experiments show the interactions of phosphate groups with Lys88 and Arg91, with adenosine ring interaction with Glu62 explaining the unfavorable binding of GTP. The finding that ATP and dATP are unique among the nucleotides in being able to turn non-native states of cyt c back to native conformation is discussed in the light of cyt c involvement in cell apoptosis.     

Temperature and driving force dependence of the folding rate of reduced horse heart cytochrome c

Utilizing the stability difference between the ferro and ferri forms of horse heart cytochrome c (cyt c), folding of reduced cyt c was triggered by laser-induced reduction of unfolded oxidized cyt c. Measurements were made of the kinetics of the main folding phase (1 ms-10 s) in which collapsed reduced cyt c transforms to the native conformation. The folding rates were studied extensively as a function of temperature (5-75 degrees C) and guanidine hydrochloride (GdnHCl) concentration (1.6-4.9 M). At constant [GdnHCl], the Arrhenius plot of the folding rate constant (k) is nonlinear. At temperatures above 40 degrees C, the decrease in protein stability counteracts the expected increase in folding rate. Introducing free energy (DeltaG), derived from protein stability data, into the Eyring and Arrhenius equations leads to: ln k = ln(k(b)T/h) + DeltaS()/R - DeltaH()/RT - theta(m)DeltaG/RT = ln A - E(a)/RT - theta(m)DeltaG/RT, where theta(m) is the ratio between the denaturant dependence of the folding rate and the stability. By using this equation at constant DeltaG [or constant equilibrium constant (K)], linear Arrhenius plots are obtained. For the main folding phase of reduced cyt c, a positive DeltaS() is obtained indicating that the transition state is less ordered than the reactant. A model is proposed in which reduced cyt c first collapses into a compact intermediate, which needs to expand to reach the transition state of the rate-limiting folding reaction.     

Highly nonexponential kinetics in the early-phase refolding of proteins at low temperatures

Theory and simulations predict that the folding kinetics of protein-like heteropolymers become nonexponential and glassy (i.e., controlled by escape from different low-energy misfolded states) at low temperatures, but there was little experimental evidence for such behavior of proteins. We have developed a stopped-flow instrument working reliably down to -40 degrees C with high mixing capability and applied it to study the refolding kinetics of horse cytochrome c (cyt c) and hen egg white lysozyme at temperatures below 0 degrees C in the presence of antifreeze NaCl, LiCl, or ethylene glycol and above 0 degrees C in the presence and absence of antifreeze. The refolding was initiated by rapid dilution of the guanidine hydrochloride unfolded proteins, and the kinetics were monitored by intrinsic tryptophan fluorescence. Highly nonexponential kinetics extended over 3 decades in time (0.01-10 s) were observed in the early phases of the refolding of cyt c and lysozyme in the temperature range of -35 to 5 degrees C. These results are in agreement with the theoretical prediction, suggesting that the folding energy landscapes of these proteins are rugged in the upper portions.     

Cytochrome c(553), a small heme protein that lacks misligation in its unfolded state, folds with rapid two-state kinetics

Cytochrome c(553) (cyt c(553)) from Desulfovibrio vulgaris is a small helical heme protein that displays apparent two-state equilibrium-unfolding behavior. The covalently attached heme is low-spin, ligated by Met and His residues, in the native state but becomes high-spin upon unfolding at pH 7. Here, we show that in contrast to other c-type heme proteins, where misligations in the unfolded states are prominent, cyt c(553) refolding kinetics at pH 7 proceeds rapidly without detectable intermediates. The extrapolated folding rate constant in water for oxidized cyt c(553) matches exactly that predicted from the cyt c(553) native-state topology: 5300 s(-1 )(experimental) versus 5020 s(-1) (predicted). We therefore conclude that the presence of the oxidized cofactor does not affect the intrinsic formation speed of the cyt c(553 )structural motif.     

A novel method for study of protein folding kinetics by monitoring diffusion coefficient in time domain

Molecular diffusion process after the photo-induced electron injection to ferric cytochrome c (Fe(III) cyt c) in guanidine hydrochloride (GdnHCl) 3.5 M buffer solution is studied by the time-resolved transient grating technique. Circular dichroism studies have revealed that Fe(III) cyt c is unfolded under this condition but the reduced form, Fe(II) cyt c, is folded. Hence, this pulsed laser-induced reduction should initiate the folding process of cyt c. The observed transient grating signal shows prominent features, which have never been observed before. Based on several characteristic points, we concluded that the apparent diffusion coefficient (D) of Fe(II) cyt c after the reduction is time dependent, which must be associated with the protein folding dynamics. This time-dependent apparent D should reflect either the continuous time development of the hydrodynamic radius or population change of the unfolded and folded states during the folding dynamics. This is the first observation of the time-dependent apparent D during any chemical reaction, and this time-dependent measurement of D should be a unique and powerful way to study the protein folding kinetics from a viewpoint of the protein's shape or the protein-water intermolecular interaction.     

Brownian dynamics simulations of the interaction of Chlamydomonas cytochrome f with plastocyanin and cytochrome c6

The interaction of Chlamydomonas cytochrome f (cyt f) with either Chlamydomonas plastocyanin (PC) or Chlamydomonas cytochrome c(6) (cyt c(6)) was studied using Brownian dynamics simulations. The two electron acceptors (PC and cyt c(6)) were found to be essentially interchangeable despite a lack of sequence homology and different secondary structures (beta-sheet for PC and alpha-helix for cyt c(6)). Simulations using PC and cyt c(6) interacting with cyt f showed approximately equal numbers of successful complexes and calculated rates of electron transfer. Cyt f-PC and cyt f-cyt c(6) showed the same types of interactions. Hydrophobic residues surrounding the Y1 ligand to the heme on cyt f interacted with hydrophobic residues on PC (surrounding the H87 ligand to the Cu) or cyt c(6) (surrounding the heme). Both types of complexes were stabilized by electrostatic interactions between K65, K188, and K189 on cyt f and conserved anionic residues on PC (E43, D44, D53, and E85) or cyt c(6) (E2, E70, and E71). Mutations on cyt f had identical effects on its interaction with either PC or cyt c(6). K65A, K188A, and K189A showed the largest effects whereas residues such as K217A, R88A, and K110A, which are located far from the positive patch on cyt f, showed very little inhibition. The effect of mutations observed in Brownian dynamics simulations paralleled those observed in experiments.     

Substitution of the sixth axial ligand of Rhodobacter capsulatus cytochrome c1 heme yields novel cytochrome c1 variants with unusual properties.

The cytochrome (cyt) c1 heme of the ubihydroquinone:cytochrome c oxidoreductase (bc1 complex) is covalently attached to two cysteine residues of the cyt c1 polypeptide chain via two thioether bonds, and the fifth and sixth axial ligands of its iron atom are histidine (H) and methionine (M), respectively. The latter residue is M183 in Rhodobacter capsulatus cyt c1, and previous mutagenesis studies revealed its critical role for the physicochemical properties of cyt c1 [Gray, K. A., Davidson, E., and Daldal, F. (1992) Biochemistry 31, 11864-11873]. In the homologous chloroplast b6f complex, the sixth axial ligand is provided by the amino group of the amino terminal tyrosine residue. To further pursue our investigation on the role played by the sixth axial ligand in heme-protein interactions, novel cyt c1 variants with histidine-lysine (K) and histidine-histidine axial coordination were sought. Using a R. capsulatus genetic system, the cyt c1 mutants M183K and M183H were constructed by site-directed mutagenesis, and chromatophore membranes as well as purified bc1 complexes obtained from these mutants were characterized in detail. The studies revealed that these mutants incorporated the heme group into the mature cyt c1 polypeptides, but yielded nonfunctional bc1 complexes with unusual spectroscopic and thermodynamic properties, including shifted optical absorption maxima (lambdamax) and decreased redox midpoint potential values (Em7). The availability and future detailed studies of these stable cyt c1 mutants should contribute to our understanding of how different factors influence the physicochemical and folding properties of membrane-bound c-type cytochromes in general.     

Cytochrome c is released in a single step during apoptosis

Release of cytochrome c from mitochondria is a central event in apoptotic signaling. In this study, we utilized a cytochrome c fusion that binds fluorescent biarsenical ligands (cytochrome c-4CYS (cyt. c-4CYS)) as well as cytochrome c-green fluorescent protein (cyt. c-GFP) to measure its release from mitochondria in different cell types during apoptosis. In single cells, the kinetics of cyt. c-4CYS release was indistinguishable from that of cyt. c-GFP in apoptotic cells expressing both molecules. Lowering the temperature by 7 degrees C did not affect this corelease, but further separated cytochrome c release from the subsequent decrease in mitochondrial membrane potential (DeltaPsi(m)). Cyt. c-GFP rescued respiration in cells lacking endogenous cytochrome c, and the duration of cytochrome c release was approximately 5 min in a variety of cell types induced to die by various forms of cellular stress. In addition, we could observe no evidence of caspase-dependent amplification of cytochrome c release or changes in DeltaPsi(m) preceding the release of cyt. c-GFP. We conclude that there is a general mechanism responsible for cytochrome c release that proceeds in a single step that is independent of changes in DeltaPsi(m).     

In situ kinetics of cytochromes c1 and c2

In Rhodobacter sphaeroides chromatophores, cytochromes (cyt) c(1) and c(2) have closely overlapping spectra, and their spectral deconvolution provides a challenging task. As a result, analyses of the kinetics of different cytochrome components of the bc(1) complex in purple bacteria usually report only the sum cyt c(1) + cyt c(2) kinetics. Here we used newly determined difference spectra of individual components to resolve the kinetics of cyt c(1) and c(2) in situ via a least-squares (LS) deconvolution. We found that the kinetics of cyt c(1) and c(2) are significantly different from those measured using the traditional difference wavelength (DW) approach, based on the difference in the absorbance at two different wavelengths specific for each component. In particular, with the wavelength pairs previously recommended, differences in instrumental calibration led to kinetics of flash-induced cyt c(1) oxidation measured with the DW method which were faster than those determined by the LS method (half-time of approximately 120 micros vs half-time of approximately 235 micros, in the presence of antimycin). In addition, the LS approach revealed a delay of approximately 50 micros in the kinetics of cyt c(1) oxidation, which was masked when the DW approach was used. We attribute this delay to all processes leading to the oxidation of cyt c(1) after light activation of the photosynthetic reaction center, especially the dissociation of cyt c(2) from the reaction center. We also found that kinetics of both cyt c(1) and c(2) measured by the DW approach were significantly distorted at times longer than 1 ms, due to spectral contamination from changes in the b hemes. The successful spectral deconvolution of cyt c(1) and c(2), and inclusion of both cytochromes in the kinetic analysis, significantly increase the data available for mechanistic understanding of bc(1) turnover in situ.     

Kinetics and mechanism of the reduction of horse heart ferricytochrome c by glutathione.

A detailed investigation of the reduction of cytochrome c by glutathione has shown that the reaction proceeds through several steps. A rapid combination of the reducing agent with the cytochrome leads to the formation of a glutathione-cytochrome intermediate in which the glutathione most likely interacts with the edge of the heme moiety. The electron transfer takes place in a subsequent slower step. Since cytochrome c(III) exists in two conformational forms at neutral pH [Kujundzic, N., & Everse, J. (1978) Biochem. Biophys. Res. Commun. 82, 1211], the reduction of cytochrome c by glutathione may be represented by cyt c(III) + GS- reversible K1 cyt c(III) ... GS- reversible k1 products cyt c*(III) + GS- reversible K2 cyt c*(III) ... GS- reversible k2 products At 25 degrees C, pH 7.5, and an ionic strength of 1.0 (NaCl), k1 = 1.2 X 10(-3) S-1, k2 = 2.0 X 10(-3) S-1, k1 = 2.9 X 10(3) M-1, and K2 = 5.3 X 10(3) M-1. The reaction is catalyzed by trisulfides, and second-order rate constants of 4.55 X 10(3) and 7.14 X 10(3) M-1 S-1 were obtained for methyl trisulfide and cysteine trisulfide, respectively.     

Yeast cytochrome c is a sequence-specific DNA-binding protein

A protein binding to the alcohol oxidase 2 upstream activation sequence (AOX2UAS) of the methylotropic yeast, Pichia pastoris, has been purified and identified as cytochrome c (cyt c). Cyt c purified from P. pastoris or Saccharomyces cerevisiae binds to AOX2UAS. Specific point mutations in AOX2UAS abolish cyt c binding. We conclude that yeast cyt c is a sequence-specific DNA-binding protein and may have a regulatory role in the nucleus.     

Nature of the fast and slow refolding reactions of iron(III) cytochrome c

The fast and slow refolding reactions of iron(III) cytochrome c (Fe(III) cyt c), previously studied by Ikai et al. (Ikai, A., Fish, W. W., & Tanford, C. (1973) J. Mol. Biol. 73, 165--184), have been reinvestigated. The fast reaction has the major amplitude (78%) and is 100-fold faster than the slow reaction in these conditions (pH 7.2, 25 degrees C, 1.75 M guanidine hydrochloride). We show here that native cyt c is the product formed in the fast reaction as well as in the slow reaction. Two probes have been used to test for formation of native cyt c. absorbance in the 695-nm band and rate of reduction of by L-ascorbate. Different unfolded species (UF, US) give rise to the fast and slow refolding reactions, as shown both by refolding assays at different times after unfolding ("double-jump" experiments) and by the formation of native cyt c in each of the fast and slow refolding reactions. Thus the fast refolding reaction is UF leads to N and the slow refolding reaction is Us leads to N, where N is native cyt c, and there is a US in equilibrium UF equilibrium in unfolded cyt c. The results are consistent with the UF in equilibrium US reaction being proline isomerization, but this has not yet been tested in detail. Folding intermediates have been detected in both reactions. In the UF leads to N reaction, the Soret absorbance change precedes the recovery of the native 695-nm band spectrum, showing that Soret absorbance monitors the formation of a folding intermediate. In the US leads to N reaction an ascorbate-reducible intermediate has been found at an early stage in folding and the Soret absorbance change occurs together with the change at 695 nm as N is formed in the final stage of folding.     

Characterization of the interaction of Rhodobacter capsulatus cytochrome c peroxidase with charge reversal mutants of cytochrome c(2)

Steady-state kinetics for the reaction of Rhodobacter capsulatus bacterial cytochrome c peroxidase (BCCP) with its substrate cytochrome c(2) were investigated. The Rb. capsulatus BCCP is dependent on calcium for activation as previously shown for the Pseudomonas aeruginosa BCCP and Paracoccus denitrificans enzymes. Furthermore, the activity shows a bell-shaped pH dependence with optimum at pH 7.0. Enzyme activity is greatest at low ionic strength and drops off steeply as ionic strength increases, resulting in an apparent interaction domain charge product of -13. All cytochromes c(2) show an asymmetric distribution of surface charge, with a concentration of 14 positive charges near the exposed heme edge of Rb. capsulatus c(2) which potentially may interact with approximately 6 negative charges, localized near the edge of the high-potential heme of the Rb. capsulatus BCCP. To test this proposal, we constructed charge reversal mutants of the 14 positively charged residues located on the front face of Rb. capsulatus cytochrome c(2) and examined their effect on steady-state kinetics with BCCP. Mutated residues in Rb. capsulatus cytochrome c(2) that showed the greatest effects on binding and enzyme activity are K12E, K14E, K54E, K84E, K93E, and K99E, which is consistent with the site of electron transfer being located at the heme edge. We conclude that a combination of long-range, nonspecific electrostatic interactions as well as localized salt bridges between, e.g., cytochrome c(2) K12, K14, K54, and K99 with BCCP D194, D241, and D6, account for the observed kinetics.     

Characterization of a new membrane-bound cytochrome c of Rhodopseudomonas capsulata

A cytochrome c (cyt. c) was solubilized with Triton-X-100 and co-purified with cytochrome c oxidase from membranes of chemotrophically grown cells of Rhodopseudomonas capsulata. Cyt. c and cytochrome oxidase were separated on Sephadex G-50 columns. Antibodies against cytochrome c2 from the same bacterium did not cross react with the membrane-bound cyt. c. The IEP of the membrane-bound cyt. c was found to be pH 8.2, the midpoint potential was 234 +/- 11 mV at pH 7.0. This cyt. c binds CO. The native cyt. c is a dimer with an apparent Mr of 25000 containing 2 mol heme per mol dimer, which is believed to function as an electron donor for the high-potential cytochrome c oxidase.     

Brownian dynamics study of cytochrome f interactions with cytochrome c6 and plastocyanin in Chlamydomonas reinhardtii plastocyanin, and cytochrome c6 mutants

Using Brownian dynamics simulations, all of the charged residues in Chlamydomonas reinhardtii cytochrome c(6) (cyt c(6)) and plastocyanin (PC) were mutated to alanine and their interactions with cytochrome f (cyt f) were modeled. Systematic mutation of charged residues on both PC and cyt c(6) confirmed that electrostatic interactions (at least in vitro) play an important role in bringing these proteins sufficiently close to cyt f to allow hydrophobic and van der Waals interactions to form the final electron transfer-active complex. The charged residue mutants on PC and cyt c(6) displayed similar inhibition classes. Our results indicate a difference between the two acidic clusters on PC. Mutations D44A and E43A of the lower cluster showed greater inhibition than do any of the mutations of the upper cluster residues. Replacement of acidic residues on cyt c(6) that correspond to the PC's lower cluster, particularly E70 and E69, was observed to be more inhibitory than those corresponding to the upper cluster. In PC residues D42, E43, D44, D53, D59, D61, and E85, and in cyt c(6) residues D2, E54, K57, D65, R66, E70, E71, and the heme had significant electrostatic contacts with cyt f charged residues. PC and cyt c(6) showed different binding sites and orientations on cyt f. As there are no experimental cyt c(6) mutation data available for algae, our results could serve as a good guide for future experimental work on this protein. The comparison between computational values and the available experimental data (for PC-cyt f interactions) showed overall good agreement, which supports the predictive power of Brownian dynamics simulations in mutagenesis studies.     

Electrostatic properties of cytochrome f: implications for docking with plastocyanin.

The electrostatic properties of cytochrome f (cyt f), a member of the cytochrome b6f complex and reaction partner with plastocyanin (PC) in photosynthetic electron transport, are qualitatively studied with the goal of determining the mechanism of electron transfer between cyt f and PC. A crystal structure for cyt f was analyzed with the software package GRASP, revealing a large region of positive potential generated by a patch of positively charged residues (including K58, K65, K66, K122, K185, K187, and R209) and reinforced by the iron center of the heme. This positive field attracts the negative charges of the two acidic patches on the mobile electron carrier PC. Three docked complexes are obtained for the two proteins, based on electrostatic or hydrophobic interactions or both and on steric fits by manual docking methods. The first of these three complexes shows strong electrostatic interactions between K187 on cyt f and D44 on PC and between E59 on PC and K58 on cyt f. Two other manually docked complexes are proposed, implicating H87 on PC as the electron-accepting site from the iron center of cyt f through Y1. The second complex maintains the D44/K187 cross-link (but not the E59/K58 link) while increasing hydrophobic interactions between PC and cyt f. Hydrophobic interactions are increased still further in the third complex, whereas the link between K187 on cyt f and D44 on PC is broken. The proposed reaction mechanism, therefore, involves an initial electrostatic docking complex that gives rise to a nonpolar attraction between the regions surrounding H87 on PC and Y1 on cyt f, providing for an electron-transfer active complex.     

Membrane-associated cytochrome cy of Rhodobacter capsulatus is an electron carrier from the cytochrome bc1 complex to the cytochrome c oxidase during respiration.

We have recently established that the facultative phototrophic bacterium Rhodobacter capsulatus has two different pathways for reduction of the photooxidized reaction center during photosynthesis (F.E. Jenney and F. Daldal, EMBO J. 12:1283-1292, 1993; F.E. Jenney, R.C. Prince, and F. Daldal, Biochemistry 33:2496-2502, 1994). One pathway is via the well-characterized, water-soluble cytochrome c2 (cyt c2), and the other is via a novel membrane-associated c-type cytochrome named cyt cy. In this work, we probed the role of cyt cy in respiratory electron transport by isolating a set of R. capsulatus mutants lacking either cyt c2 or cyt cy, in the presence or in the absence of a functional quinol oxidase-dependent alternate respiratory pathway. The growth and inhibitor sensitivity patterns of these mutants, their respiratory rates in the presence of specific inhibitors, and the oxidation-reduction kinetics of c-type cytochromes monitored under appropriate conditions demonstrated that cyt cy, like cyt c2, connects the bc1 complex and the cyt c oxidase during respiratory electron transport. Whether cyt c2 and cyt cy are the only electron carriers between these two energy-transducing membrane complexes of R. capsulatus is unknown.     

A Brownian dynamics study of the effects of cytochrome f structure and deletion of its small domain in interactions with cytochrome c6 and plastocyanin in Chlamydomonas reinhardtii

The availability of seven different structures of cytochrome f (cyt f) from Chlamydomonas reinhardtii allowed us, using Brownian dynamics simulations, to model interactions between these molecules and their redox partners, plastocyanin (PC) and cytochrome c6 (cyt c6) in the same species to study the effect of cyt f structure on its function. Our results showed that different cyt f structures, which are very similar, produced different reaction rates in interactions with PC and cyt c6. We were able to attribute this to structural differences among these molecules, particularly to a small flexible loop between A-184 and G-191 (which has some of the highest crystallographic temperature factors in all of the cyt f structures) on the cyt f small domain. We also showed that deletion of the cyt f small domain affected cyt c6 more than PC, due to their different binding positions on cyt f. One function of the small domain in cyt f may be to guide PC or cyt c6 to a uniform dock with cyt f, especially due to electrostatic interactions with K-188 and K-189 on this domain. Our results could serve as a good guide for future experimental work on these proteins to understand better the electron transfer process between them. Also, these results demonstrated the sensitivity and the power of the Brownian dynamics simulations in the study of molecular interactions.     

Spectral and kinetic resolution of the bc1 complex components in situ: a simple and robust alternative to the traditional difference wavelength approach

The kinetics of the cytochrome (cyt) components of the bc(1) complex (ubiquinol: cytochrome c oxidoreductase, Complex III) are traditionally followed by using the difference of absorbance changes at two or more different wavelengths. However, this difference-wavelength (DW) approach is of limited accuracy in the separation of absorbance changes of components with overlapping spectral bands. To resolve the kinetics of individual components in Rhodobacter sphaeroides chromatophores, we have tested a simplified version of a least squares (LS) analysis, based on measurement at a minimal number of different wavelengths. The success of the simplified LS analysis depended significantly on the wavelengths used in the set. The "traditional" set of 6 wavelengths (542, 551, 561, 566, 569 and 575 nm), normally used in the DW approach to characterize kinetics of cyt c(tot) (cyt c(1)+cyt c(2)), cyt b(L), cyt b(H), and P870 in chromatophores, could also be used to determine these components via the simplified LS analysis, with improved resolution of the individual components. However, this set is not sufficient when information about cyts c(1) and c(2) is needed. We identified multiple alternative sets of 5 and 6 wavelengths that could be used to determine the kinetics of all 5 components (P870 and cyts c(1), c(2), b(L), and b(H)) simultaneously, with an accuracy comparable to that of the LS analysis based on a full set of wavelengths (1 nm intervals). We conclude that a simplified version of LS deconvolution based on a small number of carefully selected wavelengths provides a robust and significant improvement over the traditional DW approach, since it accounts for spectral interference of the different components, and uses fewer measurements when information about all five individual components is needed. Using the simplified and complete LS analyses, we measured the simultaneous kinetics of all cytochrome components of bc(1) complex in the absence and presence of specific inhibitors and found that they correspond well to those expected from the modified Q-cycle. This is the first study in which the kinetics of all cytochrome and reaction center components of the bc(1) complex functioning in situ have been measured simultaneously, with full deconvolution over an extended time range.     

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.