Tuning heme redox Potentials in the cytochrome c subunit of photosynthetic reaction centers

The photosynthetic reaction center ( RC) from Rhodopseudomonas viridis contains four cytochrome c hemes. They establish the initial part of the electron transfer ( ET) chain through the RC. Despite their chemical identity, their midpoint potentials cover an interval of 440 mV. The individual heme midpoint potentials determine the ET kinetics and are therefore tuned by specific interactions with the protein environment. Here, we use an electrostatic approach based on the solution of the linearized Poisson-Boltzmann equation to evaluate the determinants of individual heme redox potentials. Our calculated redox potentials agree within 25 meV with the experimentally measured values. The heme redox potentials are mainly governed by solvent accessibility of the hemes and propionic acids, by neutralization of the negative charges at the propionates through either protonation or formation of salt bridges, by interactions with other hemes, and to a lesser extent, with other titratable protein side chains. In contrast to earlier computations on this system, we used quantum chemically derived atomic charges, considered an equilibrium-distributed protonation pattern, and accounted for interdependencies of site-site interactions. We provide values for the working potentials of all hemes as a function of the solution redox potential, which are crucial for calculations of ET rates. We identify residues whose site-directed mutation might significantly influence ET processes in the cytochrome c part of the RC. Redox potentials measured on a previously generated mutant could be reproduced by calculations based on a model structure of the mutant generated from the wild type RC.