The role of amino-acid residues in the hydrophobic patch surrounding the haem group of cytochrome f in the interaction with plastocyanin.

Soluble turnip cytochrome f has been purified from the periplasmic fraction of Escherichia coli expressing a truncated petA gene encoding the precursor protein lacking the C-terminal 33 amino-acid residues. The protein is identical [as judged by 1H-NMR spectroscopy, midpoint redox potential (+ 365 mV) and electron transfer reactions with plastocyanin] to cytochrome f purified from turnip leaves. Several residues in the hydrophobic patch surrounding the haem group have been changed by site-directed mutagenesis, and the proteins purified from E. coli. The Y1F and Q7N mutants showed only minor changes in the plastocyanin-binding constant Ka and the second-order rate constant for electron transfer to plastocyanin, whereas the Y160S mutant showed a 30% decrease in the overall rate of electron transfer caused in part by a 60% decrease in binding constant and partially compensated by an increased driving force due to a 27-mV decrease in redox potential. In contrast, the F4Y mutant showed increased rates of electron transfer which may be ascribed to an increased binding constant and a 14-mV decrease in midpoint redox potential. This indicates that subtle changes in the hydrophobic patch can influence rates of electron transfer to plastocyanin by changing the binding constants and altering the midpoint redox potential of the cytochrome haem group.     

Electron-transfer reactions of cytochrome f with flavin semiquinones and with plastocyanin. Importance of protein-protein electrostatic interactions and of donor-acceptor coupling.

Reduction of turnip ferricytochrome f by flavin semiquinones and oxidation of this ferrocytochrome f by French bean cupriplastocyanin are studied by laser flash photolysis over a wide range of ionic strengths. Second-order rate constants (+/- 15%) at extreme values of ionic strength, all at pH 7.0 and 22 degrees C, are as follows: with FMN semiquinone at 1.00 and 0.0040 M, 5.0 x 10(7) and 3.9 x 10(8) M-1 s-1; with riboflavin semiquinone at 1.00 and 0.0040 m, 1.7 x 10(8) and 1.9 x 10(8) M-1 s-1; with lumiflavin semiquinone at 1.00 and 0.0045 M, 1.8 x 10(8) and 4.5 x 10(8) M-1 s-1; with cupriplastocyanin at 1.00 and 0.100 M, 1.4 x 10(6) and 2.0 x 10(8) M-1 s-1. These reactions of cytochrome f are governed by the local positive charge of the interaction domain (the exposed heme edge), not by the overall negative charge of the protein. Lumiflavin semiquinone behaves as if it carried a small negative charge, probably because partial localization of the odd electron gives this electroneutral molecule some polarity; local charge seems to be more important than overall charge even for relatively small redox agents. The dependence of the rate constants on ionic strength was fitted to the equation of Watkins; this model recognizes the importance of local charges of the domains through which redox partners interact. There is kinetic evidence that a noncovalent complex between cytochrome f and plastocyanin exists at low ionic strength.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

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.     

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

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

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

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

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.     

The structural and functional role of lysine residues in the binding domain of cytochrome c in the electron transfer to cytochrome c oxidase.

The interactions of yeast iso-1 cytochrome c with bovine cytochrome c oxidase were studied using cytochrome c variants in which lysines of the binding domain were substituted by alanines. Resonance Raman spectra of the fully oxidized complexes of both proteins reveal structural changes of both the heme c and the hemes a and a3. The structural changes in cytochrome c are the same as those observed upon binding to phospholipid vesicles where the bound protein exists in two conformers, B1 and B2. Whereas the structure of B1 is the same as that of the unbound cytochrome c, the formation of B2 is associated with substantial alterations of the heme pocket. In cytochrome c oxidase, the structural changes in both hemes refer to more subtle perturbations of the immediate protein environment and may be a result of a conformational equilibrium involving two states. These changes are qualitatively different to those observed for cytochrome c oxidase upon poly-l-lysine binding. The resonance Raman spectra of the various cytochrome c/cytochrome c oxidase complexes were analyzed quantitatively. The spectroscopic studies were paralleled by steady-state kinetic measurements of the same protein combinations. The results of the spectra analysis and the kinetic studies were used to determine the stability of the complexes and the conformational equilibria B2/B1 for all cytochrome c variants. The complex stability decreases in the order: wild-type WT > J72K > K79A > K73A > K87A > J72A > K86A > K73A/K79A (where J is the natural trimethyl lysine). This order is not exhibited by the conformational equilibria. The electrostatic control of state B2 formation does not depend on individual intermolecular salt bridges, but on the charge distribution in a specific region of the front surface of cytochrome c that is defined by the lysyl residues at positions 72, 73 and 79. On the other hand, the conformational changes in cytochrome c oxidase were found to be independent of the identity of the bound cytochrome c variant. The maximum rate constants determined from steady-state kinetic measurements could be related to the conformational equilibria of the bound cytochrome c using a simple model that assumes that the conformational transitions are faster than product formation. Within this model, the data analysis leads to the conclusion that the interprotein electron transfer rate constant is around two times higher in state B2 than in B1. These results can be interpreted in terms of an increase of the driving force in state B2 as a result of the large negative shift of the reduction potential.     

The role of cytochrome P450 lysine residues in the interaction between cytochrome P450IA1 and NADPH-cytochrome P450 reductase.

Cytochrome P450IA1 (purified from hepatic microsomes of beta-naphthoflavone-treated rats) has been covalently modified with the lysine-modifying reagent acetic anhydride. Different levels of lysine residue modification in cytochrome P450IA1 can be achieved by varying the concentration of acetic anhydride. Modification of lysine residues in P450IA1 greatly inhibits the interaction of P450IA1 with NADPH-cytochrome P450 reductase. Modification of 1.0 and 3.3 mol lysine residues per mole P450IA1 resulted in 30 and 95% decreases, respectively, in 7-ethoxycoumarin hydroxylation by a reconstituted P450IA1/reductase complex. However, modification of 3.3 mol lysine residues per mole P450IA1 decreased only cumene hydroperoxide-supported P450-dependent 7-ethoxycoumarin hydroxylation by 30%. Spectral and fluorescence studies showed no indication of global conformational change of P450IA1 even with up to 8.8 mol lysine residues modified per mole P450IA1. These data suggest that at least three lysine residues in P450IA1 may be involved in the interaction with reductase. Identification of lysine residues in P450IA1 possibly involved in this interaction was carried out by [14C]acetic anhydride modification, trypsin digestion, HPLC separation, and amino acid sequencing. The lysine residue candidates identified in this manner were K97, K271, K279, and K407.     

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.     

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.     

Photoinduced electron-transfer reaction in a ternary system involving zinc cytochrome c and plastocyanin. Interplay of monopolar and dipolar electrostatic interactions between metalloproteins.

A carbodiimide promotes noninvasive cross-linking between amino groups surrounding the exposed heme edge in zinc cytochrome c and carboxylic groups in the acidic patch in plastocyanin. Eight derivatives of the covalent complex Zncyt/pc(I), which have similar structures but different overall charges because of different numbers and locations of N-acylurea groups, are separated by cation-exchange chromatography. Kinetics of electron transfer from the diprotein complex in the triplet excited state, 3Zncyt/pc(I), to free cupriplastocyanin at pH 7.0 and various ionic strengths is studied by laser flash spectroscopy. This reaction is purely bimolecular for all eight N-acylurea derivatives of the diprotein complex. The overall charges of the derivatives 1 and 2 at pH 7.0 are -2 and 0, respectively; both of them, however, have very large dipole moments of 410-480 D. The rate constants for their reactions with cupriplastocyanin, whose charge at pH 7.0 is -8 and whose dipole moment is 362 D, are determined over the range of ionic strengths from 2.5 mM to 3.00 M. The observed dependence of the rate constants on ionic strength cannot be explained in terms of net charges (monopole-monopole interactions) alone, but it can be fitted quantitatively with a theory that recognizes also monopole-dipole and dipole-dipole interactions [van Leeuwen, J. W. (1983) Biochim. Biophys. Acta 743, 408]. At ionic strengths up to ca. 10 mM monopole-monopole interactions predominate and Br?nsted-Debye-Hückel theory applies.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

Specificity of DNA basic polypeptide interactions. II+ Influence of aromatic amino acid residues investigated with agarose bound lysine copolypeptides.

Binding affinities towards DNA and base pair specificities of lysine copolymers, containing different amounts of Phe, Tyr or Trp residues, were estimated using a previously described chromatographic method. Incorporation of few aromatic residues into polylysine causes a decrease in the binding affinity, however, further raising the aromatic residue - lysine ratio results in a continous increase of affinity, which is most pronounced with the Tyr copolymers and not observed with polymers containing neutral aliphatic amino acid residues. AT-specificity increases concomitant with binding affinity in the case of the Tyr copolymers but not with the Phe copolymers. The interaction of DNA with the alternating Phe-Lys polymer is significantly stronger than with the random copolymer of equal residue composition. The molecular and conformational reasons determining specificity are discussed.     

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.     

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.     

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.