Calculated protein and proton motions coupled to electron transfer: electron transfer from QA- to QB in bacterial photosynthetic reaction centers.

Reaction centers from Rhodobacter sphaeroides were subjected to Monte Carlo sampling to determine the Boltzmann distribution of side-chain ionization states and positions and buried water orientation and site occupancy. Changing the oxidation states of the bacteriochlorophyll dimer electron donor (P) and primary (QA) and secondary (QB) quinone electron acceptors allows preparation of the ground (all neutral), P+QA-, P+QB-, P0QA-, and P0QB- states. The calculated proton binding going from ground to other oxidation states and the free energy of electron transfer from QA-QB to form QAQB- (DeltaGAB) compare well with experiment from pH 5 to pH 11. At pH 7 DeltaGAB is measured as -65 meV and calculated to be -80 meV. With fixed protein positions as in standard electrostatic calculations, DeltaGAB is +170 meV. At pH 7 approximately 0.2 H+/protein is bound on QA reduction. On electron transfer to QB there is little additional proton uptake, but shifts in side chain protonation and position occur throughout the protein. Waters in channels leading from QB to the surface change site occupancy and orientation. A cluster of acids (GluL212, AspL210, and L213) and SerL223 near QB play important roles. A simplified view shows this cluster with a single negative charge (on AspL213 with a hydrogen bond to SerL233) in the ground state. In the QB- state the cluster still has one negative charge, now on the more distant AspL210. AspL213 and SerL223 move so SerL223 can hydrogen bond to QB-. These rearrangements plus other changes throughout the protein make the reaction energetically favorable.