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.

Physical quantity of residue electrostatic energy in flavin mononucleotide binding protein dimer.

The electrostatic (ES) energy of each residue was for the first time quantitatively evaluated in a flavin mononucleotide binding protein (FBP). A residue electrostatic energy (RES) was obtained as the sum of the ES energies between atoms in each residue and all other atoms in the FBP dimer using atomic coordinates obtained by a molecular dynamics (MD) simulation. ES is one of the most important energies among the interaction energies in a protein. It is determined from the RES, the residues which mainly contribute to stabilize the structure of each subunit, and the binding energy between two subunits can be estimated. The RES of all residues in subunit A (Sub A) and subunit B (Sub B) were attractive forces, even though the residues contain net negative or positive charges. This reveals that the ES energies of any of the residues can contribute to stabilize the protein structure. The total binding ES energy over all residues among the subunits was distributed between -0.2 to -1.2 eV (mean = -0.67 eV) from the MD simulation time.

Dynamics of the protein structure of T169S pyranose 2-oxidase in solution: Molecular dynamics simulation.

Pyranose 2-oxidase (P2O) from Trametes multicolor contains FAD as cofactor, and forms a tetramer. The protein structure of a mutated P2O, T169S (Thr169 is replaced by Ser), in solution was studied by means of molecular dynamics simulation and analyses of photoinduced electron transfer (ET) from Trp168 to excited isoalloxazine (Iso*), and was compared with wild type (WT) P2O. Hydrogen bonding between Iso and nearby amino acids was very similar as between T169S and WT protein. Distances between Iso and Tyr456 were extremely heterogeneous among the subunits, 1.7 (1.5 in WT) in subunit A (Sub A), 0.97 (2.2 in WT) in Sub B, 1.3 (2.1 in WT) in Sub C, 1.3 nm (2.0 in WT) in Sub D. Mean values of root of mean square fluctuation over all residues were greater by four times than those in WT. This suggests that the protein structure of T169S is much more flexible than that of WT. Electrostatic (ES) energies between Iso anion in one subunit and ionic groups in the entire protein were evaluated. It was found that more than 50% of the total ES energy in each subunit is contributed from other subunits. Reported fluorescence decays were analyzed by a method as WT, previously reported. Electron affinities of Iso* in T169S were appreciably higher than those in WT. Static dielectric constants near Iso and Trp168 were also quite higher in T169S than those in WT.

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

In many flavoproteins photoinduced electron transfer (ET) efficiently takes place from aromatic amino acids such as tryptophan or tyrosine to the excited isoalloxazine, so that the fluorescence lifetimes of isoalloxazine in some flavoproteins become ultrashort. The mechanism of ET in the flavoproteins was classified into four classes from the relationship between logarithmic ET rates (ln Rate) and the donor-acceptor distances (Rc), using reported data. The physical quantity, GT, is defined as the sum of solvent reorganization energy, electrostatic energy between a donor cation and an Iso anion, the standard free energy gap between the photoproducts and reactants, and net electrostatic energy between the photoproducts and other ionic groups in the flavoproteins (NetES). When GT fluctuates around zero with Rc, the ET rate becomes fastest (faster than 1 ps(-1)) in Kakitani and Mataga rates. In the ultrafast ET processes, the ln Rate becomes a parabolic function (category 1) of Rc as in FMN binding proteins and pyranose 2-oxidase at the shorter emission wavelengths, when NetES is negligible compared to the other quantities in the GT function. In the ultrafast ET processes, the ln Rate does not display any clear function of Rc (category 2) when NetES is dominant in the GT function, because of no direct relation between NetES and Rc. ET in flavodoxin from Helicobacter pylori may be classified into category 2. When GT linearly varies with Rc around a certain positive value, the ET rates become much slower (<1 ps(-1)). In this case the ln Rate linearly decreases with Rc (category 3), as Tyr224 in d-amino acid oxidase dimers. It is also conceivable that the ln Rate decreases with much scattered function of Rc (category 4), when NetES is dominant in the GT function, as Tyr314 in d-amino acid oxidase dimers. In ET processes of category 1, ET rates decrease as Rc becomes shorter than the distance at the maximum values of ln Rates, where GT is negative. Conditions and physical meanings were discussed for the GT-negative region.

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.

Non-equivalent conformations of D-amino acid oxidase dimer from porcine kidney between the two subunits. Molecular dynamics simulation and photoinduced electron transfer.

The structural difference between two subunits of D-amino acid oxidase dimer from porcine kidney was studied by molecular dynamics simulation (MDS) and rate of photoinduced electron transfer (ET) from aromatic amino acids as tyrosines (Tyr) and tryptophanes (Trp) to the excited isoalloxazine (Iso*). The donor-acceptor distances (Rc) between isoalloxazine (Iso) and the donors were shortest in Tyr224 (0.74 nm) in Sub A at 10 °C (Sub A10), in Tyr224 (0.79 nm) in Sub B at 10 °C (Sub B10), in Tyr228 (0.85 nm) in Sub A at 30 °C (Sub A30), and in Tyr224 (0.72 nm) in Sub B at 30 °C (Sub B30). The Rcs were mostly shorter in the dimer than those in the monomer. Hydrogen bonding (H-bond) pairs between Iso and surrounding amino acids varied with the subunit and temperature. O2 of the Iso ring formed an H-bond exclusively with Thr317OG1 (side chain) in both Sub A10 and Sub A30, while it formed with Gly315N (peptide), Leu316N and Thr317N in Sub B10 and Sub B30. N3H of Iso formed an H-bond with Leu51O (peptide) in Sub A10 and Sub A30, but not in Sub B10 and Sub B30. Electron affinity of Iso* was appreciably lower in Sub A10 compared to Sub B10, while it was opposite at 30 °C. ET rate to Iso* was fastest from Tyr224 in Sub A10, while it was fastest from Tyr314 in Sub B10. The ET rate was fastest from Tyr314 in Sub A30, while it was fastest from Tyr224 in Sub B30. The greater ET rates in the dimer as compared to those in the monomer were elucidated with shorter Rc in the dimer as compared to the monomer. The static dielectric constants inside the subunits and the static dielectric constant between Iso and Tyr224 or Tyr228 were not different appreciably. A few water molecules and sometimes an amino acid were located between Iso and Tyr224, which may be the reason why the dielectric constant of the entire subunits did not differ from that between Iso and Tyr224.

Structural basis for the temperature-induced transition of D-amino acid oxidase from pig kidney revealed by molecular dynamic simulation and photo-induced electron transfer.

The structural basis for the temperature-induced transition in the D-amino acid oxidase (DAAO) monomer from pig kidney was studied by means of molecular dynamic simulations (MDS). The center to center (Rc) distances between the isoalloxazine ring (Iso) and all aromatic amino acids (Trp and Tyr) were calculated at 10 °C and 30 °C. Rc was shortest in Tyr224 (0.82 and 0.88 nm at 10 and 30 °C, respectively), and then in Tyr228. Hydrogen bonding (H-bond) formed between the Iso N1 and Gly315 N (peptide), between the Iso N3H and Leu51 O (peptide) and between the Iso N5 and Ala49 N (peptide) at 10 °C, whilst no H-bond was formed at the Iso N1 and Iso N3H at 30 °C. The H-bond of Iso O4 with Leu51 N (peptide) at 10 °C switched to that with Ala49 N (peptide) at 30 °C. The reported fluorescence lifetimes (228 and 182 ps at 10 and 30 °C, respectively) of DAAO were analyzed with Kakitani and Mataga (KM) ET theory. The calculated fluorescence lifetimes displayed an excellent agreement with the observed lifetimes. The ET rate was fastest from Tyr224 to the excited Iso (Iso*) at 10 °C and from Tyr314 at 30 °C, despite the fact that the Rc was shortest between Iso and Tyr224 at both temperatures. This was explained by the electrostatic energy in the protein. The differences in the observed fluorescence lifetimes at 10 and 30 °C were ascribed to the differences in electron affinity of the Iso* at both temperatures, in which the free energies of the electron affinity of Iso* at 10 and 30 °C were -8.69 eV and -8.51 eV respectively. The other physical quantities related to ET did not differ appreciably at both temperatures. The electron affinities at both temperatures were calculated with a semi-empirical molecular orbital method (MO) of PM6. Mean calculated electron affinities over 100 snapshots with 0.1 ps intervals were -7.69 eV at 10 °C and -7.59 eV at 30 °C. The difference in the calculated electron affinities, -0.11 eV, was close to the observed difference in the free energies, -0.18 eV. The present quantitative analysis predicts that the highest ET rate can occur from a donor with longer donor-acceptor distance, which was explained by differences in electrostatic energy.