Computer simulation of plastocyanin interaction with cytochrome f and photosystem I in cyanobacterium Phormidium laminosum

The paper is devoted to computer simulation of complex formation of protein plastocyanin with transmembrane pigment-protein complex photosystem I and subunit f of cytochrome b ₆ f complex in the cyanobacterium Phormidium laminosum. The computer algorithm considers diffusion and electrostatic interactions of protein molecules. The computer models have shown that electrostatic interactions in the cyanobacterium play a less important role than in higher plants because of different electrostatic potentials created by charged amino acid residues on the protein surfaces.     

Role of charges on cytochrome f from the cyanobacterium Phormidium laminosum in its interaction with plastocyanin

The role of charge on the surface of cytochrome f from the cyanobacterium Phormidium laminosum in the reaction with plastocyanin was investigated in vitro using site-directed mutagenesis. Charge was neutralized at five acidic residues individually and introduced at a residue close to the interface between the two proteins. The effects on the kinetics of the reaction were measured using stopped-flow spectrophotometry, and the midpoint potentials of the mutant proteins were determined. The dependence of the bimolecular rate constant of reaction, k(2), on ionic strength was determined for the reactions of the cytochrome f mutants with wild-type and mutant forms of plastocyanin. Double mutant cycle analysis was carried out to probe for the presence of specific electrostatic interactions. The effects of mutations on Cyt f were smaller than those seen previously for mutants of plastocyanin [Schlarb-Ridley, B. G. et al. (2002) Biochemistry 41, 3279-3285]. One specific short-range interaction between charged residues of wild-type plastocyanin (Arg93) and wild-type cytochrome f (Asp63) was identified. The kinetic evidence from this study and that of Schlarb-Ridley et al., 2002, appears to conflict with the NMR structure of the P. laminosum complex, which suggests the absence of electrostatic interactions in the final complex [Crowley, P. et al. (2001) J. Am. Chem. Soc. 123, 10444-10453]. The most likely explanation of the apparent paradox is that the overall rate is diffusion controlled and that electrostatics specifically influence the encounter complex and not the reaction complex.     

A Brownian dynamics study of the interaction of Phormidium cytochrome f with various cyanobacterial plastocyanins

Brownian dynamics simulations were used to study the role of electrostatic forces in the interactions of cytochrome f from the cyanobacterium Phormidium laminosum with various cyanobacterial plastocyanins. Both the net charge on the plastocyanin molecule and the charge configuration around H92 (H87 in higher plants) are important in determining the interactions. Those plastocyanins (PCs) with a net charge more negative than -2.0, including those from Synechococcus sp. PCC7942, Synechocystis sp. 6803, and P. laminosum showed very little complex formation. On the other hand, complex formation for those with a net charge more positive than -2.0 (including Nostoc sp. PCC7119 and Prochlorothrix hollandica) as well as Nostoc plastocyanin mutants showed a linear dependence of complex formation upon the net charge on the plastocyanin molecule. Mutation of charged residues on the surface of the PC molecules also affected complex formation. Simulations involving plastocyanin mutants K35A, R93A, and K11A (when present) showed inhibition of complex formation. In contrast, D10A and E17A mutants showed an increase in complex formation. All of these residues surround the H92 (H87 in higher plant plastocyanins) ligand to the copper. An examination of the closest electrostatic contacts shows that these residues interact with D63, E123, R157, D188, and the heme on Phormidium cytochrome f. In the complexes formed, the long axis of the PC molecule lies perpendicular to the long axis of cytochrome f. There is considerable heterogeneity in the orientation of plastocyanin in the complexes formed.     

Role of electrostatics in the interaction between cytochrome f and plastocyanin of the cyanobacterium Phormidium laminosum

The role of charged residues on the surface of plastocyanin from the cyanobacterium Phormidium laminosum in the reaction with soluble cytochrome f in vitro was studied using site-directed mutagenesis. The charge on each of five residues on the eastern face of plastocyanin was neutralized and/or inverted, and the effect of the mutation on midpoint potentials was determined. The dependence of the overall rate constant of reaction, k(2), on ionic strength was investigated using stopped-flow spectrophotometry. Removing negative charges (D44A or D45A) accelerated the reaction and increased the dependence on ionic strength, whereas removing positive charges slowed it down. Two mutations (K46A, K53A) each almost completely abolished any influence of ionic strength on k(2), and three mutations (R93A, R93Q, R93E) each converted electrostatic attraction into repulsion. At low ionic strength, wild type and all mutants showed an inhibition which might be due to changes in the interaction radius as a consequence of ionic strength dependence of the Debye length or to effects on the rate constant of electron transfer, k(et). The study shows that the electrostatics of the interaction between plastocyanin and cytochrome f of P. laminosum in vitro are not optimized for k(2). Whereas electrostatics are the major contributor to k(2) in plants [Kannt, A., et al. (1996) Biochim. Biophys. Acta 1277, 115-126], this role is taken by nonpolar interactions in the cyanobacterium, leading to a remarkably high rate at infinite ionic strength (3.2 x 10(7) M(-1) s(-1)).     

Brownian dynamics study of the interaction between plastocyanin and cytochrome f.

The electrostatic interaction between plastocyanin (PC) and cytochrome f (cyt f), electron transfer partners in photosynthesis was studied using Brownian dynamics (BD) simulations. By using the software package MacroDox, which implements the BD algorithm of Northrup et al. (Northrup, S. H., J. O. Boles, and J. C. L. Reynolds. 1987. J. Phys. Chem. 91:5991-5998), we have modeled the interaction of the two proteins based on crystal structures of poplar PC and turnip cyt f at pH 7 and a variety of ionic strengths. We find that the electrostatic attraction between positively charged residues (K58, K65, K187, and R209, among others) on cyt f and negatively charged residues (E43, D44, E59, and E60, among others) on PC steers PC into a single dominant orientation with respect to cyt f, and furthermore, that the single dominant orientation that we observe is one that we had predicted in our previous work (Pearson, D. C., E. L. Gross, and E. S. David. 1996. Biophys. J. 71:64-76). This dominant orientation permits the formation of hydrophobic interactions, which are not implemented in the MacroDox algorithm. This proposed complex between PC and cyt f implicates H87, a copper ligand on PC, as the residue that accepts electrons from the heme on cyt f (and possibly through Y1 as we proposed previously). We argue for the existence of this single dominant complex on the basis of observations that the most favorable orientations of the interaction between PC and cyt f, as determined by grouping successful BD trajectories on the basis of closest contacts of charged residues, tend to overlap one another and have very close distances between the metal centers on the two proteins (copper on PC, iron on cyt f). We use this knowledge to develop a model for PC/cyt f interaction that places a reaction between the two proteins occurring when the copper-to-iron distance is between 16 and 17 A. This reaction distance gives a good estimate of the experimentally observed rate constant for PC-cyt f interaction. Analysis of BD results as a function of ionic strength predicts an interaction that happens less frequently and becomes less specific as ionic strength increases.     

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.     

A Brownian dynamics computational study of the interaction of spinach plastocyanin with turnip cytochrome f: the importance of plastocyanin conformational changes

Brownian Dynamics (BD) computer simulations were used to study electrostatic interactions between turnip cytochrome f (cyt f) and spinach plastocyanin (PC). Three different spinach PC structures were studied: The X-ray crystal structure of Xue and coworkers [(1998) Protein Sci 7:2099–2105] and the NMR structure of Musiani et al. [(2005) J Biol Chem 280:18833–18841] and Ubbink and co-workers [(1998) Structure 6:323–335]. Significant differences exist in the backbone conformation between the PC taken from Ubbink and coworkers and the other two PC structures particularly the regions surrounding G10, E59–E60, and D51. Complexes formed in BD simulations using the PC of Ubbink and colleagues had a smaller Cu–Fe distance than the other two. These results suggest that different PC conformations may exist in solution with different capabilities of forming electron-transfer-active docks. All three types of complexes show electrostatic contacts between D42, E43, and D44 on PC and K187 on cyt f as well as between E59 on PC and K58 on cyt f. However, the PC of Ubbink and coworkers reveals additional contacts between D51 and cyt f as a result of the difference in backbone configuration. A second minor complex component was observed for the PC of Ubbink and co-workers and Xue and co-workers which had contacts between K187 on cyt f and E59 and E60 on PC rather than between K187 on cyt f and D42-D44 on PC as observed for the major components. This second type of complex may represent an earlier complex which rearranges to form a final complex capable of electron transfer. Professor Elizabeth L. Gross, Professor Emeritus of Biochemistry, The Ohio State University, passed away on June 27, 2007.     

A thermal unfolding study of plastocyanin from the thermophilic cyanobacterium Phormidium laminosum

The thermal unfolding of the plastocyanin from Phormidium laminosum, a thermophilic cyanobacterium, is herein described. The main objective of this work is to identify structural factors responsible for the higher stability observed in proteins from thermophilic organisms. With the aid of fluorescence spectroscopy, EPR, and NMR, the factors influencing the unfolding process of the protein were investigated, and procedures for its study have been standardized. The different spectroscopic techniques used provided consistent results showing that the thermal unfolding of plastocyanin is irreversible under all the conditions investigated and that this irreversibility does not appear to be related to the presence of oxygen. The oxidized plastocyanin species has proven to be more stable than the reduced one, with respect to both the required temperature for protein unfolding (up to a 9 degrees C difference between the two forms) and the kinetics of the process. The behavior of this plastocyanin contrasts with that of other cupredoxins whose unfolding had previously been studied. The unfolding pH dependence and kinetic studies indicate a process with a tight control around the physiological pH in which plastocyanin plays its redox role and the protein's isoelectric point (5.2), suggesting a close compromise between function and stability.     

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.     

Relation between interface properties and kinetics of electron transfer in the interaction of cytochrome f and plastocyanin from plants and the cyanobacterium Phormidium laminosum

Cytochrome f and plastocyanin from the cyanobacterium Phormidium laminosum react an order of magnitude faster than their counterparts from chloroplasts when long-range electrostatic interactions have been screened out by high salt concentration [Schlarb-Ridley, B. G., et al. (2002) Biochemistry 41, 3279-3285]. To investigate the relative contributions of the reaction partners to these differences, the reactions of turnip cytochrome f with P. laminosum plastocyanin and P. laminosum cytochrome f with pea plastocyanin were examined. Exchanging one of the plant reaction partners with the corresponding cyanobacterial protein nearly abolished electron transfer at low ionic strength but increased the rate at high ionic strength. This increase was larger for P. laminosum cytochrome f than for P. laminosumplastocyanin. To identify molecular features of P. laminosum cytochrome f that contribute to the increase, the effect of mutations in the N-terminal heme-shielding peptide on the reaction with P. laminosum plastocyanin was determined. Phenylalanine-3 was converted to valine and tryptophan-4 to phenylalanine or leucine. The mutations lowered the rate constant at 0.1 M ionic strength by factors of 0.71 for F4V, 0.42 for W4F, and 0.63 for W4L while introducing little change in the shape of the ionic strength dependence curve. When the N-terminal tetrapeptide (sequence YPFW) was converted into that found in the chloroplast of Chlamydomonas reinhardtii (YPVF), the reaction was slowed further (factor of 0.26). The N-terminal heme-shielding peptide was found to be responsible for 75% of the kinetic differences between cytochrome f from chloroplasts and the cyanobacterium when electrostatic interactions were eliminated.     

Implications of the effects of viscosity, macromolecular crowding, and temperature for the transient interaction between cytochrome f and plastocyanin from the cyanobacterium Phormidium laminosum

The reaction between cytochrome f and plastocyanin is a central feature of the photosynthetic electron-transport system of all oxygenic organisms. We have studied the reaction in solution to understand how the very weak binding between the two proteins from Phormidium laminosum can nevertheless lead to fast rates of electron transfer. In a previous publication [Schlarb-Ridley, B. G., et al. (2003) Biochemistry 42, 4057-4063], we suggested that the reaction is diffusion-controlled because of a strong effect of viscosity of the medium. The effects of viscosity and temperature have now been examined in detail. High molecular mass viscogens (Ficoll 70 and Dextran 70), which might mimic in vivo conditions, had little effect up to a relative viscosity of 4. Low molecular mass viscogens (ethane diol, glycerol, and sucrose) strongly decreased the bimolecular rate constant (k(2)) over a similar viscosity range. The effects correlated well with the viscosities of the solutions of the three reagents but not with their dielectric constants or molalities. A power law dependence of k(2) on viscosity suggested that k(2) depends on two viscosity-sensitive reactions in series, while the reverse reactions are little affected by viscosity. The results were incompatible with diffusion control of the overall reaction. Determination of the effect of temperature on k(2) gave an activation enthalpy, DeltaH(++) = 45 kJ mol(-)(1), which is also incompatible with diffusion control. The results were interpreted in terms of a model in which the stable form of the protein-protein complex requires further thermal activation to be competent for electron transfer.     

Structure of the soluble domain of cytochrome f from the cyanobacterium Phormidium laminosum.

Cytochrome f from the photosynthetic cytochrome b(6)f complex is unique among c-type cytochromes in its fold and heme ligation. The 1. 9-A crystal structure of the functional, extrinsic portion of cytochrome f from the thermophilic cyanobacterium Phormidium laminosum demonstrates that an unusual buried chain of five water molecules is remarkably conserved throughout the biological range of cytochrome f from cyanobacteria to plants [Martinez et al. (1994) Structure 2, 95-105]. Structure and sequence conservation of the cytochrome f extrinsic portion is concentrated at the heme, in the buried water chain, and in the vicinity of the transmembrane helix anchor. The electrostatic surface potential is variable, so that the surface of P. laminosum cytochrome f is much more acidic than that from turnip. Cytochrome f is unrelated to cytochrome c(1), its functional analogue in the mitochondrial respiratory cytochrome bc(1) complex, although other components of the b(6)f and bc(1) complexes are homologous. Identical function of the two complexes is inferred for events taking place at sites of strong sequence conservation. Conserved sites throughout the entire cytochrome b(6)f/bc(1) family include the cluster-binding domain of the Rieske protein and the heme b and quinone-binding sites on the electrochemically positive side of the membrane within the b cytochrome, but not the putative quinone-binding site on the electrochemically negative side.     

Role of electrostatics in the interaction between plastocyanin and photosystem I of the cyanobacterium Phormidium laminosum.

The interactions between photosystem I and five charge mutants of plastocyanin from the cyanobacterium Phormidium laminosum were investigated in vitro. The dependence of the overall rate constant of reaction, k2, on ionic strength was investigated using laser flash photolysis. The rate constant of the wild-type reaction increased with ionic strength, indicating repulsion between the reaction partners. Removing a negative charge on plastocyanin (D44A) accelerated the reaction and made it independent of ionic strength; removing a positive charge adjacent to D44 (K53A) had little effect. Neutralizing and inverting the charge on R93 slowed the reaction down and increased the repulsion. Specific effects of MgCl2 were observed for mutants K53A, R93Q and R93E. Thermodynamic analysis of the transition state revealed positive activation entropies, suggesting partial desolvation of the interface in the transition state. In comparison with plants, plastocyanin and photosystem I of Phormidium laminosum react slowly at low ionic strength, whereas the two systems have similar rates in the range of physiological salt concentrations. We conclude that in P. laminosum, in contrast with plants in vitro, hydrophobic interactions are more important than electrostatics for the reactions of plastocyanin, both with photosystem I (this paper) and with cytochrome f[Schlarb-Ridley, B.G., Bendall, D.S. & Howe, C.J. (2002) Biochemistry41, 3279-3285]. We discuss the implications of this conclusion for the divergent evolution of cyanobacterial and plant plastocyanins.     

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f.

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.     

Changes in non-core regions stabilise plastocyanin from the thermophilic cyanobacterium Phormidium laminosum

We report a theoretical investigation on the different stabilities of two plastocyanins. The first one belongs to the thermophilic cyanobacterium Phormidium laminosum and the second one belongs to its mesophilic relative Synechocystis sp. These proteins share the same topology and secondary-structure elements; however, the melting temperatures of their oxidised species differ by approximately 15 K. Long-time-scale molecular dynamics simulations, performed at different temperatures, show that the thermophilic protein optimises a set of intramolecular interactions (interstrand hydrogen bonding, salt bridging and hydrophobic clustering) within the region that comprises the strands β₅ and β₆, loop L₅ and the helix. This region exhibits most of the differences in the primary sequence between the two proteins and, in addition, it is involved in the interaction with known physiological partners. Further work is in progress to unveil the specific structural features responsible for the different thermal stability of the two proteins.     

Expression of plastocyanin and cytochrome f of the cyanobacterium Phormidium laminosum in Escherichia coli and Paracoccus denitrificans and the role of leader peptides.

The gene for plastocyanin from the cyanobacterium Phormidium laminosum was successfully expressed in Escherichia coli. Expression of the gene for cytochrome f resulted in the production of holocytochrome f in the periplasmic space of E. coli, but the yield was low. Expression in Paracoccus denitrificans yielded no holoprotein. When the region encoding the cytochrome f leader sequence was replaced with more typical bacterial leader sequences (those from the P. laminosum plastocyanin gene and the Paracoccus versutus cytochrome c-550 gene), much higher yields were consistently obtained in both species. Overexpressed proteins were compared to those isolated from P. laminosum and found to be identical in mass, isoelectric point, redox midpoint potential and (for plastocyanin) 1H-NMR spectrum.     

On the interaction between cytochrome f and plastocyanin

The interaction between cytochrome f and its electron acceptor plastocyanin (PC) was studied. To address the question of which specific regions and which of the positively charged residues of cytochrome f are important for the interaction with the negatively charged residues of PC we have used two different experimental approaches. Cytochrome f was proteolytically cleaved and fragments that could bind to a PC-affinity column were isolated. The smallest of these fragments was analysed to give information on the minimum structural requirement for binding to PC. By this procedure, we identified a peptide of approx. 11 kDa, containing the heme binding site, and having an N-terminal sequence identical to that of the mature cytochrome f. This finding suggests that the first 90 amino acids of cytochrome f contain at least some of the residues interacting with PC. The second approach involved modification of Arg residues of cytochrome f with the specific chemical modifier, hydroxyphenylglyoxal (HPG). Cytochrome f modification was performed in the absence of PC to enable identification of residues that are protected from modification when PC is bound to cytochrome f. Two peptides containing Arg residues which are modified in the absence of PC, but are not modified when PC is present, were isolated. Sequence analysis of these two peptides revealed that Arg residues no. 88 and 154 of cytochrome f are the residues that are protected from modification when cytochrome f is bound to PC, suggesting a role for these residues in the binding of cytochrome f to PC.     

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.     

A Brownian dynamics study of the interactions of the luminal domains of the cytochrome b6f complex with plastocyanin and cytochrome c6: the effects of the Rieske FeS protein on the interactions

The availability of the structures of the cytochrome b6f complex (cyt b6f), plastocyanin (PC), and cytochrome c6 (cyt c6) from Chlamydomonas reinhardtii allowed us, for the first time, to model electron transfer interactions between the luminal domains of this complex (including cyt f and the Rieske FeS protein) and its redox partners in the same species. We also generated a model structure in which the FeS center of the Rieske protein was positioned closer to the heme of cyt f than observed in the crystal structure and studied its interactions with both PC and cyt c6. Our data showed that the Rieske protein in both the original crystal structure and in our modeled structure of the cyt b6f complex did not physically interfere with binding position or orientation of PC or cyt c6 on cyt f. PC docked on cyt f with the same orientation in the presence or the absence of the Rieske protein, which matched well with the previously reported NMR structures of complexes between cyt f and PC. When the FeS center of the Rieske protein was moved close to the heme of cyt f, it even enhanced the interaction rates. Studies using a cyt f modified in the 184-191 loop showed that the cyt f structure is a more important factor in determining the rate of complex formations than is the presence or the absence of the Rieske protein or its position with respect to cyt f.