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.