Classical and quantum simulation of electron transfer through a polypeptide

Quantum rate theory, molecular dynamics simulations, and semiempirical electronic structure calculations are used to fully investigate electron transfer mediated by a solvated polypeptide for the first time. Using a stationary-phase approximation, the nonadiabatic electron-transfer rate constant is calculated from the nuclear free energies and the electronic coupling between the initial and final states. The former are obtained from quantum path integral and classical molecular dynamics simulations; the latter are calculated using semiempirical electronic structure calculations and the generalized Mulliken-Hush method. Importantly, no parameters are fit to kinetic data. The simulated system consists of a solvated four-proline polypeptide with a tris(bipylidine)ruthenium donor group and an oxypentamminecobalt acceptor group. From the simulation data entropy and energy contributions to the free energies are distinguished. Quantum suppression of the barrier, including important solvent contributions, is demonstrated. Although free energy profiles along the reaction coordinate are nearly parabolic, pronounced departures from harmonic behavior are found for the separate energy and entropy functions. Harmonic models of the system are compared to simulation results in order to quantify anharmonic effects. Electronic structure calculations show that electronic coupling elements vary considerably with system conformation, even when the effective donor-acceptor separation remains roughly constant. The calculations indicate that electron transfer in a significant range of conformations linking the polypeptide to the acceptor may contribute to the overall rate constant. After correction for Limitations of the solvent model, the simulations and calculations agree well with the experimental activation energy and Arrhenius prefactor.