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.     

Electrostatic orientation of the electron-transfer complex between plastocyanin and cytochrome c.

To understand the specificity and efficiency of protein-protein interactions promoting electron transfer, we evaluated the role of electrostatic forces in precollision orientation by the development of two new methods, computer graphics alignment of protein electrostatic fields and a systematic orientational search of intermolecular electrostatic energies for two proteins at present separation distances. We applied these methods to the plastocyanin/cytochrome c interaction, which is faster than random collision, but too slow for study by molecular dynamics techniques. Significant electrostatic potentials were concentrated on one-fourth (969 A2) of the plastocyanin surface, with the greatest negative potential centered on the Tyr-83 hydroxyl within the acidic patch, and on one-eighth (632 A2) of the cytochrome c surface, with the greatest positive potential centered near the exposed heme edge. Coherent electrostatic fields occurred only over these regions, suggesting that local, rather than global, charge complementarity controls productive recognition. The three energetically favored families of pre-collision orientations all directed the positive region surrounding the heme edge of cytochrome c toward the acidic patch of plastocyanin but differed in heme plane orientation. Analysis of electrostatic fields, electrostatic energies of precollision orientations with 12 and 6 A separation distances, and surface topographies suggested that the favored orientations should converge to productive complexes promoting a single electron-transfer pathway from the cytochrome c heme edge to Tyr-83 of plastocyanin. Direct interactions of the exposed Cu ligand in plastocyanin with the cytochrome c heme edge are not unfavorable sterically or electrostatically but should occur no faster than randomly, indicating that this is not the primary pathway for electron transfer.     

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.     

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.     

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.     

Direct simulation of plastocyanin and cytochrome f interactions in solution

Most biological functions, including photosynthetic activity, are mediated by protein interactions. The proteins plastocyanin and cytochrome f are reaction partners in a photosynthetic electron transport chain. We designed a 3D computer simulation model of diffusion and interaction of spinach plastocyanin and turnip cytochrome f in solution. It is the first step in simulating the electron transfer from cytochrome f to photosystem 1 in the lumen of thylakoid. The model is multiparticle and it can describe the interaction of several hundreds of proteins. In our model the interacting proteins are represented as rigid bodies with spatial fixed charges. Translational and rotational motion of proteins is the result of the effect of stochastic Brownian force and electrostatic force. The Poisson-Boltzmann formalism is used to determine the electrostatic potential field generated around the proteins. Using this model we studied the kinetic characteristics of plastocyanin-cytochrome f complex formation for plastocyanin mutants at pH 7 and a variety of ionic strength values.     

The purification of cytochrome f and plastocyanin using affinity chromatography

Both plastocyanin and cytochrome f were purified using a combination of affinity chromatography together with established methods. Plastocyanin was partially purified using the method of Davis and San Pietro (Anal. Biochem. 95 (1979) 254-259), after which it was further purified using a column of cytochrome c covalently attached to Sepharose 4B. The affinity column was prepared using the method of Godinot and Gautheron (Methods Enzymol. 54 (1979) 112-114). The final purity index ratio (A278/A597) was less than 1.2, which is equal to that obtained using the more expensive FPLC procedure (Anderson, G.P., Sanderson, D.G., Lee, C.H., Durell, S., Anderson, L.B. and Gross, E.L. (1987) Biochim. Biophys. Acta 894, issue 3). Cytochrome f was partially purified using a modification of the method of Matazaki et al. (Plant Cell. Physiol. 16 (1975) 237-246) and bound to an affinity column of plastocyanin covalently attached to Sepharose 4B. Cytochrome f purified using this procedure had a purity index ratio (A554.5/A277) of 1.2. Both proteins are tyrosine proteins containing no tryptophan residues. After the affinity chromatography step, the fluorescence emission spectrum of either plastocyanin or cytochrome f was typical of a tyrosine protein free from tryptophan contamination.     

Kinetic studies on a cross-linked complex between plastocyanin cytochrome f

A cross-linked complex between plastocyanin and cytochrome f was prepared by incubation in the presence of a water soluble carbodiimide and its kinetic properties were studied. The optical spectra, oxidation-reduction potentials and isoelectric pH of plastocyanin and cytochrome f did not change upon the formation of the cross-linked complex. Studies on the ionic strength effect on the electron transfer rate from cross-linked plastocyanin to ferricyanide indicated that the negative charge on the reaction site of plastocyanin was masked upon the cross-linking. It was also suggested that the sign of the net charge near the cytochrome f heme edge changed from positive to negative upon the cross-linking. On the other hand, electrostatic interactions between cross-linked plastocyanin and P700 seemed to be essentially the same as those in the case of native plastocyanin, although the rate of electron transfer from cross-linked plastocyanin to P700 was severely reduced. We also measured the intra-complex electron transfer from cytochrome f to plastocyanin. This suggested that the covalently cross-linked complex is a valid model of the electron transfer encounter complex. Based on these results, the reaction sites of plastocyanin with P700 and cytochrome f were discussed.     

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.     

The transient complex of poplar plastocyanin with cytochrome f: effects of ionic strength and pH

The orientation of poplar plastocyanin in the complex with turnip cytochrome f has been determined by rigid-body calculations using restraints from paramagnetic NMR measurements. The results show that poplar plastocyanin interacts with cytochrome f with the hydrophobic patch of plastocyanin close to the heme region on cytochrome f and via electrostatic interactions between the charged patches on both proteins. Plastocyanin is tilted relative to the orientation reported for spinach plastocyanin, resulting in a longer distance between iron and copper (13.9 A). With increasing ionic strength, from 0.01 to 0.11 M, all observed chemical-shift changes decrease uniformly, supporting the idea that electrostatic forces contribute to complex formation. There is no indication for a rearrangement of the transient complex in this ionic strength range, contrary to what had been proposed earlier on the basis of kinetic data. By decreasing the pH from pH 7.7 to pH 5.5, the complex is destabilized. This may be attributed to the protonation of the conserved acidic patches or the copper ligand His87 in poplar plastocyanin, which are shown to have similar pK(a) values. The results are interpreted in a two-step model for complex formation.     

Preparation of a covalently linked adduct between plastocyanin and cytochrome f

Plastocyanin can be covalently cross-linked to the monomeric cytochrome f from turnip by incubation in the presence of a water-soluble carbodiimide. The adduct between the two proteins has a molecular weight of approximately 43,000 suggesting a 1:1 stoichiometry between the two proteins of the adduct. This stoichiometry has been verified by spectral characterization of the adduct. The efficiency of the cross-linking reaction is pH dependent with a higher degree of cross-linking being observed at pH 6.5 than at pH 7.0.     

Photosynthetic electron transfer at the level of cytochrome f and plastocyanin]

Light-induced redox conversions of cytochrome f and plastocyanin in situ and electron transporst from H2O to NADP+ were studied by a dual-wave differential spectrophotometry under identical conditions and subsequently compared. The analysis in red and far red light, treatment by inhibitors, e. g. diurone and dibromothymoquinone, and the analysis of photoreactions during the greening of etiolated seedlings demonstrated that cytochrome f functions only in the non-cyclic chain of electron transport, whereas plastocyanin--both in the non-cyclic and in the cyclic electron transport chains. The jositions of cytochrome f and plastocyanin in various electron-transport chains are proposed.     

Studies on photosystem I. I. Relationship of plastocyanin, cytochrome f and P700

DEAE-Sephadex column chromatography now has been used for the final step in purification of d-amino acid oxidase apoenzyme. A specific enzymatic activity of 35–37 units/mg has been obtained for the pure holoenzyme. The purity has been established by disc and SDS gel electrophoreses and by sedimentation equilibrium. The molecular weight per enzyme monomer has been found to be 38,000 ± 1000. Each enzyme monomer binds one FAD and one benzoate with dissociation constants at 23 °C and pH 8.5 of 5.35 × 10^?7m and 1.96 × 10^?6m, respectively. The holoenzyme is more negatively charged than the apoenzyme at alkaline pH. The amino acid composition and some other physicochemical properties of the oxidase are reported.     

Different modes of interaction in cyanobacterial complexes of plastocyanin and cytochrome f

The highly efficient electron-transfer chain in photosynthesis demonstrates a remarkable variation among organisms in the type of interactions between the soluble electron-transfer protein plastocyanin and it partner cytochrome f. The complex from the cyanobacterium Nostoc sp. PCC 7119 was studied using nuclear magnetic resonance spectroscopy and compared to that of the cyanobacterium Phormidium laminosum. In both systems, the main site of interaction on plastocyanin is the hydrophobic patch. However, the interaction in the Nostoc complex is highly dependent on electrostatics, contrary to that of Phormidium, resulting in a binding constant that is an order of magnitude larger at low ionic strength for the Nostoc complex. Studies of the mixed complexes show that these differences in interactions are mainly attributable to the surface properties of the plastocyanins.     

Structure of the intermolecular complex between plastocyanin and cytochrome f from spinach

In oxygenic photosynthesis, plastocyanin shuttles electrons between the membrane-bound complexes cytochrome b6f and photosystem I. The homologous complex between cytochrome f and plastocyanin, both from spinach, is the object of this study. The solution structure of the reduced spinach plastocyanin was determined using high field NMR spectroscopy, whereas the model structure of oxidized cytochrome f was obtained by homology modeling calculations and molecular dynamics. The model structure of the intermolecular complex was calculated using the program AUTODOCK, taking into account biological information obtained from mutagenesis experiments. The best electron transfer pathway from the heme group of cytochrome f to the copper ion of plastocyanin was calculated using the program HARLEM, obtaining a coupling decay value of 1.8 x 10(-4). Possible mechanisms of interaction and electron transfer between plastocyanin and cytochrome f were discussed considering the possible formation of a supercomplex that associates one cytochrome b6f, one photosystem I, and one plastocyanin.     

Electron transfer reactions of chemically modified plastocyanin with P700 and cytochrome f. Importance of local charges

Chemically modified spinach plastocyanin, in which negatively charged carboxyl residues are replaced with positively charged amino residues, has been prepared. Four distinct species of chemically modified plastocyanin, having 1 to 4 mol of modified carboxyl residue per mol of plastocyanin, could be separated by ion-exchange chromatography on DEAE-Sephacel. The rate of electron transfer from reduced cytochrome f to oxidized singly substituted plastocyanin was 30% of that of the native unmodified plastocyanin, and the reaction rate decreased further with increasing number of modified carboxyl residues. These results indicate the importance of electrostatic interactions between the negative charges on plastocyanin and the positive charges on cytochrome f in this reaction. Since the overall net charge of cytochrome f is negative at neutral pH, the positive charges on cytochrome f involved in the reaction should be localized ones. On the other hand, the rates of electron transfer from reduced singly and doubly substituted plastocyanin to photooxidized P700 in the P700-chlorophyll alpha protein complex were similar to that of native plastocyanin, which suggests that these carboxyl residues have only a minor role in the electron transfer to P700. Although divalent cation is essential for the electron transfer from native plastocyanin to P700 at neutral pH, the triply substituted plastocyanin could donate electrons to P700 even without MgCl2, and the rate of this reaction reached the maximum at a low concentration of MgCl2 (less than 2.5 mM). The modification of four carboxyl residues per plastocyanin molecule activated this reaction to the maximum level without MgCl2.(ABSTRACT TRUNCATED AT 250 WORDS)     

Electron transfers amongst cytochrome f, plastocyanin and photosystem I: kinetics and mechanisms

The review covers the theory and practice of the determination of kinetic constants for the electron transfer reactions in chloroplast thylakoid membranes between plastocyanin and cytochrome f in cytochrome bf complexes, and between plastocyanin and the reaction centre of photosystem I. Effects of ionic strength and pH are featured. The contribution of mutant studies is included. It is concluded that nearly all data from in vitro experiments can be interpreted with a reaction scheme in which an encounter complex between donor and acceptor is formed by long-range electrostatic attraction, followed by rearrangement during which metal centres become close enough for rapid intra-complex electron transfer. In vivo experiments so far cast doubt on this particular sequence, but their interpretation is not straightforward. Means of modelling the bimolecular complex between cytochrome f and plastocyanin are outlined, and two likely structures are illustrated. The complex formed by plastocyanin and photosystem I in higher plants involves the PsaF subunit, but its structure has not been fully determined.     

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.