Effect of microhydration on the electronic structure of the chromophores of the photoactive yellow and green fluorescent proteins

Unraveling the mystery of solvation-dependent fluorescence of fluorescein dianion using computational study

Fluorescein, one of the brightest fluorescent dye molecules, is a widely used fluorophore for various applications from biomedicine to industry. The dianionic form of fluorescein is responsible for its high fluorescence quantum yield. Interestingly, the molecule was found to be nonfluorescent in the gas phase. This characteristic is attributed to the photodetachment process, which out-competes the fluorescence emission in the gas phase. In this work, we show that the calculated vertical and adiabatic detachment energies of fluorescein dianion in the gas and solvent phases account for the drastic differences observed in their fluorescence characteristics. The functional dependence of these detachment energies on the dianion's microsolvation was systematically investigated. The performance of different solvent models was also assessed. The higher thermodynamic stability of fluorescein dianion over the monoanion doublet in the solvent phase plays a crucial role in quenching photodetachment and activating the radiative channel with a high fluorescence quantum yield.

TDDFT Investigation of the Raman and Resonant Raman Spectra of Fluorescent Protein Chromophore Models

The off-resonance and resonant Raman spectra have been simulated for models of fluorescent protein chromophores, those of the green fluorescent protein (GFP, called FP1) and of DsRed (called FP2), which presents a longer π-conjugated path, with the aim of providing a systematic investigation of structural but also computational aspects. These were performed at the (time-dependent) density functional theory [(TD)DFT] level. The off-resonance intensities have been calculated from the derivatives of the frequency-dependent polarizability with respect to the normal coordinates while the resonant ones have been evaluated using Huang-Rhys factors determined from the gradients of the excitation energies with respect to the normal coordinates. When applied with the M05 meta-GGA exchange-correlation functional, this simple computational scheme can reproduce most of the experimental Raman signatures of FP1 in its protonated and deprotonated forms, the differences of vibrational signatures of the cis (Z) and trans (E) isomers, as well as their changes as a function of the excitation wavelength. On the other hand, testing the predictions made for FP2 would require new experimental work. It was also observed that simulations with methods that inadequately predict the resonant Raman spectra could nevertheless produce a UV-vis absorption spectrum that is quite similar to the one obtained with better methods, once realistic peak broadening has been applied.

Complexation of Green and Red Kaede Fluorescent Protein Chromophores by a Zwitterion to Probe Electrostatic and Induction Field Effects

The photophysics of green fluorescent protein (GFP) and red Kaede fluorescent protein (rKFP) are defined by the intrinsic properties of the light-absorbing chromophore and its interaction with the protein binding pocket. This work deploys photodissociation action spectroscopy to probe the absorption profiles for a series of synthetic GFP and rKFP chromophores as the bare anions and as complexes with the betaine zwitterion, which is assumed as a model for dipole microsolvation. Electronic structure calculations and energy decomposition analysis using Symmetry-Adapted Perturbation Theory are used to characterize gas-phase structures and complex cohesion forces. The calculations reveal a preponderance for coordination of betaine to the phenoxide deprotonation site predominantly through electrostatic forces. Calculations using the STEOM-DLPNO-CCSD method are able to reproduce absolute and relative vertical excitation energies for the bare anions and anion-betaine complexes. On the other hand, treatment of the betaine molecule with a point-charge model, in which the charges are computed from some common electron density population analysis schemes, show that just electrostatic and point-charge induction interactions are unable to account for the betaine-induced spectral shift. The present methodology could be applied to investigate cluster forces and optical properties in other gas-phase ion-zwitterion complexes.

Photoelectron photofragment coincidence spectroscopy of aromatic carboxylates: benzoate and p-coumarate

Photoelectron-photofragment coincidence spectroscopy was used to study the dissociation dynamics of the conjugate bases of benzoic acid and p-coumaric acid. Upon photodetachment at 266 nm (4.66 eV) both aromatic carboxylates undergo decarboxylation, as well as the formation of stable carboxyl radicals. The key energetics are computed using high-level electronic structure methods. The dissociation dynamics of benzoate were dominated by a two-body DPD channel resulting in CO₂ + C₆H₅ + e^-, with a very small amount of stable C₆H₅CO₂ showing that the radical ground state is stable and the excited states are dissociative. For p-coumarate (p-CA^-) the dominant channel is photodetachment resulting in a stable radical and a photoelectron with electron kinetic energy (eKE) <2 eV. We also observed a minor two-body dissociative photodetachment (DPD) channel resulting in CO₂ + HOC₆H₄CHCH + e^-, characterized by eKE <0.8 eV. Evidence was also found for a three-body ionic photodissociation channel producing HOC₆H₅ + HCC^- + CO₂. The ion beam contained both the phenolate and carboxylate isomers of p-CA^-, but DPD only occurred from the carboxylate form. For both species DPD is seen from the first and second excited states of the radical, where vibrational excitation is required for decarboxylation from the first excited radical state.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Probing the Microsolvation Environment of the Green Fluorescent Protein Chromophore In Vacuo

We present vibrational and electronic photodissociation spectra of a model chromophore of the green fluorescent protein in complexes with up to two water molecules, prepared in a cryogenic ion trap at 160-180 K. We find the band origin of the singly hydrated chromophore at 20 985 cm^-1 (476.5 nm) and observe partially resolved vibrational signatures. While a single water molecule induces only a small shift of the S₁ electronic band of the chromophore, without significant change of the Franck-Condon envelope, the spectrum of the dihydrate shows significant broadening and a greater blue shift of the band edge. Comparison of the vibrational spectra with predicted infrared spectra from density functional theory indicates that water molecules can interact with the oxygen atom on the phenolate group or on the imidazole moiety, respectively.

ortho and para chromophores of green fluorescent protein: controlling electron emission and internal conversion

Green fluorescent protein (GFP) continues to play an important role in the biological and biochemical sciences as an efficient fluorescent probe and is also known to undergo light-induced redox transformations. Here, we employ photoelectron spectroscopy and quantum chemistry calculations to investigate how the phenoxide moiety controls the competition between electron emission and internal conversion in the isolated GFP chromophore anion, following photoexcitation with ultraviolet light in the range 400-230 nm. We find that moving the phenoxide group from the para position to the ortho position enhances internal conversion back to the ground electronic state but that adding an additional OH group to the para chromophore, at the ortho position, impedes internal conversion. Guided by quantum chemistry calculations, we interpret these observations in terms of torsions around the C-C-C bridge being enhanced by electrostatic repulsions or impeded by the formation of a hydrogen-bonded seven-membered ring. We also find that moving the phenoxide group from the para position to the ortho position reduces the energy required for detachment processes, whereas adding an additional OH group to the para chromophore at the ortho position increases the energy required for detachment processes. These results have potential applications in tuning light-induced redox processes of this biologically and technologically important fluorescent protein.

Probing microhydration effect on the electronic structure of the GFP chromophore anion: Photoelectron spectroscopy and theoretical investigations

The photophysics of the Green Fluorescent Protein (GFP) chromophore is critically dependent on its local structure and on its environment. Despite extensive experimental and computational studies, there remain many open questions regarding the key fundamental variables that govern this process. One outstanding problem is the role of autoionization as a possible relaxation pathway of the excited state under different environmental conditions. This issue is considered in our work through combined experimental and theoretical studies of microsolvated clusters of the deprotonated p-hydroxybenzylidene-2,3-dimethylimidazolinone anion (HBDI(-)), an analog of the GFP chromophore. Through selective generation of microsolvated structures of predetermined size and subsequent analysis of experimental photoelectron spectra by high level ab initio methods, we are able to precisely identify the structure of the system, establish the accuracy of theoretical data, and provide reliable description of auto-ionization process as a function of hydrogen-bonding environment. Our study clearly illustrates the first few water molecules progressively stabilize the excited state of the chromophore anion against the autodetached neutral state, which should be an important trait for crystallographic water molecules in GFPs that has not been fully explored to date.

Controlling electron emission from the photoactive yellow protein chromophore by substitution at the coumaric acid group

Understanding how the interactions between a chromophore and its surrounding protein control the function of a photoactive protein remains a challenge. Here, we present the results of photoelectron spectroscopy measurements and quantum chemistry calculations aimed at investigating how substitution at the coumaryl tail of the photoactive yellow protein chromophore controls competing relaxation pathways following photoexcitation of isolated chromophores in the gas phase with ultraviolet light in the range 350-315 nm. The photoelectron spectra are dominated by electrons resulting from direct detachment and fast detachment from the 2(1)ππ* state but also have a low electron kinetic energy component arising from autodetachment from lower lying electronically excited states or thermionic emission from the electronic ground state. We find that substituting the hydrogen atom of the carboxylic acid group with a methyl group lowers the threshold for electron detachment but has very little effect on the competition between the different relaxation pathways, whereas substituting with a thioester group raises the threshold for electron detachment and appears to 'turn off' the competing electron emission processes from lower lying electronically excited states. This has potential implications in terms of tuning the light-induced electron donor properties of photoactive yellow protein.

Effect of a Single Water Molecule on the Electronic Absorption by o- and p-Nitrophenolate: A Shift to the Red or to the Blue?

Many photoactive biomolecules are anions and exhibit ππ* optical transitions but with a degree of charge transfer (CT) character determined by the local environment. The phenolate moiety is a common structural motif among biochromophores and luminophores, and nitrophenolates are good model systems because the nitro substituent allows for CT-like transitions. Here we report gas-phase absorption spectra of o- and p-nitrophenolate·H2O complexes to decipher the effect of just one H2O and compare them with ab initio calculations of vertical excitation energies. The experimental band maximum is at 3.01 and 3.00 eV for ortho and para isomers, respectively, and is red-shifted by 0.10 and 0.13 eV relative to the bare ions, respectively. These shifts indicate that the transition has become more CT-like because of localization of negative charge on the phenolate oxygen, i.e., diminished delocalization of the negative excess charge. However, the transition bears less CT than that of m-nitrophenolate·H2O because this complex absorbs further to the red (2.56 eV). Our work emphasizes the importance of local perturbations: one water causes a larger shift than experienced in bulk for para isomer and almost the full shift for ortho isomer. Predicting microenvironmental effects in the boundary between CT and non-CT with high accuracy is nontrivial. However, in agreement with experiment, our calculations show a competition between the effects of electronic delocalization and electrostatic interaction with the solvent molecule. As a result, the excitation energy of ortho and para isomers is less sensitive to hydration than that of the meta isomer because donor and acceptor orbitals are only weakly coupled in the meta isomer.

Spectroscopy of nitrophenolates in vacuo: effect of spacer, configuration, and microsolvation on the charge-transfer excitation energy

In a charge-transfer (CT) transition, electron density moves from one end of the molecule (donor) to the other end (acceptor). This type of transition is of paramount importance in nature, for example, in photosynthesis, and it governs the excitation of several protein biochromophores and luminophores such as the oxyluciferin anion that accounts for light emission from fireflies. Both transition energy and oscillator strength are linked to the coupling between the donor and acceptor groups: The weaker the coupling, the smaller the excitation energy. But a weak coupling necessarily also causes a low oscillator strength possibly preventing direct excitation (basically zero probability in the noncoupling case). The coupling is determined by the actual spacer between the two groups, and whether the spacer acts as an insulator or a conductor. However, it can be difficult or even impossible to distinguish the effect of the spacer from that of local solvent molecules that often cause large solvent shifts due to different ground-state and excited-state stabilization. This calls for gas-phase spectroscopy experiments where absorption by the isolated molecule is identified to unequivocally establish the intrinsic molecular properties with no perturbations from a microenvironment. From such insight, the effect of a protein microenvironment on the CT excited state can be deduced. In this Account, we review our results over the last 5 years from mass spectroscopy experiments using specially designed apparatus on several charged donor-acceptor ions that are based on the nitrophenolate moiety and π-extended derivatives, which are textbook examples of donor-acceptor chromophores. The phenolate oxygen is the donor, and the nitro group is the acceptor. The choice of this system is also based on the fact that phenolate is a common structural motif of biochromophores and luminophores, for example, it is a constituent of the oxyluciferin anion. A presentation of the setups used for gas-phase ion spectroscopy in Aarhus is given, and we address issues of whether double bonds or triple bonds best convey electronic coupling between the phenolate oxygen and the nitro group, the significance of separating the donor and acceptor spatially, the influence of cross-conjugation versus linear conjugation, and along this line ortho versus meta versus para configuration, and not least the effect of a single solvent molecule (water, methanol, or acetonitrile). From systematic studies, a clear picture has emerged that has been supported by high-level calculations of electronically excited states. Our work shows that CC2 coupled-cluster calculations of vertical excitation energies are within 0.2 eV of experimental band maxima, and importantly, that the theoretical method is excellent in predicting the relative order of excitation energies of a series of nitrophenolates. Finally, we discuss future challenges such as following the change in absorption as a function of the number of solvent molecules and when gradually approaching the bulk limit.

Microhydration effects on the electronic properties of protonated phenol: a theoretical study

The CC2 (second-order approximate coupled cluster method) has been employed to investigate microhydration effect on electronic properties of protonated phenol (PhH(+)) According to the CC2 calculation results on electronic excited states of microhydrated PhH(+), for the S1 and S2 electronic states, which are of (1)ππ* nature and belong to the A' representation of molecular Cs point group, a significant blue shift effect on the S1 and S2 electronic states, which are of 1ππ* nature and belong to the A' representation of molecular Cs point group, in comparison to corresponding transitions on bare cation (PhH(+)), has been predicted. Nevertheless, for the S3-S0 (1A'', 1σπ*) transition, a large red shift effect has been predicted. Furthermore, it has been found that the lowest (1)σπ* state plays a prominent role in the photochemistry of these systems. In the bare protonated phenol, the (1)σπ* state is a bound state with a broad potential curve along the OH stretching coordinate, while it is dissociative in microhydrated species. This indicates to a predissociation of the S1((1)ππ*) state by a low-lying (1)σπ* state, which leads the excited system to a concerted proton-transfer reaction from protonated chromophore to the solvent. The dissociative (1)σπ* state in monohydrated PhH(+) has small barrier, while increasing the solvent molecules up to three removes the barrier and consequently expedites the proton-transfer reaction dynamics.

On the optical absorption of the anionic GFP chromophore in vacuum, solution, and protein

In spite of the large number of experimental and theoretical studies, the optical absorption trend of the green fluorescent protein (GFP) chromophore in several environments has not been fully understood. We calculated at the same level of time dependent density functional theory the vertical excitation energy of the anionic GFP chromophore in the protein and in ethanol, dioxane, methanol and water solutions. As result, we reproduced for the first time the experimental trend of the absorption peaks with 0.015 eV as the standard deviation of the accuracy. This systematic error allowed us to analyze with confidence the relative weight of several solvation effects on the vertical excitation energy. Experimental trends not correlated with the solvent polarity were therefore explained with a fine balance of different steric and electronic effects on the photophysics of the chromophore. As an indirect and remarkable result, the present analysis confirms that the optical absorption of the chromophore in the gas phase is close to the value of 2.84 eV extrapolated by Dong et al. (J. Am. Chem. Soc., 2006, 128, 12038), and, as a consequence, that the protein environment induces a red shift of 0.23 eV.

Exploring structural and optical properties of fluorescent proteins by squeezing: modeling high-pressure effects on the mStrawberry and mCherry red fluorescent proteins

Molecular dynamics calculations of pressure effects on mStrawberry and mCherry fluorescent proteins are reported. The simulations reveal that mStrawberry has much floppier structure at atmospheric pressure, as evidenced by larger backbone fluctuations and the coexistence of two conformers that differ by Ser146 orientation. Consequently, pressure increase has a larger effect on mStrawberry, making its structure more rigid and reducing the population of one of the conformers. The most significant effect of pressure increase is in the hydrogen-bonding network between the chromophore and the nearby residues. The quantum-mechanics/molecular mechanics calculations of excitation energies in mStrawberry explain the observed blue shift and identify Lys70 as the residue that has the most pronounced effect on the spectra. The results suggest that pressure increase causes an initial increase of fluorescence yield only for relatively floppy fluorescent proteins, whereas the fluorescent proteins that have more rigid structures have quantum yields close to their maximum. The results suggest that a low quantum yield in fluorescent proteins is dynamic in nature and depends on the range of thermal motions of the chromophore and fluctuations in the H-bonding network rather than on their average structure.

Toward understanding the redox properties of model chromophores from the green fluorescent protein family: an interplay between conjugation, resonance stabilization, and solvent effects

The redox properties of model chromophores from the green fluorescent protein family are characterized computationally using density functional theory with a long-range corrected functional, the equation-of-motion coupled-cluster method, and implicit solvation models. The analysis of electron-donating abilities of the chromophores reveals an intricate interplay between the size of the chromophore, conjugation, resonance stabilization, presence of heteroatoms, and solvent effects. Our best estimates of the gas-phase vertical/adiabatic detachment energies of the deprotonated (i.e., anionic) model red, green, and blue chromophores are 3.27/3.15, 2.79/2.67, and 2.75/2.35 eV, respectively. Vertical/adiabatic ionization energies of the respective protonated (i.e., neutral) species are 7.64/7.35, 7.38/7.15, and 7.70/7.32 eV, respectively. The standard reduction potentials (E(red)(0)) of the anionic (Chr•/Chr–) and neutral (Chr+•/Chr) model chromophores in acetonitrile are 0.34/1.40 V (red), 0.22/1.24 V (green), and −0.12/1.02 V (blue), suggesting, counterintuitively, that the red chromophore is more difficult to oxidize than the green and blue ones (in both neutral and deprotonated forms). The respective redox potentials in water follow a similar trend but are more positive than the acetonitrile values.

Non-Equilibrium Isomer Distribution of the Gas-Phase Photoactive Yellow Protein Chromophore

The conjugate base of para-coumaric acid, which can be conveniently generated in the gas phase by electrospray ionization (ESI), is a commonly used model system for the chromophore of the photoactive yellow protein. Here we report its gas-phase IR spectrum, which shows that the anion easily adopts a carboxylate structure lying 60 kJ/mol higher in energy than the global minimum phenoxide structure. Generation of the biologically more relevant phenoxide isomer by ESI can be achieved using dry acetonitrile as solvent.