Charge transfer through a cytochrome multiheme chain: Theory and simulation

We study sequential charge transfer within a chain of four heme cofactors located in the c-type cytochrome subunit of the photoreaction center of Rhodopseudomonas viridis from a theoretical perspective. Molecular dynamics simulations of the thermodynamic integration type are used to compute two key energies of Marcus' theory of charge transfer, the driving force ΔG and the reorganization energy λ. Due to the small exposure of the cofactors to the solvent and to charged amino acids, the outer sphere contribution to the reorganization energy almost vanishes. Interheme effective electronic couplings are estimated using ab initio wave functions and a well-parameterized semiempirical scheme for long-range interactions. From the resulting charge transfer rates, we conclude that at most the two heme molecules closest to themembrane participate in a fast recharging of the photoreaction center, whereas the remaining hemes are likely to have a different function, such as intermediate electron storage. Finally, we suggest means to verify or falsify this hypothesis