Tautomer-Selective Fluorescence Spectroscopy of Oxyluciferin Anions

Bioluminescence in fireflies and related insects arises as emission from the fluorophore oxyluciferin, yet the color of the emission in these insects can range from red to green. The chromophore's microenvironment or multiple tautomeric forms may be responsible for the color tuning; however, these effects are difficult to separate in condensed phases. To investigate the role of oxyluciferin tautomerization in the color tuning mechanism, gas-phase spectroscopy eliminates solvent effects and allows us to study the fluorescence from individual tautomers. Using a home-built mass-spectrometry setup with a cylindrical ion trap cooled with liquid nitrogen, we measure fluorescence from the enol-locked form of oxyluciferin in the gas phase and characterize the photophysics of both keto and enol forms. At 100 K, the enol-locked form has an emission maximum of 564 ± 1 nm, coinciding with a previously reported assignment in oxyluciferin. We measure the absorption spectrum and find a maximum at 560.5 ± 0.5 nm, which implies a Stokes shift of 110 cm-1. The absorption spectrum is compared to Franck-Condon simulated spectra that identify one dominant vibrational mode in the transition. Additionally, we ultimately separated the emission by the enol and keto forms present in the trap by selectively exciting each form. We demonstrate that fluorescence measured close to the 0-0 transition limits the reheating of the ions, thereby providing the coldest ions and therefore the narrowest emission spectra. These experimental data are also crucial benchmarks for computational studies, offering actual emission spectra in the gas phase for both tautomeric forms. Thus, our findings serve as essential reference points for excited-state calculations aimed at understanding the color tuning mechanism of bioluminescence computationally.

Hydration effects on the vibrational properties of carboxylates: From continuum models to QM/MM simulations

The presence of carboxyl groups in a molecule delivers an affinity to metal cations and a sensitivity to the chemical environment, especially for an environment that can give rise to intermolecular hydrogen bonds. Carboxylate groups can also induce intramolecular interactions, such as the formation of hydrogen bonds with donor groups, leading to an impact on the conformational space of biomolecules. In the latter case, the protonation state of the amino groups plays an important role. In order to provide an accurate description of the modifications induced in a carboxylated molecule by the formation of hydrogen bonds, one needs a compromise between a quantum chemical description of the system and the necessity to take into account explicit solvent molecules. In this work, we propose a bottom-up approach to study the conformational space and the carboxylate stretching band of (bio)organic anions. Starting from the anions in a continuum solvent, we then move to calculations using a microsolvation approach including one explicit water molecule per polar group, immersed in a continuum. Finally, we run QM/MM molecular dynamics simulations to analyze the solvation properties and to explore the anions conformational space. The results thus obtained are in good agreement with the description given by the microsolvation approach and they bring a more detailed description of the solvation shell and of the intermolecular hydrogen bonds.

Unraveling the contributions to the spectral shape of flexible dyes in solution: insights on the absorption spectrum of an oxyluciferin analogue

We present a computational investigation of the absorption spectrum in water of 5,5-spirocyclopropyl-oxyluciferin (5,5-CprOxyLH), an analogue of the emitter compound responsible for the bioluminescence in fireflies. Several factors participate in determining the 5,5-CprOxyLH's spectral shape: (i) the contribution of the four close-energy excited states, which show significant non-adiabatic couplings, (ii) the flexible molecular structure and (iii) the specific interactions established with the surrounding environment, which strongly couple the protic solvent dynamics with the dye's spectral response. To tackle the challenge to capture and dissect the role of all these effects we preliminarily investigate the role of non-adiabatic couplings with quantum dynamics simulations and a linear vibronic coupling model in the gas phase. Then, we account for both the molecular flexibility and solvent interactions by resorting to a mixed quantum classical protocol, named Adiabatic Molecular Dynamics generalized Vertical Gradient (Ad-MD|gVG), which is built on a method recently proposed by some of us. It is rooted in the partition between stiff degrees of freedom of the dye, accounted for at the vibronic level within the harmonic approximation, and flexible degrees of freedom of the solute (and of the solvent), described classically through a sampling based on Molecular Dynamics (MD). Ad-MD|gVG avoids spurious effects arising in the excited state Hessians due to non-adiabatic couplings, and can therefore be applied to account for the contributions of the first four excited states to the 5,5-CprOxyLH absorption spectrum. The final simulated spectrum is in very good agreement with the experiment, especially when the MD is driven by a refined quantum-mechanically derived force-field. More importantly, the origin of each separate contribution to the spectral shape is appropriately accounted for, paving the way to future applications of the method to more complex systems or alternative spectroscopies, as emission or circular dichroism.

Computational Investigation of Substituent Effects on the Fluorescence Wavelengths of Oxyluciferin Analogs

Oxyluciferin, which is the light emitter for firefly bioluminescence, has been subjected to extensive chemical modifications to tune its emission wavelength and quantum yield. However, the exact mechanisms for various electron-donating and withdrawing groups to perturb the photophysical properties of oxyluciferin analogs are still not fully understood. To elucidate the substituent effects on the fluorescence wavelength of oxyluciferin analogs, we applied the absolutely localized molecular orbitals (ALMO)-based frontier orbital analysis to assess various types of interactions (i.e. permanent electrostatics/exchange repulsion, polarization, occupied-occupied orbital mixing, virtual-virtual orbital mixing, and charge-transfer) between the oxyluciferin and substituent orbitals. We suggested two distinct mechanisms that can lead to red-shifted oxyluciferin emission wavelength, a design objective that can help increase the tissue penetration of bioluminescence emission. Within the first mechanism, an electron-donating group (such as an amino or dimethylamino group) can contribute its highest occupied molecular orbital (HOMO) to an out-of-phase combination with oxyluciferin's HOMO, thus raising the HOMO energy of the substituted analog and narrowing its HOMO-LUMO gap. Alternatively, an electron-withdrawing group (such as a nitro or cyano group) can participate in an in-phase virtual-virtual orbital mixing of fragment LUMOs, thus lowering the LUMO energy of the substituted analog. Such an ALMO-based frontier orbital analysis is expected to lead to intuitive principles for designing analogs of not only the oxyluciferin molecule, but also many other functional dyes.

Spectrochemistry of Firefly Bioluminescence

The chemical reactions underlying the emission of light in fireflies and other bioluminescent beetles are some of the most thoroughly studied processes by scientists worldwide. Despite these remarkable efforts, fierce academic arguments continue around even some of the most fundamental aspects of the reaction mechanism behind the beetle bioluminescence. In an attempt to reach a consensus, we made an exhaustive search of the available literature and compiled the key discoveries on the fluorescence and chemiluminescence spectrochemistry of the emitting molecule, the firefly oxyluciferin, and its chemical analogues reported over the past 50+ years. The factors that affect the light emission, including intermolecular interactions, solvent polarity, and electronic effects, were analyzed in the context of both the reaction mechanism and the different colors of light emitted by different luciferases. The collective data points toward a combined emission of multiple coexistent forms of oxyluciferin as the most probable explanation for the variation in color of the emitted light. We also highlight realistic research directions to eventually address some of the remaining questions related to firefly bioluminescence. It is our hope that this extensive compilation of data and detailed analysis will not only consolidate the existing body of knowledge on this important phenomenon but will also aid in reaching a wider consensus on some of the mechanistic details of firefly bioluminescence.

Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface

Theoretical determinations of absorption cross sections (σ) in the gas phase and molar extinction coefficients (ε) in condensed phases (water solution, interfaces or surfaces, protein or nucleic acids embeddings, etc.) are of interest when rates of photochemical processes, J = ∫ ϕ(λ) σ(λ) I(λ) dλ, are needed, where ϕ(λ) and I(λ) are the quantum yield of the process and the irradiance of the light source, respectively, as functions of the wavelength λ. Efficient computational strategies based on single-reference quantum-chemistry methods have been developed enabling determinations of line shapes or, in some cases, achieving rovibrational resolution. Developments are however lacking for strongly correlated problems, with many excited states, high-order excitations, and/or near degeneracies between states of the same and different spin multiplicities. In this work, we define and compare the performance of distinct computational strategies using multiconfigurational quantum chemistry, nuclear sampling of the chromophore (by means of molecular dynamics, ab initio molecular dynamics, or Wigner sampling), and conformational and statistical sampling of the environment (by means of molecular dynamics). A new mathematical approach revisiting previous absolute orientation algorithms is also developed to improve alignments of geometries. These approaches are benchmarked through the nπ* band of acrolein not only in the gas phase and water solution but also in a gas-phase/water interface, a common situation for instance in atmospheric chemistry. Subsequently, the best strategy is used to compute the absorption band for the adduct formed upon addition of an OH radical to the C6 position of uracil and compared with the available experimental data. Overall, quantum Wigner sampling of the chromophore with molecular dynamics sampling of the environment with CASPT2 electronic-structure determinations arise as a powerful methodology to predict meaningful σ(λ) and ε(λ) band line shapes with accurate absolute intensities.

Quantum-mechanical hydration plays critical role in the stability of firefly oxyluciferin isomers: State-of-the-art calculations of the excited states

Stabilizing mechanisms of three possible isomers (phenolate-keto, phenolate-enol, and phenol-enolate) of the oxyluciferin anion hydrated with quantum explicit water molecules in the first singlet excited state were investigated using first-principles Born-Oppenheimer molecular dynamics simulations for up to 1.8 ns (or 3.7 × 10⁶ MD steps), revealing that the surrounding water molecules were distributed to form clear single-layered structures for phenolate-keto and multi-layered structures for phenolate-enol and phenol-enolate isomers. The isomers employed different stabilizing mechanisms compared to the ground state. Only the phenolate-keto isomer became attracted to the water molecules in its excited state and was stabilized by increasing the number of hydrogen bonds with nearby water molecules. The most stable isomer in the excited state was the phenolate-keto, and the phenolate-enol and phenol-enolate isomers were higher in energy by ∼0.38 eV and 0.57 eV, respectively, than the phenolate-keto. This was in contrast to the case of ground state in which the phenolate-enol was the most stable isomer.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

pH-Dependent absorption spectrum of oxyluciferin analogues in the active site of firefly luciferase

In the quest for the identification of the light emitter(s) responsible for the firefly bioluminescence, the study of oxyluciferin analogues with controlled chemical and electronic structures is of particular importance. In this article, we report the results of our experimental and computational investigation of the pH-dependent absorption spectra characterizing three analogues bound into the luciferase cavity, together with adenosine-monophosphate (AMP). While the analogue microscopic pKa values do not differ much from their reference values, it turns out that the AMP protonation state is analogue-dependent and never doubly-deprotonated. A careful analysis of the interactions evidences the main role of E344 glutamic acid, as well as the flexibility of the cavity which can accommodate any oxyluciferin analogue. The consideration of the absorption spectra suggests that the oxyluciferin enolate form has to be excluded from the list of the bioluminescence reaction products.

Fungal Light Emitter: Understanding Its Chemical Nature and pH-Dependent Emission in Water Solution

Fungal bioluminescence is a fascinating natural process, standing out for the continuous conversion of chemical energy into light. The structure of fungal oxyluciferin (light emitter) was proposed in 2017, being different and more complex than other oxyluciferins. The complexity of fungal oxyluciferin arises from diverse equilibria such as keto/enol tautomerization or deprotonation equilibria of four titratable groups. For this reason, still some crucial details of its structure remain unexplored. To obtain further structural information, a combined experimental and computational study of natural and three synthetic fungal oxyluciferin analogues has been performed. Here, we state the most stable chemical form of fungal oxyluciferin regarding its keto and enol tautomers, in the ground and excited states. We propose the (3Z,5E)-6-(3,4-dihydroxyphenyl)-4-hydroxy-2-oxohexa-3,5-dienoic acid form as the light emitter (fluorescent state) in water solution. Moreover, we show that chemical modifications on fungal oxyluciferin can affect the relative stability of the conformers. Furthermore, we show the clear effect of pH on emission. General conclusions about the role of these titratable groups in emission modulation have been drawn, such as the key role of dihydroxyphenyl deprotonation. This study is key to further analyze the properties of fungal bioluminescence and propose novel synthetic analogues.

Light emission colour modulation study of oxyluciferin synthetic analogues via QM and QM/MM approaches

Study of emission spectra of three analogues of firefly oxyluciferin in gas phase, in water (PCM) and in protein.

Effect of Protein Conformation and AMP Protonation State on Fireflies' Bioluminescent Emission

The emitted color in fireflies’ bioluminescent systems depends on the beetle species the system is extracted from and on different external factors (pH, temperature…) among others. Controlling the energy of the emitted light (i.e., color) is of crucial interest for the use of such bioluminescent systems. For instance, in the biomedical field, red emitted light is desirable because of its larger tissue penetration and lower energies. In order to investigate the influence of the protein environment and the AMP protonation state on the emitted color, the emission spectra of the phenolate-keto and phenolate-enol oxyluciferin forms have been simulated by means of MD simulations and QM/MM calculations, considering: two different protein conformations (with an open or closed C-terminal domain with respect to the N-terminal) and two protonation states of AMP. The results show that the emission spectra when considering the protein characterized by a closed conformation are blue-shifted compared to the open conformation. Moreover, the complete deprotonation of AMP phosphate group (AMP²⁻) can also lead to a blue-shift of the emission spectra but only when considering the closed protein conformation (open form is not sensitive to changes of AMP protonation state). These findings can be reasoned by the different interactions (hydrogen-bonds) found between oxyluciferin and the surrounding (protein, AMP and water molecules). This study gets partial insight into the possible origin of the emitted color modulation by changes of the pH or luciferase conformations.