Cardioprotection by resveratrol: A human clinical trial in patients with stable coronary artery disease

Several beneficial effects of resveratrol (RES), a natural antioxidant present in red wine have already been described. The aim of our study was to investigate if RES had a clinically measurable cardioprotective effect in patients after myocardial infarction. In this double-blind, placebo controlled trial 40 post-infarction Caucasian patients were randomized into two groups. One group received 10 mg RES capsule daily for 3 months. Systolic and diastolic left ventricular function, flow-mediated vasodilation (FMD), several laboratory and hemorheological parameters were measured before and after the treatment. Left ventricular ejection fraction showed an increasing tendency (ns) by RES treatment. However, left ventricular diastolic function was improved significantly (p < 0.01) by RES. A significant improvement in endothelial function measured by FMD was also observed (p < 0.05). Low-density lipoprotein (LDL) level significantly decreased (p < 0.05) in the RES treated group. Red blood cell deformability decreased and platelet aggregation increased significantly in the placebo group (p < 0.05), while resveratrol treatment has prevented these unfavourable changes. Concerning other measured parameters no significant changes were observed neither in placebo nor in RES group. Our results show that resveratrol improved left ventricle diastolic function, endothelial function, lowered LDL-cholesterol level and protected against unfavourable hemorheological changes measured in patients with coronary artery disease (CAD).

Isotope effect on electron transfer in reaction centers from Rhodopseudomonas sphaeroides

Previous ENDOR studies on reaction centers from Rhodopseudomonas sphaeroides have shown the presence of two hydrogen-bonded protons associated with the primary, ubiquinone, acceptor Q(A). These protons exchange with deuterons from solvent (2)H(2)O. The effect of this deuterium substitution on the charge-recombination kinetics (BChl)(2) (+)Q(A) (-) --> (BChl)(2)Q(A) has been studied with a sensitive kinetic difference technique. The electron-transfer rate was found to increase with deuterium exchange up to a maximum Deltak/k of 5.7 +/- 0.3%. The change in rate was found to have an exchange time of 2 hr, which matched the disappearance of the ENDOR lines due to the exchangeable protons. These results indicate that these protons play a role in the vibronic coupling associated with electron transfer. A simple model for the isotope effect on electron transfer predicts a maximum rate increase of 20%, which is consistent with the experimental results.

Primary photochemistry of iron-depleted and zinc-reconstituted reaction centers from Rhodopseudomonas sphaeroides

The primary photochemistry of Fe-depleted and Zn-reconstituted reaction centers from Rhodopseudomonas sphaeroides R-26.1 was studied by transient absorption spectroscopy and compared with native, Fe(2+)-containing reaction centers. Excitation of metal-free reaction centers with 30-ps flashes produced the initial charge-separated state P(+)I(-) (P(+)BPh(-), where P is the primary donor and BPh is bacteriopheophytin) with a yield and visible/near-infrared absorption difference spectrum indistinguishable from that observed in native reaction centers. However, the lifetime of P(+)I(-) was found to increase approximately 20-fold to 4.2 +/- 0.3 ns (compared to 205 ps in native reaction centers), and the yield of formation of the subsequent state P(+)Q(A) (-) (Q(A) is the primary quinone acceptor) was reduced to 47 +/- 5% (compared to essentially 100% in native reaction centers). The remaining 53% of the metal-free reaction centers were found to undergo charge recombination during the P(+)I(-) lifetime to yield both the ground state (28 +/- 5%) and the triplet state P(R) (25 +/- 5%). Reconstitution of Fe-depleted reaction centers with Zn(2+) restored the "native" photochemistry. Possible mechanisms responsible for the reduced decay rate of P(+)I(-) in metal-free reaction centers are discussed.

Primary structure of the M subunit of the reaction center from Rhodopseudomonas sphaeroides

The reaction center is a membrane-bound bacteriochlorophyll-protein complex that mediates the primary photochemical events in the photosynthetic bacterium Rhodopseudomonas sphaeroides. The previously determined amino-terminal sequences of the three subunits of the reaction center protein were used to design synthetic mixed oligonucleotide probes for the structural genes encoding the subunits. One of these probes was used to isolate and clone a fragment of DNA from R. sphaeroides that contained the gene encoding the M subunit. The nucleotide sequence of this gene was determined by the dideoxy method. In addition, a number of tryptic and chymotryptic peptides from the M protein were isolated and subjected to sequence analysis, and the sequence of the carboxyl terminus was determined. Together with the amino-terminal sequence, the data establish the primary structure of the M protein. The distribution of hydrophobic residues in the amino acid sequence suggests the presence of five membrane-spanning segments. A significant homology was found between the amino acid sequence of the M subunit and a thylakoid membrane protein (M(r) 32,000) from spinach that has been implicated in herbicide and quinone binding.

Functional reconstitution of photosynthetic reaction centers in planar lipid bilayers

Planar lipid bilayers containing reaction centers from Rhodopseudomonas sphaeroides R-26 were formed by apposing two reaction center-lipid monolayers formed from a reaction center-lipid complex in hexane. Secondary donors (cytochrome c) and acceptors (ubiquinone-0) were added on opposite sides of the membrane. Upon illumination, this system generated transient and steady-state voltages and currents. The wavelength dependence of the photoresponse matched the absorption spectrum of reaction centers. A simple model based on the transfer of charges across the membrane that explains the salient features of the photoresponse is presented.

Fluctuation spectroscopy: determination of chemical reaction kinetics from the frequency spectrum of fluctuations

The kinetic parameters of a chemical reaction were obtained from analysis of the frequency spectrum of the fluctuations (i.e., "noise") in the concentrations of the reactants. In "fluctuation spectroscopy," no external perturbation is applied and the system remains in macroscopic chemical equilibrium during the experiment. Results obtained by this method for the dissociation reaction of beryllium sulfate agree well with those obtained by relaxation methods in which the approach to equilibrium is analyzed. Other noise sources not originating from a chemical reaction were observed and analyzed. The most prominent of these arose from the flow of an electrolyte through a capillary. The method of fluctuation spectroscopy should be applicable to problems of physical, chemical, and biological interest.

The effect of carotid stenting on rheological parameters, free radical production and platelet aggregation

INTRODUCTION

Carotid artery stenting has become a possible treatment of significant carotid stenosis. The risk of stent occlusion and restenosis might be increased by abnormal rheological conditions amplified platelet aggregation and free radical production during the operation.

AIMS

The aim of our study was to assess the changes in hemorheological parameters, platelet aggregation, and catalase activity after endovascular treatment of carotid stenosis.

METHODS

18 patients (11 men, ages 68 +/- 9 years and 7 women, ages 62 +/- 8 years) suffering from significant carotid stenosis and treated with carotid endovascular intervention were examined. Alteration in hemorheological parameters as well as epinephrine-, ADP-, and collagen-induced platelet aggregation were evaluated. Antioxidant reserve was characterized by the determination of catalase activity. The measurements were carried out directly before and after the procedure and 1, 2, 5 days and 1 month following the intervention. Preceding the operation the patients were administered a maximum dose (300 mg) of clopidogrel.

RESULTS

Hematocrit, plasma fibrinogen concentration (PFC) and whole blood-, and plasma viscosity values (WBV and PV) significantly decreased immediately after stenting (p<0.001). By the fifth day following the intervention the PFC, WBV, PV, red blood cell (RBC) aggregation and ADP-induced platelet aggregation significantly increased (p<0.0001) compared to values measured postprocedurally. At 1 month follow-up these parameters, except whole blood viscosity, decreased significantly compared to measurements made on the 5th day. On the other hand, catalase activity showed significant elevation by the end of the first month.

CONCLUSION

Hemorheological parameters and platelet aggregation showed specific changes following carotid stenting. Abnormal changes of the rheological conditions and increasing platelet activation are the most pronounced in the first week following stenting, which may lead to early stent occlusion. Oxidative stress production returned to baseline levels only by the end of the first month.

Protein-cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: II. Geometry of the hydrogen bonds to the primary quinone formula by 1H and 2H ENDOR spectroscopy

The geometry of the hydrogen bonds to the two carbonyl oxygens of the semiquinone Q(A)(. -) in the reaction center (RC) from the photosynthetic purple bacterium Rhodobacter sphaeroides R-26 were determined by fitting a spin Hamiltonian to the data derived from (1)H and (2)H ENDOR spectroscopies at 35 GHz and 80 K. The experiments were performed on RCs in which the native Fe(2+) (high spin) was replaced by diamagnetic Zn(2+) to prevent spectral line broadening of the Q(A)(. -) due to magnetic coupling with the iron. The principal components of the hyperfine coupling and nuclear quadrupolar coupling tensors of the hydrogen-bonded protons (deuterons) and their principal directions with respect to the quinone axes were obtained by spectral simulations of ENDOR spectra at different magnetic fields on frozen solutions of deuterated Q(A)(. -) in H(2)O buffer and protonated Q(A)(. -) in D(2)O buffer. Hydrogen-bond lengths were obtained from the nuclear quadrupolar couplings. The two hydrogen bonds were found to be nonequivalent, having different directions and different bond lengths. The H-bond lengths r(OH) are 1.73 +/- 0.03 Angstrom and 1.60 +/- 0.04 Angstrom, from the carbonyl oxygens O(1) and O(4) to the NH group of Ala M260 and the imidazole nitrogen N(delta) of His M219, respectively. The asymmetric hydrogen bonds of Q(A)(. -) affect the spin density distribution in the quinone radical and its electronic structure. It is proposed that the H-bonds play an important role in defining the physical properties of the primary quinone, which affect the electron transfer processes in the RC.

Protein-cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: I. Identification of the ENDOR lines associated with the hydrogen bonds to the primary quinone QA*-

Hydrogen bonds are important in determining the structure and function of biomolecules. Of particular interest are hydrogen bonds to quinones, which play an important role in the bioenergetics of respiration and photosynthesis. In this work we investigated the hydrogen bonds to the two carbonyl oxygens of the semiquinone QA*- in the well-characterized reaction center from the photosynthetic bacterium Rhodobacter sphaeroides R-26. We used electron paramagnetic resonance and electron nuclear double resonance techniques at 35 GHz at a temperature of 80 K. The goal of this study was to identify and assign sets of 1H-ENDOR lines to protons hydrogen bonded to each of the two oxygens. This was accomplished by preferentially exchanging the hydrogen bond on one of the oxygens with deuterium while concomitantly monitoring the changes in the amplitudes of the 1H-ENDOR lines. The preferential deuteration of one of the oxygens was made possible by the different 1H --> 2H exchange times of the protons bonded to the two oxygens. The assignment of the 1H-ENDOR lines sets the stage for the determination of the geometries of the H-bonds by a detailed field selection ENDOR study to be presented in a future article.

Gender differences in hemorheological parameters of coronary artery disease patients

Plasma fibrinogen concentration, plasma and whole blood viscosity (WBV) are independent risk factors of coronary artery disease (CAD). Fibrinogen seems to be a relatively stronger risk factor for women than for men, but men are more endangered by higher hematocrit (Hct) and WBV than women are. We have previously reported that a theoretically optimal Hct value can be determined using Hct/WBV ratio in healthy subjects, hyperlipidemic and Raynaud's disease patients. Our aim was to examine whether Hct/WBV ratio is differently correlated with Hct in men and women with proven CAD. In a retrospective study we analysed the hemorheological data of 162 CAD outpatients (107 men and 55 women). Coronary angiography, echocardiography and impedance cardiography were performed. Hemorheological parameters (Hct, fibrinogen level, plasma viscosity, WBV), blood picture, serum lipid concentrations were determined and Hct/WBV ratio was calculated. Mean ages of male and female patients were similar (54.9 and 55.4 years, respectively), but men had significantly higher coronary angiography score than women. Mean left ventricular ejection fraction, stroke volume index and cardiac index showed no significant differences in men and women. Similarly, lipid concentrations, fibrinogen levels and plasma viscosities demonstrated no statistical differences. However, Hct, WBV and Hct/WBV ratios were significantly higher in male than in female patients (p < 0.00001; p < 0.00001 and p < 0.005, respectively). The most striking gender difference was found in the correlation between Hct/WBV ratio and cardiac index. Men older than 56 years showed negative, women positive correlation (r = -0.485, p = 0.01; r = 0.468, p = 0.006, respectively). This study demonstrates that Hct/WBV ratio as a rheological oxygen carrying capacity parameter is positively correlated with the cardiac index as it can be expected. However, the correlation is negative in elder men indicating an unhealthy relation between hemodynamic and hemorheologic parameters.

Quinone (QB) Reduction by B-Branch Electron Transfer in Mutant Bacterial Reaction Centers from Rhodobacter sphaeroides: Quantum Efficiency and X-ray Structure,

The photosynthetic reaction center (RC) from purple bacteria converts light into chemical energy. Although the RC shows two nearly structurally symmetric branches, A and B, light-induced electron transfer in the native RC occurs almost exclusively along the A-branch to a primary quinone electron acceptor Q(A). Subsequent electron and proton transfer to a mobile quinone molecule Q(B) converts it to a quinol, Q(B)H(2). We report the construction and characterization of a series of mutants in Rhodobacter sphaeroides designed to reduce Q(B) via the B-branch. The quantum efficiency to Q(B) via the B-branch Phi(B) ranged from 0.4% in an RC containing the single mutation Ala-M260 --> Trp to 5% in a quintuple mutant which includes in addition three mutations to inhibit transfer along the A-branch (Gly-M203 --> Asp, Tyr-M210 --> Phe, Leu-M214 --> His) and one to promote transfer along the B-branch (Phe-L181 --> Tyr). Comparing the value of 0.4% for Phi(B) obtained in the AW(M260) mutant, which lacks Q(A), to the 100% quantum efficiency for Phi(A) along the A-branch in the native RC, we obtain a ratio for A-branch to B-branch electron transfer of 250:1. We determined the structure of the most effective (quintuple) mutant RC at 2.25 A (R-factor = 19.6%). The Q(A) site did not contain a quinone but was occupied by the side chain of Trp-M260 and a Cl(-). In this structure a nonfunctional quinone was found to occupy a new site near M258 and M268. The implications of this work to trap intermediate states are discussed.

X-Ray Structure Determination of Three Mutants of the Bacterial Photosynthetic Reaction Centers from Rb. sphaeroides

In the photosynthetic reaction center (RC) from Rhodobacter sphaeroides, the reduction of a bound quinone molecule Q(B) is coupled with proton uptake. When Asp-L213 is replaced by Asn, proton transfer is inhibited. Proton transfer was restored by two second-site revertant mutations, Arg-M233-->Cys and Arg-H177-->His. Kinetic effects of Cd(2+) on proton transfer showed that the entry point in revertant RCs to be the same as in the native RC. The structures of the parental and two revertant RCs were determined at resolutions of 2.10, 1.80, and 2.75 A. From the structures, we were able to delineate alternate proton transfer pathways in the revertants. The main changes occur near Glu-H173, which allow it to substitute for the missing Asp-L213. The electrostatic changes near Glu-H173 cause it to be a good proton donor and acceptor, and the structural changes create a cavity which accommodates water molecules that connect Glu-H173 to other proton transfer components.

Proton transfer pathways and mechanism in bacterial reaction centers

The focus of this minireview is to discuss the state of knowledge of the pathways and rates of proton transfer in the bacterial reaction center (RC) from Rhodobacter sphaeroides. Protons involved in the light driven catalytic reduction of a quinone molecule QB to quinol QBH2 travel from the aqueous solution through well defined proton transfer pathways to the oxygen atoms of the quinone. Three main topics are discussed: (1) the pathways for proton transfer involving the residues: His-H126, His-H128, Asp-L210, Asp-M17, Asp-L213, Ser-L223 and Glu-L212, which were determined by a variety of methods including the use of proton uptake inhibiting metal ions (e.g. Zn2+ and Cd2+); (2) the rate constants for proton transfer, obtained from a 'chemical rescue' study was determined to be 2 x 10(5) s(-1) and 2 x 10(4) s(-1) for the proton uptake to Glu-L212 and QB-*, respectively; (3) structural studies of altered proton transfer pathways in revertant RCs that lack the key amino acid Asp-L213 show a series of structural changes that propagate toward L213 potentially allowing Glu-H173 to participate in the proton transfer processes.

Mechanism of proton transfer inhibition by Cd(2+) binding to bacterial reaction centers: determination of the pK(A) of functionally important histidine residues

The bacterial photosynthetic reaction center (RC) uses light energy to catalyze the reduction of a bound quinone molecule Q(B) to quinol Q(B)H(2). In RCs from Rhodobacter sphaeroides the protons involved in this process come from the cytoplasm and travel through pathways that involve His-H126 and His-H128 located near the proton entry point. In this study, we measured the pH dependence from 4.5 to 8.5 of the binding of the proton transfer inhibitor Cd(2+), which ligates to these surface His in the RC and inhibits proton-coupled electron transfer. At pH <6, the negative slope of the logarithm of the dissociation constant, K(D), versus pH approaches 2, indicating that, upon binding of Cd(2+), two protons are displaced; i.e., the binding is electrostatically compensated. At pH >7, K(D) becomes essentially independent of pH. A theoretical fit to the data over the entire pH range required two protons with pK(A) values of 6.8 and 6.3 (+/-0.5). To assess the contribution of His-H126 and His-H128 to the observed pH dependence, K(D) was measured in mutant RCs that lack the imidazole group of His-H126 or His-H128 (His --> Ala). In both mutant RCs, K(D) was approximately pH independent, showing that Cd(2+) does not displace protons upon binding in the mutant RCs, in contrast to the native RC in which His-H126 and His-H128 are the predominant contributors to the observed pH dependence of K(D). Thus, Cd(2+) inhibits RC function by binding to functionally important histidines.

Determination of proton transfer rates by chemical rescue: application to bacterial reaction centers

The bacterial reaction center (RC) converts light into chemical energy through the reduction of an internal quinone molecule Q(B) to Q(B)H(2). In the native RC, proton transfer is coupled to electron transfer and is not rate-controlling. Consequently, proton transfer is not directly observable, and its rate was unknown. In this work, we present a method for making proton transfer rate-controlling, which enabled us to determine its rate. The imidazole groups of the His-H126 and His-H128 proton donors, located at the entrance of the transfer pathways, were removed by site-directed mutagenesis (His --> Ala). This resulted in a reduction in the observed proton-coupled electron transfer rate [(Q(A)(-)(*)Q(B))Glu(-) + H(+) --> (Q(A)Q(B)(-)(*))GluH], which became rate-controlled by proton uptake to Glu-L212 [Adelroth, P., et al. (2001) Biochemistry 40, 14538-14546]. The proton uptake rate was enhanced (rescued) in a controlled fashion by the addition of imidazole or other amine-containing acids. From the dependence of the observed rate on acid concentration, an apparent second-order rate constant k((2)) for the "rescue" of the rate was determined. k((2)) is a function of the proton transfer rate and the binding of the acid. The dependence of k((2)) on the acid pK(a) (i.e., the proton driving force) was measured over 9 pK(a) units, resulting in a Brönsted plot that was characteristic of general acid catalysis. The results were fitted to a model that includes the binding (facilitated by electrostatic attraction) of the cationic acid to the RC surface, proton transfer to an intermediate proton acceptor group, and subsequent proton transfer to Glu-L212. A proton transfer rate constant of approximately 10(5) s(-)(1) was determined for transfer from the bound imidazole group to Glu-L212 (over a distance of approximately 20 A). The same method was used to determine a proton transfer rate constant of 2 x 10(4) s(-)(1) for transfer to Q(B)(-)(*). The relatively fast proton transfer rates are explained by the presence of an intermediate acceptor group that breaks the process into sequential proton transfer steps over shorter distances. This study illustrates an approach that could be generally applied to obtain information about the individual rates and energies for proton transfer processes, as well as the pK(a)s of transfer components, in a variety of proton translocating systems.

Spin-lattice relaxation of coupled metal-radical spin-dimers in proteins: application to Fe(2+)-cofactor (Q(A)(-.), Q(B)(-.), phi(-.)) dimers in reaction centers from photosynthetic bacteria

The spin-lattice relaxation times (T(1)) for the reduced quinone acceptors Q(A)(-.) and Q(B)(-.), and the intermediate pheophytin acceptor phi(-.), were measured in native photosynthetic reaction centers (RC) containing a high spin Fe(2+) (S = 2) and in RCs in which Fe(2+) was replaced by diamagnetic Zn(2+). From these data, the contribution of the Fe(2+) to the spin-lattice relaxation of the cofactors was determined. To relate the spin-lattice relaxation rate to the spin-spin interaction between the Fe(2+) and the cofactors, we developed a spin-dimer model that takes into account the zero field splitting and the rhombicity of the Fe(2+) ion. The relaxation mechanism of the spin-dimer involves a two-phonon process that couples the fast relaxing Fe(2+) spin to the cofactor spin. The process is analogous to the one proposed by R. Orbach (Proc. R. Soc. A. (Lond.). 264:458-484) for rare earth ions. The spin-spin interactions are, in general, composed of exchange and dipolar contributions. For the spin dimers studied in this work the exchange interaction, J(o), is predominant. The values of J(o) for Q(A)(-.)Fe(2+), Q(B)(-.)Fe(2+), and phi(-.)Fe(2+) were determined to be (in kelvin) -0.58, -0.92, and -1.3 x 10(-3), respectively. The |J(o)| of the various cofactors (obtained in this work and those of others) could be fitted with the relation exp(-beta(J)d), where d is the distance between cofactor spins and beta(J) had a value of (0.66-0.86) A(-1). The relation between J(o) and the matrix element |V(ij)|(2) involved in electron transfer rates is discussed.

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.

Identification of the proton pathway in bacterial reaction centers: decrease of proton transfer rate by mutation of surface histidines at H126 and H128 and chemical rescue by imidazole identifies the initial proton donors

The pathway for proton transfer to Q(B) was studied in the reaction center (RC) from Rhodobacter sphaeroides. The binding of Zn(2+) or Cd(2+) to the RC surface at His-H126, His-H128, and Asp-H124 inhibits the rate of proton transfer to Q(B), suggesting that the His may be important for proton transfer [Paddock, M. L., Graige, M. S., Feher, G. and Okamura, M. Y. (1999) Proc. Natl. Acad. Sci. U.S.A. 96, 6183-6188]. To assess directly the role of the histidines, mutant RCs were constructed in which either one or both His were replaced with Ala. In the single His mutant RCs, no significant effects were observed. In contrast, in the double mutant RC at pH 8.5, the observed rates of proton uptake associated with both the first and the second proton-coupled electron-transfer reactions k(AB)(()(1)()) [Q(A)(-)(*)Q(B)-Glu(-) + H(+) --> Q(A)(-)(*)Q(B)-GluH --> Q(A)Q(B)(-)(*)-GluH] and k(AB)(()(2)()) [Q(A)(-)(*)Q(B)(-)(*) + H(+) --> Q(A)(-)(*)(Q(B)H)(*) --> Q(A)(Q(B)H)(-)], were found to be slowed by factors of approximately 10 and approximately 4, respectively. Evidence that the observed changes in the double mutant RC are due to a reduction in the proton-transfer rate constants are provided by the observations: (i) k(AB)(1) at pH approximately pK(a) of GluH became biphasic, indicating that proton transfer is slower than electron transfer and (ii) k(AB)(2) became independent of the driving force for electron transfer, indicating that proton transfer is the rate-limiting step. These changes were overcome by the addition of exogenous imidazole which acts as a proton donor in place of the imidazole groups of His that were removed in the double mutant RC. Thus, we conclude that His-H126 and His-H128 facilitate proton transfer into the RC, acting as RC-bound proton donors at the entrance of the proton-transfer pathways.

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

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

Identification of the proton pathway in bacterial reaction centers: cooperation between Asp-M17 and Asp-L210 facilitates proton transfer to the secondary quinone (QB)

The reaction center (RC) from Rhodobacter sphaeroides uses light energy to reduce and protonate a quinone molecule, Q(B) (the secondary quinone electron acceptor), to form quinol, Q(B)H2. Asp-L210 and Asp-M17 have been proposed to be components of the pathway for proton transfer [Axelrod, H. L., Abresch, E. C., Paddock, M. L., Okamura, M. Y., and Feher, G. (2000) Proc. Natl. Acad. Sci. U.S.A. 97, 1542-1547]. To test the importance of these residues for efficient proton transfer, the rates of the proton-coupled electron-transfer reaction k(AB)(2) (Q(A-*)Q(B-*) + H+ <==>Q(A-*)Q(B)H* --> Q(A)Q(B)H-) and its associated proton uptake were measured in native and mutant RCs, lacking one or both Asp residues. In the double mutant RCs, the k(AB)(2) reaction and its associated proton uptake were approximately 300-fold slower than in native RCs (pH 8). In contrast, single mutant RCs displayed reaction rates that were < or =3-fold slower than native (pH 8). In addition, the rate-limiting step of k(AB)(2) was changed from electron transfer (native and single mutants) to proton transfer (double mutant) as shown from the lack of a dependence of the observed rate on the driving force for electron transfer in the double mutant RCs compared to the native or single mutants. This implies that the rate of the proton-transfer step was reduced (> or =10(3)-fold) upon replacement of both Asp-L210 and Asp-M17 with Asn. Similar, but less drastic, differences were observed for k(AB)(1), which at pH > or =8 is coupled to the protonation of Glu-L212 [(Q(A-*)Q(B))-Glu- + H+ --> (Q(A)Q(B-*)-GluH]. These results show that the pathway for proton transfer from solution to reduced Q(B) involves both Asp-L210 and Asp-M17, which provide parallel branches to the proton-transfer pathway and through their electrostatic interaction have a cooperative effect on the proton-transfer rate. A possible mechanism for the cooperativity is discussed.

Identification of the proton pathway in bacterial reaction centers: both protons associated with reduction of QB to QBH2 share a common entry point

The reaction center from Rhodobacter sphaeroides uses light energy for the reduction and protonation of a quinone molecule, Q(B). This process involves the transfer of two protons from the aqueous solution to the protein-bound Q(B) molecule. The second proton, H(+)(2), is supplied to Q(B) by Glu-L212, an internal residue protonated in response to formation of Q(A)(-) and Q(B)(-). In this work, the pathway for H(+)(2) to Glu-L212 was studied by measuring the effects of divalent metal ion binding on the protonation of Glu-L212, which was assayed by two types of processes. One was proton uptake from solution after the one-electron reduction of Q(A) (DQ(A)-->D(+)Q(A)(-)) and Q(B) (DQ(B)-->D(+)Q(B)(-)), studied by using pH-sensitive dyes. The other was the electron transfer k(AB)((1)) (Q(A)(-)Q(B)-->Q(A)Q(B)(-)). At pH 8.5, binding of Zn(2+), Cd(2+), or Ni(2+) reduced the rates of proton uptake upon Q(A)(-) and Q(B)(-) formation as well as k(AB)((1)) by approximately an order of magnitude, resulting in similar final values, indicating that there is a common rate-limiting step. Because D(+)Q(A)(-) is formed 10(5)-fold faster than the induced proton uptake, the observed rate decrease must be caused by an inhibition of the proton transfer. The Glu-L212-->Gln mutant reaction centers displayed greatly reduced amplitudes of proton uptake and exhibited no changes in rates of proton uptake or electron transfer upon Zn(2+) binding. Therefore, metal binding specifically decreased the rate of proton transfer to Glu-L212, because the observed rates were decreased only when proton uptake by Glu-L212 was required. The entry point for the second proton H(+)(2) was thus identified to be the same as for the first proton H(+)(1), close to the metal binding region Asp-H124, His-H126, and His-H128.

Proton and electron transfer in bacterial reaction centers

The bacterial reaction center couples light-induced electron transfer to proton pumping across the membrane by reactions of a quinone molecule Q(B) that binds two electrons and two protons at the active site. This article reviews recent experimental work on the mechanism of the proton-coupled electron transfer and the pathways for proton transfer to the Q(B) site. The mechanism of the first electron transfer, k((1))(AB), Q(-)(A)Q(B)-->Q(A)Q(-)(B), was shown to be rate limited by conformational gating. The mechanism of the second electron transfer, k((2))(AB), was shown to involve rapid reversible proton transfer to the semiquinone followed by rate-limiting electron transfer, H(+)+Q(-)(A)Q(-)(B) ifQ(-)(A)Q(B)H-->Q(A)(Q(B)H)(-). The pathways for transfer of the first and second protons were elucidated by high-resolution X-ray crystallography as well as kinetic studies showing changes in the rate of proton transfer due to site directed mutations and metal ion binding.

Identification of the proton pathway in bacterial reaction centers: replacement of Asp-M17 and Asp-L210 with asn reduces the proton transfer rate in the presence of Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the reduction and protonation of a bound quinone molecule Q(B) (the secondary quinone electron acceptor). We investigated the proton transfer pathway by measuring the proton-coupled electron transfer, k(AB)((2)) [Q(A)Q(B) + H(+) --> Q(A)(Q(B)H)(-)] in native and mutant RCs in the absence and presence of Cd(2+). Previous work has shown that the binding of Cd(2+) decreases k(AB)((2)) in native RCs approximately 100-fold. The preceding paper shows that bound Cd(2+) binds to Asp-H124, His-H126, and His-H128. This region represents the entry point for protons. In this work we investigated the proton transfer pathway connecting the entry point with Q(B) by searching for mutations that greatly affect k(AB)((2)) ( greater, similar10-fold) in the presence of Cd(2+), where k(AB)((2)) is limited by the proton transfer rate (k(H)). Upon mutation of Asp-L210 or Asp-M17 to Asn, k(H) decreased from approximately 60 s(-1) to approximately 7 s(-1), which shows the important role that Asp-L210 and Asp-M17 play in the proton transfer chain. By comparing the rate of proton transfer in the mutants (k(H) approximately 7 s(-1)) with that in native RCs in the absence of Cd(2+) (k(H) >/= 10(4) s(-1)), we conclude that alternate proton transfer pathways, which have been postulated, are at least 10(3)-fold less effective.

Determination of the binding sites of the proton transfer inhibitors Cd2+ and Zn2+ in bacterial reaction centers

The reaction center (RC) from Rhodobacter sphaeroides couples light-driven electron transfer to protonation of a bound quinone acceptor molecule, Q(B), within the RC. The binding of Cd(2+) or Zn(2+) has been previously shown to inhibit the rate of reduction and protonation of Q(B). We report here on the metal binding site, determined by x-ray diffraction at 2.5-A resolution, obtained from RC crystals that were soaked in the presence of the metal. The structures were refined to R factors of 23% and 24% for the Cd(2+) and Zn(2+) complexes, respectively. Both metals bind to the same location, coordinating to Asp-H124, His-H126, and His-H128. The rate of electron transfer from Q(A)(-) to Q(B) was measured in the Cd(2+)-soaked crystal and found to be the same as in solution in the presence of Cd(2+). In addition to the changes in the kinetics, a structural effect of Cd(2+) on Glu-H173 was observed. This residue was well resolved in the x-ray structure-i.e., ordered-with Cd(2+) bound to the RC, in contrast to its disordered state in the absence of Cd(2+), which suggests that the mobility of Glu-H173 plays an important role in the rate of reduction of Q(B). The position of the Cd(2+) and Zn(2+) localizes the proton entry into the RC near Asp-H124, His-H126, and His-H128. Based on the location of the metal, likely pathways of proton transfer from the aqueous surface to Q(B) are proposed.

Observation of the protonated semiquinone intermediate in isolated reaction centers from Rhodobacter sphaeroides: implications for the mechanism of electron and proton transfer in proteins

A proton-activated electron transfer (PAET) mechanism, involving a protonated semiquinone intermediate state, had been proposed for the electron-transfer reaction k(2)AB [Q(A)(-)(*)Q(B)(-)(*) + H(+) <--> Q(A)(-)(*)(Q(B)H)(*) --> Q(A)(Q(B)H)(-)] in reaction centers (RCs) from Rhodobacter sphaeroides [Graige, M. S., Paddock, M. L., Bruce, M. L., Feher, G., and Okamura, M. Y. (1996) J. Am. Chem. Soc. 118, 9005-9016]. Confirmation of this mechanism by observing the protonated semiquinone (Q(B)H)(*) had not been possible, presumably because of its low pK(a). By replacing the native Q(10) in the Q(B) site with rhodoquinone (RQ), which has a higher pK(a), we were able to observe the (Q(B)H)(*) state. The pH dependence of the semiquinone optical spectrum gave a pK(a) = 7.3 +/- 0.2. At pH < pK(a), the observed rate for the reaction was constant and attributed to the intrinsic electron-transfer rate from Q(A)(-)(*) to the protonated semiquinone (i.e., k(2)AB = k(ET)(RQ) = 2 x 10(4) s(-)(1)). The rate decreased at pH > pK(a) as predicted by the PAET mechanism in which fast reversible proton transfer precedes rate-limiting electron transfer. Consequently, near pH 7, the proton-transfer rate k(H) >> 10(4) s(-)(1). Applying the two step mechanism to RCs containing native Q(10) and taking into account the change in redox potential, we find reasonable values for the fraction of (Q(B)H)(*) congruent with 0.1% (consistent with a pK(a)(Q(10)) of approximately 4.5) and k(ET)(Q(10)) congruent with 10(6) s(-)(1). These results confirm the PAET mechanism in RCs with RQ and give strong support that this mechanism is active in RCs with Q(10) as well.

Identification of the proton pathway in bacterial reaction centers: inhibition of proton transfer by binding of Zn2+ or Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the light induced two-electron, two-proton reduction of a bound quinone molecule QB (the secondary quinone acceptor). A unique pathway for proton transfer to the QB site had so far not been determined. To study the molecular basis for proton transfer, we investigated the effects of exogenous metal ion binding on the kinetics of the proton-assisted electron transfer kAB(2) (QA-*QB-* + H+ --> QA(QBH)-, where QA is the primary quinone acceptor). Zn2+ and Cd2+ bound stoichiometrically to the RC (KD /= 10(2)-fold) and has become the rate-limiting step. The lack of an effect of the metal binding on the charge recombination reaction D+*QAQB-* --> DQAQB suggests that the binding site is located far (>10 A) from QB. This hypothesis is confirmed by preliminary x-ray structure analysis. The large change in the rate of proton transfer caused by the stoichiometric binding of the metal ion shows that there is one dominant site of proton entry into the RC from which proton transfer to QB-* occurs.

Conformational gating of the electron transfer reaction QA-.QB --> QAQB-. in bacterial reaction centers of Rhodobacter sphaeroides determined by a driving force assay

The mechanism of the electron transfer reaction, QA-.QB --> QAQB-., was studied in isolated reaction centers from the photosynthetic bacterium Rhodobacter sphaeroides by replacing the native Q10 in the QA binding site with quinones having different redox potentials. These substitutions are expected to change the intrinsic electron transfer rate by changing the redox free energy (i.e., driving force) for electron transfer without affecting other events that may be associated with the electron transfer (e.g., protein dynamics or protonation). The electron transfer from QA-. to QB was measured by three independent methods: a functional assay involving cytochrome c2 to measure the rate of QA-. oxidation, optical kinetic spectroscopy to measure changes in semiquinone absorption, and kinetic near-IR spectroscopy to measure electrochromic shifts that occur in response to electron transfer. The results show that the rate of the observed electron transfer from QA-. to QB does not change as the redox free energy for electron transfer is varied over a range of 150 meV. The strong temperature dependence of the observed rate rules out the possibility that the reaction is activationless. We conclude, therefore, that the independence of the observed rate on the driving force for electron transfer is due to conformational gating, that is, the rate limiting step is a conformational change required before electron transfer. This change is proposed to be the movement, controlled kinetically either by protein dynamics or intermolecular interactions, of QB by approximately 5 A as observed in the x-ray studies of Stowell et al. [Stowell, M. H. B., McPhillips, T. M., Rees, D. C., Soltis, S. M., Abresch, E. & Feher, G. (1997) Science 276, 812-816].

Proton and electron transfer to the secondary quinone (QB) in bacterial reaction centers: the effect of changing the electrostatics in the vicinity of QB by interchanging asp and glu at the L212 and L213 sites

The bacterial reaction center (RC) plays a central role in photosynthetic energy conversion by facilitating the light induced double reduction and protonation of a bound quinone molecule, QB. Two carboxylic acid residues, Asp-L213 and Glu-L212, located near QB, were previously shown to be important for proton transfer to QB. In this work, the ability of Glu to substitute for Asp at L213 and Asp to substitute for Glu at L212 was tested by site-directed mutagenesis. Both single mutants and a double mutant in which Asp and Glu were exchanged between the two sites were constructed. The electron transfer rate constants kBD (D+QAQB- --> DQAQB), and kAB(2) (DQA-QB- + H+ --> DQA(QBH)-), that are known to be sensitive to the energy of the QB- state, were found to be altered by Asp/Glu substitutions. Both rates were fastest ( approximately 10-fold) in RCs with Asp at both sites, slowest with Glu at both sites ( approximately 50-fold) and relatively unchanged by the caboxylic acid exchange. These changes could be explained if Asp was predominantly ionized and Glu was predominantly protonated at both sites (pH 7.5). The charge recombination kBD suggests an observed approximately 5 pKa unit difference of Glu over Asp. Modeling of kBD by strong electrostatic interactions ( approximately 3-4 pKa units) among negatively charged acids and QB- indicated a lower intrinsic pKa for Asp compared to Glu at either site of approximately 2-3 units. The mechanism of the kAB(2) reaction was determined to be the same in all mutant RCs as for native RCs. A quantitative explanation of the effect of the electrostatic environment on kAB(2) was obtained using the two-step model proposed for native RCs [Graige, M. S., Paddock, M. L., Bruce, J. M., Feher, G., & Okamura, M. Y. (1996) J. Am. Chem. Soc. 118, 9005-9016] which involves fast protonation of the semiquinone followed by rate-limiting electron transfer. Using simple models for the quinone/quinol conversion rate, it is shown that the optimal electrostatic potential for the QB site is close to that found in native RCs.

Light-induced electrogenic events associated with proton uptake upon forming QB- in bacterial wild-type and mutant reaction centers

Light-induced voltage changes (electrogenic events) were measured in wild-type and site-directed mutants of reaction centers (RCs) from Rhodobacter sphaeroides oriented in a lipid monolayer adsorbed to a Teflon film. A rapid increase in voltage associated with charge separation was followed by a slower increase attributed to proton transfer from solution to protonatable amino-acid residues in the vicinity of the QB site. In native reaction centers the proton-transfer voltage had a pH-dependent amplitude with two peaks at pH 4.5 and pH 9.7, respectively. In the Glu-L212-->Gln RCs the high-pH peak was absent, whereas in the Asp-L213-->Asn RCs the low-pH peak was absent and the high-pH peak was shifted to lower pH by about 1.3 pH units. The amplitudes of the electrogenic phases as a function of pH follow approximately the measured proton uptake from solution (P.H. McPherson, M.Y. Okamura, G. Feher, Biochim. Biophys. Acta, vol. 934, 1988, pp. 348-368) and are ascribed to proton transfer to amino acid residues upon QB- formation. The peak around pH 9.7 is ascribed to proton uptake predominantly by Glu-L212 and the peak around pH 4.5 to proton uptake predominantly by Asp-L213 or a residue strongly interacting with Asp-L213.

Light-induced structural changes in photosynthetic reaction center: implications for mechanism of electron-proton transfer

High resolution x-ray diffraction data from crystals of the Rhodobacter sphaeroides photosynthetic reaction center (RC) have been collected at cryogenic temperature in the dark and under illumination, and the structures were refined at 2.2 and 2.6 angstrom resolution, respectively. In the charge-separated D+QAQB- state (where D is the primary electron donor (a bacteriochlorophyll dimer), and QA and QB are the primary and secondary quinone acceptors, respectively), QB- is located approximately 5 angstroms from the QB position in the charge-neutral (DQAQB) state, and has undergone a 180 degrees propeller twist around the isoprene chain. A model based on the difference between the two structures is proposed to explain the observed kinetics of electron transfer from QA-QB to QAQB- and the relative binding affinities of the different ubiquinone species in the QB pocket. In addition, several water channels (putative proton pathways) leading from the QB pocket to the surface of the RC were delineated, one of which leads directly to the membrane surface.

Protein crystallization

Crystallization is necessary to obtain the three-dimensional structure of proteins and nucleic acids; it often represents the bottleneck in structure determination. Our understanding of crystallization mechanisms is still incomplete. In this review, we emphasize fundamental aspects of the crystallization process. Protein-protein contacts in crystals are complex, involving a delicate balance of specific and nonspecific interactions. Depending on solution conditions, these interactions can lead to nucleation of crystals or to amorphous aggregation; this stage of crystallization has been successfully studied by light scattering. Post-nucleation crystal growth may proceed by mechanisms involving crystal defects or two-dimensional nucleation, as observed by atomic force and interference microscopy. Cessation of growth has been observed but remains incompletely understood. Impurities may play important roles during all stages of crystallization. Phase diagrams can guide optimization of conditions for nucleation and subsequent crystal growth; a theoretical understanding relating these to the intermolecular interactions is beginning to develop.

Co-crystallization and characterization of the photosynthetic reaction center-cytochrome c2 complex from Rhodobacter sphaeroides

The photosynthetic reaction center (RC) of Rhodobacter sphaeroides and cytochrome c2 (cyt c2), its physiological secondary electron donor, have been co-crystallized. The molar ratio of RC/cyt c2 was found by SDS-PAGE and optical absorbance changes in the co-crystals to be 4. The crystals diffracted X-rays to 3.5 angstroms. However, the resolution degraded during data collection. A data set, 82.5% complete, was collected to 4.5 angstroms. The crystals belong to the tetragonal space group P4(3)2(1)2, with unit cell dimensions of a = b = 142.7 angstroms and c = 254.8 angstroms. The positions of the RCs in the unit cell were determined by molecular replacement. A comparable search for the cyt c2 by this method was unsuccessful because of the small contribution of the cytochrome to the total scattering and because of its low occupancy. The cyt c2 was positioned manually into patches of difference electron density, adjacent to the periplasmic surface of the M polypeptide subunit of the RC. The difference electron density was not sufficient for precise positioning of the cyt c2, and its orientation was modeled by placing the exposed edge of the heme toward the primary donor of the reaction center D and by forming pairs for electrostatically interacting RC and cyt c2 amino acid residues. The RC-cyt c2 structure derived from the co-crystal data was supported by use of omit maps and structure refinement analyses. Cyt c2 reduces the photooxidized primary donor D+ in 0.9 +/- 0.1 micros in the co-crystals, which is the same as the fast electron transfer rate in vivo and in solution. This result provides strong evidence that the structure of the complex in the co-crystal is the same as in solution. Two additional methods were used to investigate the structure of the RC-cyt c2 complex: (i) Docking calculations based on interprotein electrostatic interactions identified possible binding positions of the cyt c2 on the RC. The cyt c2 position with the lowest electrostatic energy is very similar to that of the cyt c2 in the proposed co-crystal structure. (ii) Site-directed mutagenesis was used to modify two aspartic acid residues (M184 and L155) on the periplasmic surface of the RC. Cyt c2 binding affinity to these RCs and electron transfer rates to D+ in these RCs support the co-crystal structure of th RC-cyt c2 complex.

Pathway of proton transfer in bacterial reaction centers: further investigations on the role of Ser-L223 studied by site-directed mutagenesis

The role of Ser-L223 in proton transfer to reduced QB in the reaction center (RC) from Rhodobacter sphaeroides was studied by site-directed replacement of Ser with residues having different proton donor properties, e.g., the aliphatic residues Ala and Gly, the hydroxyl residue Thr, the amide residue Asn, the sulhydryl residue Cys, the imidazole residue His, and the carboxylic acid residue Asp. Compared to native reaction centers, RCs with Ala or Asn at L223 had greatly reduced (approximately 300-fold) proton-coupled electron transfer rates, kAB(2), associated with the second electron reduction of QB (QA(-)QB(-) + H+ --> QAQBH-). In contrast, RCs containing Thr, Asp, or Gly at L223 retained fast proton-coupled electron transfer rates. RCs with His or Cys at L223 did not bind the secondary quinone QB. These results show that kAB(2) is larger when a good proton transfer group, e.g., a hydroxyl residue (Ser, Thr) or a carboxylic acid (Asp), occupies the L223 site, supporting the proposal that Ser-L223 is a component of a proton transfer chain [Paddock, M. L., McPherson, P. H., Feher, G., & Okamura, M. Y. (1990) Proc. Natl. Acad. Sci. U.S.A. 87, 6803--6807]. The surprising result that kAB(2) is not significantly reduced in RCs with Gly at L223 suggests that a water molecule functionally replaces the missing Ser hydroxyl group in the mutant RCs. The importance of Ser-L223 in internal proton transfer reactions within the RC is discussed.

Electronic structure of Q-A in reaction centers from Rhodobacter sphaeroides. I. Electron paramagnetic resonance in single crystals

The magnitude and orientation of the electronic g-tensor of the primary electron acceptor quinone radical anion, Q-A, has been determined in single crystals of zinc-substituted reaction centers of Rhodobacter sphaeroides R-26 at 275 K and at 80 K. To obtain high spectral resolution, EPR experiments were performed at 35 GHz and the native ubiquinone-10 (UQ10) in the reaction center was replaced by fully deuterated UQ10. The principal values and the direction cosines of the g-tensor axes with respect to the crystal axes a, b, c were determined. Freezing of the single crystals resulted in only minor changes in magnitude and orientation of the g-tensor. The orientation of Q-A as determined by the g-tensor axes deviates only by a few degrees (< or = 8 degrees) from the orientation of the neutral QA obtained from an average of four different x-ray structures of Rb. sphaeroides reaction centers. This deviation lies within the accuracy of the x-ray structure determinations. The g-tensor values measured in single crystals agree well with those in frozen solutions. Variations in g-values between Q-A, Q-B, and UQ10 radical ion in frozen solutions were observed and attributed to different environments.

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.