Understanding Short-Range Electron-Transfer Dynamics in Proteins

Elementary Reactions in the Functional Triads of the Blue-Light Photoreceptor BLUF Domain

The blue light using the flavin (BLUF) domain is one of the smallest photoreceptors in nature, which consists of a unique bidirectional electron-coupled proton relay process in its photoactivation reaction cycle. This perspective summarizes our recent efforts in dissecting the photocycle into three elementary processes, including proton-coupled electron transfer (PCET), proton rocking, and proton relay. Using ultrafast spectroscopy, we have determined the temporal sequence, rates, kinetic isotope effects (KIEs), and concertedness of these elementary steps. Our findings provide important implications for illuminating the photoactivation mechanism of the BLUF domain and suggest an engineering platform to characterize intricate reactions involving proton motions that are ubiquitous in nonphotosensitive protein machines.

Origin of the multi-phasic quenching dynamics in the BLUF domains across the species

Blue light using flavin (BLUF) photoreceptors respond to light via one of nature's smallest photo-switching domains. Upon photo-activation, the flavin cofactor in the BLUF domain exhibits multi-phasic dynamics, quenched by a proton-coupled electron transfer reaction involving the conserved Tyr and Gln. The dynamic behavior varies drastically across different species, the origin of which remains controversial. Here, we incorporate site-specific fluorinated Trp into three BLUF proteins, i.e., AppA, OaPAC and SyPixD, and characterize the percentages for the Wout, WinNHin and WinNHout conformations using 19F nuclear magnetic resonance spectroscopy. Using femtosecond spectroscopy, we identify that one key WinNHin conformation can introduce a branching one-step proton transfer in AppA and a two-step proton transfer in OaPAC and SyPixD. Correlating the flavin quenching dynamics with the active-site structural heterogeneity, we conclude that the quenching rate is determined by the percentage of WinNHin, which encodes a Tyr-Gln configuration that is not conducive to proton transfer.

Tryptophan to Tryptophan Hole Hopping in an Azurin Construct

Electron transfer (ET) between neutral and cationic tryptophan residues in the azurin construct [ReI(H126)(CO)3(dmp)](W124)(W122)CuI (dmp = 4,7-Me2-1,10-phenanthroline) was investigated by Born-Oppenheimer quantum-mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) simulations. We focused on W124•+ ← W122 ET, which is the middle step of the photochemical hole-hopping process *ReII(CO)3(dmp•-) ← W124 ← W122 ← CuI, where sequential hopping amounts to nearly 10,000-fold acceleration over single-step tunneling (ACS Cent. Sci. 2019, 5, 192-200). In accordance with experiments, UKS-DFT QM/MM/MD simulations identified forward and reverse steps of W124•+ ↔ W122 ET equilibrium, as well as back ET ReI(CO)3(dmp•-) → W124•+ that restores *ReII(CO)3(dmp•-). Strong electronic coupling between the two indoles (≥40 meV in the crossing region) makes the productive W124•+ ← W122 ET adiabatic. Energies of the two redox states are driven to degeneracy by fluctuations of the electrostatic potential at the two indoles, mainly caused by water solvation, with contributions from the protein dynamics in the W122 vicinity. ET probability depends on the orientation of Re(CO)3(dmp) relative to W124 and its rotation diminishes the hopping yield. Comparison with hole hopping in natural systems reveals structural and dynamics factors that are important for designing efficient hole-hopping processes.

Ultrafast Nonequilibrium Dynamics of Vibrationally Hot Electron Transfer in Flavodoxin

The understanding of ultrafast short-range electron transfer (ET) in proteins remains challenging, and thorough studies on well-defined biological systems are demanding. Here, we utilized two types of flavodoxins and designed a series of mutants on two positions to systematically characterize the complete photoinduced redox cycles. We identified one position with a favorable orientation and distance for ultrafast ET in a few femtoseconds and the other position is relatively flexible with a longer ET time scale. We found that all forward and back ET dynamics are ultrafast nonequilibrium processes, occurring through highly vibronic states and ending in vibrationally hot ground states with subsequent cooling relaxation to efficiently dissipate photon energy into the protein environment.

Ultrafast Dynamics of Nonequilibrium Short-Range Electron Transfer in Semiquinone Flavodoxin

Short-range protein electron transfer (ET) is crucially important in light-induced biological processes such as in photoenzymes and photoreceptors and often occurs on time scales similar to those of environment fluctuations, leading to a coupled dynamic process. Herein, we use semiquinone Anabaena flavodoxin to characterize the ultrafast photoinduced redox cycle of the wild type and seven mutants by ultrafast spectroscopy. We have found that the forward and backward ET dynamics show stretched behaviors in a few picoseconds (1-5 ps), indicating a coupling with the local protein fluctuations. By comparison with the results from semiquinone D. vulgaris flavodoxin, we find that the electronic coupling is crucial to the ET rates. With our new nonergodic model, we obtain smaller values of the outer reorganization energy (λ_o^γ) of environment fluctuations and the reaction free energy force (ΔG^γ), a signature of nonequilibrium ET dynamics.

Ultrafast nonequilibrium dynamics of short-range protein electron transfer in flavodoxin

Short-range protein electron transfer (ET) is ubiquitous in biology and is often observed in photosynthesis, photoreceptors and photoenzymes. These ET processes occur on an ultrafast timescale from femtoseconds to picoseconds at a short donor-acceptor distance within 10 Å, and thus couple with local environmental fluctuations. Here, we use oxidized Anabaena flavodoxin as a model system and have systematically studied the photoinduced redox cycle of the wild type and seven mutant proteins by femtosecond spectroscopy. We observed a series of ultrafast dynamics from the initial charge separation in 100-200 fs, subsequent charge recombination in 1-2 ps and final vibrational cooling process of the products in 3-6 ps. We further characterized the active-site solvation and observed the relaxations in 1-200 ps, indicating a nonergodic ET dynamics. With our new ET model, we uncovered a minor outer (solvent) reorganization energy and a large inner (donor and acceptor) reorganization energy, suggesting a frozen active site in the initial ultrafast ET while the back ET couples with the environment relaxations. The vibronically coupled back ET dynamics was first reported in D. vulgaris flavodoxin and here is observed in Anabaena flavodoxin again, completely due to the faster ET dynamics than the cooling relaxations. We also compared the two flavodoxin structures, revealing a stronger coupling with the donor tyrosine in Anabaena. All ultrafast ET dynamics are from the large donor-acceptor couplings and the minor activation barriers due to the reaction free energies being close to the inner reorganization energies. These observations should be general to many redox reactions in flavoproteins.

Electron Transfer-Induced Active Site Structural Relaxation in 64-Photolyase of Drosophila melanogaster

While catalytic electron flow and photoreactivation of CPD-photolyases are increasingly understood, the microscopic details of the 64-photolyase repair mechanism are perpetually debated. Here, we investigate in long-time (μs) molecular dynamics simulations combined with extensive quantum mechanical/molecular mechanical (QM/MM) simulations the primary electron transfer (ET) reactions in 64-photolyase of Drosophila melanogaster (D. melanogaster). The characterization of the relative energetics of locally excited and charge separated states in the (6-4) photoproduct enzyme repair complex reveals a charge-separated state involving the adenine moiety of the FADH^- cofactor that facilitates reduction of the photoproduct. Microscopic details of the collective reaction coordinate of ET reactions are identified that involve the reorganization of the hydrogen bond network and structural relaxation of the active site. The simulations reveal complex active site relaxation dynamics involving distinguished amino acids (Lys246, His365, and His369), conformational reorganization of the hydroxyl group of the (6-4) photoproduct, and a strengthening of hydrogen bonds with immobilized water molecules. In particular, rotation of the Lys246 side chain is found to impose a double-well character along the reaction coordinate of the ET reaction. Our findings suggest that the primary ET reactions in the (6-4) photoproduct enzyme repair complex of D. melanogaster are governed by a complex multi-minima active site relaxation dynamics and potentially precede the equilibration of the protein. ET pathways mediated by the adenine moiety and the 5' side of the photoproduct are proposed to be relevant for triggering the catalytic (6-4) photoproduct reactivation.

Exact eigenenergies of a model of vibronically coupled electron transfer reactions

We present the first exact solution to the time-independent Schrödinger equation of a model Hamiltonian consisting of a vibrational mode coupled to three electronic states. This Hamiltonian serves as a generic model for photo-induced electronic transfer reactions. The solution is non-perturbative and can be applied to ET reactions with weak and strong electronic and vibrational coupling strengths. This work suggests a new direction towards understanding the vibronic effects in ET dynamics beyond the non-adiabatic limit and Condon approximation.

Electron transfer pathways from quantum dynamics simulations

This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, -1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant.

Effects of nonequilibrium fluctuations on ultrafast short-range electron transfer dynamics

A variety of electron transfer (ET) reactions in biological systems occurs at short distances and is ultrafast. Many of them show behaviors that deviate from the predictions of the classic Marcus theory. Here, we show that these ultrafast ET dynamics highly depend on the coupling between environmental fluctuations and ET reactions. We introduce a dynamic factor, γ (0 ≤ γ ≤ 1), to describe such coupling, with 0 referring to the system without coupling to a “frozen” environment, and 1 referring to the system’s complete coupling with the environment. Significantly, this system’s coupling with the environment modifies the reaction free energy, ΔG^γ, and the reorganization energy, λ^γ, both of which become smaller. This new model explains the recent ultrafast dynamics in flavodoxin and elucidates the fundamental mechanism of nonequilibrium ET dynamics, which is critical to uncovering the molecular nature of many biological functions.

Ultrafast electron-transfer reactions are fundamental to protein functions. Here the authors show that these reaction dynamics are affected by the ruggedness of protein energy landscape, which even modifies the reaction free energy and reorganization energy.

Nonequilibrium dynamics of photoinduced forward and backward electron transfer reactions

The biological functions of photoenzymes are often triggered by photoinduced electron transfer (ET) reactions. An ultrafast backward ET (BET) reaction follows the initial photoinduced forward ET (FET), which dissipates the energy of absorbed photons and terminates the biological function in vain. Based upon our previous works, we reasoned that the dynamics of the BET is coupled with that of the FET and other local motions. In this work, the dynamics of the FET and BET is modeled as the master equation of the reduced density operator of a three-state system coupled with a classical harmonic reservoir. The coupling of the FET and BET is reflected in the time-evolution of the charge-transfer state's population, which is generated by a source, the reaction flux for the FET, and annihilated by a sink, the reaction flux for the BET. Surprisingly, numerical simulations show that when the BET is in the Marcus normal region, the BET can be accelerated by nonequilibrium local motions and becomes faster than what is predicted by the Marcus theory. The experimental confirmation of this novel dynamics would provide qualitative evidence for nonequilibrium effects on ultrafast ET dynamics. Additionally, the effects of quantum vibrational modes on the dynamics are discussed. This work can help understand the dynamical interactions between the chain of ultrafast reactions and the complex local environmental motions, revealing the physical nature underlying biological functions.

Light-Induced Nanosecond Relaxation Dynamics of Rhenium-Labeled Pseudomonas aeruginosa Azurins

Time-resolved phosphorescence spectra of Re(CO)₃(dmp)⁺ and Re(CO)₃(phen)⁺ chromophores (dmp = 4,7-dimethyl-1,10-phenanthroline, phen = 1,10-phenanthroline) bound to surface histidines (H83, H124, and H126) of Pseudomonas aeruginosa azurin mutants exhibit dynamic band maxima shifts to lower wavenumbers following 3-exponential kinetics with 1-5 and 20-100 ns major phases and a 1.1-2.5 μs minor (5-16%) phase. Observation of slow relaxation components was made possible by using an organometallic Re chromophore as a probe whose long phosphorescence lifetime extends the observation window up to ∼3 μs. Integrated emission-band areas also decay with 2- or 3-exponential kinetics; the faster decay phase(s) is relaxation-related, whereas the slowest one [360-680 ns (dmp); 90-140 ns (phen)] arises mainly from population decay. As a result of shifting bands, the emission intensity decay kinetics depend on the detection wavelength. Detailed kinetics analyses and comparisons with band-shift dynamics are needed to disentangle relaxation and population decay kinetics if they occur on comparable timescales. The dynamic phosphorescence Stokes shift in Re-azurins is caused by relaxation motions of the solvent, the protein, and solvated amino acid side chains at the Re binding site in response to chromophore electronic excitation. Comparing relaxation and decay kinetics of Re(dmp)124K122Cu^II and Re(dmp)124W122Cu^II suggests that electron transfer (ET) and relaxation motions in the W122 mutant are coupled. It follows that nanosecond and faster photo-induced ET steps in azurins (and likely other redox proteins) occur from unrelaxed systems; importantly, these reactions can be driven (or hindered) by structural and solvational dynamics.