Next Article in Journal
In Vitro α-Amylase and α-Glucosidase Inhibitory Effects, Antioxidant Activities, and Lutein Content of Nine Different Cultivars of Marigold Flowers (Tagetes spp.)
Next Article in Special Issue
Electrodeposition of High-Quality Ni/SiC Composite Coatings by Using Binary Non-Ionic Surfactants
Previous Article in Journal
Liquid-Phase Dehydration of Glycerol to Acrolein with ZSM-5-Based Catalysts in the Presence of a Dispersing Agent
Previous Article in Special Issue
Intramolecular and Intermolecular Interaction Switching in the Aggregates of Perylene Diimide Trimer: Effect of Hydrophobicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systematic Theoretical Study on the pH-Dependent Absorption and Fluorescence Spectra of Flavins

1
Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
2
Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(8), 3315; https://doi.org/10.3390/molecules28083315
Submission received: 6 February 2023 / Revised: 29 March 2023 / Accepted: 7 April 2023 / Published: 8 April 2023
(This article belongs to the Special Issue Exclusive Feature Papers in Physical Chemistry)

Abstract

:
Flavins are a class of organic compounds with the basic structure of 7,8-dimethy-10-alkyl isoalloxazine. They are ubiquitous in nature and participate in many biochemical reactions. Due to various existing forms, there is a lack of systematic research on the absorption and fluorescence spectra of flavins. In this study, employing the density functional theory (DFT) and time-dependent (TD) DFT, we calculated the pH-dependent absorption and fluorescence spectra of flavin of three redox states (quinone, semiquinone, and hydroquinone) in solvents. The chemical equilibrium of three redox states of flavins and the pH effect on the absorption spectra and fluorescence spectra of flavins were carefully discussed. The conclusion helps with identifying the existing forms of flavins in solvent with different pH values.
Keywords:
flavin; pH; spectra; DFT

Graphical Abstract

1. Introduction

Flavins are highly active molecules used by a variety of enzymes to perform a number of biological functions, including biocatalysis [1,2,3], biosynthesis [4], bioluminescence [5,6,7], and so on. Among flavins, flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are two common flavin molecules with enzymatic activity [8]. Their basic structures are 7, 8-dimethyl isoalloxazine, which has the potential to transfer single and double electrons. Therefore, flavins exist in three different redox states: an oxidized state (quinone), one-electron-reduced state (semiquinone), and two-electron-reduced state (hydroquinone). For quinone, semiquinone, or hydroquinone, acid-base equilibrium involves a cationic form, neutral form, and anionic form. Notably, the cationic quinone form in the ground (S0) and the first singlet excited (S1) states have different protonated species, which are cationic oxidized form F L H o x + _ N 1 and F L H o x + _ N 5 , respectively [9]. The structures and acid-base equilibrium of flavins are shown in Figure 1.
Flavins have been widely studied in various forms due to their versatility [10,11,12,13]. Different forms of flavins have characteristic UV-visible absorption spectra, which can be inversely used to identify flavin. The UV-visible absorption spectra of flavins have been studied in gas [14], solutions [14,15,16,17], and protein [12,18]. However, due to the instability of flavin semiquinone, generally the absorption spectra of only flavin quinone or flavin hydroquinone have been studied in experimental studies [15,16]. The absorption spectra of five flavin forms existing under the physiological conditions shown in Figure 1 were studied theoretically [17]. The fluorescence spectra of flavins have also been studied as another tool to characterize their photophysical properties [19,20,21,22]. Flavin quinone is highly fluorescent, and the fluorescence of FAD can be used for autofluorescence imaging to monitor for subcellular activity [22]. The anionic semiquinone F M N in the nitronate monooxygenase exhibits no fluorescence, whereas the neutral semiquinone radical F M N H shows a relatively strong fluorescence [23]. Flavin hydroquinone shows weak fluorescence in solution, but is enhanced in a rigid protein environment [12]. To sum up, previous studies on the absorption spectra and fluorescence spectra of flavin have involved only limited forms of flavin. It is necessary to study the absorption and fluorescence spectra of other forms of flavin. Therefore, there is still a lack of systematic research on the absorption spectra and fluorescence spectra of various forms of flavins. In this study, we systematically studied the absorption and fluorescence spectra of the ten forms of flavins shown in Figure 1 in aqueous solution. In addition, the chromophore of the reduced flavin was found to be very sensitive to variations in the pH of the solvent [12]. To consider the pH effect on spectra of flavins, the p K a of the S0-state and p K a * of the S1-state flavin were calculated, which was helpful to confirm whether the flavin was a photoacid. For photoacids, the p K a * was considerably less than the pKa. Moreover, based on the p K a and p K a * of each flavin, the chemical equilibrium of the three-redox state of the flavin was studied, and pH-dependent absorption spectra and fluorescence spectra were simulated.

2. Computational Details

Quantum mechanics (QM) calculations. The density functional theory (DFT) [24,25] and the time-dependent (TD) DFT [26,27] were employed. Functionals, basis sets, and solvent models were tested based on the reaction ( H 2 F L r e d H F L r e d + H + ) with an experimental p K a of 6.5 [28,29]. As the test results in Table S1 show, the M06-2X [30] /6-31 + G** level with solvation model density (SMD) approach [31] provided sufficient precision for calculating p K a . Therefore, all calculations in the S0 state were carried out at the M06-2X/6-31 + G** computational level in this paper. In addition, the SMD approach considering the polarity and nonpolarity was adapted to simulate the solution conditions. The geometries of the ten forms of flavin in Figure 1 were optimized at the M06-2X/6-31 + G** level for the S0 state in gas and water. Gibbs free energies were obtained from the vibrational analysis. For the photophysical properties of flavin, the B3LYP is usually employed for calculating the spectra and behaves well [17,18]. For radicals and anions, a relatively big basis set should be used. Therefore, the vertical absorption spectra of nine forms of flavin in water in the S0 state were predicted at the TD B3LYP/6-311 ++G** level, and 20 excited states were involved. Solvent relaxation was not included during excitation. The geometries of nine forms of flavin in the S1 state in water were optimized, and their fluorescence spectra were predicted at the TD B3LYP/6-311 ++G** level. All the above calculations were performed by the Gaussian09 program package [32].
Calculation of the  p K a  and  p K a * . The ability of AH to deprotonate to A can be expressed by the negative logarithm of the acid constant of AH, which is p K a for the S0 state and p K a * for the S1 state. p K a is related to the dissolution free energy G a q , which can be computed by the Born–Haber thermodynamic cycle [33] (see Figure 2). The specific steps are as follows [34]:
G g a s = G g a s A H + G g a s A + G g a s H +
where G g a s refers to the variation of the Gibbs free energy of the gas phase. G g a s A H and G g a s A refer to the Gibbs free energy of AH and A, respectively, which can be obtained by vibration analysis of the optimized optimal geometric configuration of AH and A. G g a s H + refers to the free energy of the proton in gas phase (−6.28 kcal/mol) [35].
G a q = G g a s G s o l v A H + G s o l v A + G s o l v H + + G a d d
where G s o l v A H refers to the solvation free energy of AH, calculated from the variation of the single-point energy between the gas and water. Similarly, G s o l v A refers to the solvation free energy of A, calculated from the variation of the single-point energy between the gas and water. G s o l v H + is the solvation free energy of the hydrogen ion (−265.9 kcal/mol) [36], and G a d d represents the transfer of a solute molecule from the 1 atm gas phase to the 1 M solvent standard state (1.89 kcal/mol).
p K a = G a q 2.303 R T
where R refers to the molar constant of gas, which is valued at 8.314 J/(mol/K). T is the temperature, which is valued at 298.15 K.
p K a * can be obtained by the following equation [37]:
p K a * = p K a + Δ E a b s 2.303 R T
where Δ E a b s is the absorption energy variation between protonated ( A H ) and anionic ( A ) species.
Calculation of the relative concentration and the pH-dependent spectra. The relative concentration of each form of flavins in the S0 state at a specific pH can be calculated with the p K a . The calculation process is as follows:
A H 2 + A H + H + K a 1 = A H H + A H 2 +
A H A + H + K a 2 = A H + A H
where A H 2 + , A H , and A represent the cation, neutral molecule, and anion in the same redox state, respectively. K a 1 and K a 2 are the acid constant of A H 2 + and A H , respectively.
Then the concentration ratio of three forms of each redox state can be presented as follows:
[ A H 2 + ] : [ A H ] : [ A ] = 1 : [ H + ] 1 K a 1 : [ H + ] 2 K a 1 K a 2
Then the relative concentrations of the three chemical forms for each redox state at a specific pH value can be obtained, and the pH-dependent absorption spectra can be simulated. At a specific pH value, multiply the relative concentrations of the three forms by the wavelengths of the three forms of absorption spectrum to obtain the wavelength at that pH value, and multiply the relative concentrations of the three forms by the intensity of the three forms of absorption spectrum to obtain the intensity at the pH value. The same method is used for obtaining absorption spectra at other pH values. These spectra at different pH values can be integrated to obtain pH-dependent absorption spectra.
Similarly, the relative concentration of each form of flavins in the S1 state at a specific pH can be calculated with the p K a * . In addition, the pH-dependent fluorescence spectra can be simulated by the same method with the pH-dependent absorption spectra.

3. Results and Discussion

3.1. Absorption and Fluorescence Spectra of Flavins in Solution

The calculated vertical absorption spectra and fluorescence spectra of the ten forms of flavins in solution are shown in Figure 3. The TD state order, transition, absorption maximum ( λ m a x ), and oscillator strength (f) of the absorption spectra of the flavins are listed in Table S2. The full width at half-maximum (FWHM) of the spectra were determined empirically and are listed in Table S3. The spectra were normalized, no scaling factor was used, and the electronic excited states and their corresponding transition analysis next to them are shown in Figure 3. The frontier molecular orbital π 3 plays a central role for all redox forms of flavin; it mainly delocalizes over the isoalloxazine ring and exhibits a bonding character between C4a and C10a and an anti-bonding character between N10 and C10a. The π orbitals with energy higher than π 3 are π 4 , π 5 , π 6 , and those with energy lower than π 3 are π 2 , π 1 , π 0 . These definitions and sequences are consistent with previous studies [18,38].
For flavin quinone, π 3 is the lowest unoccupied molecular orbital (LUMO). In Figure 3A, the absorption spectrum of F L H o x + _ N 1 has an absorption peak, which comes from the contribution of the π 2 π 3 transition from the S0 state to the S1 state and the π 1 π 3 transition from the S0 state to the S2 state. The π 1 π 3 transition corresponds to the λ m a x at 366.0 nm, corresponding to the experimental maximum ( λ e x p ) around 394.0 nm [15]. The absorption spectrum of F L o x has two absorption peaks: the lowest-energy peak is represented by the π 2 π 3 transition from S0 to S1, and the highest-energy peak is represented by the π 1 π 3 transition from S0 to S2. The π 2 π 3 and the π 1 π 3 transitions correspond to the λ m a x values at 427.5 and 358.7 nm, corresponding to the λ e x p values around 445.0 and 370.0 nm [15,39], respectively. The absorption spectrum of F L o x also has two peaks: the lowest-energy peak is mainly represented by the π 2 π 3 transition from S0 to S2, and the highest-energy peak involves the π 1 π 3 transition from S0 to S4 and the π 0 π 3 transition from S0 to S5. The π 1 π 3 transition and π 0 π 3 transition correspond to the λ m a x values at 424.0 and 354.2 nm, corresponding to the λ e x p values of 450.0 and 350.0 nm [15], respectively. The absorption spectra of F L H o x + _ N 1 , F L o x and F L o x all reproduce the spectral characteristics (single-peaked spectrum for F L H o x + _ N 1 , double-peaked spectra for F L o x and F L o x ). In Figure 3B, the transitions from S1 to S0 for F L H o x + _ N 5 , F L o x , and F L o x are π 3 π 2 transitions, and the π 3 π 2 transition for F L o x corresponds to the fluorescence wavelength ( λ F ) at 526.0 nm, closely reproducing the experimental value around 530.0 nm [15]. The large f of flavin quinone is consistent with the experimental high bright fluorescence of oxidized flavin [39].
For flavin semiquinone, π 3 is a singly occupied molecular orbital (SOMO). In Figure 3C, the lowest-energy peak of the absorption spectrum of H 2 F L s q + mainly involves the contributions of the π 2 β π 3 β , π 3 α π 4 α , and π 1 β π 3 β transitions from the S0 state to the first three excited singlet states. The highest-energy peak mainly comes from the contributions of the π 3 α π 5 α and π 3 α π 6 α transitions from the S0 state to the S5 state and the S6 state. The triple-peaked absorption spectra of H F L s q involves many electronic transitions, where the bright π 2 β π 3 β , π 1 β π 3 β , and π 3 α π 6 α transitions are found at 568.6, 460.3, and 340.2 nm, respectively, corresponding to the λ e x p around 571.0, 481.0, and 340.0 nm, respectively [40]. For F L s q , the π 2 β π 3 β transition and the π 3 α π 6 α transition correspond to the λ m a x values at 450.8 and 363.1 nm, corresponding to the λ e x p around 480.0 and 370.0 nm, respectively [40]. There is no experimental report on the absorption spectra of H 2 F L s q + . The absorption spectra of H F L s q and of F L s q well reproduce the experimental spectral characteristics (triple-peaked spectrum for H F L s q and double-peaked spectrum for F L s q ). In Figure 3D, the transitions from S1 to S0 are the π 3 β π 2 β transitions for H 2 F L s q + and H F L s q and the π 4 α π 3 α transition for F L s q . The f of H F L s q is large, and it is much larger than that of F L s q . This is consistent with the results that flavin semiquinone yielded in the nitronate monooxygenase [23]. Due to flavin semiquinone not being stable in solution and being a transient species in the excited state, there is no experimental study on the fluorescence spectra of flavin semiquinone in aqueous.
For flavin hydroquinone, π 3 is the highest occupied molecular orbital (HOMO). For H 2 F L H r e d + , the λ m a x corresponding to the π 3 π 4 transition from S0 to S1 is 295.6 nm. Therefore, there is no absorption peak for H 2 F L H r e d + in Figure 3E. In Figure 3E, the absorption spectra of both H 2 F L r e d and H F L r e d have a peak [39]. For H 2 F L r e d , the absorption peak comes from the contribution of the π 3 π 4 transition, which corresponds to the λ m a x at 392.4 nm, matching the λ e x p around 395.0 nm. For H F L r e d , the π 3 π 5 transition corresponds to the λ m a x at 330.4 nm, corresponding to the λ e x p around 342.0 nm [39]. In Figure 3F, the transitions from S1 to S0 for H 2 F L H r e d + , H 2 F L r e d , and H F L r e d are π 4 π 3 transitions, which correspond to the λ F values at 496.9, 626.8, and 563.7 nm, respectively. H 2 F L H r e d + does not exist in the pH range of 1~14. In addition, the f of H 2 F L r e d and H F L r e d is low. All these are consistent with the experimental result of flavin hydroquinone that λ F is not below 500 nm and essentially nonfluorescent in aqueous solutions [39].
In summary, our calculation results can almost reproduce the existing experimental values. As listed in Rajiv K. Kar’s review [41], some theoretical studies have calculated the vertical excitation energy ( E e x ) of flavins [7,8,42,43,44]. The E e x and f of the five flavin forms ( F L o x , H F L s q , F L s q , H 2 F L r e d , and H F L r e d ) in aqueous solution are calculated at the level of TD B3LYP/def2-TZVPP [38] and TD B3LYP/cc-pVTZ [17]. The comparison of our calculation results with their reported results and the experimental values are listed in Table 1. For all five forms of flavin, the λ m a x values we predicted were closer to the λ e x p than those predicted in Ref. [17] and Ref [38]. For F L o x , H 2 F L r e d , and H F L r e d , the transitions corresponding to the λ e x p were completely consistent with Ref. [38]. For H F L s q and F L s q , the transitions corresponding to the λ e x p of 340.0 and 370.0 nm involved in the highest peak of spectra were π 3 α π 6 α transitions, which were different from the π 0 β π 3 β transition and the π 1 β π 3 β transition reported in Ref. [38]. The π 0 β π 3 β transition for H F L s q and the π 1 β π 3 β transition for F L s q in Table S2 were found to have lower f than those of the π 3 α π 6 α transition. The λ m a x and f predicted in this paper, in most cases, were similar to those predicted in Ref. [17], but the corresponding excited state orders were different. All these differences may be due to the computational level and solvation model. In addition, the lowest-energy peak of the absorption spectra of F L o x , H F L s q   F L s q , a n d H F L r e d involved the contributions of the same electronic transition as those in Ref. [18].

3.2. Chemical Equilibrium of Flavins in Solution

The proton affinity of flavins is easily affected by the pH of the solution. To confirm the existence forms of flavin at different pH values, the calculated (calc.) and experimental (exp.) p K a and p K a * values are listed in Table 1. All calculated p K a and p K a * values of flavins are in good agreement with the reported experimental values. In addition, the p K a * value of F L H o x + _ N 5 in the S1 state is greater than the p K a value of F L H o x + _ N 1 in the S0 state. This shows that the acidity of the oxidized cationic form in the S1 state is weaker than that in the S0 state, which is consistent with the previous research results [9]. This can indicate that the cationic quinone is not a photoacid but a photobase. From Table 2, we can see that the p K a * values of H 2 F L s q + and H 2 F L H r e d + are less than their p K a values, which indicates that they are photoacids. By contrast, the p K a * values of H F L s q , H 2 F L r e d , and cationic quinone are greater than their p K a values, which indicates that they are photobases. There is little difference between the p K a * value and p K a value of F L o x . Based on the p K a and p K a * , the relative concentrations of each form of flavins in the S0 state and S1 state at different values are listed in Tables S4 and S5, respectively. In addition, the diagrams of the relative concentrations of different forms of flavins in the S0 state and S1 state with various pH values are shown in Figure 4. For the flavin quinone (Figure 4A,B), in the S0 state, F L H o x + _ N 1 almost does not exist in the entire pH range. When the pH value is less than 8, F L o x is the dominant form. When the pH is larger than 8, the relative concentration of F L o x decreases gradually with the increase of pH, and the relative concentration of F L o x increases gradually, and they are equal at the pH value of 10.8. Later, with the pH increasing, F L o x gradually becomes the dominant form. Different from the S0 state, F L H o x + _ N 5 in the S1 state exists in the pH range of 1.0~2.1. For the flavin semiquinone (Figure 4C,D), H 2 F L s q + exists when the pH is less than 2 for the S0 state but does not exist in the whole range of pH for the S1 state. For the S0 and S1 states, H F L s q are the dominant forms when the pH is less than 6 and 7, and when the pH is larger than 6 and 7, the relative concentration of H F L s q gradually decreases, and the relative concentration of F L s q gradually increases. They reach equilibrium when the pH is 8.6 and 9.2 for the S0 state and the S1 state, respectively. Then, with the increase of the pH, F L s q gradually becomes the dominant form for both the S0 and S1 states. For the flavin hydroquinone (Figure 4E,F), the cationic form H 2 F L H r e d + does not exist in the whole pH range for both the S0 and S1 states. The neutral form H 2 F L r e d and the anionic form H F L r e d are the dominant forms before and after they reach equilibrium, respectively.

3.3. The pH-Dependent Absorption and Fluorescence Spectra of Flavins in Solution

We systematically studied the dependence of the absorption and fluorescence spectra on the pH values for both the S0 and S1 states of the ten forms of flavins. The pH-dependent absorption and fluorescence spectra of flavins in solution are shown in Figure 5. In Figure 5A, the absorption spectra of flavin quinone in the pH range of 1~8 and 13~14 are consistent with the absorption spectra of F L o x and F L o x in Figure 2A, respectively. In addition, the peak positions or intensities of the spectra in the pH range of 9~12 are different from those of F L o x and F L o x , which is attributed to the coexistence of F L o x and F L o x in this pH range (see Figure 4A). In Figure 5B, the fluorescence spectra of flavin quinone in the pH range of 2~14 only have a shark peak when the pH equals to 1; the fluorescence spectrum of flavin quinone has a hump, which can be attributed to the fact that F L H o x + _ N 5 and F L o x coexist when the pH equals to 1 (see Figure 4B); and the fluorescence intensity of F L o x is stronger than that of F L H o x + _ N 5 (see Figure 3B). In the range of 8~13, the intensity of spectra gradually increases, which can be attributed to the fact that F L o x gradually becomes the dominant form (see Figure 4B). In Figure 5C, the absorption spectra of flavin semiquinone in the pH range of 2~6 and 11~14 are consistent with the absorption spectra of H F L s q and F L s q in Figure 3C, respectively. In the pH range of 1~6, the spectra present a short and wide triple-peak feature. In the pH range of 7~10, the triple-peak feature of the spectra gradually disappears with the increase of pH, because the relative concentration of H F L s q decreases, and the relative concentration of F L s q gradually increases in this pH range (see Figure 4C).The spectrum at pH equal to 1 is different from others, which can be attributed to the coexistence of H 2 F L s q + and H F L s q at pH equal to 1. In Figure 5D, the fluorescence intensity of flavin semiquinone gradually decreases as the pH increases to 7, which is attributed to F L s q gradually becoming the dominant existing form (see Figure 4D). In Figure 5E, the spectra in the pH range of 1~4 and 9~14 are consistent with the absorption spectra of H 2 F L r e d and H F L r e d in Figure 2E, respectively. In the pH range of 5~8, the shape of the spectrum gradually changes from that of H 2 F L r e d to that of H F L r e d , because the relative concentration of H 2 F L r e d decreases in this pH range, while the relative concentration of H F L r e d increases in this pH range (see Figure 4E). The fluorescence spectra of flavin hydroquinone are short and wide in Figure 5F, which indicates that there is no H 2 F L H r e d + in the pH range of 1~14 (see Figure 3F and Figure 4F). The λ F values and the fluorescence intensity of the fluorescence spectra at the pH range of 1~9 are larger and stronger than those at the pH range of 11~14, and the fluorescence spectrum have two peaks when pH is 10.

4. Conclusions

The versatility of flavins is reflected in their various existence forms. Since the available experimental spectra of flavin are limited, we systematically studied the absorption and fluorescence spectra of ten forms of flavin in solution and provided the precise assignment of the excited states in this paper. Every form of flavin has unique spectral characteristics. Most absorption peaks of the absorption spectra of flavin contain the contribution of multiple electronic transitions. In addition, the predicted λ m a x values almost reproduce the existing experimental values. The p K a and p K a * of flavin were calculated, and we found that the p K a values of H 2 F L s q + and H 2 F L H r e d + were greater than their p K a * values, indicating that they are photoacids. The chemical equilibrium between different protonated forms of flavins was obtained, and the existing forms of flavins at different pH were confirmed. The cationic flavin quinone and hydroquinone in the S0 state and the cationic flavin semiquinone and hydroquinone in the S1 state do not exist throughout the pH range. The cationic flavin quinone in the S0 state and the cationic flavin semiquinone in the S1 state exist when the pH is lower than 2. For all the redox states of flavins in both the S0 and S1 states, their dominant forms are their neutral forms before reaching chemical equilibrium, and their anionic forms become the dominant forms after reaching chemical equilibrium when the pH is in the range of 2~14. Moreover, the pH-dependent absorption and fluorescence spectra of flavin were simulated, which provided a spectral basis for determining the presence form of each flavin. For flavin quinone, the pH-dependent absorption spectra retained the double-peaked characteristics in the pH range of 1~14, and the fluorescence intensity of the pH-dependent fluorescence spectra increased as the pH increased. For flavin semiquinone, the triple-peaked absorption spectrum significantly changed to a double-peaked spectrum in the pH range of 7~9, and the fluorescence intensity of the pH-dependent fluorescence spectra decreased as the pH increased. For flavin hydroquinone, the λ m a x values of the absorption spectra underwent a blue shift in the pH range of 4~9. In addition, the fluorescence spectra were short and wide in the pH range of 1~14, and the spectrum had two peaks when the pH was 10.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/molecules28083315/s1: Table S1: The p K a values under different functional, basis set, and solvation models; Table S2: The TD state order, transition, absorption maximum ( λ m a x ), and oscillator strength (f) of absorption spectra of flavins; Table S3: FWHM of absorption and fluorescence spectra of ten forms of flavins; Table S4: The relative concentrations of each form of flavins in S0 state at different pH values; Table S5: The relative concentrations of each form of flavins in S1 state at different pH values; Table S6: The Cartesian coordinates was optimized equilibrium structures reported in this paper.

Author Contributions

Manuscript conception, Y.L.; writing—original draft preparation, J.W.; writing—review and editing, Y.L.; calculation, J.W.; methodology, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 21973005).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are available from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are available from the authors.

References

  1. Huijbers, M.M.; Montersino, S.; Westphal, A.H.; Tischler, D.; van Berkel, W.J. Flavin Dependent Monoxygenases. Arch. Biochem. Biophys. 2014, 544, 2–17. [Google Scholar] [CrossRef] [PubMed]
  2. van Berkel, W.J.; Kamerbeek, N.M.; Fraaije, M.W. Flavoprotein Monooxygenases, A Diverse Class of Oxidative Biocatalysts. J. Biotechnol. 2006, 124, 670–689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Koetter, J.W.A.; Schulz, G.E. Crystal Structure of 6-Hydroxy-d-nicotine Oxidase from Arthrobacter nicotinovorans. J. Mol. Biol. 2005, 352, 418–428. [Google Scholar] [CrossRef] [PubMed]
  4. Hassan-Abdallah, A.; Bruckner, R.C.; Zhao, G.; Jorns, M.S. Biosynthesis of Covalently Bound Flavin: Isolation and In vitro Flavinylation of the Monomeric Sarcosine Oxidase Apoprotein. Biochemistry 2005, 44, 6452–6462. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, Y.-J. Understanding the Complete Bioluminescence Cycle from a Multiscale Computational Perspective: A review. J. Photochem. Photobiol. C 2022, 52, 1. [Google Scholar] [CrossRef]
  6. O’Kane, D.J.; Prasher, D.C. Evolutionary Origins of Bacterial Bioluminescence. Mol. Microbiol. 1992, 6, 443–449. [Google Scholar] [CrossRef]
  7. Mager, H.I.; Tu, S.C. Chemical Aspects of Bioluminescence. Photochem. Photobiol. 1995, 62, 607–614. [Google Scholar] [CrossRef]
  8. Colombo, B.; Saraceno, L.; Comi, G. Riboflavin and Migraine: The Bridge Over Troubled Mitochondria. Neurol. Sci. 2014, 35, 141–144. [Google Scholar] [CrossRef]
  9. Schulman, S.G. pH Dependence of Fluorescence of Riboflavin and Related Isoalloxazine Derivatives. J. Pharm. Sci. 1971, 60, 628–631. [Google Scholar] [CrossRef]
  10. Ahmad, I.; Anwar, Z.; Iqbal, K.; Ali, S.A.; Mirza, T.; Khurshid, A.; Khurshid, A.; Arsalan, A. Effect of Acetate and Carbonate buffers on the Photolysis of Riboflavin in Aqueous Solution: A Kinetic Study. AAPS PharmSciTech 2014, 15, 550–559. [Google Scholar] [CrossRef] [Green Version]
  11. Villabona-Monsalve, J.P.; Varnavski, O.; Palfey, B.A.; Goodson, T. Two-Photon Excitation of Flavins and Flavoproteins with Classical and Quantum Light. J. Am. Chem. Soc. 2018, 140, 14562–14566. [Google Scholar] [CrossRef] [PubMed]
  12. Ghisla, S.; Massey, V.; Lhoste, J.M.; Mayhew, S.G. Fluorescence and Optical Characteristics of Reduced Flavines and Flavoproteins. Biochemistry 1974, 13, 589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Van den Berg, P.A.; Widengren, J.; Hink, M.A.; Rigler, R.; Visser, A.J. Fluorescence Correlation Spectroscopy of Flavins and Flavoenzymes: Photochemical and Photophysical Aspects. Spectrochim. Acta Part A 2001, 57, 2135–2144. [Google Scholar] [CrossRef]
  14. Ai, Y.; Zhao, C.; Xing, J.; Liu, Y.; Wang, Z.; Jin, J.; Xia, S.; Cui, G.; Wang, X. Excited-State Decay Pathways of Flavin Molecules in Five Redox Forms: The Role of Conical Intersections. J. Phys. Chem. A. 2018, 122, 7954–7961. [Google Scholar] [CrossRef] [PubMed]
  15. Tyagi, A.; Penzkofer, A. pH Dependence of the Absorption and Emission Behaviour of Lumiflavin in Aqueous Solution. J. Photochem. Photobiol. A 2010, 215, 108–117. [Google Scholar] [CrossRef]
  16. Song, S.H.; Dick, B.; Penzkofer, A.; Pokorny, R.; Batschauer, A.; Essen, L.O. Absorption and Fluorescence Spectroscopic Characterization of Cryptochrome 3 from Arabidopsis Thaliana. J. Photochem. Photobiol. B 2006, 85, 1–16. [Google Scholar] [CrossRef]
  17. Kabir, M.P.; Orozco-Gonzalez, Y.; Gozem, S. Electronic Spectra of Flavin in Different Redox and Protonation States: A Computational Perspective on the Effect of the Electrostatic Environment. Phys. Chem. Chem. Phys. 2019, 21, 16526–16537. [Google Scholar] [CrossRef]
  18. Schwinn, K.; Ferre, N.; Huix-Rotllant, M. UV-visible Absorption Spectrum of FAD and Its Reduced Forms Embedded in A Cryptochrome Protein. Phys. Chem. Chem. Phys. 2020, 22, 12447–12455. [Google Scholar] [CrossRef]
  19. Sato, K.; Nishina, Y.; Shiga, K. Decomposition of the Fluorescence Spectra of Two FAD Molecules in Electron-transferring Flavoprotein from Megasphaera elsdenii. J. Biochem. 2013, 154, 61–66. [Google Scholar] [CrossRef]
  20. Islam, M.S.; Honma, M.; Nakabayashi, T.; Kinjo, M.; Ohta, N. pH Dependence of the Fluorescence Lifetime of FAD in Solution and in Cells. Int. J. Mol. Sci. 2013, 14, 1952–1963. [Google Scholar] [CrossRef] [Green Version]
  21. Nakabayashi, T.; Islam, M.S.; Ohta, N. Fluorescence Decay Dynamics of Flavin Adenine Dinucleotide in a Mixture of Alcohol and Water in the Femtosecond and Nanosecond Time Range. J. Phys. Chem. B 2010, 114, 15254–15260. [Google Scholar] [CrossRef] [PubMed]
  22. Galban, J.; Sanz-Vicente, I.; Navarro, J.; de Marcos, S. The Intrinsic Fluorescence of FAD and Its Application in Analytical Chemistry: A Review. Methods Appl. Fluoresc. 2016, 4, 042005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Su, D.; Kabir, M.P.; Orozco-Gonzalez, Y.; Gozem, S.; Gadda, G. Fluorescence Properties of Flavin Semiquinone Radicals in Nitronate Monooxygenase. Chembiochem 2019, 20, 1646–1652. [Google Scholar] [CrossRef]
  24. Perdew, J.P.; Ruzsinszky, A.; Tao, J.; Staroverov, V.N.; Scuseria, G.E.; Csonka, G.I. Prescription for the Design and Selection of Density Functional Approximations: More Constraint Satisfaction with Fewer Fits. J. Chem. Phys. 2005, 123, 062201. [Google Scholar] [CrossRef] [PubMed]
  25. Dreuw, A.; Head-Gordon, M. Single-reference ab Initio Methods for the Calculation of Excited States of Large Molecules. Chem. Rev. 2005, 105, 4009–4037. [Google Scholar] [CrossRef] [PubMed]
  26. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Runge, E.; Gross, E.K.U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997–1000. [Google Scholar] [CrossRef]
  28. Ghisla, S.; Massey, V. New Flavins for Old: Artificial Flavins as Active Site Probes of Flavoproteins. Biochem. J. 1986, 1, 1–12. [Google Scholar] [CrossRef] [Green Version]
  29. Kılıç, M.; Ensing, B. Acidity Constants of Lumiflavin from First Principles Molecular Dynamics Simulations. Phys. Chem. Chem. Phys. 2014, 16, 18993–19000. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Y.; Truhlar, D.G. The M06 suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, noncovalent Interactions, Excited states, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2007, 120, 215–241. [Google Scholar]
  31. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal Solvation Model Based on Solute Electron Density and on A Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B. 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
  32. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2013. [Google Scholar]
  33. Haber, F. Betrachtungen zur theorie der warmetonung. Ber. Dtsch. Phys. Ges. 1919, 21, 750. [Google Scholar]
  34. Liptak, M.D.; Shields, G.C. Accurate pK(a) Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods. J. Am. Chem. Soc. 2001, 123, 7314–7319. [Google Scholar] [CrossRef] [PubMed]
  35. Rashin, A.A.; Topol, I.A.; Tawa, G.J.; Burt, S.K. Charge Distributions in Water and Ion-water Clusters. Chem. Phys. Lett. 2001, 335, 327–333. [Google Scholar] [CrossRef]
  36. Kelly, C.P.; Cramer, C.J.; Truhlar, D.G. Aqueous Solvation Free Energies of Ions and Ion-water Clusters Based on An Accurate Value for the Absolute Aqueous Solvation Free Energy of the Proton. J. Phys. Chem. B 2006, 110, 16066–16081. [Google Scholar] [CrossRef] [Green Version]
  37. Houari, Y.; Jacquemin, D.; Laurent, A.D. Methodological Keys for Accurate pKa* Simulations. Phys. Chem. Chem. Phys. 2013, 15, 11875–11882. [Google Scholar] [CrossRef]
  38. Mondal, P.; Schwinn, K.; Huix-Rotllant, M. Impact of The Redox State of Flavin Chromophores on the UV–vis Spectra, Redox and Acidity Constants and Electron Affinities. J. Photochem. Photobiol. A 2020, 387, 1. [Google Scholar] [CrossRef]
  39. McBride, R.A.; Barnard, D.T.; Jacoby-Morris, K.; Harun-Or-Rashid, M.; Stanley, R.J. Reduced Flavin in Aqueous Solution Is Nonfluorescent. Biochemistry 2023, 62, 759–769. [Google Scholar] [CrossRef]
  40. Sakai, M. One-electron Photoreduction of Flavin Mononucleotide: Time-resolved Resonance Raman and Absorption Study. J. Mol. Struct. 1996, 379, 9–18. [Google Scholar] [CrossRef]
  41. Kar, R.K.; Miller, A.F.; Mroginski, M.A. Understanding Flavin Electronic Structure and Spectra. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 12, 2. [Google Scholar] [CrossRef]
  42. Choe, Y.K.; Nagase, S.; Nishimoto, K. Theoretical Study of the Electronic Spectra of Oxidized and Reduced States of Lumiflavin and Its derivative. J. Comput. Chem. 2007, 28, 727–739. [Google Scholar] [CrossRef] [PubMed]
  43. Kar, R.K.; Borin, V.A.; Ding, Y.; Matysik, J.; Schapiro, I. Spectroscopic Properties of Lumiflavin: A Quantum Chemical Study. Photochem. Photobiol. 2019, 95, 662–674. [Google Scholar] [CrossRef] [PubMed]
  44. Brisker-Klaiman, D.; Dreuw, A. On the Influence of Dimerisation of Lumiflavin in Aqueous Solution on Its Optical Spectra—A Quantum Chemical Study. Mol. Phys. 2019, 117, 2167–2178. [Google Scholar] [CrossRef]
  45. Land, E.J.; Swallow, A.J. One-electron Reactions in Biochemical Systems as Studied by Pulse Radiolysis. II. Riboflavin. Biochemistry 1969, 5, 2117. [Google Scholar] [CrossRef]
Figure 1. Structures and acid-base equilibrium of flavins and labels. The molecules in the black dotted box are the forms that exist under physiological conditions.
Figure 1. Structures and acid-base equilibrium of flavins and labels. The molecules in the black dotted box are the forms that exist under physiological conditions.
Molecules 28 03315 g001
Figure 2. Born–Haber thermodynamic cycle.
Figure 2. Born–Haber thermodynamic cycle.
Molecules 28 03315 g002
Figure 3. The computational absorption spectra (A,C,E) and fluorescence spectra (B,D,F) of flavin quinone, flavin semiquinone, and flavin hydroquinone in aqueous solution. Thin vertical lines represent the electronic excited states, with the corresponding transition analysis next to them. Values in parentheses represent the f of fluorescence spectra.
Figure 3. The computational absorption spectra (A,C,E) and fluorescence spectra (B,D,F) of flavin quinone, flavin semiquinone, and flavin hydroquinone in aqueous solution. Thin vertical lines represent the electronic excited states, with the corresponding transition analysis next to them. Values in parentheses represent the f of fluorescence spectra.
Molecules 28 03315 g003
Figure 4. The diagrams of relative concentrations for the S0 state (A,C,E) and S1 state (B,D,F) of different forms of flavins in aqueous solution at different pH values.
Figure 4. The diagrams of relative concentrations for the S0 state (A,C,E) and S1 state (B,D,F) of different forms of flavins in aqueous solution at different pH values.
Molecules 28 03315 g004
Figure 5. The pH-dependent absorption (A,C,E) and fluorescence (B,D,F) spectra of flavins.
Figure 5. The pH-dependent absorption (A,C,E) and fluorescence (B,D,F) spectra of flavins.
Molecules 28 03315 g005
Table 1. Comparison of absorption spectra of flavin between this paper and previous theoretical studies. Experimental results for F L o x were taken from Refs. [15,39]; for H F L s q and F L s q from Ref. [40]; and for H 2 F L r e d and H F L r e d from Ref. [39].
Table 1. Comparison of absorption spectra of flavin between this paper and previous theoretical studies. Experimental results for F L o x were taken from Refs. [15,39]; for H F L s q and F L s q from Ref. [40]; and for H 2 F L r e d and H F L r e d from Ref. [39].
TD State OrderTransition λ m a x f λ e x p
F L o x This paper1 π 2 π 3 427.50.22445.0
Ref. [38] π 2 π 3 391.00.25
Ref. [17]1 422.70.24
This paper2 π 1 π 3 358.70.31370.0
Ref. [38] π 1 π 3 326.00.26
Ref. [17]4 345.30.25
H F L s q This paper1 π 2 β π 3 β 568.60.11571.0
Ref. [38] π 2 β π 3 β 535.00.13
Ref. [17]1 581.00.13
This paper2 π 1 β π 3 β 460.30.06485.0
Ref. [38] π 1 β π 3 β 406.00.06
Ref. [17]3 431.30.06
This paper6 π 3 α π 6 α 340.20.08340.0
Ref. [38] π 0 β π 3 β 296.00.11
Ref. [17]5 360.30.09
F L s q This paper2 π 2 β π 3 β 450.80.04480.0
Ref. [38] π 2 β π 3 β 423.00.13
Ref. [17]3 437.50.14
This paper7 π 3 α π 6 α 363.10.37370.0
Ref. [38] π 1 β π 3 β 359.00.101
Ref. [17]6 357.60.297
H 2 F L r e d This paper1 π 3 π 4 392.40.06395.0
Ref. [38] π 3 π 4 400.00.03
Ref. [17]1 411.40.03
H F L r e d This paper2 π 3 π 5 330.40.08342.0
Ref. [38] π 3 π 5 347.00.12
Ref. [17]2 345.40.13
Table 2. The p K a and p K a * of the flavins.
Table 2. The p K a and p K a * of the flavins.
Form p K a p K a *
calc.exp.calc.exp.
F L H o x + _ N 1 −2.80.0 [9]−4.9-
F L H o x + _ N 5 −10.1-0.21.7 [9]
F L o x 10.810.8 [15]10.510.8 [15]
H 2 F L s q + 0.12.3 [45]−3.1-
H F L s q 8.68.5 [28,29]9.2-
H 2 F L H r e d + −3.6-−21.1-
H 2 F L r e d 6.66.5 [28,29]9.7-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, J.; Liu, Y. Systematic Theoretical Study on the pH-Dependent Absorption and Fluorescence Spectra of Flavins. Molecules 2023, 28, 3315. https://doi.org/10.3390/molecules28083315

AMA Style

Wang J, Liu Y. Systematic Theoretical Study on the pH-Dependent Absorption and Fluorescence Spectra of Flavins. Molecules. 2023; 28(8):3315. https://doi.org/10.3390/molecules28083315

Chicago/Turabian Style

Wang, Jinyu, and Yajun Liu. 2023. "Systematic Theoretical Study on the pH-Dependent Absorption and Fluorescence Spectra of Flavins" Molecules 28, no. 8: 3315. https://doi.org/10.3390/molecules28083315

APA Style

Wang, J., & Liu, Y. (2023). Systematic Theoretical Study on the pH-Dependent Absorption and Fluorescence Spectra of Flavins. Molecules, 28(8), 3315. https://doi.org/10.3390/molecules28083315

Article Metrics

Back to TopTop