Transient kinetics of electron transfer from a variety of c-type cytochromes to plastocyanin

New Insights into the Evolution of the Electron Transfer from Cytochrome f to Photosystem I in the Green and Red Branches of Photosynthetic Eukaryotes

In cyanobacteria and most green algae of the eukaryotic green lineage, the copper-protein plastocyanin (Pc) alternatively replaces the heme-protein cytochrome c6 (Cc6) as the soluble electron carrier from cytochrome f (Cf) to photosystem I (PSI). The functional and structural equivalence of 'green' Pc and Cc6 has been well established, representing an example of convergent evolution of two unrelated proteins. However, plants only produce Pc, despite having evolved from green algae. On the other hand, Cc6 is the only soluble donor available in most species of the red lineage of photosynthetic organisms, which includes, among others, red algae and diatoms. Interestingly, Pc genes have been identified in oceanic diatoms, probably acquired by horizontal gene transfer from green algae. However, the mechanisms that regulate the expression of a functional Pc in diatoms are still unclear. In the green eukaryotic lineage, the transfer of electrons from Cf to PSI has been characterized in depth. The conclusion is that in the green lineage, this process involves strong electrostatic interactions between partners, which ensure a high affinity and an efficient electron transfer (ET) at the cost of limiting the turnover of the process. In the red lineage, recent kinetic and structural modeling data suggest a different strategy, based on weaker electrostatic interactions between partners, with lower affinity and less efficient ET, but favoring instead the protein exchange and the turnover of the process. Finally, in diatoms the interaction of the acquired green-type Pc with both Cf and PSI may not yet be optimized.

Enhancing the population of the encounter complex affects protein complex formation efficiency

Optimal charge distribution is considered to be important for efficient formation of protein complexes. Electrostatic interactions guide encounter complex formation that precedes the formation of an active protein complex. However, disturbing the optimized distribution by introduction of extra charged patches on cytochrome c peroxidase does not lead to a reduction in productive encounters with its partner cytochrome c. To test whether a complex with a high population of encounter complex is more easily affected by suboptimal charge distribution, the interactions of cytochrome c mutant R13A with wild‐type cytochrome c peroxidase and a variant with an additional negative patch were studied. The complex of the peroxidase and cytochrome c R13A was reported to have an encounter state population of 80%, compared to 30% for the wild‐type cytochrome c. NMR analysis confirms the dynamic nature of the interaction and demonstrates that the mutant cytochrome c samples the introduced negative patch. Kinetic experiments show that productive complex formation is fivefold to sevenfold slower at moderate and high ionic strength values for cytochrome c R13A but the association rate is not affected by the additional negative patch on cytochrome c peroxidase, showing that the total charge on the protein surface can compensate for less optimal charge distribution. At low ionic strength (44 mm), the association with the mutant cytochrome c reaches the same high rates as found for wild‐type cytochrome c, approaching the diffusion limit.

The Charge Distribution on a Protein Surface Determines Whether Productive or Futile Encounter Complexes Are Formed

Protein complex formation depends strongly on electrostatic interactions. The distribution of charges on the surface of redox proteins is often optimized by evolution to guide recognition and binding. To test the degree to which the electrostatic interactions between cytochrome c peroxidase (CcP) and cytochrome c (Cc) are optimized, we produced five CcP variants, each with a different charge distribution on the surface. Monte Carlo simulations show that the addition of negative charges attracts Cc to the new patches, and the neutralization of the charges in the regular, stereospecific binding site for Cc abolishes the electrostatic interactions in that region entirely. For CcP variants with the charges in the regular binding site intact, additional negative patches slightly enhance productive complex formation, despite disrupting the optimized charge distribution. Removal of the charges in the regular binding site results in a dramatic decrease in the complex formation rate, even in the presence of highly negative patches elsewhere on the surface. We conclude that additional charge patches can result in either productive or futile encounter complexes, depending on whether negative residues are located also in the regular binding site.

Photosynthetic Electron Transfer by Dint of Protein Mobile Carriers. Multi-particle Brownian and Molecular Modeling

The paper presents the review of works on modeling the interaction of photosynthetic proteins using the multiparticle Brownian dynamics method developed at the Department of Biophysics, Biological Faculty, Lomonosov Moscow State University. The method describes the displacement of individual macromolecules – mobile electron carriers, and their electrostatic interactions between each other and with pigment-protein complexes embedded in photosynthetic membrane. Three-dimensional models of the protein molecules were constructed on the basis of the data from the Protein Data Bank. We applied the Brownian methods coupled to molecular dynamic simulations to reveal the role of electrostatic interactions and conformational motions in the transfer of an electron from the cytochrome complex Cyt b6f) membrane we developed the model which combines events of proteins Pc diffusion along the thylakoid membrane, electrostatic interactions of Pc with the membrane charges, formation of Pc super-complexes with multienzyme complexes of Photosystem I and to the molecule of the mobile carrier plastocyanin (Pc) in plants, green algae and cyanic bacteria. Taking into account the interior of photosynthetic membrane we developed the model which combines events of proteins Pc diffusion along the thylakoid membrane, electrostatic interactions of Pc with the membrane charges, formation of Pc super-complexes with multienzyme complexes of Photosystem I and Cyt b6f, embedded in photosynthetic membrane, electron transfer and complex dissociation. Multiparticle Brownian simulation method can be used to consider the processes of protein interactions in subcellular systems in order to clarify the role of individual stages and the biophysical mechanisms of these processes.

Non-covalent forces tune the electron transfer complex between ferredoxin and sulfite reductase to optimize enzymatic activity

Although electrostatic interactions between negatively charged ferredoxin (Fd) and positively charged sulfite reductase (SiR) have been predominantly highlighted to characterize complex formation, the detailed nature of intermolecular forces remains to be fully elucidated. We investigated interprotein forces for the formation of an electron transfer complex between Fd and SiR and their relationship to SiR activity using various approaches over NaCl concentrations between 0 and 400 mM. Fd-dependent SiR activity assays revealed a bell-shaped activity curve with a maximum ∼40-70 mM NaCl and a reverse bell-shaped dependence of interprotein affinity. Meanwhile, intrinsic SiR activity, as measured in a methyl viologen-dependent assay, exhibited saturation above 100 mM NaCl. Thus, two assays suggested that interprotein interaction is crucial in controlling Fd-dependent SiR activity. Calorimetric analyses showed the monotonic decrease in interprotein affinity on increasing NaCl concentrations, distinguished from a reverse bell-shaped interprotein affinity observed from Fd-dependent SiR activity assay. Furthermore, Fd:SiR complex formation and interprotein affinity were thermodynamically adjusted by both enthalpy and entropy through electrostatic and non-electrostatic interactions. A residue-based NMR investigation on the addition of SiR to ¹⁵N-labeled Fd at the various NaCl concentrations also demonstrated that a combination of electrostatic and non-electrostatic forces stabilized the complex with similar interfaces and modulated the binding affinity and mode. Our findings elucidate that non-electrostatic forces are also essential for the formation and modulation of the Fd:SiR complex. We suggest that a complex configuration optimized for maximum enzymatic activity near physiological salt conditions is achieved by structural rearrangement through controlled non-covalent interprotein interactions.

Fusion of Ferredoxin and Cytochrome P450 Enables Direct Light-Driven Biosynthesis

Cytochrome P450s (P450s) are key enzymes in the synthesis of bioactive natural products in plants. Efforts to harness these enzymes for in vitro and whole-cell production of natural products have been hampered by difficulties in expressing them heterologously in their active form, and their requirement for NADPH as a source of reducing power. We recently demonstrated targeting and insertion of plant P450s into the photosynthetic membrane and photosynthesis-driven, NADPH-independent P450 catalytic activity mediated by the electron carrier protein ferredoxin. Here, we report the fusion of ferredoxin with P450 CYP79A1 from the model plant Sorghum bicolor, which catalyzes the initial step in the pathway leading to biosynthesis of the cyanogenic glucoside dhurrin. Fusion with ferredoxin allows CYP79A1 to obtain electrons for catalysis by interacting directly with photosystem I. Furthermore, electrons captured by the fused ferredoxin moiety are directed more effectively toward P450 catalytic activity, making the fusion better able to compete with endogenous electron sinks coupled to metabolic pathways. The P450-ferredoxin fusion enzyme obtains reducing power solely from its fused ferredoxin and outperforms unfused CYP79A1 in vivo. This demonstrates greatly enhanced electron transfer from photosystem I to CYP79A1 as a consequence of the fusion. The fusion strategy reported here therefore forms the basis for enhanced partitioning of photosynthetic reducing power toward P450-dependent biosynthesis of important natural products.

Nanodomains of cytochrome b6f and photosystem II complexes in spinach grana thylakoid membranes

The cytochrome b6f (cytb6f) complex plays a central role in photosynthesis, coupling electron transport between photosystem II (PSII) and photosystem I to the generation of a transmembrane proton gradient used for the biosynthesis of ATP. Photosynthesis relies on rapid shuttling of electrons by plastoquinone (PQ) molecules between PSII and cytb6f complexes in the lipid phase of the thylakoid membrane. Thus, the relative membrane location of these complexes is crucial, yet remains unknown. Here, we exploit the selective binding of the electron transfer protein plastocyanin (Pc) to the lumenal membrane surface of the cytb6f complex using a Pc-functionalized atomic force microscope (AFM) probe to identify the position of cytb6f complexes in grana thylakoid membranes from spinach (Spinacia oleracea). This affinity-mapping AFM method directly correlates membrane surface topography with Pc-cytb6f interactions, allowing us to construct a map of the grana thylakoid membrane that reveals nanodomains of colocalized PSII and cytb6f complexes. We suggest that the close proximity between PSII and cytb6f complexes integrates solar energy conversion and electron transfer by fostering short-range diffusion of PQ in the protein-crowded thylakoid membrane, thereby optimizing photosynthetic efficiency.

External loops at the ferredoxin-NADP(+) reductase protein-partner binding cavity contribute to substrates allocation

Ferredoxin-NADP(+) reductase (FNR) is the structural prototype of a family of FAD-containing reductases that catalyze electron transfer between low potential proteins and NAD(P)(+)/H, and that display a two-domain arrangement with an open cavity at their interface. The inner part of this cavity accommodates the reacting atoms during catalysis. Loops at its edge are highly conserved among plastidic FNRs, suggesting that they might contribute to both flavin stabilization and competent disposition of substrates. Here we pay attention to two of these loops in Anabaena FNR. The first is a sheet-loop-sheet motif, loop102-114, that allocates the FAD adenosine. It was thought to determine the extended FAD conformation, and, indirectly, to modulate isoalloxazine electronic properties, partners binding, catalytic efficiency and even coenzyme specificity. The second, loop261-269, contains key residues for the allocation of partners and coenzyme, including two glutamates, Glu267 and Glu268, proposed as candidates to facilitate the key displacement of the C-terminal tyrosine (Tyr303) from its stacking against the isoalloxazine ring during the catalytic cycle. Our data indicate that the main function of loop102-114 is to provide the inter-domain cavity with flexibility to accommodate protein partners and to guide the coenzyme to the catalytic site, while the extended conformation of FAD must be induced by other protein determinants. Glu267 and Glu268 appear to assist the conformational changes that occur in the loop261-269 during productive coenzyme binding, but their contribution to Tyr303 displacement is minor than expected. Additionally, loop261-269 appears a determinant to ensure reversibility in photosynthetic FNRs.

Electron transfer pathways and dynamics of chloroplast NADPH-dependent thioredoxin reductase C (NTRC)

NADPH-dependent thioredoxin reductases (NTRs) contain a flavin cofactor and a disulfide as redox-active groups. The catalytic mechanism of standard NTR involves a large conformational change between two configurations. Oxygenic photosynthetic organisms possess a plastid-localized NTR, called NTRC, with a thioredoxin module fused at the C terminus. NTRC is an efficient reductant of 2-Cys peroxiredoxins (2-Cys Prxs) and thus is involved in the protection against oxidative stress, among other functions. Although the mechanism of electron transfer of canonical NTRs is well established, it is not yet known in NTRC. By employing stopped-flow spectroscopy, we have carried out a comparative kinetic study of the electron transfer reactions involving NTRC, the truncated NTR module of NTRC, and NTRB, a canonical plant NTR. Whereas the three NTRs maintain the conformational change associated with the reductive cycle of catalysis, NTRC intramolecular electron transfer to the thioredoxin module presents two kinetic components (k(ET) of ~2 and 0.1 s(-1)), indicating the occurrence of additional dynamic motions. Moreover, the dynamic features associated with the electron transfer to the thioredoxin module are altered in the presence of 2-Cys Prx. NTRC shows structural constraints that may locate the thioredoxin module in positions with different efficiencies for electron transfer, the presence of 2-Cys Prx shifting the conformational equilibrium of the thioredoxin module to a specific position, which is not the most efficient.

ArsH from the cyanobacterium Synechocystis sp. PCC 6803 is an efficient NADPH-dependent quinone reductase

The cyanobacterium Synechocystis sp. PCC 6803 possesses an arsenic resistance operon that encodes, among others, an ArsH protein. ArsH is a flavin mononucleotide (FMN)-containing protein of unknown function and a member of the family of NADPH-dependent FMN reductases. The nature of its final electron acceptor and the role of ArsH in the resistance to arsenic remained to be clarified. Here we have expressed and purified Synechocystis ArsH and conducted an intensive biochemical study. We present kinetic evidence supporting a quinone reductase activity for ArsH, with a preference for quinones with hydrophobic substituents. By using steady-state activity measurements, as well as stopped-flow and laser-flash photolysis kinetic analyses, it has been possible to establish the mechanism of the process and estimate the values of the kinetic constants. Although the enzyme is able to stabilize the anionic semiquinone form of the FMN, reduction of quinones involves the hydroquinone form of the flavin cofactor, and the enzymatic reaction occurs through a ping-pong-type mechanism. ArsH is able to catalyze one-electron reactions (oxygen and cytocrome c reduction), involving the FMN semiquinone form, but with lower efficiency. In addition, arsH mutants are sensitive to the oxidizing agent menadione, suggesting that ArsH plays a role in the response to oxidative stress caused by arsenite.

Dual role of FMN in flavodoxin function: electron transfer cofactor and modulation of the protein-protein interaction surface

Flavodoxin (Fld) replaces Ferredoxin (Fd) as electron carrier from Photosystem I (PSI) to Ferredoxin-NADP(+) reductase (FNR). A number of Anabaena Fld (AnFld) variants with replacements at the interaction surface with FNR and PSI indicated that neither polar nor hydrophobic residues resulted critical for the interactions, particularly with FNR. This suggests that the solvent exposed benzenoid surface of the Fld FMN cofactor might contribute to it. FMN has been replaced with analogues in which its 7- and/or 8-methyl groups have been replaced by chlorine and/or hydrogen. The oxidised Fld variants accept electrons from reduced FNR more efficiently than Fld, as expected from their less negative midpoint potential. However, processes with PSI (including reduction of Fld semiquinone by PSI, described here for the first time) are impeded at the steps that involve complex re-arrangement and electron transfer (ET). The groups introduced, particularly chlorine, have an electron withdrawal effect on the pyrazine and pyrimidine rings of FMN. These changes are reflected in the magnitude and orientation of the molecular dipole moment of the variants, both factors appearing critical for the re-arrangement of the finely tuned PSI:Fld complex. Processes with FNR are also slightly modulated. Despite the displacements observed, the negative end of the dipole moment points towards the surface that contains the FMN, still allowing formation of complexes competent for efficient ET. This agrees with several alternative binding modes in the FNR:Fld interaction. In conclusion, the FMN in Fld not only contributes to the redox process, but also to attain the competent interaction of Fld with FNR and PSI.

Characterisation of the electron transfer and complex formation between flavodoxin from D. vulgaris and the haem domain of cytochrome P450 BM3 from B. megaterium

Investigation of the complex formation and electron transfer kinetics between P450 BMP and flavodoxin was carried out following the suggested involvement of flavodoxin in modulating the electron transfer to BMP in artificial redox chains bound to an electrode surface. While electron transfer measurements show the formation of a tightly bound complex, the NMR data indicate the formation of shortly lived complexes. The measured k(obs) ranged from 24.2 s(-1) to 44.1 s(-1) with k(on) ranging from 0.07 x 10(6) to 1.1 x 10(6) s(-1) M(-1) and K(d) ranging from 300 microM to 24 microM in buffers of different ionic strength. This apparent contradiction is due to the existence of two events in the complex formation prior to electron transfer. A stable complex is initially formed. Within such tightly bound complex, flavodoxin rocks rapidly between different positions. The rocking of the bound flavodoxin between several different orientations gives rise to the transient complexes in fast exchange as observed in the NMR experiments. Docking simulations with two different approaches support the theory that there is no highly specific orientation in the complex, but instead one side of the flavodoxin binds the P450 with high overall affinity but with a number of different orientations. The level of functionality of each orientation is dependent on the distance between cofactors, which can vary between 8 and 25 A, with some of the transient complexes showing distances compatible with the measured electron transfer rate constants.

Demonstration of short-lived complexes of cytochrome c with cytochrome bc1 by EPR spectroscopy: implications for the mechanism of interprotein electron transfer

One of the steps of a common pathway for biological energy conversion involves electron transfer between cytochrome c and cytochrome bc1. To clarify the mechanism of this reaction, we examined the structural association of those two proteins using the electron transfer-independent electron paramagnetic resonance (EPR) techniques. Drawing on the differences in the continuous wave EPR spectra and saturation recoveries of spin-labeled bacterial and mitochondrial cytochromes c recorded in the absence and presence of bacterial cytochrome bc1, we have exposed a time scale of dynamic equilibrium between the bound and the free state of cytochrome c at various ionic strengths. Our data show a successive decrease of the bound cytochrome c fraction as the ionic strength increases, with a limit of approximately 120 mm NaCl above which essentially no bound cytochrome c can be detected by EPR. This limit does not apply to all of the interactions of cytochrome c with cytochrome bc1 because the cytochrome bc1 enzymatic activity remained high over a much wider range of ionic strengths. We concluded that EPR monitors just the tightly bound state of the association and that an averaged lifetime of this state decreases from over 100 micros at low ionic strength to less than 400 ns at an ionic strength above 120 mm. This suggests that at physiological ionic strength, the tightly bound complex on average lasts less than the time needed for a single electron exchange between hemes c and c1, indicating that productive electron transfer requires several collisions of the two molecules. This is consistent with an early idea of diffusion-coupled reactions that link the soluble electron carriers with the membranous complexes, which, we believe, provides a robust means of regulating electron flow through these complexes.

Reduction of plastocyanin by tyrosine-containing oligopeptides

Oxidized plastocyanin (PC) was reduced with TyrTyrTyr and LysLysLysLysTyrTyrTyr (KKKKYYY) oligopeptides at neutral pH. The TyrTyrTyr site of the peptides provided an electron to the copper active site of PC, whereas the tetralysine site of KKKKYYY functioned as the recognition site for the negative patch of PC. The reciprocal initial rate constant (1/k(int)) increased linearly with the reciprocal TyrTyrTyr concentration and proton concentration, although the electron transfer rate decreased gradually with time. The results showed that PC was reduced by the deprotonated species of TyrTyrTyr. A linear increase of log k(int) with increase in the ionic strength was observed due to decrease in the electrostatic repulsion between negatively charged PC and deprotonated (TyrTyrTyr)(-). PC was reduced faster by an addition of KKKKYYY to the PC-TyrTyrTyr solution, although KKKKYYY could not reduce PC without TyrTyrTyr. The ESI-LCMS spectrum of the products from the reaction between PC and TyrTyrTyr showed molecular ion peaks at m/z 1015.7 and 1037.7, which suggested formation of a dimerized peptide that may be produced from the reaction of a tyrosyl radical. The results indicate that PC and the tyrosine-containing oligopeptides form an equilibrium, PC(ox)/(oligopeptide)(-)-->/<--PC(red)/(oligopeptide)(*). The equilibrium is usually shifted to the left, but could shift to the right when the produced oligopeptide radical reacts with unreacted peptides. For the reaction of PC with KKKKYYY in the absence of TyrTyrTyr, the produced KKKK(YYY)(*) radical peptide could not react with other KKKKYYY peptides, since they were positively charged. In the presence of both KKKKYYY and TyrTyrTyr, PC may interact effectively with KKKKYYY through its tetralysine site and receive an electron from its TyrTyrTyr site, where the produced KKKK(YYY)(*) may interact with TyrTyrTyr peptides.

A comparative kinetic analysis of the reactivity of plant, horse, and human respiratory cytochrome c towards cytochrome c oxidase

Two synthetic genes coding for human and Arabidopsis cytochrome c, respectively, have been designed and constructed, and the recombinant proteins have been over-expressed in Escherichia coli cells. Thus a comparative analysis of the two heme proteins, including horse cytochrome c as a reference, has been performed. In addition to their physico-chemical properties, the redox behavior of the three proteins has been analyzed by following the kinetics of both their reduction by flavin semiquinones (lumiflavin, riboflavin, and FMN) and oxidation by cytochrome c oxidase. The resulting data indicate that the accessibility and electrostatic charge of the active site do not differ in a significant way among the three proteins, but human cytochrome c exhibits some intriguing differences when interacting with cytochrome c oxidase that could be related to the amino acid changes underwent by the latter along evolution.

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.

Respiratory cytochromecoxidase can be efficiently reduced by the photosynthetic redox proteins cytochromec6and plastocyanin in cyanobacteria

Plastocyanin and cytochrome c6 are two small soluble electron carriers located in the intrathylacoidal space of cyanobacteria. Although their role as electron shuttle between the cytochrome b6f and photosystem I complexes in the photosynthetic pathway is well established, their participation in the respiratory electron transport chain as donors to the terminal oxidase is still under debate. Here, we present the first time-resolved analysis showing that both cytochrome c6 and plastocyanin can be efficiently oxidized by the aa3 type cytochrome c oxidase in Nostoc sp. PCC 7119. The apparent electron transfer rate constants are ca. 250 and 300 s(-1) for cytochrome c6 and plastocyanin, respectively. These constants are 10 times higher than those obtained for the oxidation of horse cytochrome c by the oxidase, in spite of being a reaction thermodynamically more favourable.

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.

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.

Functional characterization of the evolutionarily divergent fern plastocyanin

Plastocyanin (Pc) is a soluble copper protein that transfers electrons from cytochrome b(6)f to photosystem I (PSI), two protein complexes that are localized in the thylakoid membranes in chloroplasts. The surface electrostatic potential distribution of Pc plays a key role in complex formation with the membrane-bound partners. It is practically identical for Pcs from plants and green algae, but is quite different for Pc from ferns. Here we report on a laser flash kinetic analysis of PSI reduction by Pc from various eukaryotic and prokaryotic organisms. The reaction of fern Pc with fern PSI fits a two-step kinetic model, consisting of complex formation and electron transfer, whereas other plant systems exhibit a mechanism that requires an additional intracomplex rearrangement step. The fern Pc interacts inefficiently with spinach PSI, showing no detectable complex formation. This can be explained by assuming that the unusual surface charge distribution of fern Pc impairs the interaction. Fern PSI behaves in a similar way as spinach PSI in reaction with other Pcs. The reactivity of fern Pc towards several soluble c-type cytochromes, including cytochrome f, has been analysed by flavin-photosensitized laser flash photolysis, demonstrating that the specific surface motifs for the interaction with cytochrome f are conserved in fern Pc.

Interaction of plastocyanin with oligopeptides: effect of lysine distribution within the peptide

We synthesized and purified four oligopeptides containing four lysines (KKKK, GKKGGKK, KKGGGKK, and KGKGKGK) as models for the plastocyanin (PC) interacting site of cytochrome f. These peptides competitively inhibited electron transfer between cytochrome c and PC. The inhibitory effect increased as the peptide concentrations were increased. The association constants between PC and the peptides did not differ significantly (3500-5100 M(-1)), although the association constant of PC-KGKGKGK was a little larger than the constants between PC and other peptides. Changes in the absorption spectrum of PC were observed when the peptides were added to the PC solution: peaks and troughs were detected at about 460 and 630 nm and at about 560 and 700 nm, respectively, in the difference absorption spectra between the spectra with and without peptides. These changes were attributed to the structural change at the copper site of PC by interaction with the peptides. The structural change was most significant when tetralysine was used. These results show that binding of the oligopeptide to PC is slightly more efficient when lysines are distributed uniformly within the peptide, whereas the structural change of PC becomes larger when the lysines are close to each other within the peptide.

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.

Conformational reorganisation in interfacial protein electron transfer

Protein-protein electron transfer (ET) plays an essential role in all redox chains. Earlier studies which used cross-linking and increased solution viscosity indicated that the rate of many ET reactions is limited (i.e., gated) by conformational reorientations at the surface interface. These results are later supported by structural studies using NMR and molecular modelling. New insights into conformational gating have also come from electrochemical experiments in which proteins are noncovalently adsorbed on the electrode surface. These systems have the advantage that it is relatively easy to vary systematically the driving force and electronic coupling. In this review we summarize the current knowledge obtained from these electrochemical experiments and compare it with some of the results obtained for protein-protein ET.

Docking and electron transfer studies between rubredoxin and rubredoxin:oxygen oxidoreductase

The interaction and electron transfer (ET) between rubredoxin (Rd) and rubredoxin:oxygen oxidoreductase (ROO) from Desulfovibrio gigas is studied by molecular modelling techniques. Experimental kinetic assays using recombinant proteins show that the Rd reoxidation by ROO displays a bell-shaped dependence on ionic strength, suggesting a non-trivial electrostatic dependence of the interaction between these two proteins. Rigid docking studies reveal a prevalence for Rd to interact, in a very specific way, with the surface of the ROO dimer near its FMN cofactors. The optimization of the lowest energy complexes, using molecular dynamics simulation, shows a very tight interaction between the surface of the two proteins, with a high probability for Rd residues (but not the iron centre directly) to be in direct contact with the FMN cofactors of ROO. Both electrostatics and van der Waals interactions contribute to the final energy of the complex. In these complexes, the major contributions for complex formation are polar interactions between acidic residues of Rd and basic residues of ROO, plus substantial non-polar interactions between different groups. Important residues for this process are identified. ET estimates (using the Pathways model), in the optimized lowest energy complexes, suggest that these configurations are efficient for transferring electrons. The experimental bell-shaped dependence of kinetics on ionic strength is analysed in view of the molecular modelling results, and hypotheses for the molecular basis of this phenomenon are discussed.

Mutagenesis of prochlorothrix plastocyanin reveals additional features in photosystem I interactions

Three surface residues of plastocyanin from Prochlorothrix hollandica have been modified by site-directed mutagenesis. Changes have been made in methionine 33, located in the hydrophobic patch of the copper protein, and in arginine 86 and proline 53, both located in the eastern hydrophilic area. The reactivity toward photosystem I of single mutants M33N, P53A, P53E, R86Q, R86E, and the double mutant M33N/P14L has been studied by laser flash absorption spectroscopy. All the mutations yield increased reactivity of plastocyanin toward photosystem I as compared with wild type plastocyanin, thus indicating that in Prochlorothrix electron donation to photosystem I is not optimized. The most drastic increases in the intracomplex electron transfer rate are obtained with mutants in methionine 33, whereas replacing arginine 86 only modestly affects the plastocyanin-photosystem I equilibrium constant for complex formation. Mutations at position 53 also promote major changes in the association of plastocyanin with photosystem I, yielding a change from a mechanism involving complex formation to a simpler collisional interaction. Molecular dynamics calculations indicate that mutations at position 33 promote changes in the H-bond network around the copper center. The comparative kinetic analysis of the reactivity of Prochlorothrix plastocyanin mutants toward photosystem I from other cyanobacteria reveals that mutations M33N, P53A, and P53E result in enhanced general reactivity.

Electron transfer from the Rieske iron-sulfur protein (ISP) to cytochrome f in vitro. Is a guided trajectory of the ISP necessary for competent docking?

The time course of electron transfer in vitro between soluble domains of the Rieske iron-sulfur protein (ISP) and cytochrome f subunits of the cytochrome b(6)f complex of oxygenic photosynthesis was measured by stopped-flow mixing. The domains were derived from Chlamydomonas reinhardtii and expressed in Escherichia coli. The expressed 142-residue soluble ISP apoprotein was reconstituted with the [2Fe-2S] cluster. The second-order rate constant, k(2)((ISP-f)) = 1.5 x 10(6) m(-1) s(-1), for ISP to cytochrome f electron transfer was <10(-2) of the rate constant at low ionic strength, k(2)((f-PC))(> 200 x 10(6) m(-1) s(-1)), for the reduction of plastocyanin by cytochrome f, and approximately 1/30 of k(2)((f-PC)) at the ionic strength estimated for the thylakoid interior. In contrast to k(2)((f-PC)), k(2)((ISP-f)) was independent of pH and ionic strength, implying no significant role of electrostatic interactions. Effective pK values of 6.2 and 8.3, respectively, of oxidized and reduced ISP were derived from the pH dependence of the amplitude of cytochrome f reduction. The first-order rate constant, k(1)((ISP-f)), predicted from k(2)((ISP-f)) is approximately 10 and approximately 150 times smaller than the millisecond and microsecond phases of cytochrome f reduction observed in vivo. It is proposed that in the absence of electrostatic guidance, a productive docking geometry for fast electron transfer is imposed by the guided trajectory of the ISP extrinsic domain. The requirement of a specific electrically neutral docking configuration for ISP electron transfer is consistent with structure data for the related cytochrome bc(1) complex.

Electron transfer between yeast cytochrome bc(1) complex and cytochrome c: a structural analysis

The structure of the complex between cytochrome c (CYC) and the cytochrome bc(1) complex (QCR) from yeast crystallized with an antibody fragment has been recently determined at 2.97 A resolution [Proc. Natl. Acad. Sci. U. S. A. 99 (2002) 2800]. CYC binds to subunit cytochrome c(1) of the enzyme stabilized by hydrophobic interactions surrounding the heme crevices creating a small, compact contact site. A central cation-pi interaction is an important and conserved feature of CYC binding. Peripheral patches with highly conserved complementary charges further stabilize the enzyme-substrate complex by long-range electrostatic forces and may affect the orientation of the substrate. Size and characteristics of the contact site are optimal for a transient electron transfer complex. Kinetic data show a bell-shaped ionic strength dependence of the cytochrome c reduction with a maximum activity near physiological ionic strength. The dependence is less pronounced in yeast compared to horse heart CYC indicating less impact of electrostatic interactions in the yeast system. Interestingly, a local QCR activity minimum is found for both substrates at 120-140 mM ionic strength. The architecture of the complex results in close distance of both c-type heme groups allowing the rapid reduction of cytochrome c by QCR via direct heme-to-heme electron transfer. Remarkably, CYC binds only to one of the two possible binding sites of the homodimeric complex and binding appears to be coordinated with the presence of ubiquinone at the Q(i) site. Regulatory aspects of CYC reduction are discussed.

Anabaena sp. PCC 7119 flavodoxin as electron carrier from photosystem I to ferredoxin-NADP+ reductase. Role of Trp(57) and Tyr(94)

The influence of the amino acid residues sandwiching the flavin ring in flavodoxin (Fld) from the cyanobacterium Anabaena sp. PCC 7119 in complex formation and electron transfer (ET) with its natural partners, photosystem I (PSI) and ferredoxin-NADP(+) reductase (FNR), was examined in mutants of the key residues Trp(57) and Tyr(94). The mutants' ability to form complexes with either FNR or PSI is similar to that of wild-type Fld. However, some of the mutants exhibit altered kinetic properties in their ET processes that can be explained in terms of altered flavin accessibility and/or thermodynamic parameters. The most noticeable alteration is produced upon replacement of Tyr(94) by alanine. In this mutant, the processes that involve the transfer of one electron from either PSI or FNR are clearly accelerated, which might be attributable to a larger accessibility of the flavin to the reductant. However, when the opposite ET flow is analyzed with FNR, the reduced Y94A mutant transfers electrons to FNR slightly more slowly than wild type. This can be explained thermodynamically from a decrease in driving force due to the significant shift of 137 mV in the reduction potential value for the semiquinone/hydroquinone couple (E(1)) of Y94A, relative to wild type (Lostao, A., Gómez-Moreno, C., Mayhew, S. G., and Sancho, J. (1997) Biochemistry 36, 14334-14344). The behavior of the rest of the mutants can be explained in the same way. Overall, our data indicate that Trp(57) and Tyr(94) do not play any active role in flavodoxin redox reactions providing a path for the electrons but are rather involved in setting an appropriate structural and electronic environment that modulates in vivo ET from PSI to FNR while providing a tight FMN binding.

Surface interactions in the complex between cytochrome f and the E43Q/D44N and E59K/E60Q plastocyanin double mutants as determined by (1)H-NMR chemical shift analysis

A combination of site-directed mutagenesis and NMR chemical shift perturbation analysis of backbone and side-chain protons has been used to characterize the transient complex of the photosynthetic redox proteins plastocyanin and cytochrome f. To elucidate the importance of charged residues on complex formation, the complex of cytochrome f and E43Q/D44N or E59K/E60Q spinach plastocyanin double mutants was studied by full analysis of the (1)H chemical shifts by use of two-dimensional homonuclear NMR spectra. Both mutants show a significant overall decrease in chemical shift perturbations compared with wild-type plastocyanin, in agreement with a large decrease in binding affinity. Qualitatively, the E43Q/D44N mutant showed a similar interaction surface as wild-type plastocyanin. The interaction surface in the E59K/E60Q mutant was distinctly different from wild type. It is concluded that all four charged residues contribute to the affinity and that residues E59 and E60 have an additional role in fine tuning the orientation of the proteins in the complex.

The unique proline of the Prochlorothrix hollandica plastocyanin hydrophobic patch impairs electron transfer to photosystem I

A number of surface residues of plastocyanin from Prochlorothrix hollandica have been modified by site-directed mutagenesis. Changes have been made in amino acids located in the amino-terminal hydrophobic patch of the copper protein, which presents a variant structure as compared with other plastocyanins. The single mutants Y12G, Y12F, Y12W, P14L, and double mutant Y12G/P14L have been produced. Their reactivity toward photosystem I has been analyzed by laser flash absorption spectroscopy. Plots of the observed rate constant with all mutants versus plastocyanin concentration show a saturation profile similar to that with wild-type plastocyanin, thus suggesting the formation of a plastocyanin-photosystem I transient complex. The mutations do not induce relevant changes in the equilibrium constant for complex formation but induce significant variations in the electron transfer rate constant, mainly with the two mutants at proline 14. Additionally, molecular dynamics calculations indicate that mutations at position 14 yield small changes in the geometry of the copper center. The comparative kinetic analysis of the reactivity of plastocyanin mutants toward photosystem I from different organisms (plants and cyanobacteria) reveals that reversion of the unique proline of Prochlorothrix plastocyanin to the conserved leucine of all other plastocyanins at this position enhances the reactivity of the Prochlorothrix protein.

Dynamic interaction of plastocyanin with the cytochrome bf complex

The interaction between plastocyanin and the intact cytochrome bf complex, both from spinach, has been studied by stopped-flow kinetics with mutant plastocyanin to elucidate the site of electron transfer and the docking regions of the molecule. Mutation of Tyr-83 to Arg or Leu provides no evidence for a second electron transfer path via Tyr-83 of plastocyanin, which has been proposed to be the site of electron transfer from cytochrome f. The data found with mutations of acidic residues indicate that both conserved negative patches are essential for the binding of plastocyanin to the intact cytochrome bf complex. Replacing Ala-90 and Gly-10 at the flat hydrophobic surface of plastocyanin by larger residues slowed down and accelerated, respectively, the rate of electron transfer as compared with wild-type plastocyanin. These opposing effects reveal that the hydrophobic region around the electron transfer site at His-87 is divided up into two regions, of which only that with Ala-90 contributes to the attachment to the cytochrome bf complex. These binding sites of plastocyanin are substantially different from those interacting with photosystem I. It appears that each of the two binding regions of plastocyanin is split into halves, which are used in different combinations in the molecular recognition at the two membrane complexes.

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.

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

The FMN to heme electron transfer in cytochrome P450BM-3. Effect of chemical modification of cysteines engineered at the FMN-heme domain interaction site

The crystal structure of the complex between the heme and FMN-containing domains of Bacillus megaterium cytochrome P450BM-3 (Sevrioukova, I. F., Li, H., Zhang, H., Peterson, J. A., and Poulos, T. L. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 1863-1868) indicates that the proximal side of the heme domain molecule is the docking site for the FMN domain and that the Pro(382)-Gln(387) peptide may provide an electron transfer (ET) path from the FMN to the heme iron. In order to evaluate whether ET complexes formed in solution by the heme and FMN domains are structurally relevant to that seen in the crystal structure, we utilized site-directed mutagenesis to introduce Cys residues at positions 104 and 387, which are sites of close contact between the domains in the crystal structure and at position 372 as a control. Cys residues were modified with a bulky sulfhydryl reagent, 1-dimethylaminonaphthalene-5-sulfonate-L-cystine (dansylcystine (DC)), to prevent the FMN domain from binding at the site seen in the crystal structure. The DC modification of Cys(372) and Cys(387) resulted in a 2-fold decrease in the rates of interdomain ET in the reconstituted system consisting of the separate heme and FMN domains and had no effect on heme reduction in the intact heme/FMN-binding fragment of P450BM-3. DC modification of Cys(104) caused a 10-20-fold decrease in the interdomain ET reaction rate in both the reconstituted system and the intact heme/FMN domain. This indicates that the proximal side of the heme domain molecule represents the FMN domain binding site in both the crystallized and solution complexes, with the area around residue 104 being the most critical for the redox partner docking.

Site-directed mutagenesis of cytochrome c(6) from Anabaena species PCC 7119. Identification of surface residues of the hemeprotein involved in photosystem I reduction

A number of surface residues of cytochrome c(6) from the cyanobacterium Anabaena sp. PCC 7119 have been modified by site-directed mutagenesis. Changes were made in six amino acids, two near the heme group (Val-25 and Lys-29) and four in the positively charged patch (Lys-62, Arg-64, Lys-66, and Asp-72). The reactivity of mutants toward the membrane-anchored complex photosystem I was analyzed by laser flash absorption spectroscopy. The experimental results indicate that cytochrome c(6) possesses two areas involved in the redox interaction with photosystem I: 1) a positively charged patch that may drive its electrostatic attractive movement toward photosystem I to form a transient complex and 2) a hydrophobic region at the edge of the heme pocket that may provide the contact surface for the transfer of electrons to P(700). The isofunctionality of these two areas with those found in plastocyanin (which acts as an alternative electron carrier playing the same role as cytochrome c(6)) are evident.

Auracyanin A from the thermophilic green gliding photosynthetic bacterium Chloroflexus aurantiacus represents an unusual class of small blue copper proteins

The amino acid sequence of the small copper protein auracyanin A isolated from the thermophilic photosynthetic green bacterium Chloroflexus aurantiacus has been determined to be a polypeptide of 139 residues. His58, Cys123, His128, and Met132 are spaced in a way to be expected if they are the evolutionary conserved metal ligands as in the known small copper proteins plastocyanin and azurin. Secondary structure prediction also indicates that auracyanin has a general beta-barrel structure similar to that of azurin from Pseudomonas aeruginosa and plastocyanin from poplar leaves. However, auracyanin appears to have sequence characteristics of both small copper protein sequence classes. The overall similarity with a consensus sequence of azurin is roughly the same as that with a consensus sequence of plastocyanin, namely 30.5%. We suggest that auracyanin A, together with the B forms, is the first example of a new class of small copper proteins that may be descendants of an ancestral sequence to both the azurin proteins occurring in prokaryotic nonphotosynthetic bacteria and the plastocyanin proteins occurring in both prokaryotic cyanobacteria and eukaryotic algae and plants. The N-terminal sequence region 1-18 of auracyanin is remarkably rich in glycine and hydroxy amino acids, and required mass spectrometric analysis to be determined. The nature of the blocking group X is not yet known, although its mass has been determined to be 220 Da. The auracyanins are the first small blue copper proteins found and studied in anoxygenic photosynthetic bacteria and are likely to mediate electron transfer between the cytochrome bc1 complex and the photosynthetic reaction center.

The effect of multiple binding modes on empirical modeling of ligand docking to proteins

A popular approach to the computational modeling of ligand/receptor interactions is to use an empirical free energy like model with adjustable parameters. Parameters are learned from one set of complexes, then used to predict another set. To improve these empirical methods requires an independent way to study their inherent errors. We introduce a toy model of ligand/receptor binding as a workbench for testing such errors. We study the errors incurred from the two state binding assumption--the assumption that a ligand is either bound in one orientation, or unbound. We find that the two state assumption can cause large errors in free energy predictions, but it does not affect rank order predictions significantly. We show that fitting parameters using data from high affinity ligands can reduce two state errors; so can using more physical models that do not use the two state assumption. We also find that when using two state models to predict free energies, errors are more severe on high affinity ligands than low affinity ligands. And we show that two state errors can be diagnosed by systematically adding new binding modes when predicting free energies: if predictions worsen as the modes are added, then the two state assumption in the fitting step may be at fault.

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.

Coupled plasmon-waveguide resonance spectroscopy studies of the cytochrome b6f/plastocyanin system in supported lipid bilayer membranes

The incorporation of cytochrome (cyt) b6f into a solid-supported planar egg phosphatidylcholine (PC) bilayer membrane and complex formation with plastocyanin have been studied by a variant of surface plasmon resonance called coupled plasmon-waveguide resonance (CPWR) spectroscopy, developed in our laboratory. CPWR combines greatly enhanced sensitivity and spectral resolution with direct measurement of anisotropies in refractive index and optical extinction coefficient, and can therefore probe structural properties of lipid-protein and protein-protein interactions. Cyt b6f incorporation into the membrane proceeds in two stages. The first occurs at low protein concentration and is characterized by an increase in total proteolipid mass without significant changes in the molecular order of the system, as demonstrated by shifts of the resonance position to larger incident angles without changing the refractive index anisotropy. The second stage, occurring at higher protein concentrations, results in a decrease in both the mass density and the molecular order of the system, evidenced by shifts of the resonance position to smaller incident angles and a large decrease in the membrane refractive index anisotropy. Plastocyanin can bind to such a proteolipid system in three different ways. First, the addition of plastocyanin before the second stage of b6f incorporation begins results in complex formation between the two proteins with a KD of approximately 10 microM and induces structural changes in the membrane that are similar to those occurring during the second stage of complex incorporation. The addition of larger amounts of plastocyanin under these conditions leads to nonspecific binding to the lipid phase with a KD of approximately 180 microM. Finally, the addition of plastocyanin after the completion of the second phase of b6f incorporation results in tighter binding between the two proteins (KD approximately 1 microM). Quantitation of the binding stoichiometry indicates that two plastocyanin molecules bind tightly to the dimeric form of the cyt b6f complex, assuming random insertion of the cytochrome into the bilayer. The structural basis for these results and formation of the proteolipid membrane are discussed.

Electron entry in a CuA mutant of cytochrome c oxidase from Paracoccus denitrificans. Conclusive evidence on the initial electron entry metal center

A cytochrome c oxidase subunit II C216S mutant from Paracoccus denitrificans in which the CuA site was changed by site-directed mutagenesis to a mononuclear copper site [Zickermann, V., Wittershagen, A., Kolbesen, B.O. and Ludwig, B. Biochemistry 36 (1997) 3232-3236] was investigated by stopped-flow spectroscopy. Contrary to the behavior of the wild type enzyme, in this mutant cytochrome a cannot be reduced by excess cytochrome c in the millisecond time scale in which cytochrome c oxidation is observed. The results conclusively identify and establish CuA as the initial electron entry site in cytochrome c oxidase. Partial rapid reduction (ca. 20%) of the modified CuA site suggests that the mononuclear copper ion has a redox potential ca. 100 mV lower than the wild type, and that internal electron transfer to cytochrome a is > or = 10(3)-fold slower than with the wild type enzyme.

Photoinduced electron transfer from carboxymethylated cytochrome c to plastocyanin

Photoinduced electron transfer from cytochrome c to plastocyanin was investigated using a novel method. Reduced carboxymethylated cytochrome c (CmCyt c), with carbon monoxide bound to the heme iron, and oxidized plastocyanin were mixed. At 1 mM CO the reduced state of CmCyt c is stabilized by about 350 meV. After flash photolysis of CO the apparent redox potential of CmCyt c drops resulting in electron transfer to plastocyanin. The electron transfer characteristics were investigated at approximately 30 different wavelengths in the range 390-460 nm. A global fit of the data showed that the electron transfer rate is 960+/-30 s-1 at pH 7.