Role of the Molecular Environment in Flavoprotein Color and Redox Tuning: QM Cluster versus QM/MM Modeling

Spectroscopic and Computational Observation of Glutamine Tautomerization in the Blue Light Sensing Using Flavin Domain Photoreaction

Blue light sensing using flavin (BLUF) domains constitute a family of flavin-binding photoreceptors of bacteria and eukaryotic algae. BLUF photoactivation proceeds via a light-driven hydrogen-bond switch among flavin adenine dinucleotide (FAD) and glutamine and tyrosine side chains, whereby FAD undergoes electron and proton transfer with tyrosine and is subsequently re-oxidized by a hydrogen back-shuttle in picoseconds, constituting an important model system to understand proton-coupled electron transfer in biology. The specific structure of the hydrogen-bond patterns and the prevalence of glutamine tautomeric states in dark-adapted (DA) and light-activated (LA) states have remained controversial. Here, we present a combined femtosecond stimulated Raman spectroscopy (FSRS), computational chemistry, and site-selective isotope labeling Fourier-transform infrared spectroscopy (FTIR) study of the Slr1694 BLUF domain. FSRS showed distinct vibrational bands from the FADS₁ singlet excited state. We observed small but significant shifts in the excited-state vibrational frequency patterns of the DA and LA states, indicating that these frequencies constitute a sensitive probe for the hydrogen-bond arrangement around FAD. Excited-state model calculations utilizing four different realizations of hydrogen bond patterns and glutamine tautomeric states were consistent with a BLUF reaction model that involved glutamine tautomerization to imidic acid, accompanied by a rotation of its side chain. A combined FTIR and double-isotope labeling study, with ¹³C labeling of FAD and ¹⁵N labeling of glutamine, identified the glutamine imidic acid C═N stretch vibration in the LA state and the Gln C═O in the DA state. Hence, our study provides support for glutamine tautomerization and side-chain rotation in the BLUF photoreaction.

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.

Non-statistical fragmentation in photo-activated flavin mononucleotide anions

The spectroscopy and photo-induced dissociation of flavin mononucleotide anions in vacuo are investigated over the 300-500 nm wavelength range. Comparison of the dependence of fragment ion yields as a function of deposited photon energy with calculated dissociation energies and collision-induced dissociation measurements performed under single-collision conditions suggests that a substantial fraction of photo-activated ions decompose through non-statistical fragmentation pathways. Among these pathways is the dominant photo-induced fragmentation channel, the loss of a fragment identified as formylmethylflavin. The fragment ion specific action spectra reveal electronic transition energies close to those for flavins in solution and previously published gas-phase measurements, although the photo-fragment yield upon excitation of the S₂ ← S₀ transition appears to be suppressed.

Mechanism and dynamics of fatty acid photodecarboxylase

Fatty acid photodecarboxylase (FAP) is a photoenzyme with potential green chemistry applications. By combining static, time-resolved, and cryotrapping spectroscopy and crystallography as well as computation, we characterized Chlorella variabilis FAP reaction intermediates on time scales from subpicoseconds to milliseconds. High-resolution crystal structures from synchrotron and free electron laser x-ray sources highlighted an unusual bent shape of the oxidized flavin chromophore. We demonstrate that decarboxylation occurs directly upon reduction of the excited flavin by the fatty acid substrate. Along with flavin reoxidation by the alkyl radical intermediate, a major fraction of the cleaved carbon dioxide unexpectedly transformed in 100 nanoseconds, most likely into bicarbonate. This reaction is orders of magnitude faster than in solution. Two strictly conserved residues, R451 and C432, are essential for substrate stabilization and functional charge transfer.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

Electronic spectra of flavin in different redox and protonation states: a computational perspective on the effect of the electrostatic environment

This study discusses how UV/vis absorption spectra of flavin in different redox and protonation states are shifted by the nearby electrostatic microenvironment.

Novel flavin-based fluorescent proteins with red-shifted emission bands: a computational study

The iLOV protein is a promising member of the class of flavin mononucleotide (FMN) based fluorescent proteins (FbFPs). It is becoming a popular tool for bioanalytical applications and bioimaging as a competitor of the well-known green fluorescent protein and its analogues. The main limitation of FbFPs is that all the members have close values of their absorption and emission band maxima. Therefore the upcoming challenge is to introduce novel variants of FbFPs to extend their color palette. We report the results of computational studies of iLOV variants, introducing point mutations and chromophore analogues. We found that point mutations of the apoprotein and substitution of FMN with either 8-amino-FMN or 8-methylamino-FMN lead to the red shift of emission bands up to 100 nm. Substitution with 1-deaza-FMN and the point mutations of the apoprotein result in a set of novel fluorescent proteins with emission bands in the "transparent" window where light readily penetrates through mammalian tissues. Newly suggested FbFPs can be used for multicolor imaging and also as components of FRET pairs.

Ion mobility action spectroscopy of flavin dianions reveals deprotomer-dependent photochemistry

Photo-induced proton transfer, deprotomer-dependent photochemistry, and intramolecular charge transfer in flavin anions are investigated using action spectroscopy.

The intrinsic optical properties and photochemistry of flavin adenine dinucleotide (FAD) dianions are investigated using a combination of tandem ion mobility spectrometry and action spectroscopy. Two principal isomers are observed, the more stable form being deprotonated on the isoalloxazine group and a phosphate (N-3,PO₄ deprotomer), and the other on the two phosphates (PO₄,PO₄ deprotomer). Ion mobility data and electronic action spectra suggest that photo-induced proton transfer occurs from the isoalloxazine group to a phosphate group, converting the PO₄,PO₄ deprotomer to the N-3,PO₄ deprotomer. Comparisons of the isomer selective action spectra of FAD dianions and flavin monoanions with solution spectra and gas-phase photodissociation action spectra suggests that solvation shifts the electronic absorption of the deprotonated isoalloxazine group to higher energy. This is interpreted as evidence for significant charge transfer in the lowest optical transition of deprotonated isoalloxazine. Overall, this work demonstrates that the site of deprotonation of flavin anions strongly affects their electronic absorptions and photochemistry.

Absorption and luminescence spectroscopy of mass-selected flavin adenine dinucleotide mono-anions

We report the absorption profile of isolated Flavin Adenine Dinucleotide (FAD) mono-anions recorded using photo-induced dissociation action spectroscopy. In this charge state, one of the phosphoric acid groups is deprotonated and the chromophore itself is in its neutral oxidized state. These measurements cover the first four optical transitions of FAD with excitation energies from 2.3 to 6.0 eV (210-550 nm). The S₀ → S₂ transition is strongly blue shifted relative to aqueous solution, supporting the view that this transition has a significant charge-transfer character. The remaining bands are close to their solution-phase positions. This confirms that the large discrepancy between quantum chemical calculations of vertical transition energies and solution-phase band maxima cannot be explained by solvent effects. We also report the luminescence spectrum of FAD mono-anions in vacuo. The gas-phase Stokes shift for S₁ is 3000 cm^-1, which is considerably larger than any previously reported for other molecular ions and consistent with a significant displacement of the ground and excited state potential energy surfaces. Consideration of the vibronic structure is thus essential for simulating the absorption and luminescence spectra of flavins.

The Dynamic Origin of Color Tuning in Proteins Revealed by a Carotenoid Pigment

Understanding the microscopic origin of the color tuning in pigment-protein complexes is a challenging yet fundamental issue in photoactive biological systems. Here, we propose a possible interpretation by using a state-of-the-art multiscale strategy based on the integration of quantum chemistry and polarizable atomistic embeddings into a dynamic description. By means of such a strategy we are able to resolve the long-standing dispute over the coloration mechanism in the crustacyanin protein. It is shown that the combination of the dynamical flexibility of the carotenoid pigments (astaxanthin) with the responsive protein environment is essential to obtain quantitative predictions of the spectral tuning. The strong linear correlation between the excitation energies and the bond length alternation in the long-chain carotenoids modulated by the dynamical protein environment is a novel finding explaining the high color tunability in crustacyanin.

Role of active site conformational changes in photocycle activation of the AppA BLUF photoreceptor

Blue light using flavin adenine dinucleotide (BLUF) proteins are essential for the light regulation of a variety of physiologically important processes and serve as a prototype for photoinduced proton-coupled electron transfer (PCET). Free-energy simulations elucidate the active site conformations in the AppA (activation of photopigment and puc expression) BLUF domain before and following photoexcitation. The free-energy profile for interconversion between conformations with either Trp104 or Met106 closer to the flavin, denoted Trp_in/Met_out and Trp_out/Met_in, reveals that both conformations are sampled on the ground state, with the former thermodynamically favorable by ∼3 kcal/mol. These results are consistent with the experimental observation of both conformations. To analyze the proton relay from Tyr21 to the flavin via Gln63, the free-energy profiles for Gln63 rotation were calculated on the ground state, the locally excited state of the flavin, and the charge-transfer state associated with electron transfer from Tyr21 to the flavin. For the Trp_in/Met_out conformation, the hydrogen-bonding pattern conducive to the proton relay is not thermodynamically favorable on the ground state but becomes more favorable, corresponding to approximately half of the configurations sampled, on the locally excited state. The calculated energy gaps between the locally excited and charge-transfer states suggest that electron transfer from Tyr21 to the flavin is more facile for configurations conducive to proton transfer. When the active site conformation is not conducive to PCET from Tyr21, Trp104 can directly compete with Tyr21 for electron transfer to the flavin through a nonproductive pathway, impeding the signaling efficiency.

Point charge embedding for ONIOM excited states calculations

Hybrid quantum mechanical methods can assist in the interpretation and prediction of the electronic spectra of large molecular structures. In this work, we study the performance of the ONIOM (Our own N-layered Integrated molecular Orbital molecular Mechanics) hybrid method for the calculation of transition energies and oscillator strengths by embedding the core region in a field of fixed point charges. These charges introduce polarization effects from the substituent groups to the core region. We test various charge definitions, with particular attention to the issue of overpolarization near the boundary between layers. To minimize this issue, we fit the charges on the electrostatic potential of the entire structure in the presence of the link atoms used to cap dangling bonds. We propose two constrained fitting strategies: one that produces an average set of charges common to both model system calculations, EE(L1), and one that produces two separate sets of embedding charges, EE(L2). The results from our tests show that indeed electronic embedding with constrained-fitted charges tends to improve the performance of ONIOM compared to non-embedded calculations. However, the EE(L2) charges work best for transition energies, and the EE(L1) charges work best for oscillator strengths. This may be an indication that fixed point charges do not have enough flexibility to adapt to each system, and other effects (e.g., polarization of the embedding field) may be necessary.

Evidence for Tautomerisation of Glutamine in BLUF Blue Light Receptors by Vibrational Spectroscopy and Computational Chemistry

BLUF (blue light sensor using flavin) domains regulate the activity of various enzymatic effector domains in bacteria and euglenids. BLUF features a unique photoactivation through restructuring of the hydrogen-bonding network as opposed to a redox reaction or an isomerization of the chromophore. A conserved glutamine residue close to the flavin chromophore plays a central role in the light response, but the underlying modification is still unclear. We labelled this glutamine with (15)N in two representative BLUF domains and performed time-resolved infrared double difference spectroscopy. The assignment of the signals was conducted by extensive quantum chemical calculations on large models with 187 atoms reproducing the UV-vis and infrared signatures of BLUF photoactivation. In the dark state, the comparatively low frequency of 1,667 cm(-1) is assigned to the glutamine C=O accepting a hydrogen bond from tyrosine. In the light state, the signature of a tautomerised glutamine was extracted with the C=N stretch at ~1,691 cm(-1) exhibiting the characteristic strong downshift by (15)N labelling. Moreover, an indirect isotope effect on the flavin C4=O stretch was found. We conclude that photoactivation of the BLUF receptor does not only involve a rearrangement of hydrogen bonds but includes a change in covalent bonds of the protein.