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.     

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.     

A Brownian dynamics study of the interaction of Phormidium laminosum plastocyanin with Phormidium laminosum cytochrome f

The interaction of Phormidium laminosum plastocyanin (PC) with P. laminosum cytochrome f (cyt f) was studied using Brownian dynamics (BD) simulations. Few complexes and a low rate of electron transfer were observed for wild-type PC. Increasing the positive electrostatic field on PC by the addition of a Zn(2+) ion in the neighborhood of D44 and D45 on PC (as found in crystal structure of plastocyanin) increased the number of complexes formed and the calculated rates of electron transfer as did PC mutations D44A, D45A, E54A, and E57A. Mutations of charged residues on Phormidium PC and Phormidium cyt f were used to map binding sites on both proteins. In both the presence and absence of the Zn(2+) ion, the following residues on PC interact with cyt f: D44, D45, K6, D79, R93, and K100 that lie in a patch just below H92 and Y88 and D10, E17, and E70 located on the upper portion of the PC molecule. In the absence of the Zn(2+) ion, K6 and K35 on the top of the PC molecule also interact with cyt f. Cyt f residues involved in binding PC, in the absence of the Zn(2+) ion, include E165, D187, and D188 that are located on the small domain of cyt f. The orientation of PC in the complexes was quite random in accordance with NMR results. In the presence of the Zn(2+) ion, K53 and E54 in the lower patch of the PC molecule also interact with cyt f and PC interacts with E86, E95, and E123 on the large domain of cyt f. Also, the orientation of PC in the complexes was much more uniform than in the absence of the Zn(2+) ion. The difference may be due to both the larger electrostatic field and the greater asymmetry of the charge distribution on PC observed in the presence of the Zn(2+) ion. Hydrophobic interactions were also observed suggesting a model of cyt f-PC interactions in which electrostatic forces bring the two molecules together but hydrophobic interactions participate in stabilizing the final electron-transfer-active dock.     

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.     

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.     

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 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.     

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.     

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 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.     

Effects of charged amino-acid mutation on the solution structure of cytochrome b5 and binding between cytochrome b5 and cytochrome c

The solution structure of oxidized bovine microsomal cytochrome b(5) mutant (E48, E56/A, D60/A) has been determined through 1524 meaningful nuclear Overhauser effect constraints together with 190 pseudocontact shift constraints. The final family of 35 conformers has rmsd values with respect to the mean structure of 0.045+/-0.009 nm and 0.088+/-0.011 nm for backbone and heavy atoms, respectively. A characteristic of this mutant is that of having no significant changes in the whole folding and secondary structure compared with the X-ray and solution structures of wild-type cytochrome b(5). The binding of different surface mutants of cytochrome b(5) with cytochrome c shows that electrostatic interactions play an important role in maintaining the stability and specificity of the protein complex formed. The differences in association constants demonstrate the electrostatic contributions of cytochrome b(5) surface negatively charged residues, which were suggested to be involved in complex formation in the Northrup and Salemme models, have cumulative effect on the stability of cyt c-cyt b(5) complex, and the contribution of Glu48 is a little higher than that of Glu44. Moreover, our result suggests that the docking geometry proposed by Northrup, which is involved in the participation of Glu48, Glu56, Asp60, and heme propionate of cytochrome b(5), do occur in the association between cytochrome b(5) and cytochrome c.     

Continuum electrostatic model for the binding of cytochrome c2 to the photosynthetic reaction center from Rhodobacter sphaeroides

Electrostatic interactions are important for protein-protein association. In this study, we examined the electrostatic interactions between two proteins, cytochrome c(2) (cyt c(2)) and the reaction center (RC) from the photosynthetic bacterium Rhodobacter sphaeroides, that function in intermolecular electron transfer in photosynthesis. Electrostatic contributions to the binding energy for the cyt c(2)-RC complex were calculated using continuum electrostatic methods based on the recent cocrystal structure [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. Calculated changes in binding energy due to mutations of charged interface residues agreed with experimental results for a protein dielectric constant epsilon(in) of 10. However, the electrostatic contribution to the binding energy for the complex was close to zero due to unfavorable desolvation energies that compensate for the favorable Coulomb attraction. The electrostatic energy calculated as a function of displacement of the cyt c(2) from the bound position showed a shallow minimum at a position near but displaced from the cocrystal configuration. These results show that although electrostatic steering is present, other short-range interactions must be present to contribute to the binding energy and to determine the structure of the complex. Calculations made to model the experimental data on association rates indicate a solvent-separated transition state for binding in which the cyt c(2) is displaced approximately 8 A above its position in the bound complex. These results are consistent with a two-step model for protein association: electrostatic docking of the cyt c(2) followed by desolvation to form short-range van der Waals contacts for rapid electron transfer.     

Structural transformation of cytochrome c and apo cytochrome c induced by sulfonated polystyrene

The structural transformation of cytochrome c (cyt c) and its heme-free precursor, apo cyt c, induced by negatively charged sulfonated polystyrene (SPS) with different charge density (degree of sulfonation) and chain length was studied to understand the factors that influence the folding and unfolding of the protein. SPS forms stable transparent nanoparticles in aqueous solution. The hydrophobic association of the backbone chain and phenyl groups is balanced by the electrostatic repulsion of the sulfonate groups on the particle surface. The binding of cyt c to negatively charged SPS particles causes an extensive disruption of the native compact structure of cyt c: the cleavage of Fe-Met80 ligand, about 40% loss of the helical structure, and the disruption of the asymmetry environment of Trp59. On the other hand, SPS particle-bound apo cyt c undergoes a conformational change from the random coil to alpha-helical structure. The folding of apo cyt c in SPS particles was influenced by pH and ionic strength of the solution, SPS concentration, and the degree of sulfonation and chain length of SPS. The folding can reach more than 90% of the alpha-helix content of native cyt c in solution. Poly(sodium 4-styrenesulfonate) (PSS), which is 100% sulfonated polystyrene and cannot form hydrophobic cores in the solution, induces only two-thirds of the alpha-helix content compared with SPS. It appears that the electrostatic interaction between PSS/SPS and apo cyt c induces an early partially folded state of apo cyt c. The hydrophobic interaction between nonpolar residues in apo cyt c and the hydrophobic cores in SPS particles extends the alpha-helical structure of apo cyt c.     

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.     

The role of individual lysine residues in the basic patch on turnip cytochrome f for electrostatic interactions with plastocyanin in vitro.

The role of electrostatic interactions in determining the rate of electron transfer between cytochrome f and plastocyanin has been examined in vitro with mutants of turnip cytochrome f and mutants of pea and spinach plastocyanins. Mutation of lysine residues Lys58, Lys65 and Lys187 of cytochrome f to neutral or acidic residues resulted in decreased binding constants and decreased rates of electron transfer to wild-type pea plastocyanin. Interaction of the cytochrome f mutant K187E with the pea plastocyanin mutant D51K gave a further decrease in electron transfer rate, indicating that a complementary charge pair at these positions could not compensate for the decreased overall charge on the proteins. Similar results were obtained with the interaction of the cytochrome f mutant K187E with single, double and triple mutants of residues in the acidic patches of spinach plastocyanin. These results suggest that the lysine residues of the basic patch on cytochrome f are predominantly involved in long-range electrostatic interactions with plastocyanin. However, analysis of the data using thermodynamic cycles provided evidence for the interaction of Lys187 of cytochrome f with Asp51, Asp42 and Glu43 of plastocyanin in the complex, in agreement with a structural model of a cytochrome f-plastocyanin complex determined by NMR.     

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.     

Interaction between cytochrome c2 and the photosynthetic reaction center from Rhodobacter sphaeroides: effects of charge-modifying mutations on binding and electron transfer.

The electrostatic interactions governing binding and electron transfer from cytochrome c(2) (cyt c(2)) to the reaction center (RC) from the photosynthetic bacteria Rhodobacter sphaeroides were studied by using site-directed mutagenesis to change the charges of residues on the RC surface. Charge-reversing mutations (acid --> Lys) decreased the binding affinity for cyt c(2). Dissociation constants, K(D) (0.3--250 microM), were largest for mutations of Asp M184 and nearby acid residues, identifying the main region for electrostatic interaction with cyt c(2). The second-order rate constants, k(2) (1--17 x 10(8) M(-1) s(-1)), increased with increasing binding affinity (log k(2) vs log 1/K(D) had a slope of approximately 0.4), indicating a transition state structurally related to the final complex. In contrast, first-order electron transfer rates, k(e), for the bound cyt did not change significantly (<3-fold), indicating that electron tunneling pathways were unchanged by mutation. Charge-neutralizing mutations (acid --> amide) showed changes in binding free energies of approximately 1/2 the free energy changes due to the corresponding charge-reversing mutations, suggesting that the charges in the docked complex remain well solvated. Charge-enhancing mutations (amide --> acid) produced free energy changes of the same magnitude (but opposite sign) as changes due to the charge-neutralizing mutations in the same region, indicating a diffuse electrostatic potential due to cyt c(2). A two-domain model is proposed, consisting of an electrostatic docking domain with charged surfaces separated by a water layer and a hydrophobic tunneling domain with atomic contacts that provide an efficient pathway for electron transfer.     

Identification of precise electrostatic recognition sites between cytochrome c6 and the photosystem I subunit PsaF using mass spectrometry

The reduction of the photo-oxidized special chlorophyll pair P700 of photosystem I (PSI) in the photosynthetic electron transport chain of eukaryotic organisms is facilitated by the soluble copper-containing protein plastocyanin (pc). In the absence of copper, pc is functionally replaced by the heme-containing protein cytochrome c6 (cyt c6) in the green alga Chlamydomonas reinhardtii. Binding and electron transfer between both donors and PSI follows a two-step mechanism that depends on electrostatic and hydrophobic recognition between the partners. Although the electrostatic and hydrophobic recognition sites on pc and PSI are well known, the precise electrostatic recognition site on cyt c6 is unknown. To specify the interaction sites on a molecular level, we cross-linked cyt c6 and PSI using a zero-length cross-linker and obtained a cross-linked complex competent in fast and efficient electron transfer. As shown previously, cyt c6 cross-links specifically with the PsaF subunit of PSI. Mass spectrometric analysis of tryptic peptides from the cross-linked product revealed specific interaction sites between residues Lys27 of PsaF and Glu69 of cyt c6 and between Lys23 of PsaF and Glu69/Glu70 of cyt c6. Using these new data, we present a molecular model of the intermolecular electron transfer complex between eukaryotic cyt c6 and PSI.     

Metalloprotein complexes for the study of electron-transfer reactions. Characterization of diprotein complexes obtained by covalent cross-linking of cytochrome c and plastocyanin with a carbodiimide.

Cytochrome c (cyt) and zinc cytochrome c (Zncyt) are separately cross-linked to plastocyanin (pc) by the carbodiimide EDC according to a published method. The changes in the protein reduction potentials indicate the presence of approximately two amide cross-links. Chromatography of the diprotein complexes cyt/pc and Zncyt/pc on CM-52 resin yields multiple fractions, whose numbers depend on the eluent. UV-vis, EPR, CD, MCD, resonance Raman, and surface-enhanced resonance Raman spectra show that cross-linking does not significantly perturb the heme and blue copper active sites. Degrees of heme exposure show that plastocyanin covers most of the accessible heme edge in cytochrome c. Impossibility of cross-linking cytochrome c to a plastocyanin derivative whose acidic patch had been blocked by chemical modification shows that it is the acidic patch that abuts the heme edge in the covalent complex. The chromatographic fractions of the covalent diprotein complex are structurally similar to one another and to the electrostatic diprotein complex. Isoelectric points show that the fractions differ from one another in the number and distribution of N-acylurea groups, byproducts of the reaction with the carbodiimide. Cytochrome c and plastocyanin are also tethered to each other via lysine residues by N-hydroxysuccinimide diesters. Tethers, unlike direct amide bonds, allow mobility of the cross-linked molecules. Laser-flash-photolysis experiments show that, nonetheless, the intracomplex electron-transfer reaction cyt(II)/pc(II)----cyt(III)/pc(I) is undetectable in complexes of either type. Only the electrostatic diprotein complex, in which protein rearrangement from the docking configuration to the reactive configuration is unrestricted, undergoes this intracomplex reaction at a measurable rate.     

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.