Quantum and dynamical effects of proton donor-acceptor vibrational motion in nonadiabatic proton-coupled electron transfer reactions

This paper presents a general theoretical formulation for proton-coupled electron transfer (PCET) reactions. The solute is represented by a multistate valence bond model, and the active electrons and transferring proton(s) are treated quantum mechanically. This formulation enables the classical or quantum mechanical treatment of the proton donor-acceptor vibrational mode, as well as the dynamical treatment of the proton donor-acceptor mode and the solvent. Nonadiabatic rate expressions are presented for PCET reactions in a number of well-defined limits for both dielectric continuum and molecular representations of the environment. The dynamical rate expressions account for correlations between the fluctuations of the proton donor-acceptor distance and the nonadiabatic PCET coupling. The quantities in the rate expressions can be calculated with a dielectric continuum model or a molecular dynamics simulation of the full system. The significance of the quantum and dynamical effects of the proton donor-acceptor mode is illustrated with applications to model PCET systems.     

Proton-coupled electron transfer in a model for tyrosine oxidation in photosystem II

Theoretical calculations of a model for tyrosine oxidation in photosystem II are presented. In this model system, an electron is transferred to ruthenium from tyrosine, which is concurrently deprotonated. This investigation is motivated by experimental measurements of the dependence of the rates on pH and temperature (Sj?din et al. J. Am. Chem. Soc. 2000, 122, 3932). The mechanism is proton-coupled electron transfer (PCET) at pH < 10 when the tyrosine is initially protonated and is single electron transfer (ET) for pH > 10 when the tyrosine is initially deprotonated. The PCET rate increases monotonically with pH, whereas the single ET rate is independent of pH and is 2 orders of magnitude faster than the PCET rate. The calculations reproduce these experimentally observed trends. The pH dependence for the PCET reaction arises from the decrease in the reaction free energies with pH. The calculations indicate that the larger rate for single ET arises from a combination of factors, including the smaller solvent reorganization energy for ET and the averaging of the coupling for PCET over the reactant and product hydrogen vibrational wave functions (i.e., a vibrational overlap factor in the PCET rate expression). The temperature dependence of the rates, the solvent reorganization energies, and the deuterium kinetic isotope effects determined from the calculations are also consistent with the experimental results.     

Proton-coupled electron transfer in solution, proteins, and electrochemistry

Recent advances in the theoretical treatment of proton-coupled electron transfer (PCET) reactions are reviewed. These reactions play an important role in a wide range of biological processes, as well as in fuel cells, solar cells, chemical sensors, and electrochemical devices. A unified theoretical framework has been developed to describe both sequential and concerted PCET, as well as hydrogen atom transfer (HAT). A quantitative diagnostic has been proposed to differentiate between HAT and PCET in terms of the degree of electronic nonadiabaticity, where HAT corresponds to electronically adiabatic proton transfer and PCET corresponds to electronically nonadiabatic proton transfer. In both cases, the overall reaction is typically vibronically nonadiabatic. A series of rate constant expressions have been derived in various limits by describing the PCET reactions in terms of nonadiabatic transitions between electron-proton vibronic states. These expressions account for the solvent response to both electron and proton transfer and the effects of the proton donor-acceptor vibrational motion. The solvent and protein environment can be represented by a dielectric continuum or described with explicit molecular dynamics. These theoretical treatments have been applied to numerous PCET reactions in solution and proteins. Expressions for heterogeneous rate constants and current densities for electrochemical PCET have also been derived and applied to model systems.     

Nonadiabatic proton-coupled electron transfer reactions: impact of donor-acceptor vibrations, reorganization energies, and couplings on dynamics and rates

Fundamental aspects of proton-coupled electron transfer (PCET) reactions in solution are analyzed with molecular dynamics simulations for a series of model systems. The analysis addresses the impact of the solvent reorganization energy, the proton donor-acceptor mode vibrational frequency, and the distance dependence of the nonadiabatic coupling on the dynamics of the reaction and the magnitude of the rate. The rate for nonadiabatic PCET is expressed in terms of a time-dependent probability flux correlation function. The time dependence of the probability flux correlation function is determined mainly by the solvent reorganization energy and is not significantly influenced by the proton donor-acceptor frequency or the distance dependence of the nonadiabatic coupling. The magnitude of the PCET rate becomes greater as the solvent reorganization energy decreases, the proton donor-acceptor frequency decreases, and the distance dependence of the nonadiabatic coupling increases. The approximations underlying a previously derived analytical PCET rate expression are also investigated. The short-time approximation for the solvent is valid for these types of systems. In addition, solvent damping effects on the proton donor-acceptor motion are not significant on the time scale of the probability flux. The rates calculated from the molecular dynamics simulations agree well with those calculated from the analytical rate expression.     

Proton-coupled electron transfer in soybean lipoxygenase

The proton-coupled electron transfer reaction catalyzed by soybean lipoxygenase-1 is studied with a multistate continuum theory that represents the transferring hydrogen nucleus as a quantum mechanical wave function. The inner-sphere reorganization energy of the iron cofactor is calculated with density functional theory, and the outer-sphere reorganization energy of the protein is calculated with the frequency-resolved cavity model for conformations obtained with docking simulations. Both classical and quantum mechanical treatments of the proton donor-acceptor vibrational motion are presented. The temperature dependence of the calculated rates and kinetic isotope effects is in agreement with the experimental data. The weak temperature dependence of the rates is due to the relatively small free energy barrier arising from a balance between the reorganization energy and the reaction free energy. The unusually high deuterium kinetic isotope effect of 81 is due to the small overlap of the reactant and product proton vibrational wave functions and the dominance of the lowest energy reactant and product vibronic states in the tunneling process. The temperature dependence of the kinetic isotope effect is strongly influenced by the proton donor-acceptor distance with the dominant contribution to the overall rate. This dominant proton donor-acceptor distance is significantly smaller than the equilibrium donor-acceptor distance and is determined by a balance between the larger coupling and the smaller Boltzmann probability as the distance decreases. Thus, the proton donor-acceptor vibrational motion plays a vital role in decreasing the dominant donor-acceptor distance relative to its equilibrium value to facilitate the proton-coupled electron transfer reaction.     

Buffer-assisted proton-coupled electron transfer in a model rhenium-tyrosine complex

The mechanism for tyrosyl radical generation in the [Re(P-Y)(phen)(CO)3]PF6 complex is investigated with a multistate continuum theory for proton-coupled electron transfer (PCET) reactions. Both water and the phosphate buffer are considered as potential proton acceptors. The calculations indicate that the model in which the proton acceptor is the phosphate buffer species HPO(4)2- can successfully reproduce the experimentally observed pH dependence of the overall rate and H/D kinetic isotope effect, whereas the model in which the proton acceptor is water is not physically reasonable for this system. The phosphate buffer species HPO4(2-) is favored over water as the proton acceptor in part because the proton donor-acceptor distance is approximately 0.2 A smaller for the phosphate acceptor due to its negative charge. The physical quantities impacting the overall rate constant, including the reorganization energies, reaction free energies, activation free energies, and vibronic couplings for the various pairs of reactant/product vibronic states, are analyzed for both hydrogen and deuterium transfer. The dominant contribution to the rate arises from nonadiabatic transitions between the ground reactant vibronic state and the third product vibronic state for hydrogen transfer and the fourth product vibronic state for deuterium transfer. These contributions dominate over contributions from lower product states because of the larger vibronic coupling, which arises from the greater overlap between the reactant and product vibrational wave functions. These calculations provide insight into the fundamental mechanism of tyrosyl radical generation, which plays an important role in a wide range of biologically important processes.     

Cation-modulated electron-transfer channel: H-atom transfer vs proton-coupled electron transfer with a variable electron-transfer channel in acylamide units

The mechanism of proton transfer (PT)/electron transfer (ET) in acylamide units was explored theoretically using density functional theory in a representative model (a cyclic coupling mode between formamide and the N-dehydrogenated formamidic radical, FF). In FF, PT/ET normally occurs via a seven-center cyclic proton-coupled electron transfer (PCET) mechanism with a N-->N PT and an O-->O ET. However, when different hydrated metal ions are bound to the two oxygen sites of FF, the PT/ET mechanism may significantly change. In addition to their inhibition of PT/ET rate, the hydrated metal ions can effectively regulate the FF PT/ET cooperative mechanism to produce a single pathway hydrogen atom transfer (HAT) or a flexible proton coupled electron transfer (PCET) mechanism by changing the ET channel. The regulation essentially originates from the change in the O...O bond strength in the transition state, subject to the binding ability of the hydrated metal ions. In general, the high valent metal ions and those with large binding energies can promote HAT, and the low valent metal ions and those with small binding energies favor PCET. Hydration may reduce the Lewis acidity of cations, and thus favor PCET. Good correlations among the binding energies, barrier heights, spin density distributions, O...O contacts, and hydrated metal ion properties have been found, which can be used to interpret the transition in the PT/ET mechanism. These findings regarding the modulation of the PT/ET pathway via hydrated metal ions may provide useful information for a greater understanding of PT/ET cooperative mechanisms, and a possible method for switching conductance in nanoelectronic devices.     

Model system-bath Hamiltonian and nonadiabatic rate constants for proton-coupled electron transfer at electrode-solution interfaces

An extension of the Anderson-Newns-Schmickler model for electrochemical proton-coupled electron transfer (PCET) is presented. This model describes reactions in which electron transfer between a solute complex in solution and an electrode is coupled to proton transfer within the solute complex. The model Hamiltonian is derived in a basis of electron-proton vibronic states defined within a double adiabatic approximation for the electrons, transferring proton, and bath modes. The interaction term responsible for electronic transitions between the solute complex and the electrode depends on the proton donor-acceptor vibrational mode within the solute complex. This model Hamiltonian is used to derive the anodic and cathodic rate constants for nonadiabatic electrochemical PCET. The derivation is based on the master equations for the reduced density matrix of the electron-proton subsystem, which includes the electrons of the solute complex and the electrode, as well as the transferring proton. The rate constant expressions differ from analogous expressions for electrochemical electron transfer because of the summation over electron-proton vibronic states and the dependence of the couplings on the proton donor-acceptor vibrational motion. These differences lead to additional contributions to the total reorganization energy, an additional exponential temperature-dependent prefactor, and a temperature-dependent term in the effective activation energy that has different signs for the anodic and cathodic processes. This model can be generalized to describe both nonadiabatic and adiabatic electrochemical PCET reactions and provides the framework for the inclusion of additional effects, such as the breaking and forming of other chemical bonds.     

Dynamics of proton transfer and vibrational relaxation in dilute hydrofluoric acid

The molecular mechanisms in both vibrational relaxation and proton transfer (PT) associated with infrared (IR)-induced PT in a dilute hydrofluoric acid solution at ambient temperature are studied by molecular dynamics (MD) simulations with the multistate empirical valence bond model. To investigate the solvation dynamics, a collective solvent coordinate and its perpendicular bath modes are defined from the diabatic energy gap and their motions are examined by the generalized Langevin equation (GLE) formalism. The GLE analysis using the equilibrium MD simulation shows that the major solvent reorganizations in the PT are represented by the libration and hindered translation. In particular, the libration gives the stronger coupling to the solvent reorganization and the faster relaxation. The nonequilibrium MD simulation demonstrated that both the HF stretching vibration and the solvent reorganization relax on a similar time scale and thus compete in the PT. It also supported the "presolvation mechanism" for the PT in this system.     

Proton-coupled electron transfer from tryptophan: a concerted mechanism with water as proton acceptor

The mechanism of proton-coupled electron transfer (PCET) from tyrosine in enzymes and synthetic model complexes is under intense discussion, in particular the pH dependence of the PCET rate with water as proton acceptor. Here we report on the intramolecular oxidation kinetics of tryptophan derivatives linked to [Ru(bpy)(3)](2+) units with water as proton acceptor, using laser flash-quench methods. It is shown that tryptophan oxidation can proceed not only via a stepwise electron-proton transfer (ETPT) mechanism that naturally shows a pH-independent rate, but also via another mechanism with a pH-dependent rate and higher kinetic isotope effect that is assigned to concerted electron-proton transfer (CEP). This is in contrast to current theoretical models, which predict that CEP from tryptophan with water as proton acceptor can never compete with ETPT because of the energetically unfavorable PT part (pK(a)(Trp(?)H(+)) = 4.7 ? pK(a)(H(3)O(+)) ≈ -1.5). The moderate pH dependence we observe for CEP cannot be explained by first-order reactions with OH(-) or the buffers and is similar to what has been demonstrated for intramolecular PCET in [Ru(bpy)(3)](3+)-tyrosine complexes (Sjo?din, M.; et al. J. Am. Chem. Soc.2000, 122, 3932. Irebo, T.; et al. J. Am. Chem. Soc.2007, 129, 15462). Our results suggest that CEP with water as the proton acceptor proves a general feature of amino acid oxidation, and provide further experimental support for understanding of the PCET process in detail.     

Switching the redox mechanism: models for proton-coupled electron transfer from tyrosine and tryptophan

The coupling of electron and proton transfer is an important controlling factor in radical proteins, such as photosystem II, ribinucleotide reductase, cytochrome oxidases, and DNA photolyase. This was investigated in model complexes in which a tyrosine or tryptophan residue was oxidized by a laser-flash generated trisbipyridine-Ru(III) moiety in an intramolecular, proton-coupled electron transfer (PCET) reaction. The PCET was found to proceed in a competition between a stepwise reaction, in which electron transfer is followed by deprotonation of the amino acid radical (ETPT), and a concerted reaction, in which both the electron and proton are transferred in a single reaction step (CEP). Moreover, we found that we could analyze the kinetic data for PCET by Marcus' theory for electron transfer. By altering the solution pH, the strength of the Ru(III) oxidant, or the identity of the amino acid, we could induce a switch between the two mechanisms and obtain quantitative data for the parameters that control which one will dominate. The characteristic pH-dependence of the CEP rate (M. Sjodin et al. J. Am. Chem. Soc. 2000, 122, 3932) reflects the pH-dependence of the driving force caused by proton release to the bulk. For the pH-independent ETPT on the other hand, the driving force of the rate-determining ET step is pH-independent and smaller. On the other hand, temperature-dependent data showed that the reorganization energy was higher for CEP, while the pre-exponential factors showed no significant difference between the mechanisms. Thus, the opposing effect of the differences in driving force and reorganization energy determines which of the mechanisms will dominate. Our results show that a concerted mechanism is in general quite likely and provides a low-barrier reaction pathway for weakly exoergonic reactions. In addition, the kinetic isotope effect was much higher for CEP (kH/kD > 10) than for ETPT (kH/kD = 2), consistent with significant changes along the proton reaction coordinate in the rate-determining step of CEP.     

Theoretical analysis of proton relays in electrochemical proton-coupled electron transfer

The coupling of long-range electron transfer to proton transport over multiple sites plays a vital role in many biological and chemical processes. Recently the concerted proton-coupled electron transfer (PCET) reaction in a molecule with a hydrogen-bond relay inserted between the proton donor and acceptor sites was studied electrochemically. The standard rate constants and kinetic isotope effects (KIEs) were measured experimentally for this double proton transfer system and a related single proton transfer system. In the present paper, these systems are studied theoretically using vibronically nonadiabatic rate constant expressions for electrochemical PCET. Application of this approach to proton relays requires the calculation of multidimensional proton vibrational wave functions and the incorporation of multiple proton donor-acceptor motions. The decrease in proton donor-acceptor distances due to thermal fluctuations and the contributions from excited electron-proton vibronic states play important roles in these systems. The calculated KIEs and the ratio of the standard rate constants for the single and double proton transfer systems are in agreement with the experimental data. The calculations indicate that the standard PCET rate constant is lower for the double proton transfer system because of the smaller overlap integral between the ground state reduced and oxidized proton vibrational wave functions, resulting in greater contributions from excited electron-proton vibronic states with higher free energy barriers. The theory predicts that this rate constant may be increased by modifying the molecule in a manner that decreases the equilibrium proton donor-acceptor distances or alters the molecular thermal motions to facilitate the concurrent decrease of these distances. These insights may guide the design of more efficient catalysts for energy conversion devices.     

Anomalous temperature-isotope dependence in proton-coupled electron transfer

Motivated by the experiments of Hodgkiss et al. [J. Phys. Chem. (submitted)] on electron transfer (ET) through a H-bonding interface, we present a new theoretical model for proton-coupled electron transfer (PCET) in the condensed phase, that does not involve real proton transfer. These experiments, which directly probe the joint T-isotope effects in coupled charge transfer reactions, show anomalous T dependence in k(H)k(D), where k(H) and k(D) are the ET rates through the H-bonding interface with H-bonded protons and deuterons, respectively. We address the anomalous T dependence of the k(H)k(D) in our model by attributing the modulation of the electron tunneling dynamics to bath-induced fluctuations in the proton coordinate, so that the mechanism for coupled charge transfer might be better termed vibrationally assisted ET rather than PCET. We argue that such a mechanism may be relevant to understanding traditional PCET processes, i.e., those in which protons undergo a transfer from donor to acceptor during the course of ET, provided there is an appropriate time scale separating both coupled charge transfers. Likewise, it may also be useful in understanding long-range ET in proteins, where tunneling pathways between redox cofactors often pass through H-bonded amino acid residues, or other systems with sufficiently decoupled proton and electron donating functionalities.     

Prediction of reorganization free energies for biological electron transfer: a comparative study of Ru-modified cytochromes and a 4-helix bundle protein

The acceleration of electron transfer (ET) rates in redox proteins relative to aqueous solutes can be attributed to the protein's ability to reduce the nuclear response or reorganization upon ET, while maintaining sufficiently high electronic coupling. Quantitative predictions of reorganization free energy remain a challenge, both experimentally and computationally. Using density functional calculations and molecular dynamics simulation with an electronically polarizable force field, we report reorganization free energies for intraprotein ET in four heme-containing ET proteins that differ in their protein fold, hydrophilicity, and solvent accessibility of the electron-accepting group. The reorganization free energies for ET from the heme cofactors of cytochrome c and b(5) to solvent exposed Ru-complexes docked to histidine residues at the surface of these proteins fall within a narrow range of 1.2-1.3 eV. Reorganization free energy is significantly lowered in a designed 4-helix bundle protein where both redox active cofactors are protected from the solvent. For all ET reactions investigated, the major components of reorganization are the solvent and the protein, with the solvent contributing close to or more than 50% of the total. In three out of four proteins, the protein reorganization free energy can be viewed as a collective effect including many residues, each of which contributing a small fraction. These results have important implications for the design of artificial electron transport proteins. They suggest that reorganization free energy may in general not be effectively controlled by single point mutations, but to a large extent by the degree of solvent exposure of the ionizable cofactors.     

The reorganization energy of electron transfer in nonpolar solvents: molecular level treatment of the solvent

The intermolecular electron transfer in a solute pair consisting of pyrene and dimethylaniline is investigated in a nonpolar solvent, n-hexane. The earlier elaborated approach [M. Tachiya, J. Phys Chem. 97, 5911 (1993)] is used; this method provides a physically relevant background for separating inertial and inertialess polarization responses for both nonpolarizable and polarizable molecular level simulations. The molecular-dynamics technique was implemented for obtaining the equilibrium ensemble of solvent configurations. The nonpolar solvent, n-hexane, was treated in terms of OPLS-AA parametrization. Solute Lennard-Jones parameters were taken from the same parametrization. Solute charge distributions of the initial and final states were determined using ab initio level [HF/6-31G(d,p)] quantum-chemical calculations. Configuration analysis was performed explicitly taking into account the anisotropic polarizability of n-hexane. It is shown that the Gaussian law well describes calculated distribution functions of the solvent coordinate, therefore, the rate constant of the ET reaction can be characterized by the reorganization energy. Evaluated values of the reorganization energies are in a range of 0.03-0.11 eV and significant contribution (more then 40% of magnitude) comes from anisotropic polarizability. Investigation of the reorganization energy lambda dependence on the solute pair separation distance d revealed unexpected behavior. The dependence has a very sharp peak at the distance d=7 A where solvent molecules are able to penetrate into the intermediate space between the solute pair. The reason for such behavior is clarified. This new effect has a purely molecular origin and cannot be described within conventional continuum solvent models.     

Activation entropy of electron transfer reactions

We report microscopic calculations of free energies and entropies for intramolecular electron transfer reactions. The calculation algorithm combines the atomistic geometry and charge distribution of a molecular solute obtained from quantum calculations with the microscopic polarization response of a polar solvent expressed in terms of its polarization structure factors. The procedure is tested on a donor–acceptor complex in which ruthenium donor and cobalt acceptor sites are linked by a four-proline polypeptide. The reorganization energies and reaction energy gaps are calculated as a function of temperature by using structure factors obtained from our analytical procedure and from computer simulations. Good agreement between two procedures and with direct computer simulations of the reorganization energy is achieved. The microscopic algorithm is compared to the dielectric continuum calculations. We found that the strong dependence of the reorganization energy on the solvent refractive index predicted by continuum models is not supported by the microscopic theory. Also, the reorganization and overall solvation entropies are substantially larger in the microscopic theory compared to continuum models.     

Long-range proton transfer in aqueous Acid-base reactions

We study the mechanism of proton transfer (PT) in the aqueous acid-base reaction between the photoacid 8-hydroxy-1,3,6-pyrenetrisulfonic acid (HPTS) and acetate by probing the vibrational resonances of HPTS, acetate, and the hydrated proton with femtosecond mid-infrared laser pulses. We find that PT takes place in a distribution of hydrogen-bound reaction complexes that differ in the number of water molecules separating the acid and the base. The number of intervening water molecules ranges from 0 to 5, which, together with a strongly distance-dependent PT rate, explains the observed highly nonexponential reaction kinetics. The kinetic isotope effect for the reaction is determined to be 1.5, indicating that tunneling does not play a significant role in the transfer of the proton. Rather, the transfer mechanism is best described in terms of the adiabatic PT picture as it has been formulated by Hynes and co-workers [Staib, A.; Borgis, D.; Hynes, J. T. J. Chem. Phys. 1995, 102, 2487. Ando, K.; Hynes, J. T. J. Phys. Chem. B 1997, 101, 10464.], where solvent fluctuations play an essential role in forming the correct hydrogen-bond configuration and solvent polarization to facilitate PT.     

Kinetic effects of hydrogen bonds on proton-coupled electron transfer from phenols

The kinetics and mechanism of proton-coupled electron transfer (PCET) from a series of phenols to a laser flash generated [Ru(bpy)(3)](3+) oxidant in aqueous solution was investigated. The reaction followed a concerted electron-proton transfer mechanism (CEP), both for the substituted phenols with an intramolecular hydrogen bond to a carboxylate group and for those where the proton was directly transferred to water. Without internal hydrogen bonds the concerted mechanism gave a characteristic pH-dependent rate for the phenol form that followed a Marcus free energy dependence, first reported for an intramolecular PCET in Sj?din, M. et al. J. Am. Chem. Soc. 2000, 122, 3932-3962 and now demonstrated also for a bimolecular oxidation of unsubstituted phenol. With internal hydrogen bonds instead, the rate was no longer pH-dependent, because the proton was transferred to the carboxylate base. The results suggest that while a concerted reaction has a relatively high reorganization energy (lambda), this may be significantly reduced by the hydrogen bonds, allowing for a lower barrier reaction path. It is further suggested that this is a general mechanism by which proton-coupled electron transfer in radical enzymes and model complexes may be promoted by hydrogen bonding. This is different from, and possibly in addition to, the generally suggested effect of hydrogen bonds on PCET in enhancing the proton vibrational wave function overlap between the reactant and donor states. In addition we demonstrate how the mechanism for phenol oxidation changes from a stepwise electron transfer-proton transfer with a stronger oxidant to a CEP with a weaker oxidant, for the same series of phenols. The hydrogen bonded CEP reaction may thus allow for a low energy barrier path that can operate efficiently at low driving forces, which is ideal for PCET reactions in biological systems.