Simultaneous analysis of ultrafast fluorescence decays of FMN binding protein and its mutated proteins by molecular dynamic simulation and electron transfer theory

Ultrafast photooxidation of protein-bound anionic flavin radicals

Flavoproteins are colored proteins involved in a large variety of biochemical reactions. They can perform photochemical reactions, which are increasingly exploited for bioengineering new protein-derived photocatalysts. In particular, light-induced reduction of the resting oxidized state of the flavin by close-lying amino acids or substrates is extensively studied. Here, we demonstrate that the reverse and previously unknown reaction photooxidation of the anionic semireduced flavin radical, a short-lived reaction intermediate in many biochemical reactions, efficiently occurs in flavoprotein oxidases. We anticipate that this finding will allow photoreduction of external reactants and lead to exploration of novel photocatalytic pathways.

The photophysical properties of anionic semireduced flavin radicals are largely unknown despite their importance in numerous biochemical reactions. Here, we studied the photoproducts of these intrinsically unstable species in five different flavoprotein oxidases where they can be stabilized, including the well-characterized glucose oxidase. Using ultrafast absorption and fluorescence spectroscopy, we unexpectedly found that photoexcitation systematically results in the oxidation of protein-bound anionic flavin radicals on a time scale of less than ∼100 fs. The thus generated photoproducts decay back in the remarkably narrow 10- to 20-ps time range. Based on molecular dynamics and quantum mechanics computations, positively charged active-site histidine and arginine residues are proposed to be the electron acceptor candidates. Altogether, we established that, in addition to the commonly known and extensively studied photoreduction of oxidized flavins in flavoproteins, the reverse process (i.e., the photooxidation of anionic flavin radicals) can also occur. We propose that this process may constitute an excited-state deactivation pathway for protein-bound anionic flavin radicals in general. This hitherto undocumented photochemical reaction in flavoproteins further extends the family of flavin photocycles.

Characterization of Light-Induced, Short-Lived Interacting Radicals in the Active Site of Flavoprotein Ferredoxin-NADP+ Oxidoreductase

Radicals of flavin adenine dinucleotide (FAD), as well as tyrosine and tryptophan, are widely involved as key reactive intermediates during electron-transfer (ET) reactions in flavoproteins. Due to the high reactivity of these species and their corresponding short lifetime, characterization of these intermediates in functional processes of flavoproteins is usually challenging but can be achieved by ultrafast spectroscopic studies of light-activatable flavoproteins. In ferredoxin-NADP⁺ oxidoreductase from Bacillus subtilis (BsFNR), fluorescence of the FAD cofactor that very closely interacts with a neighboring tyrosine residue (Tyr50) is strongly quenched. Here we study short-lived photoproducts of this enzyme and its variants, with Tyr50 replaced by tryptophan or glycine. Using time-resolved fluorescence and absorption spectroscopies, we show that, upon the excitation of WT BsFNR, ultrafast ET from Tyr50 to the excited FAD cofactor occurs in ∼260 fs, an order of magnitude faster than the decay by charge recombination, facilitating the characterization of the reaction intermediates in the charge-separated state with respect to other recently studied systems. These studies are corroborated by experiments on the Y50W mutant protein, which yield photoproducts qualitatively similar to those observed in other tryptophan-bearing flavoproteins. By combining the experimental results with molecular dynamics simulations and quantum mechanics calculations, we investigate in detail the effects of protein environment and relaxations on the spectral properties of those radical intermediates and demonstrate that the spectral features of radical anionic FAD are highly sensitive to its environment, and in particular to the dynamics and nature of the counterions formed in the photoproducts. Altogether, comprehensive characterizations are provided for important radical intermediates that are generally involved in functional processes of flavoproteins.

Comparative studies on picosecond-resolved fluorescence of d-amino acid oxidases from human with one from porcine kidney. Photoinduced electron transfer from aromatic amino acids to the excited flavin

Fluorescence dynamics of human d-amino acid oxidase (hDAAO) and its five inhibitors have been studied in the picoseconds time domain, and compared with one in d-amino acid oxidase from porcine kidney (pkDAAO) reported. The fluorescence lifetimes were identified as 47 ps in the dimer, 235 ps in the monomer, which are compared with those of pkDAAO (45 ps-185 ps). The fluorescence lifetimes of the hDAAO did not change upon the inhibitor bindings despite of modifications in the absorption spectra. This indicates that the lifetimes of the complexes are too short to detect with the picosecond lifetime instrument. Numbers of the aromatic amino acids are similar between the both DAAOs. The fluorescence lifetimes of hDAAO were analysed with an ET theory using the crystal structure. The difference in the lifetimes of the dimer and monomer was well described in terms of difference in the electron affinity of the excited isoalloxazine (Iso*) between the two forms of the protein, though it is not known whether the structure of the monomer is different from the dimer. Three fastest ET donors were Tyr314, Trp52 and Tyr224 in the dimer, while Tyr314, Tyr224 and Tyr55 in the monomer, which are compared to those in pkDAAO, Tyr314, Tyr224 and Tyr228 in the dimer, and Tyr224, Tyr314 and Tyr228 in the monomer. The ET rate from Trp55 in hDAAO dimer was much faster compared to the rate in pkDAAO dimer. A rise component with negative pre-exponential factor was not observed in hDAAO, which are found in pkDAAO.

Conformational difference between two subunits in flavin mononucleotide binding protein dimers from Desulfovibrio vulgaris (MF): molecular dynamics simulation

The structural and dynamical properties of five FMN binding protein (FBP) dimers, WT (wild type), E13K (Glu13 replaced by Lys), E13R (Glu13 replaced by Arg), E13T (Glu13 replaced by Thr) and E13Q (Glu13 replaced by Gln), were investigated using a method of molecular dynamics simulation (MDS). In crystal structures, subunit A (Sub A) and subunit B (Sub B) were almost completely equivalent in all of the five FBP dimers. However, the predicted MDS structures of the two subunits were not equivalent in solution, revealed by the distances and inter-planar angles between isoalloxazine (Iso) and aromatic amino acids (Trp32, Tyr35 and Trp106) as well as the hydrogen bonding pairs between Iso and nearby amino acids. Residue root of mean square fluctuations (RMSF) also displayed considerable differences between Sub A and Sub B and in the five FBP dimers. The dynamics of the whole protein structures were examined with the distance (RNN) between the peptide N atom of the N terminal (Met1) and the peptide N atom of the C terminal (Leu122). Water molecules were rarely accessible to Iso in all FBP dimers which are in contrast with other flavoenzymes.

Classification of the mechanisms of photoinduced electron transfer from aromatic amino acids to the excited flavins in flavoproteins

Emission wavelength-dependence of the relationship between logarithmic ET rate vs. donor–acceptor distance in pyranose 2-oxidase.

Photoinduced electron transfer modeling to simulate flavoprotein fluorescence decay

A method of analysis is described on the photoinduced electron transfer (PET) from aromatic amino acids as tryptophans (Trp) and tyrosines (Tyr) to the excited isoalloxazine (Iso*) in FMN-binding proteins (FBP) from Desulfovibrio vulgaris (strain, Miyazaki F). Time-dependent geometrical factors as the donor-acceptor distances are determined by means of a molecular dynamics simulation (MDS) of the proteins. Fluorescence decays of the single mutated isoforms of FBP are used as experimental data. The electrostatic (ES) energy between the photoproducts and ionic groups in the proteins is introduced into the Kakitani and Mataga (KM) model, which is modeled for an electron transfer process in solution. The PET parameters contained in the KM rate are determined by means of a nonlinear least square method, according to the Marquardt algorithm. The agreement between the observed and calculated decays is quite good, but not optimal. Characteristics on PET in flavoproteins, obtained by the present method, are described. Possible improvements of the method are discussed.

Theoretical analyses of the fluorescence lifetimes of the d-amino acid oxidase–benzoate complex dimer from porcine kidney: molecular dynamics simulation and photoinduced electron transfer

The mechanism of photoinduced electron transfer from benzoate and aromatic amino acids to the excited isoalloxazine in the d-amino acid oxidase–benzoate complex dimer was studied using molecular dynamics simulation and an electron transfer theory.

Femtosecond dynamics of short-range protein electron transfer in flavodoxin

Intraprotein electron transfer (ET) in flavoproteins is important for understanding the correlation of their redox, configuration, and reactivity at the active site. Here, we used oxidized flavodoxin as a model system and report our complete characterization of a photoinduced redox cycle from the initial charge separation in 135-340 fs to subsequent charge recombination in 0.95-1.6 ps and to the final cooling relaxation of the product(s) in 2.5-4.3 ps. With 11 mutations at the active site, we observed that these ultrafast ET dynamics, much faster than active-site relaxation, mainly depend on the reduction potentials of the electron donors with minor changes caused by mutations, reflecting a highly localized ET reaction between the stacked donor and acceptor at a van der Waals distance and leading to a gas-phase type of bimolecular ET reaction confined in the active-site nanospace. Significantly, these ultrafast ET reactions ensure our direct observation of vibrationally excited reaction product(s), suggesting that the back ET barrier is effectively reduced because of the decrease in the total free energy in the Marcus inverted region, leading to the accelerated charge recombination. Such vibrationally coupled charge recombination should be a general feature of flavoproteins with similar configurations and interactions between the cofactor flavin and neighboring aromatic residues.

Effects of the disappearance of one charge on ultrafast fluorescence dynamics of the FMN binding protein

Crystal structures of E13T (Glu13 was replaced by Thr13) and E13Q (Glu13 was replaced by Gln13) FMN binding proteins (FMN-bp) from Desulfovibrio vulgaris, strain Miyazaki F, were determined by the X-ray diffraction method. Geometrical factors related to photoinduced electron transfer from Trp32, Tyr35, and Trp106 to the excited isoalloxazine (Iso*) were compared among the three forms of FMN-bp. The rate of ET is considered to be fastest from Trp32 to Iso* in FMN-bp and then from Tyr35 and Trp106. The distances between Iso and Trp32 did not change appreciably (0.705-0.712 nm) among WT, E13T, and E13Q FMN-bps, though the distances between Iso and Tyr35 or Trp106 became a little shorter by ca. 0.01 nm in both mutated FMN-bps. The distances between the residue at 13 and the ET donors or acceptor in the mutated proteins, however, changed markedly, compared to WT. Hydrogen bonding pairs and distances between Iso and surrounding amino acids were not modified when Glu13 was replaced by Thr13 or Gln13. Effects of elimination of ionic charge at Glu13 on the ultrafast fluorescence dynamics in E13T and E13Q were investigated comparing to WT, by means of a fluorescence up-conversion method. Fluorescence lifetimes were tau(1) = 107 fs (alpha(1) = 0.86), tau(2) = 475 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.02) in E13T and tau(1) = 134 fs (alpha(1) = 0.85), alpha(2) = 746 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.03) in E13Q, which are compared to the reported lifetimes in WT, tau(1) = 168 fs (alpha(1) = 0.95) and alpha(2) = 1.4 ps (alpha(2) = 0.05). Average lifetimes (tau(AV) = Sigma(i=1)(2or3)alpha(i)tau(i)) were 0.75 ps in E13T, 1.10 ps in E13Q, and 0.23 ps in WT, which implies that tau(AV) was 3.3 times longer in E13T and 4.8 times longer in E13Q, compared to WT. The ultrafast fluorescence dynamics of WT did not change when solvent changed from H(2)O to D(2)O. Static ET rates (inverse of average lifetimes) were analyzed with static structures of the three systems of FMN-bp. Net electrostatic (ES) energies of Iso and Trp32, on which ET rates depend, were 0.0263 eV in WT, 0.322 eV in E13T, and 0.412 eV in E13Q. The calculated ET rates were in excellent agreement with the observed ones in all systems.