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

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

Interactions between cytochrome c2 and the photosynthetic reaction center from Rhodobacter sphaeroides: the cation-pi interaction

The cation-pi interaction between positively charged and aromatic groups is a common feature of many proteins and protein complexes. The structure of the complex between cytochrome c(2) (cyt c(2)) and the photosynthetic reaction center (RC) from Rhodobacter sphaeroides exhibits a cation-pi complex formed between Arg-C32 on cyt c(2) and Tyr-M295 on the RC [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. The importance of the cation-pi interaction for binding and electron transfer was studied by mutating Tyr-M295 and Arg-C32. The first- and second-order rates for electron transfer were not affected by mutating Tyr-M295 to Ala, indicating that the cation-pi complex does not greatly affect the association process or structure of the state active in electron transfer. The dissociation constant K(D) showed a greater increase when Try-M295 was replaced with nonaromatic Ala (3-fold) as opposed to aromatic Phe (1.2-fold), which is characteristic of a cation-pi interaction. Replacement of Arg-C32 with Ala increased K(D) (80-fold) largely due to removal of electrostatic interactions with negatively charged residues on the RC. Replacement with Lys increased K(D) (6-fold), indicating that Lys does not form a cation-pi complex. This specificity for Arg may be due to a solvation effect. Double mutant analysis indicates an interaction energy between Tyr-M295 and Arg-C32 of approximately -24 meV (-0.6 kcal/mol). This energy is surprisingly small considering the widespread occurrence of cation-pi complexes and may be due to the tradeoff between the favorable cation-pi binding energy and the unfavorable desolvation energy needed to bury Arg-C32 in the short-range contact region between the two proteins.     

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

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

The structure and function of the cytochrome <Emphasis Type="Italic">c</Emphasis><Subscript>2</Subscript>: reaction center electron transfer complex from <Emphasis Type="Italic">Rhodobacter sphaeroides</Emphasis>

In the photosynthetic bacterium, Rhodobacter sphaeroides, the mobile electron carrier, cytochrome c₂ (cyt c₂) transfers an electron from reduced heme to the photooxidized bacteriochlorophyll dimer in the membrane bound reaction center (RC) as part of the light induced cyclic electron transfer chain. A complex between these two proteins that is active in electron transfer has been crystallized and its structure determined by X-ray diffraction. The structure of the cyt:RC complex shows the cyt c₂ (cyt c₂) positioned at the center of the periplasmic surface of the RC. The exposed heme edge from cyt c₂ is in close tunneling contact with the electron acceptor through an intervening bridging residue, Tyr L162 located on the RC surface directly above the bacteriochlorophyll dimer. The binding interface between the two proteins can be divided into two regions: a short-range interaction domain and a long-range interaction domain. The short-range domain includes residues immediately surrounding the tunneling contact region around the heme and Tyr L162 that display close intermolecular contacts optimized for electron transfer. These include a small number of hydrophobic interactions, hydrogen bonds and a pi-cation interaction. The long-range interaction domain consists of solvated complementary charged residues; positively charged residues from the cyt and negatively charged residues from the RC that provide long range electrostatic interactions that can steer the two proteins into position for rapid association.     

Electrostatic interactions in an integral membrane protein

The effects of charge-charge interactions on the midpoint reduction potential (E(m)()) of the primary electron donor (P) in the photosynthetic reaction center of Rhodobacter sphaeroides were investigated by introducing mutations of ionizable amino acids at selected sites. The mutations were designed to alter the electrostatic environment of P, a bacteriochlorophyll dimer, without greatly affecting its structure or molecular orbitals. Two arginine residues at homologous positions in the L and M subunits [residues (L135) and (M164)], Asp (L155), Tyr (L164), and Cys (L247) were changed independently. Arginine (L135) was replaced by Lys, Leu, Gln, or Glu; Arg (M164), by Leu or Glu; Asp (L155), by Asn; Tyr (L164), by Phe; and Cys (L247), by Lys or Asp. The R(L135)E/C(L247)K double mutant also was made. The shift in the E(m)() of P/P(+) was measured in each mutant and was compared with the effect predicted by electrostatics calculations using several different computational approaches. A simple distance-dependent dielectric screening factor reproduced the effects remarkably well. By contrast, microscopic methods that considered the reaction field in the protein and solvent but did not include explicit counterions overestimated the changes in the E(m)() considerably. Including counterions for the charged residues reduced the calculated effects of the mutations in molecular dynamics calculations. The results show that electrostatic interactions of P with ionizable amino acid residues are strongly screened, and suggest that counterions make major contributions to this screening. The screening also could reflect penetration of water or other relaxations not taken into account because of incomplete sampling of configurational space.     

Effects of ionizable residues on the absorption spectrum and initial electron-transfer kinetics in the photosynthetic reaction center of Rhodobacter sphaeroides

Effects of ionizable amino acids on spectroscopic properties and electron-transfer kinetics in the photosynthetic reaction center (RC) of Rhodobacter sphaeroides are investigated by site-directed mutations designed to alter the electrostatic environment of the bacteriochlorophyll dimer that serves as the photochemical electron donor (P). Arginine residues at homologous positions in the L and M subunits (L135 and M164) are changed independently: Arg L135 is replaced by Lys, Leu, Glu, and Gln and Arg M164 by Leu and Glu. Asp L155 also is mutated to Asn, Tyr L164 to Phe, and Cys L247 to Lys and Asp. The mutations at L155, L164, and M164 have little effect on the absorption spectrum, whereas those at L135 and L247 shift the long-wavelength absorption band of P to higher energies. Fits to the ground-state absorption and hole-burned spectra indicate that the blue shift and increased width of the absorption band in the L135 mutants are due partly to changes in the distribution of energies for the zero-phonon absorption line and partly to stronger electron-phonon coupling. The initial electron-transfer kinetics are not changed significantly in most of the mutants, but the time constant increases from 3.0 +/- 0.2 in wild-type RCs to 4.7 +/- 0.2 in C(L247)D and 7.0 +/- 0.3 ps in C(L247)K. The effects of the mutations on the solvation free energies of the product of the initial electron-transfer reaction (P(+)) and the charge-transfer states that contribute to the absorption spectrum ( and ) were calculated by using a distance-dependent electrostatic screening factor. The results are qualitatively in accord with the view that electrostatic interactions of the bacteriochlorophylls with ionized residues of the protein are strongly screened and make only minor contributions to the energetics and dynamics of charge separation. However, the slowing of electron transfer in the Cys L247 mutants and the blue shift of the spectrum in some of the Arg L135 and Cys L247 mutants cannot be explained consistently by electrostatic interactions of the mutated residues with P and B(L); we ascribe these effects tentatively to structural changes caused by the mutations.     

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

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

X-ray structure determination of the cytochrome c2: reaction center electron transfer complex from Rhodobacter sphaeroides

In the photosynthetic bacterium Rhodobacter sphaeroides, a water soluble cytochrome c2 (cyt c2) is the electron donor to the reaction center (RC), the membrane-bound pigment-protein complex that is the site of the primary light-induced electron transfer. To determine the interactions important for docking and electron transfer within the transiently bound complex of the two proteins, RC and cyt c2 were co-crystallized in two monoclinic crystal forms. Cyt c2 reduces the photo-oxidized RC donor (D+), a bacteriochlorophyll dimer, in the co-crystals in approximately 0.9 micros, which is the same time as measured in solution. This provides strong evidence that the structure of the complex in the region of electron transfer is the same in the crystal and in solution. X-ray diffraction data were collected from co-crystals to a maximum resolution of 2.40 A and refined to an R-factor of 22% (R(free)=26%). The structure shows the cyt c2 to be positioned at the center of the periplasmic surface of the RC, with the heme edge located above the bacteriochlorophyll dimer. The distance between the closest atoms of the two cofactors is 8.4 A. The side-chain of Tyr L162 makes van der Waals contacts with both cofactors along the shortest intermolecular electron transfer pathway. The binding interface can be divided into two domains: (i) A short-range interaction domain that includes Tyr L162, and groups exhibiting non-polar interactions, hydrogen bonding, and a cation-pi interaction. This domain contributes to the strength and specificity of cyt c2 binding. (ii) A long-range, electrostatic interaction domain that contains solvated complementary charges on the RC and cyt c2. This domain, in addition to contributing to the binding, may help steer the unbound proteins toward the right conformation.     

Mapping the electron transfer interface between cytochrome b5 and cytochrome c

To characterize the cytochrome b(5) (Cyt b(5))-cytochrome c (Cyt c) interactions during electron transfer, variants of Cyt b(5) have been employed to assess the contributions of electrostatic interactions (substitution of surface charged residues Glu44, Glu48, Glu56, and Asp60 and heme propionate), hydrophobic interactions, and the thermodynamic driving forces (substitutions for hydrophobic residues in heme pocket residues Phe35, Pro40, Val45, Phe58, and Val61). The electrostatic interactions play an important role in maintaining the stability and specificity of the Cyt b(5)-Cyt c complex that is formed. There is no essential effect on the intraprotein complex electron transfer even if most of the involved negatively charged residues on the surface of Cyt b(5) have been removed. The results support a dynamic docking paradigm for Cyt b(5)-Cyt c interactions. The orientation that is optimal for binding may not be optimal form for electron transfer. Substitution of hydrophobic residues does not have a significant effect on the binding between Cyt b(5) and Cyt c; rather, it regulates the electron transfer rates via changes in the driving force. Combining the electron transfer studies of the Cyt b(5)-Cyt c system and the Cyt b(5)-Zn-Cyt c system, we obtain the reorganization energy (0.6 eV) at an ionic strength of 150 mM.     

Models for the complexes formed between cytochrome b5 and the subunits of methemoglobin

Computer graphics-generated models for the electron transfer complexes formed between cytochrome b5 and the subunits of methemoglobin are proposed. For both complexes, the orientation allowing optimal hydrogen bonding involves interaction between negatively charged residues on cytochrome b5 and positively charged residues on methemoglobin. In each complex, the heme groups of the interacting species are coplanar with the edges of the heme groups separated by 7-8 A and with the iron atoms 16 A apart. For the alpha-chain X cytochrome b5 complex, alpha-chain residues 56 (Lys), 60 (Lys), and 90 (Lys) interact with cytochrome b5 residues 44 (Glu), 43 (Glu), and 60 (Asp) respectively. A fourth hydrogen bond involves alpha-61 (Lys) bridging between a heme propionate from cytochrome b5 and a heme propionate from the alpha-chain. The contacts present in the beta-chain X cytochrome b5 complex involve hydrogen-bonding between beta-chain lysyl residues 59, 61, 65, and 95, and cytochrome b5 residues 48 (Glu), 44 (Glu), 43 (Glu), and 60 (Asp) respectively. An additional hydrogen bond can be formed by bridging of the epsilon-amino group of beta-66 (Lys) between a heme propionate from cytochrome b5 and a beta-chain heme propionate. In each complex, two nonionic interactions, one on each side of the heme groups, are also suggested. These interactions appear to effectively exclude external water molecules from the center of the protein-protein interaction domain. Comparison of the proposed binding loci for cytochrome b5 on the methemoglobin subunits with those proposed on cytochrome c reveals considerable structural homology between the cytochrome b5 binding sites.     

Stabilization of a fibronectin type III domain by the removal of unfavorable electrostatic interactions on the protein surface

It is generally considered that electrostatic interactions on the protein surface, such as ion pairs, contribute little to protein stability, although they may play important roles in conformational specificity. We found that the tenth fibronectin type III domain of human fibronectin (FNfn10) is more stable at acidic pH than neutral pH, with an apparent midpoint of transition near pH 4. Determination of pK(a)'s for all the side chain carboxyl groups of Asp and Glu residues revealed that Asp 23 and Glu 9 have an upshifted pK(a). These residues and Asp 7 form a negatively charged patch on the surface of FNfn10, with Asp 7 centrally located between Asp 23 and Glu 9, suggesting repulsive electrostatic interactions among these residues at neutral pH. Mutant proteins, D7N and D7K, in which Asp 7 was replaced with Asn and Lys, respectively, exhibited a modest but significant increase in stability at neutral pH, compared to the wild type, and they no longer showed pH dependence of stability. The pK(a)'s of Asp 23 and Glu 9 in these mutant proteins shifted closer to their respective unperturbed values, indicating that the unfavorable electrostatic interactions have been reduced in the mutant proteins. Interestingly, the wild-type and mutant proteins were all stabilized to a similar degree by the addition of 1 M sodium chloride at both neutral and acidic pH, suggesting that the repulsive interactions between the carboxyl groups cannot be effectively shielded by 1 M sodium chloride. These results indicate that repulsive interactions between like charges on the protein surface can destabilize a protein, and protein stability can be significantly improved by relieving these interactions.     

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.     

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

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

Interactions between cytochrome c2 and photosynthetic reaction center from Rhodobacter sphaeroides: changes in binding affinity and electron transfer rate due to mutation of interfacial hydrophobic residues are strongly correlated

The structure of the complex between cytochrome c(2) (cyt) and the photosynthetic reaction center (RC) from Rhodobacter sphaeroides shows contacts between hydrophobic residues Tyr L162, Leu M191, and Val M192 on the RC and the surface of the cyt [Axelrod et al. (2002) J. Mol. Biol. 319, 501-515]. The role of these hydrophobic residues in binding and electron transfer was investigated by replacing them with Ala and other residues. Mutations of the hydrophobic residues generally resulted in relatively small changes in the second-order electron-transfer rate k(2) (Br?nsted coefficient, alpha( )()= 0.15 +/- 0.05) indicating that the transition state for association occurs before short-range hydrophobic contacts are established. Larger changes in k(2), found in some cases, were attributed to a change in the second-order mechanism from a diffusion controlled regime to a rapidly reversible binding regime. The association constant, K(A), of the cyt and the rate of electron transfer from the bound cyt, k(e), were both decreased by mutation. Replacement of Tyr L162, Leu M191, or Val M192 by Ala decreased K(A) and k(e) by factors of 130, 10, 0.6, and 120, 9, 0.6, respectively. The largest changes were obtained by mutation of Tyr L162, showing that this residue plays a key role in both binding and electron transfer. The binding affinity, K(A), and electron-transfer rate, k(e) were strongly correlated, showing that changes of hydrophobic residues affect both binding and electron transfer. This correlation suggests that changes in distance across hydrophobic interprotein contacts have similar effects on both electron tunneling and binding interactions.     

Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.     

EPR, ENDOR, and Special TRIPLE measurements of P•+ in wild type and modified reaction centers from Rb. sphaeroides

The influence of the protein environment on the primary electron donor, P, a bacteriochlorophyll a dimer, of reaction centers from Rhodobacter sphaeroides, has been investigated using electron paramagnetic resonance and electron nuclear double resonance spectroscopy. These techniques were used to probe the effects on P that are due to alteration of three amino acid residues, His L168, Asn L170, and Asn M199. The introduction of Glu at L168, Asp at L170, or Asp at M199 changes the oxidation/reduction midpoint potential of P in a pH-dependent manner (Williams et al. (2001) Biochemistry 40, 15403-15407). For the double mutant His L168 to Glu and Asn at L170 to Asp, excitation results in electron transfer along the A-side branch of cofactors at pH 7.2, but at pH 9.5, a long-lived state involving B-side cofactors is produced (Haffa et al. (2004) J Phys Chem B 108, 4-7). Using electron paramagnetic resonance spectroscopy, the mutants with alterations of each of the three individual residues and a double mutant, with changes at L168 and L170, were found to have increased linewidths of 10.1-11.0 G compared to the linewidth of 9.6 G for wild type. The Special TRIPLE spectra were pH dependent, and at pH 8, the introduction of aspartate at L170 increased the spin density ratio, rho (L)/rho (M), to 6.1 while an aspartate at the symmetry related position, M199, decreased the ratio to 0.7 compared to the value of 2.1 for wild type. These results indicate that the energy of the two halves of P changes by about 100 meV due to the mutations and are consistent with the interpretation that electrostatic interactions involving these amino acid residues contribute to the switch in pathway of electron transfer.     

Spectroscopic characterization of quinone-site mutants of the bacterial photosynthetic reaction center

Site-specific mutations in the quinone binding sites of the photosynthetic reaction center (RC) protein complexes of Rhodobacter (R.) capsulatus caused pronounced effects on sequential electron transfer. Conserved residues that break the twofold symmetry in this region of the RC – M246Ala and M247Ala in the QA binding pocket, and L212Glu and L213Asp in the QB binding pocket – were targeted. We constructed a QB-site mutant, L212Glu-L213Asp Ala-Ala, and a QA-site mutant, M246Ala–M247Ala Glu-Asp, to partially balance the differences in charge distribution normally found between the two quinone binding sites. In addition, two photocompetent revertants were isolated from the photosynthetically-incompetent M246Glu-M247Asp mutant: M246Ala–M247Asp and M246Gly–M247Asp. Sequential electron transfer was investigated by continuous light excitation and time-resolved electron paramagnetic resonance (EPR), and time-resolved optical techniques. Several lines of EPR evidence suggested that the forward electron transfer rate to QA, kQ, was slowed in those strains containing altered QA sites. The slower rates of secondary electron transfer were confirmed by time-resolved optical results with the M246Glu-M247Asp mutations in the QA site resulting in a dramatically lowered secondary electron transfer efficiency [kQ < (2 ns)-1] in comparison with either the native R. capsulatus RC or the QB site mutant [kQ (200 ps)-1]. Secondary electron transfer in the two revertants was intermediate between that of the native RC and the QA mutant. The P+ QA- PQA charge recombination rates were also changed in the strains that carried altered QA sites. We show that local mutations in the QA site, presumably through local electrostatic changes, significantly alter binding and electron transfer properties of QA.     

Electrostatic calculations of amino acid titration and electron transfer, Q-AQB-->QAQ-B, in the reaction center.

The titration of amino acids and the energetics of electron transfer from the primary electron acceptor (QA) to the secondary electron acceptor (QB) in the photosynthetic reaction center of Rhodobacter sphaeroides are calculated using a continuum electrostatic model. Strong electrostatic interactions between titrating sites give rise to complex titration curves. Glu L212 is calculated to have an anomalously broad titration curve, which explains the seemingly contradictory experimental results concerning its pKa. The electrostatic field following electron transfer shifts the average protonation of amino acids near the quinones. The pH dependence of the free energy between Q-AQB and QAQ-B calculated from these shifts is in good agreement with experiment. However, the calculated absolute free energy difference is in severe disagreement (by approximately 230 meV) with the observed experimental value, i.e., electron transfer from Q-A to QB is calculated to be unfavorable. The large stabilization energy of the Q-A state arises from the predominantly positively charged residues in the vicinity of QA in contrast to the predominantly negatively charged residues near QB. The discrepancy between calculated and experimental values for delta G(Q-AQB-->QAQ-B) points to limitations of the continuum electrostatic model. Inclusion of other contributions to the energetics (e.g., protein motion following quinone reduction) that may improve the agreement between theory and experiment are discussed.     

The fourth epidermal growth factor-like domain of thrombomodulin interacts with the basic exosite of protein C

Thrombomodulin (TM) functions as a cofactor to enhance the rate of protein C activation by thrombin approximately 1000-fold. The molecular mechanism by which TM improves the catalytic efficiency of thrombin toward protein C is not known. Molecular modeling of the protein C activation based on the crystal structure of thrombin in complex with the epidermal growth factor-like domains 4, 5, and 6 of TM (TM456) predicts that the binding of TM56 to exosite 1 of thrombin positions TM4 so that a negatively charged region on this domain juxtaposes a positively charged region of protein C. It has been hypothesized that electrostatic interactions between these oppositely charged residues of TM4 and protein C facilitate a proper docking of the substrate into the catalytic pocket of thrombin. To test this hypothesis, we have constructed several mutants of TM456 and protein C in which charges of the putative interacting residues on both TM4 (Asp/Glu) and protein C (Lys/Arg) have been reversed. Results of TM-dependent protein C activation studies by such a compensatory mutagenesis approach support the molecular model that TM4 interacts with the basic exosite of protein C.     

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

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