Multidimensional Quantum Mechanical Modeling of Electron Transfer and Electronic Coherence in Plant Cryptochromes: The Role of Initial Bath Conditions.

A multidimensional quantum mechanical protocol is used to describe the photoinduced electron transfer and electronic coherence in plant cryptochromes without any semiempirical, e.g., experimentally obtained, parameters. Starting from a two-level spin-boson Hamiltonian we look at the effect that the initial photoinduced nuclear bath distribution has on an intermediate step of this biological electron transfer cascade for two idealized cases. The first assumes a slow equilibration of the nuclear bath with respect to the previous electron transfer step that leads to an ultrafast decay with little temperature dependence; while the second assumes a prior fast bath equilibration on the donor potential energy surface leading to a much slower decay, which contrarily displays a high temperature dependence and a better agreement with previous theoretical and experimental results. Beyond Marcus and semiclassical pictures these results unravel the strong impact that the presence or not of equilibrium initial conditions has on the electronic population and coherence dynamics at the quantum dynamics level in this and conceivably in other biological electron transfer cascades.

Quantum effects in ultrafast electron transfers within cryptochromes.

Cryptochromes and photolyases are flavoproteins that may undergo ultrafast charge separation upon electronic excitation of their flavin cofactors. Charge separation involves chains of three or four tryptophan residues depending on the protein of interest. The molecular mechanisms of these processes are not completely clear. In the present work we investigate the relevance of quantum effects like the occurrence of nuclear tunneling and of coherences upon charge transfer in Arabidopsis thaliana cryptochromes. The possible breakdown of the Condon approximation is also investigated. We have devised a simulation protocol based on the realization of molecular dynamics simulations on diabatic potential energy surfaces defined at the hybrid constrained density functional theory/molecular mechanics level. The outcomes of the simulations are analyzed through various dedicated kinetics schemes related to the Marcus theory that account for the aforementioned quantum effects. MD simulations also provide a basic material to define realistic model Hamiltonians for subsequent quantum dissipative dynamics. To carry out quantum simulations, we have implemented an algorithm based on the Hierarchical Equations of Motion. With this new tool in hand we have been able to model the electron transfer chain considering either two- or three-state models. Kinetic models and quantum simulations converge to the conclusion that quantum effects have a significant impact on the rate of charge separation. Nuclear tunneling involving atoms of the tryptophan redox cofactors as well as of the environment (protein atoms and water molecules) is significant. On the other hand non-Condon effects are negligible in most simulations. Taken together, the results of the present work provide new insights into the molecular mechanisms controlling charge separation in this family of flavoproteins.

Energetics of Photoinduced Charge Migration within the Tryptophan Tetrad of an Animal (6-4) Photolyase.

Cryptochromes and photolyases are flavoproteins that undergo cascades of electron/hole transfers after excitation of the flavin cofactor. It was recently discovered that animal (6-4) photolyases, as well as animal cryptochromes, feature a chain of four tryptophan residues, while other members of the family contain merely a tryptophan triad. Transient absorption spectroscopy measurements on Xenopus laevis (6-4) photolyase have shown that the fourth residue is effectively involved in photoreduction but at the same time could not unequivocally ascertain the final redox state of this residue. In this article, polarizable molecular dynamics simulations and constrained density functional theory calculations are carried out to reveal the energetics of charge migration along the tryptophan tetrad. Migration toward the fourth tryptophan is found to be thermodynamically favorable. Electron transfer mechanisms are sought either through an incoherent hopping mechanism or through a multiple sites tunneling process. The Jortner-Bixon formulation of electron transfer (ET) theory is employed to characterize the hopping mechanism. The interplay between electron transfer and relaxation of protein and solvent is analyzed in detail. Our simulations confirm that ET in (6-4) photolyase proceeds out of equilibrium. Multiple site tunneling is modeled with the recently proposed flickering resonance mechanism. Given the position of energy levels and the distribution of electronic coupling values, tunneling over three tryptophan residues may become competitive in some cases, although a hopping mechanism is likely to be the dominant channel. For both reactive channels, computed rates are very sensitive to the starting protein configuration, suggesting that both can take place and eventually be mixed, depending on the state of the system when photoexcitation takes place.

Diffusion Rates for Hydrogen on Pd(111) from Molecular Quantum Dynamics Calculations.

The van Hove formula for the dynamical structure factor (DSF) related to particle scattering at mobile adsorbates is extended to include the relaxation of the adsorbates' vibrational states. The total rate obtained from the DSF is assumed to be the sum of a diffusion and a relaxation rate. A simple kinetic model to support this assumption is presented. To illustrate its potential applicability, the formula is evaluated using wave functions, energies, and lifetimes of vibrational states obtained for H/Pd(111) from first-principle calculations. Results show that quantum effects can be expected to be important even at room temperature.