Reorganization energy of the initial electron-transfer step in photosynthetic bacterial reaction centers

The reorganization energy (lambda) for electron transfer from the primary electron donor (P*) to the adjacent bacteriochlorophyll (B) in photosynthetic bacterial reaction centers is explored by molecular-dynamics simulations. Relatively long (40 ps) molecular-dynamics trajectories are used, rather than free energy perturbation techniques. When the surroundings of the reaction center are modeled as a membrane, lambda for P* B --> P+ B- is found to be approximately 1.6 kcal/mol. The results are not sensitive to the treatment of the protein's ionizable groups, but surrounding the reaction center with water gives higher values of lambda (approximately 6.5 kcal/mol). In light of the evidence that P+ B- lies slightly below P* in energy, the small lambda obtained with the membrane model is consistent with the speed and temperature independence of photochemical charge separation. The calculated reorganization energy is smaller than would be expected if the molecular dynamics trajectories had sampled the full conformational space of the system. Because the system does not relax completely on the time scale of electron transfer, the lambda obtained here probably is more pertinent than the larger value that would be obtained for a fully equilibrated system.