Potential Energy Landscape of the Electronic States of the GFP Chromophore in Different Protonation Forms: Electronic Transition Energies and Conical Intersections

Analysis of computational models for an accurate study of electronic excitations in GFP

The accuracy of PERI-CC2 is successfully assessed against RVS-CC2 for GFP model clusters and subsequently applied to a whole protein model.

Hidden photoinduced reactivity of the blue fluorescent protein mKalama1

We report a complete photocycle of the blue fluorescent protein exhibiting two delayed branches coupled to hidden proton transfer events.

Quantum dynamics of the photostability of pyrazine

We propose a new mechanism for the radiationless decay of photoexcited pyrazine to its ground electronic state involving a conical intersection between the dark A_u(nπ) state and the ground state.

An experimental and theoretical study of the photoisomerization and thermal reversion on 5-arylmethylene-2-thioxoimidazolidin-4-one

Unraveling the photochemical behaviour of the GFP chromophore has attracted widespread attention among scientists. Results for a new chromophore analog are presented.

Resonantly Enhanced Multiphoton Ionization Spectrum of the Neutral Green Fluorescent Protein Chromophore

The photophysics of the green fluorescent protein is governed by the electronic structure of the chromophore at the heart of its β-barrel protein structure. We present the first two-color, resonance-enhanced, multiphoton ionization spectrum of the isolated neutral chromophore in vacuo with supporting electronic structure calculations. We find the absorption maximum to be 3.65 ± 0.05 eV (340 ± 5 nm), which is blue-shifted by 0.5 eV (55 nm) from the absorption maximum of the protein in its neutral form. Our results show that interactions between the chromophore and the protein have a significant influence on the electronic structure of the neutral chromophore during photoabsorption and provide a benchmark for the rational design of novel chromophores as fluorescent markers or photomanipulators.

Vibrationally Resolved Photoelectron Spectroscopy of the Model GFP Chromophore Anion Revealing the Photoexcited S1 State Being Both Vertically and Adiabatically Bound against the Photodetached D0 Continuum

The first excited state of the model green fluorescence protein (GFP) chromophore anion (S1) and its energy level against the electron-detached neutral radical D0 state are crucial in determining the photophysics and the photoinduced dynamics of GFP. Extensive experimental and theoretical studies, particularly several very recent gas-phase investigations, concluded that S1 is a bound state in the Franck-Condon vertical region with respect to D0. However, what remains unknown and challenging is if S1 is bound adiabatically, primarily due to lack of accurate experimental measurements as well as due to the close proximity in energy for these two states that even sophisticated high-level ab initio calculations cannot reliably predict. Here, we report a negative ion photoelectron spectroscopy study on the model GFP chromophore anion, the deprotonated p-hydroxybenzylidene-2,3-dimethylimidazolinone anion (HBDI(-)) taken under low-temperature conditions with improved energy resolution. Despite the considerable size and low symmetry of the molecule, resolved vibrational structures were obtained with the 0-0 transition being the most intense peak. The adiabatic (ADE) and vertical detachment (VDE) energies therefore are determined both to be 2.73 ± 0.01 eV, indicating that the detached D0 state is 0.16 eV higher in energy than the photon excited S1 state. The accurate ADE and VDE values and the well-resolved photoelectron spectra reported here provide much needed robust benchmarks for future theoretical investigations.

New insight into the potential energy landscape and relaxation pathways of photoexcited aniline from CASSCF and XMCQDPT2 electronic structure calculations

New insight into the nonadiabatic relaxation dynamics of aniline following excitation to its first three singlet excited states, 1¹ππ*, 1¹π3s/πσ* and 2¹ππ*.

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.

Thermal isomerization of the chromoprotein asFP595 and its kindling mutant A143G: QM/MM molecular dynamics simulations

Chromoprotein asFP595 and its A143G variant called kindling fluorescent protein (KFP) are among the chronologically first species for which trans-cis chromophore isomerization has been proposed as a driving force of photoswitching. In spite of long-lasting efforts to characterize the route between protein conformations referring to the trans and cis forms of the chromophore, the molecular mechanism of this transformation is still under debate. We report the results of computational studies of the trans-cis isomerization of the anionic and neutral chromophore inside the protein matrices in the ground electronic state for both variants, asFP595 and KFP. Corresponding free energy profiles (potentials of mean force) were evaluated by using molecular dynamics simulations with the quantum mechanical-molecular mechanical (QM/MM) forces. The computed free energy barrier for the cis-trans ground state (thermal) isomerization reaction is about 2 kcal/mol higher in KFP than that in asFP595. These results provide interpretation of experimental studies on thermal relaxation from the light-induced activation of fluorescence of these proteins and correctly show that the A143G mutation in asFP595 noticeably increases the lifetime of the fluorescence species.

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.

Direct and indirect electron emission from the green fluorescent protein chromophore

Photoelectron spectra of the deprotonated green fluorescent protein chromophore have been measured in the gas phase at several wavelengths within and beyond the S(0)-S(1) photoabsorption band of the molecule. The vertical detachment energy (VDE) was determined to be 2.68 ± 0.1 eV. The data show that the first electronically excited state is bound in the Franck-Condon region, and that electron emission proceeds through an indirect (resonant) electron-emission channel within the corresponding absorption band.

Ab initio study of the solvent H-bonding effect on ESIPT reaction and electronic transitions of 3-hydroxychromone derivatives

The electronic transitions occurring in 4-(N,N-dimethylamino)-3-hydroxyflavone (DMAF) and 2-furanyl-3-hydroxychromone (FHC) were investigated using the TDDFT method in aprotic and protic solvents. The solvent effect was incorporated into the calculations via the PCM formalism. The H-bonding between solute and protic solvent was taken into account by considering a molecular complex between these molecules. To examine the effect of the H-bond on the ESIPT reaction, the absorption and emission wavelengths as well as the energies of the different states that intervene during these electronic transitions were calculated in acetonitrile, ethanol and methanol. The calculated positions of the absorption and emission wavelengths in various solvents were in excellent agreement with the experimental spectra, validating our approach. We found that in DMAF, the hydrogen bonding with protic solvents makes the ESIPT reaction energetically unfavourable, which explains the absence of the ESIPT tautomer emission in protic solvents. In contrast, the excited tautomer state of FHC remains energetically favourable in both aprotic and protic solvents. Comparing our calculations with the previously reported time-resolved fluorescence data, the ESIPT reaction of DMAF in aprotic solvents is reversible because the emitting states are energetically close, whereas in FHC, ESIPT is irreversible because the tautomer state is below the corresponding normal state. Therefore, the ESIPT reaction in DMAF is controlled by the relative energies of the excited states (thermodynamic control), while in FHC the ESIPT is controlled probably by the energetic barrier (kinetic control).

Quantum chemistry behind bioimaging: insights from ab initio studies of fluorescent proteins and their chromophores

The unique properties of green fluorescent protein (GFP) have been harnessed in a variety of bioimaging techniques, revolutionizing many areas of the life sciences. Molecular-level understanding of the underlying photophysics provides an advantage in the design of new fluorescent proteins (FPs) with improved properties; however, because of its complexity, many aspects of the GFP photocycle remain unknown. In this Account, we discuss computational studies of FPs and their chromophores that provide qualitative insights into mechanistic details of their photocycle and the structural basis for their optical properties. In a reductionist framework, studies of well-defined model systems (such as isolated chromophores) help to understand their intrinsic properties, while calculations including protein matrix and/or solvent demonstrate, on the atomic level, how these properties are modulated by the environment. An interesting feature of several anionic FP chromophores in the gas phase is their low electron detachment energy. For example, the bright excited ππ* state of the model GFP chromophore (2.6 eV) lies above the electron detachment continuum (2.5 eV). Thus, the excited state is metastable with respect to electron detachment. This autoionizing character needs to be taken into account in interpreting gas-phase measurements and is very difficult to describe computationally. Solvation (and even microsolvation by a single water molecule) stabilizes the anionic states enough such that the resonance excited state becomes bound. However, even in stabilizing environments (such as protein or solution), the anionic chromophores have relatively low oxidation potentials and can act as light-induced electron donors. Protein appears to affect excitation energies very little (<0.1 eV), but alters ionization or electron detachment energies by several electron volts. Solvents (especially polar ones) have a pronounced effect on the chromophore's electronic states; for example, the absorption wavelength changes considerably, the ground-state barrier for cis-trans isomerization is reduced, and fluorescence quantum yield drops dramatically. Calculations reveal that these effects can be explained in terms of electrostatic interactions and polarization, as well as specific interactions such as hydrogen bonding. The availability of efficient computer implementations of predictive electronic structure methods is essential. Important challenges include developing faster codes (to enable better equilibrium sampling and excited-state dynamics modeling), creating algorithms for properties calculations (such as nonlinear optical properties), extending standard excited-state methods to autoionizing (resonance) states, and developing accurate QM/MM schemes. The results of sophisticated first-principle calculations can be interpreted in terms of simpler, qualitative molecular orbital models to explain general trends. In particular, an essential feature of the anionic GFP chromophore is an almost perfect resonance (mesomeric) interaction between two Lewis structures, giving rise to charge delocalization, bond-order scrambling, and, most importantly, allylic frontier molecular orbitals spanning the methine bridge. We demonstrate that a three-center Hückel-like model provides a useful framework for understanding properties of FPs. It can explain changes in absorption wavelength upon protonation or other structural modifications of the chromophore, the magnitude of transition dipole moment, barriers to isomerization, and even non-Condon effects in one- and two-photon absorption.

Insight into the common mechanism of the chromophore formation in the red fluorescent proteins: the elusive blue intermediate revealed

Understanding the chromophore maturation process in fluorescent proteins is important for the design of proteins with improved properties. Here, we present the results of electronic structure calculations identifying the nature of a blue intermediate, a key species in the process of the red chromophore formation in DsRed, TagRFP, fluorescent timers, and PAmCherry. The chromophore of the blue intermediate has a structure in which the π-system of the imidazole ring is extended by the acylimine bond, which can be represented by the model N-[(5-hydroxy-1H-imidazole-2yl)methylidene]acetamide (HIMA) compound. Ab initio and QM/MM calculations of the isolated model and protein-bound (mTagBFP) chromophores identify the anionic form of HIMA as the only structure that has absorption that is consistent with the experiment and is stable in the protein binding pocket. The anion and zwitterion are the only protonation forms of HIMA whose absorption (421 and 414 nm, or 2.95 and 3.00 eV) matches the experimental spectrum of the blue form in DsRed (the absorption maximum is 408 nm or 3.04 eV) and mTagBFP (400 nm or 3.10 eV). The QM/MM optimization of the protein-bound anionic form results in a structure that is close to the X-ray one, whereas the zwitterionic chromophore is unstable in the protein binding pocket and undergoes prompt proton transfer. The computed excitation energy of the protein-bound anionic form of the mTagBFP-like chromophore (3.04 eV) agrees with the experimental absorption spectrum of the protein. The DsRed-like chromophore formation in red fluorescent proteins is revisited on the basis of ab initio results and verified by directed mutagenesis revealing a key role of the amino acid residue 70, which is the second after the chromophore tripeptide, in the formation process.

Intramolecular hydrogen bonding plays a crucial role in the photophysics and photochemistry of the GFP chromophore

In commonly studied GFP chromophore analogues such as 4-(4-hydroxybenzylidene)-1,2-dimethyl-1H-imidazol-5(4H)-one (PHBDI), the dominant photoinduced processes are cis-trans isomerization and subsequent S(1) → S(0) decay via a conical intersection characterized by a highly twisted double bond. The recently synthesized 2-hydroxy-substituted isomer (OHBDI) shows an entirely different photochemical behavior experimentally, since it mainly undergoes ultrafast intramolecular excited-state proton transfer, followed by S(1) → S(0) decay and ground-state reverse hydrogen transfer. We have chosen 4-(2-hydroxybenzylidene)-1H-imidazol-5(4H)-one (OHBI) to model the gas-phase photodynamics of such 2-hydroxy-substituted chromophores. We first use various electronic structure methods (DFT, TDDFT, CC2, DFT/MRCI, OM2/MRCI) to explore the S(0) and S(1) potential energy surfaces of OHBI and to locate the relevant minima, transition state, and minimum-energy conical intersection. These static calculations suggest the following decay mechanism: upon photoexcitation to the S(1) state, an ultrafast adiabatic charge-transfer induced excited-state intramolecular proton transfer (ESIPT) occurs that leads to the S(1) minimum-energy structure. Nearby, there is a S(1)/S(0) minimum-energy conical intersection that allows for an efficient nonadiabatic S(1) → S(0) internal conversion, which is followed by a fast ground-state reverse hydrogen transfer (GSHT). This mechanism is verified by semiempirical OM2/MRCI surface-hopping dynamics simulations, in which the successive ESIPT-GSTH processes are observed, but without cis-trans isomerization (which is a minor path experimentally with less than 5% yield). These gas-phase simulations of OHBI give an estimated first-order decay time of 476 fs for the S(1) state, which is larger but of the same order as the experimental values measured for OHBDI in solution: 270 fs in CH(3)CN and 230 fs in CH(2)Cl(2). The differences between the photoinduced processes of the 2- and 4-hydroxy-substituted chromophores are attributed to the presence or absence of intramolecular hydrogen bonding between the two rings.

Excited-state intramolecular proton transfer molecules bearing o-hydroxy analogues of green fluorescent protein chromophore

o-Hydroxy analogues, 1a-g, of the green fluorescent protein chromophore have been synthesized. Their structures and electronic properties were investigated by X-ray single-crystal analyses, electrochemistry, and luminescence properties. In solid and nonpolar solvents 1a-g exist mainly as Z conformers that possess a seven-membered-ring hydrogen bond and undergo excited-state intramolecular proton transfer (ESIPT) reactions, resulting in a proton-transfer tautomer emission. Fluorescence upconversion dynamics have revealed a coherent type of ESIPT, followed by a fast vibrational/solvent relaxation (<1 ps) to a twisted (regarding exo-C(5)-C(4)-C(3) bonds) conformation, from which a fast population decay of a few to several tens of picoseconds was resolved in cyclohexane. Accordingly, the proton-transfer tautomer emission intensity is moderate (0.08 in 1e) to weak (∼10(-4) in 1a) in cyclohexane. The stronger intramolecular hydrogen bonding in 1g suppresses the rotation of the aryl-alkene bond, resulting in a high yield of tautomer emission (Φ(f) ≈ 0.2). In the solid state, due to the inhibition of exo-C(5)-C(4)-C(3) rotation, intense tautomer emission with a quantum yield of 0.1-0.9 was obtained for 1a-g. Depending on the electronic donor or acceptor strength of the substituent in either the HOMO or LUMO site, a broad tuning range of the emission from 560 (1g) to 670 nm (1a) has been achieved.