Electrostatic effects on electron-transfer kinetics in the cytochrome f-plastocyanin complex

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.

Structure-function of the cytochrome b6f lipoprotein complex: a scientific odyssey and personal perspective

Structure-function studies of the cytochrome b₆f complex, the central hetero-oligomeric membrane protein complex in the electron transport chain of oxygenic photosynthesis, which formed the basis for a high-resolution (2.5 Å) crystallographic solution of the complex, are described. Structure-function differences between the structure of subunits of the bc complexes, b₆f, and bc₁ from mitochondria and photosynthetic bacteria, which are often assumed to function identically, are discussed. Major differences which suggest that quinone-dependent electron transport pathways can vary in b₆f and bc₁ complexes are as follows: (a) an additional c-type heme, c_n, and bound single copies of chlorophyll a and β-carotene in the b₆f complex; and (b) a cyclic electron transport pathway that encompasses the b₆f and PSI reaction center complexes. The importance of including lipid in crystallization of the cytochrome complex, or with any hetero-oligomeric membrane protein complex, is emphasized, and consequences to structure-function of b₆f being a lipoprotein complex discussed, including intra-protein dielectric heterogeneity and resultant pathways of trans-membrane electron transport. The role of the b₆f complex in trans-membrane signal transduction from reductant generated on the p-side of the electron transport chain to the regulation of light energy to the two photosystems by trans-side phosphorylation of the light-harvesting chlorophyll protein is presented. Regarding structure aspects relevant to plastoquinol-quinone entrance-egress: (i) modification of the p-side channel for plastoquinone access to the iron-sulfur protein would change the rate-limiting step in electron transport; (ii) the narrow niche for entry of plastoquinol into b₆f from the PSII reaction center complex would seem to require close proximity between the complexes.

Modeling and computations of the intramolecular electron transfer process in the two-heme protein cytochrome c(4)

The di-heme protein Pseudomonas stutzeri cytochrome c(4) (cyt c(4)) has emerged as a useful model for studying long-range protein electron transfer (ET). Recent experimental observations have shown a dramatically different pattern of intramolecular ET between the two heme groups in different local environments. Intramolecular ET in homogeneous solution is too slow (>10 s) to be detected but fast (ms-μs) intramolecular ET in an electrochemical environment has recently been achieved by controlling the molecular orientation of the protein assembled on a gold electrode surface. In this work we have performed computational modeling of the intramolecular ET process by a combination of density functional theory (DFT) and quantum mechanical charge transfer theory to disclose reasons for this difference. We first address the electronic structures of the model heme core with histidine and methionine axial ligands in both low- and high-spin states by structure-optimized DFT. The computations enable estimating the intramolecular reorganization energy of the ET process for different combinations of low- and high-spin heme couples. Environmental reorganization free energies, work terms ("gating") and driving force were determined using dielectric continuum models. We then calculated the electronic transmission coefficient of the intramolecular ET rate using perturbation theory combined with the electronic wave functions determined by the DFT calculations for different heme group orientations and Fe-Fe separations. The reactivity of low- and high-spin heme groups was notably different. The ET rate is exceedingly low for the crystallographic equilibrium orientation but increases by several orders of magnitude for thermally accessible non-equilibrium configurations. Deprotonation of the propionate carboxyl group was also found to enhance the ET rate significantly. The results are discussed in relation to the observed surface immobilization effect and support the notion of conformationally gated ET.

Reorganization energy for internal electron transfer in multicopper oxidases

We have calculated the reorganization energy for the intramolecular electron transfer between the reduced type 1 copper site and the peroxy intermediate of the trinuclear cluster in the multicopper oxidase CueO. The calculations are performed at the combined quantum mechanics and molecular mechanics (QM/MM) level, based on molecular dynamics simulations with tailored potentials for the two copper sites. We obtain a reorganization energy of 91-133 kJ/mol, depending on the theoretical treatment. The two Cu sites contribute by 12 and 22 kJ/mol to this energy, whereas the solvent contribution is 34 kJ/mol. The rest comes from the protein, involving small contributions from many residues. We have also estimated the energy difference between the two electron-transfer states and show that the reduction of the peroxy intermediate is exergonic by 43-87 kJ/mol, depending on the theoretical method. Both the solvent and the protein contribute to this energy difference, especially charged residues close to the two Cu sites. We compare these estimates with energies obtained from QM/MM optimizations and QM calculations in a vacuum and discuss differences between the results obtained at various levels of theory.

Redox potential of the Rieske iron-sulfur protein quantum-chemical and electrostatic study

Quantum-chemical study of structures, energies, and effective partial charge distribution for several models of the Rieske protein redox center is performed in terms of the B3LYP density functional method in combination with the broken symmetry approach using three different atomic basis sets. The structure of the redox complex optimized in vacuum differs markedly from that inside the protein. This means that the protein matrix imposes some stress on the active site resulting in distortion of its structure. The redox potentials calculated for the real active site structure are in a substantially better agreement with the experiment than those calculated for the idealized structure. This shows an important role of the active site distortion in tuning its redox potential. The reference absolute electrode potential of the standard hydrogen electrode is used that accounts for the correction caused by the water surface potential. Electrostatic calculations are performed in the framework of the polarizable solute model. Two dielectric permittivities of the protein are employed: the optical permittivity for calculation of the intraprotein electric field, and the static permittivity for calculation of the dielectric response energy. Only this approach results in a reasonable agreement of the calculated and experimental redox potentials.

Semi-continuum electrostatic calculations of redox potentials in photosystem I

The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill.

Structure of the Complex between Plastocyanin and Cytochrome f from the Cyanobacterium Nostoc sp. PCC 7119 as Determined by Paramagnetic NMR

The complex between cytochrome f and plastocyanin from the cyanobacterium Nostoc has been characterized by NMR spectroscopy. The binding constant is 16 mM(-1), and the lifetime of the complex is much less than 10 ms. Intermolecular pseudo-contact shifts observed for the plastocyanin amide nuclei, caused by the heme iron, as well as the chemical-shift perturbation data were used as the sole experimental restraints to determine the orientation of plastocyanin relative to cytochrome f with a precision of 1.3 angstroms. The data show that the hydrophobic patch surrounding tyrosine 1 in cytochrome f docks the hydrophobic patch of plastocyanin. Charge complementarities are found between the rims of the respective recognition sites of cytochrome f and plastocyanin. Significant differences in the relative orientation of both proteins are found between this complex and those previously reported for plants and Phormidium, indicating that electrostatic and hydrophobic interactions are balanced differently in these complexes.

Structure of the Intermolecular Complex between Plastocyanin and Cytochrome f from Spinach

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

Modeling electron transfer thermodynamics in protein complexes: interaction between two cytochromes c(3)

Redox protein complexes between type I and type II tetraheme cytochromes c(3) from Desulfovibrio vulgaris Hildenborough are here analyzed using theoretical methodologies. Various complexes were generated using rigid-body docking techniques, and the two lowest energy complexes (1 and 2) were relaxed using molecular dynamics simulations with explicit solvent and subjected to further characterization. Complex 1 corresponds to an interaction between hemes I from both cytochromes c(3). Complex 2 corresponds to an interaction between the heme IV from type I and the heme I from type II cytochrome c(3). Binding free energy calculations using molecular mechanics, Poisson-Boltzmann, and surface accessibility methods show that complex 2 is more stable than complex 1. Thermodynamic calculations on complex 2 show that complex formation induces changes in the reduction potential of both cytochromes c(3), but the changes are larger in the type I cytochrome c(3) (the largest one occurring on heme IV, of approximately 80 mV). These changes are sufficient to invert the global titration curves of both cytochromes, generating directionally in electron transfer from type I to type II cytochrome c(3), a phenomenon of obvious thermodynamic origin and consequences, but also with kinetic implications. The existence of processes like this occurring at complex formation may constitute a natural design of efficient redox chains.

Ironies in photosynthetic electron transport: a personal perspective

In this brief personal perspective, studies, mainly in my laboratory, directed toward understanding of the structure and function of the Fe-metalloprotein electron transport carriers of oxygenic photosynthesis are discussed. This level of understanding and the description of the events that led to this level are inevitably incomplete. The Fe-proteins considered are those, cytochromes b-559 (560), f, b (6), and x, and the Rieske [2Fe-2S] iron-sulfur protein, which have now been placed in a structural context as a result of X-ray diffraction analysis of crystals of Photosystem II and the cytochrome b (6) f complex and the Photosystem II reaction center.

Metalloprotein association, self-association, and dynamics governed by hydrophobic interactions: simultaneous occurrence of gated and true electron-transfer reactions between cytochrome f and cytochrome c(6) from Chlamydomonas reinhardtii

Noninvasive reconstitution of the heme in cytochrome c(6) with zinc(II) ions allowed us to study the photoinduced electron-transfer reaction (3)Zncyt c(6) + cyt f(III) --> Zncyt c(6)(+) + cyt f(II) between physiological partners cytochrome c(6) and cytochrome f, both from Chlamydomonas reinhardtii. The reaction kinetics was analyzed in terms of protein docking and electron transfer. In contrast to various protein pairs studied before, both the unimolecular and the bimolecular reactions of this oxidative quenching take place at all ionic strengths from 2.5 through 700 mM. The respective intracomplex rate constants are k(uni) (1.2 +/- 0.1) x 10(4) s(-1) for persistent and k(bi) (9 +/- 4) x 10(2) s(-1) for the transient protein complex. The former reaction seems to be true electron transfer, and the latter seems to be electron transfer gated by a structural rearrangement. Remarkably, these reactions occur simultaneously, and both rate constants are invariant with ionic strength. The association constant K(a) for zinc cytochrome c(6) and cytochrome f(III) remains (5 +/- 3) x 10(5) M(-1) in the ionic strength range from 700 to 10 mM and then rises slightly to (7 +/- 2) x 10(6) M(-1), as ionic strength is lowered to 2.5 mM. Evidently, docking of these proteins from C. reinhardtii is due to hydrophobic interaction, slightly augmented by weak electrostatic attraction. Kinetics, chromatography, and cross-linking consistently show that cytochrome f self-dimerizes at ionic strengths of 200 mM and higher. Cytochrome f(III) quenches triplet state (3)Zncyt c(6), but its dimer does not. Formation of this unreactive dimer is an important step in the mechanism of electron transfer. Not only association between the reacting proteins, but also their self-association, should be considered when analyzing reaction mechanisms.

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.

Geometry, reduction potential, and reorganization energy of the binuclear Cu(A) site, studied by density functional theory

The dimeric Cu(A) site found in cytochrome c oxidase and nitrous oxide reductase has been studied with the density functional B3LYP method. We have optimized the structure of the realistic (Im)(S(CH(3))(2))Cu(SCH(3))(2)Cu(Im)(CH(3)CONHCH(3)) model in the fully reduced, mixed-valence, and fully oxidized states. The optimized structures are very similar to crystal structures of the protein, which shows that the protein does not strain the site significantly. Instead, inorganic model complexes of the protein site are strained by the macrocyclic connections between the ligand models. For the mixed-valence (Cu(I)+Cu(II)) state, two distinct equilibrium structures were found, one with a short Cu-Cu distance, 248 pm, similar to the protein structure, and one with a longer distance, 310 pm, similar to what is found in inorganic models. In the first state, the unpaired electron is delocalized over both copper ions, whereas in the latter, it is more localized to one of the ions. The two states are nearly degenerate. The potential energy surfaces for the Cu-Cu, Cu-S(Met), and Cu-O interactions are extremely flat. In fact, all three distances can be varied between 230 and 310 pm at an expense in energy of less than 8 kJ/mol, which explains the large variation observed in crystal structures for these interactions. Inclusion of solvation effects does not change this significantly. Therefore, we can conclude that a variation in these distances can change the reduction potential of the Cu(A) site by at most 100 mV. The model complex has a reorganization energy of 43 kJ/mol, 20 kJ/mol lower than for a monomeric blue-copper site. This lowering is caused by the delocalization of the unpaired electron in the mixed-valence state.

Electron transfer and stability of the cytochrome b6f complex in a small domain deletion mutant of cytochrome f

The lumen segment of cytochrome f consists of a small and a large domain. The role of the small domain in the biogenesis and stability of the cytochrome b(6)f complex and electron transfer through the cytochrome b(6)f complex was studied with a small domain deletion mutant in Chlamydomonas reinhardtii. The mutant is able to grow photoautotrophically but with a slower rate than the wild type strain. The heme group is covalently attached to the polypeptide, and the visible absorption spectrum of the mutant protein is identical to that of the native protein. The kinetics of electron transfer in the mutant were measured by flash kinetic spectroscopy. Our results show that the rate for the oxidation of cytochrome f was unchanged (t(12) = approximately 100 micros), but the half-time for the reduction of cytochrome f is increased (t(12) = 32 ms; for wild type, t(12) = 2.1 ms). Cytochrome b(6) reduction was slower than that of the wild type by a factor of approximately 2 (t(12) = 8.6 ms; for wild type, t(12) = 4.7 ms); the slow phase of the electrochromic band shift also displayed a slower kinetics (t(12) = 5.5 ms; for wild type, t(12) = 2.7 ms). The stability of the cytochrome b(6)f complex in the mutant was examined by following the kinetics of the degradation of the individual subunits after inhibiting protein synthesis in the chloroplast. The results indicate that the cytochrome b(6)f complex in the small domain deletion mutant is less stable than in the wild type. We conclude that the small domain is not essential for the biogenesis of cytochrome f and the cytochrome b(6)f complex. However, it does have a role in electron transfer through the cytochrome b(6)f complex and contributes to the stability of the complex.

X-ray structure of a truncated form of cytochrome f from chlamydomonas reinhardtii

A truncated form of cytochrome f from Chlamydomonas reinhardtii (an important eukaryotic model organism for photosynthetic electron transfer studies) has been crystallized (space group P2(1)2(1)2(1); three molecules/asymmetric unit) and its structure determined to 2.0 A resolution by molecular replacement using the coordinates of a truncated turnip cytochrome f as a model. The structure displays the same folding and detailed features as turnip cytochrome f, including (a) an unusual heme Fe ligation by the alpha-amino group of tyrosine 1, (b) a cluster of lysine residues (proposed docking site of plastocyanin), and (c) the presence of a chain of seven water molecules bound to conserved residues and extending between the heme pocket and K58 and K66 at the lysine cluster. For this array of waters, we propose a structural role. Two cytochrome f molecules are related by a noncrystallographic symmetry operator which is a distorted proper 2-fold rotation. This may represent the dimeric relation of the monomers in situ; however, the heme orientation suggested by this model is not consistent with previous EPR measurements on oriented membranes.

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.

Side-chain interactions in the plastocyanin-cytochrome f complex

Cytochrome f and plastocyanin are redox partners in the photosynthetic electron-transfer chain. Electron transfer from cytochrome f to plastocyanin occurs in a specific short-lived complex. To obtain detailed information about the binding interface in this transient complex, the effects of binding on the backbone and side-chain protons of plastocyanin have been analyzed by mapping NMR chemical-shift changes. Cytochrome f was added to plastocyanin up to 0.3 M equiv, and the plastocyanin proton chemical shifts were measured. Out of approximately 500 proton resonances, 86% could be observed with this method. Nineteen percent demonstrate significant chemical-shift changes and these protons are located in the hydrophobic patch (including the copper ligands) and the acidic patches of plastocyanin, demonstrating that both areas are part of the interface in the complex. This is consistent with the recently determined structure of the complex [Ubbink, M., Ejdebäck, M., Karlsson, B. G., and Bendall, D. S. (1998) Structure 6, 323-335]. The largest chemical-shift changes are found around His87 in the hydrophobic patch, which indicates tight contacts and possibly water exclusion from this part of the protein interface. These results support the idea that electron transfer occurs via His87 to the copper in plastocyanin and suggest that the hydrophobic patch determines the specificity of the binding. The chemical-shift changes in the acidic patches are significant but small, suggesting that the acidic groups are involved in electrostatic interactions but remain solvent exposed. The existence of small differences between the present data and those used for the structure may imply that the redox state of the metals in both proteins slightly affects the structure of the complex. The chemical-shift mapping is performed on unlabeled proteins, making it an efficient way to analyze effects of mutations on the structure of the complex.

Low dielectric response in enzyme active site

The kinetics of charge transfer depend crucially on the dielectric reorganization of the medium. In enzymatic reactions that involve charge transfer, atomic dielectric response of the active site and of its surroundings determines the efficiency of the protein as a catalyst. We report direct spectroscopic measurements of the reorganization energy associated with the dielectric response in the active site of alpha-chymotrypsin. A chromophoric inhibitor of the enzyme is used as a spectroscopic probe. We find that water strongly affects the dielectric reorganization in the active site of the enzyme in solution. The reorganization energy of the protein matrix in the vicinity of the active site is similar to that of low-polarity solvents. Surprisingly, water exhibits an anomalously high dielectric response that cannot be described in terms of the dielectric continuum theory. As a result, sequestering the active site from the aqueous environment inside low-dielectric enzyme body dramatically reduces the dielectric reorganization. This reduction is particularly important for controlling the rate of enzymatic reactions.

Redox- and pH-dependent association of plastocyanin with lipid bilayers: effect on protein conformation and thermal stability

The effect of electrostatic interactions on the conformation and thermal stability of plastocyanin (Pc) was studied by infrared spectroscopy. Association of any of the two redox states of the protein with positively charged membranes at neutral pH does not significantly change the secondary structure of Pc. However, upon membrane binding, the denaturation temperature decreases, regardless of the protein redox state. The extent of destabilization depends on the proportion of positively charged lipid headgroups in the membrane, becoming greater as the surface density of basic phospholipids increases. In contrast, at pH 4.8 the membrane binding-dependent conformational change becomes redox-sensitive. While the secondary structures and thermal stabilities of free and membrane-bound oxidized Pc are similar under acidic conditions, the conformation of the reduced form of the protein drastically rearranges upon membrane association. This rearrangement does not depend on electrostatic interactions to occur, since it is also observed in the presence of uncharged lipid bilayers. The conformational transition, only observed for reduced Pc, involves the exposure of hydrophobic regions that leads to intermolecular interactions at the membrane surface. Membrane-mediated partial unfolding of reduced Pc can be reversed by readjusting the pH to neutrality, in the absence of electrostatic interactions. This redox-dependent behavior might reflect specific structural requirements for the interaction of Pc with its redox partners.

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.

Free energy of charge transfer and intraprotein electric field: method of calculation depends on the charge state of protein at a given structure

Free energy of charge transfer presents a basic characteristic of reactions such as protonation, oxido-reduction and similar. Evaluation of this quantity requires calculation of charging energy. Proteins are structured dielectrics, and a consistent incorporation of their structure into calculation of intraprotein electric field results in expression for charging energy of an active group in protein, which is essentially different from that for a simple dielectric. An algorithm for semi-continuum calculation of relevant free energies is described. First of the two components of charging energy in protein, energy of the medium response to charge redistribution in reactants, should be always calculated as the charging energy by the charge redistribution using the static dielectric constant of protein. The second term is interaction energy of the charge redistribution with the 'frozen' electric field of the system before reaction. Charges of protein groups, at which the protein structure has been determined, are often different from those before reaction of charge transfer, so is the corresponding intraprotein field. The field is expressed through either both the optical and static dielectric constants of protein or only optical one depending on whether the charges of protein groups before reaction and upon structural analysis are the same or not. Proper allowance for difference in charges of reacting groups before reaction and upon structural analysis of protein is thermodynamically necessary and quantitatively important. The expression for activation free energy for charge transfer in proteins is derived in the form presenting explicitly an invariant contribution of protein structure.

Quantum chemical calculations of the reorganization energy of blue-copper proteins

The inner-sphere reorganization energy for several copper complexes related to the active site in blue-copper protein has been calculated with the density functional B3LYP method. The best model of the blue-copper proteins, Cu(Im)2(SCH3)(S(CH3)2)(0/+), has a self-exchange inner-sphere reorganization energy of 62 kJ/mol, which is at least 120 kJ/mol lower than for Cu(H2O)4(+/2+). This lowering of the reorganization energy is caused by the soft ligands in the blue-copper site, especially the cysteine thiolate and the methionine thioether groups. Soft ligands both make the potential surfaces of the complexes flatter and give rise to oxidized structures that are quite close to a tetrahedron (rather than tetragonal). Approximately half of the reorganization energy originates from changes in the copper-ligand bond lengths and half of this contribution comes from the Cu-S(Cys) bond. A tetragonal site, which is present in the rhombic type 1 blue-copper proteins, has a slightly higher (16 kJ/mol) inner-sphere reorganization energy than a trigonal site, present in the axial type 1 copper proteins. A site with the methionine ligand replaced by an amide group, as in stellacyanin, has an even higher reorganization energy, about 90 kJ/mol.

Identification of the basic residues of cytochrome f responsible for electrostatic docking interactions with plastocyanin in vitro: relevance to the electron transfer reaction in vivo

The prominent basic patch seen in the atomic structure of the lumen-side domain of turnip cytochrome f, consisting of Arg209 and Lys187, 58, 65, and 66, was proposed to be an electrostatically complementary docking site for its physiological electron acceptor, plastocyanin [Martinez, S. E., Huang, D., Szczepaniak, A., Cramer, W. A., and Smith, J. L. (1994) Structure 2, 95-105]. This proposal agrees with solution studies on the cytochrome f/plastocyanin electron-transfer reaction that showed a major contribution of electrostatic interactions to the docking, but not with studies on the oxidation rate of cyt f in vivo using mutants in which the basic patch of cyt f was neutralized. The apparent contradiction might be explained by an unknown electron acceptor protein for cyt f. However, (i) flash-induced oxidation of cyt f is absent in a PC-deficient mutant. (ii) Lys58, 65, and 66 in the large domain and Lys188 and 189 in the small domain are major contributors to the ionic strength dependence of the electron-transfer reaction in solution. Replacement of Lys58 and 65 by neutral residues and of Lys66 by the acidic residue Glu66 resulted in a >10-fold decrease in the rate of electron transfer in solution and complete loss of its ionic strength dependence. Replacement of Lys188 and Lys189 in the small domain of cyt f resulted in a 3-4-fold decrease in the second-order rate constant and a smaller dependence of the overall rate of electron transfer on ionic strength, corresponding to a loss of two positive charges. (iii) Acidification of the thylakoid lumen cannot explain the absence of electrostatic interactions. (iv) Changing the five lysines to acidic residues did not result in any significant retardation of the rate of cyt f oxidation in vivo. If the docking of cyt f and plastocyanin in vivo is mediated by basic residues of cyt f, they are different from those that mediate electron transfer in vitro or that are implicated by simulations of electrostatic interactions of the docking. Alternatively, docking of cyt f/PC in vivo is limited by spatial constraints or release of PC from P700 that precludes a rate-limiting mediation of the cyt f/PC reaction by specific electrostatic interactions. The cyt f/PC system in Chlamydomonas reinhardtii is the first electron-transfer couple for which the role of electrostatics in mediating the docking reaction has been studied both in vitro and in vivo.

Electrogenicity accompanies photoreduction of the iron-sulfur clusters F(A) and F(B) in photosystem I

Photovoltage responses accompanying electron transfer on the acceptor side of photosystem I (PS I) were investigated in proteoliposomes containing PS I complexes from the cyanobacterium Synechococcus sp. PCC 6301 using a direct electrometrical technique. The relative contributions of the F(X) --> F(B) and the F(X) --> F(A) electron transfer reactions to the overall electrogenicity were elucidated by comparing the sodium dithionite-induced decrease in the magnitude of the total photoelectric responses in control and in F(B)-less (HgCl2-treated) PS I complexes. The results obtained suggest that the electrogenesis on the acceptor side of PS I is related to electron transfers between both F(X) and F(A) and F(A) and F(B). Based on the electrogenic nature of the latter reaction in PS I complexes, we conclude that F(A) rather than F(B) is the acceptor proximal to F(X).