Understanding flavin electronic structure and spectra

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

Exploring the Reactivity of Flavins with Nucleophiles Using a Theoretical and Experimental Approach

Covalent adducts of flavin cofactors with nucleophiles play an important role in non-canonical function of flavoenzymes as well as in flavin-based catalysis. Herein, the interaction of flavin derivatives including substituted flavins (isoalloxazines), 1,10-ethylene-bridged flavinium salts, and non-substituted alloxazine and deazaflavin with selected nucleophiles was investigated using an experimental and computational approach. Triphenylphosphine or trimethylphosphine, 1-nitroethan-1-ide, and methoxide were selected as representatives of neutral soft, anionic soft, and hard nucleophiles, respectively. The interactions were investigated using UV/Vis and 1H NMR spectroscopy as well as by DFT calculations. The position of nucleophilic attack estimated using the calculated Gibbs free energy values was found to correspond with the experimental data, favouring the addition of phosphine and 1-nitroethan-1-ide into position N(5) and methoxide into position C(10a) of 1,10-ethylene-bridged flavinium salts. The calculated Gibbs free energy values were found to correlate with the experimental redox potentials of the flavin derivatives tested. These findings can be utilized as valuable tools for the design of artificial flavin-based catalytic systems or investigating the mechanism of flavoenzymes.

Oxidation of Flavin by Molecular Oxygen: Computational Insights into a Possible Radical Mechanism

As a highly electrophilic moiety capable of oxidizing a variety of small organic molecules and biomolecules, flavin is an important prosthetic group in many enzymes. Upon oxidation of the substrate, flavin is converted into its reduced (dihydrogenated) form. The catalytic cycle is completed through oxidation back to the oxidized form, thus restoring the enzyme's oxidizing capability. While it has been firmly established that oxidation of the reduced form of flavin is cast by molecular oxygen, yielding oxidized flavin and hydrogen peroxide, the mechanism of this process is still poorly understood. Herein, we investigate the radical mechanism, which is one of the possible reaction mechanisms, by quantum chemical calculations. Because molecular oxygen exists as a triplet in its electronic ground state, whereas the products are singlets, the reaction is accompanied by hopping between electronic surfaces. We find that the rate-limiting factor of flavin oxidation is likely associated with the change in the spin state of the system. By considering several possible reactions involving flavin and its derivatives in the radical form and by examining the corresponding parts of the potential energy surface in various spin states, we estimate the effective barrier of the kinetically and thermodynamically preferred variant of flavin oxidation to be about 15 kcal/mol in the gas phase and about 7 kcal/mol in a polar (aqueous) environment. This is in agreement with kinetic studies of the corresponding monoamine oxidase enzymes, confirming the radical mechanism as a viable option for flavin regeneration in enzymes.

Light-Oxygen-Voltage (LOV)-sensing Domains: Activation Mechanism and Optogenetic Stimulation

The light-oxygen-voltage (LOV) domains of phototropins emerged as essential constituents of light-sensitive proteins, helping initiate blue light-triggered responses. Moreover, these domains have been identified across all kingdoms of life. LOV domains utilize flavin nucleotides as co-factors and undergo structural rearrangements upon exposure to blue light, which activates an effector domain that executes the final output of the photoreaction. LOV domains are versatile photoreceptors that play critical roles in cellular signaling and environmental adaptation; additionally, they can noninvasively sense and control intracellular processes with high spatiotemporal precision, making them ideal candidates for use in optogenetics, where a light signal is linked to a cellular process through a photoreceptor. The ongoing development of LOV-based optogenetic tools, driven by advances in structural biology, spectroscopy, computational methods, and synthetic biology, has the potential to revolutionize the study of biological systems and enable the development of novel therapeutic strategies.

Two distinct mechanisms of flavoprotein spectral tuning revealed by low-temperature and time-dependent spectroscopy

Flavins such as flavin mononucleotide or flavin adenine dinucleotide are bound by diverse proteins, yet have very similar spectra when in the oxidized state. Recently, we developed new variants of flavin-binding protein CagFbFP exhibiting notable blue (Q148V) or red (I52V A85Q) shifts of fluorescence emission maxima. Here, we use time-resolved and low-temperature spectroscopy to show that whereas the chromophore environment is static in Q148V, an additional protein-flavin hydrogen bond is formed upon photoexcitation in the I52V A85Q variant. Consequently, in Q148V, excitation, emission, and phosphorescence spectra are shifted, whereas in I52V A85Q, excitation and low-temperature phosphorescence spectra are relatively unchanged, while emission spectrum is altered. We also determine the x-ray structures of the two variants to reveal the flavin environment and complement the spectroscopy data. Our findings illustrate two distinct color-tuning mechanisms of flavin-binding proteins and could be helpful for the engineering of new variants with improved optical properties.

Endogenous Photosensitizers in Human Skin

Endogenous photosensitizers play a critical role in both beneficial and harmful light-induced transformations in biological systems. Understanding their mode of action is essential for advancing fields such as photomedicine, photoredox catalysis, environmental science, and the development of sun care products. This review offers a comprehensive analysis of endogenous photosensitizers in human skin, investigating the connections between their electronic excitation and the subsequent activation or damage of organic biomolecules. We gather the physicochemical and photochemical properties of key endogenous photosensitizers and examine the relationships between their chemical reactivity, location within the skin, and the primary biochemical events following solar radiation exposure, along with their influence on skin physiology and pathology. An important take-home message of this review is that photosensitization allows visible light and UV-A radiation to have large effects on skin. The analysis presented here unveils potential causes for the continuous increase in global skin cancer cases and emphasizes the limitations of current sun protection approaches.

Tetramethylalloxazines as efficient singlet oxygen photosensitizers and potential redox-sensitive agents

Tetramethylalloxazines (TMeAll) have been found to have a high quantum yield of singlet oxygen generation when used as photosensitizers. Their electronic structure and transition energies (S0 → Si, S0 → Ti, T1 → Ti) were calculated using DFT and TD-DFT methods and compared to experimental absorption spectra. Generally, TMeAll display an energy diagram similar to other derivatives belonging to the alloxazine class of compounds, namely π,π* transitions are accompanied by closely located n,π* transitions. Photophysical data such as quantum yields of fluorescence, fluorescence lifetimes, and nonradiative rate constants were also studied in methanol (MeOH), acetonitrile (ACN), and 1,2-dichloroethane (DCE). The transient absorption spectra were also analyzed. To assess cytotoxicity of new compounds, a hemolytic assay was performed using human red blood cells (RBC) in vitro. Subsequently, fluorescence lifetime imaging experiments (FLIM) were performed on RBC under physiological and oxidative stress conditions alone or in the presence of TMeAll allowing for pinpointing changes caused by those compounds on the intracellular environment of these cells.

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., they use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. We show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near-reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

QM/MM Modeling of the Flavin Functionalization in the RutA Monooxygenase

Oxygenase activity of the flavin-dependent enzyme RutA is commonly associated with the formation of flavin-oxygen adducts in the enzyme active site. We report the results of quantum mechanics/molecular mechanics (QM/MM) modeling of possible reaction pathways initiated by various triplet state complexes of the molecular oxygen with the reduced flavin mononucleotide (FMN) formed in the protein cavities. According to the calculation results, these triplet-state flavin-oxygen complexes can be located at both re-side and si-side of the isoalloxazine ring of flavin. In both cases, the dioxygen moiety is activated by electron transfer from FMN, stimulating the attack of the arising reactive oxygen species at the C4a, N5, C6, and C8 positions in the isoalloxazine ring after the switch to the singlet state potential energy surface. The reaction pathways lead to the C(4a)-peroxide, N(5)-oxide, or C(6)-hydroperoxide covalent adducts or directly to the oxidized flavin, depending on the initial position of the oxygen molecule in the protein cavities.

Fine spectral tuning of a flavin-binding fluorescent protein for multicolor imaging

Flavin-binding fluorescent proteins are promising genetically encoded tags for microscopy. However, spectral properties of their chromophores (riboflavin, flavin mononucleotide, and flavin adenine dinucleotide) are notoriously similar even between different protein families, which limits applications of flavoproteins in multicolor imaging. Here, we present a palette of 22 finely tuned fluorescent tags based on the thermostable LOV domain from Chloroflexus aggregans. We performed site saturation mutagenesis of three amino acid positions in the flavin-binding pocket, including the photoactive cysteine, to obtain variants with fluorescence emission maxima uniformly covering the wavelength range from 486 to 512 nm. We demonstrate three-color imaging based on spectral separation and two-color fluorescence lifetime imaging of bacteria, as well as two-color imaging of mammalian cells (HEK293T), using the proteins from the palette. These results highlight the possibility of fine spectral tuning of flavoproteins and pave the way for further applications of flavin-binding fluorescent proteins in fluorescence microscopy.

Molecular properties and tautomeric equilibria of isolated flavins

Flavins are employed as redox cofactors and chromophores in a plethora of flavoenzymes. Their versatility is an outcome of intrinsic molecular properties of the isoalloxazine ring modulated by the protein scaffold and surrounding solvent. Thus, an investigation of isolated flavins with high-level electronic-structure methods and with error assessment of the calculated properties will contribute to building better models of flavin reactivity. Here, we benchmarked ground-state properties such as electron affinity, gas-phase basicity, dipole moment, torsion energy, and tautomer stability for lumiflavins in all biologically relevant oxidation and charge states. Overall, multiconfigurational effects are small and chemical accuracy is achieved by coupled-cluster treatments of energetic properties. Augmented basis sets and extrapolations to the complete basis-set limit are necessary for consistent agreement with experimental energetics. Among DFT functionals tested, M06-2X shows the best performance for most properties, except gas-phase basicity, in which M06 and CAM-B3LYP perform better. Moreover, dipole moments of radical flavins show large deviations for all functionals studied. Tautomers with noncanonical protonation states are significantly populated at normal temperatures, adding to the complexity of modeling flavins. These results will guide future computational studies of flavoproteins and flavin chemistry by indicating the limitations of electronic-structure methodologies and the contributions of multiple tautomeric states.

Comparing ultrafast excited state quenching of flavin 1,N6-ethenoadenine dinucleotide and flavin adenine dinucleotide by optical spectroscopy and DFT calculations

Flavins are photoenzymatic cofactors often exploiting the absorption of light to energize photoinduced redox chemistry in a variety of contexts. Both flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) are used for this function. The study of these photoenzymes has been facilitated using flavin analogs. Most of these analogs involve modification of the flavin ring, and there is recent evidence that adenine (Ade)-modified FAD can affect enzyme turnover, but so far this has only been shown for enzymes where the adenine and flavin rings are close to each other in a stacked conformation. FAD is also stacked in aqueous solution, and its photodynamics are quite different from unstacked FAD or FMN. Oxidized photoexcited FAD decays rapidly, presumably through PET with Ade as donor and Fl* as acceptor. Definitive identification of the spectral signatures of Ade^∙+ and Fl^∙- radicals is elusive. Here we use the FAD analog Flavin 1,N⁶-Ethenoadenine Dinucleotide (εFAD) to study how different photochemical outcomes depend on the identity of the Ade moiety in stacked FAD and its analog εFAD. We have used UV-Vis transient absorption spectroscopy complemented by TD-DFT calculations to investigate the excited state evolution of the flavins. In FAD*, no radicals were observed, suggesting that FAD* does not undergo PET. εFAD* kinetics showed a broad absorption band that suggests a charge transfer state exists upon photoexcitation with evidence for radical pair formation. Surprisingly, significant triplet flavin was produced from εFAD* We hypothesize that the dipolar (ε)Ade moieties differentially modulate the singlet-triplet energy gap, resulting in different intersystem crossing rates. The additional electron density on the etheno group of εFAD supplies better orbital overlap with the flavin S₁ state, accelerating charge transfer in that molecule.

Tuning the Quantum Chemical Properties of Flavins via Modification at C8

Flavins are central to countless enzymes but display different reactivities depending on their environments. This is understood to reflect modulation of the flavin electronic structure. To understand changes in orbital natures, energies, and correlation over the ring system, we begin by comparing seven flavin variants differing at C8, exploiting their different electronic spectra to validate quantum chemical calculations. Ground state calculations replicate a Hammett trend and reveal the significance of the flavin π-system. Comparison of higher-level theories establishes CC2 and ACD(2) as methods of choice for characterization of electronic transitions. Charge transfer character and electron correlation prove responsive to the identity of the substituent at C8. Indeed, bond length alternation analysis demonstrates extensive conjugation and delocalization from the C8 position throughout the ring system. Moreover, we succeed in replicating a particularly challenging UV/Vis spectrum by implementing hybrid QM/MM in explicit solvents. Our calculations reveal that the presence of nonbonding lone pairs correlates with the change in the UV/Vis spectrum observed when the 8-methyl is replaced by NH₂, OH, or SH. Thus, our computations offer routes to understanding the spectra of flavins with different modifications. This is a first step toward understanding how the same is accomplished by different binding environments.