Quantum Chemical Benchmark Studies of the Electronic Properties of the Green Fluorescent Protein Chromophore. 1. Electronically Excited and Ionized States of the Anionic Chromophore in the Gas Phase

Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins

Red fluorescent protein (RFP) variants are highly sought after for in vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pretwisted chromophore conformation which provides easier access to conical intersections. The main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

A photoelectron imaging study of the deprotonated GFP chromophore anion and RNA fluorescent tags

Green fluorescent protein (GFP), together with its family of variants, is the most widely used fluorescent protein for in vivo imaging. Numerous spectroscopic studies of the isolated GFP chromophore have been aimed at understanding the electronic properties of GFP. Here, we build on earlier work [A. V. Bochenkova, C. Mooney, M. A. Parkes, J. Woodhouse, L. Zhang, R. Lewin, J. M. Ward, H. Hailes, L. H. Andersen and H. H. Fielding, Chem. Sci., 2017, 8, 3154] investigating the impact of fluorine and methoxy substituents that have been employed to tune the electronic structure of the GFP chromophore for use as fluorescent RNA tags. We present photoelectron spectra following photoexcitation over a broad range of wavelengths (364-230 nm) together with photoelectron angular distributions following photoexcitation at 364 nm, which are interpreted with the aid of quantum chemistry calculations. The results support the earlier high-level quantum chemistry calculations that predicted how fluorine and methoxy substituents tune the electronic structure and we find evidence to suggest that the methoxy substituents enhance internal conversion, most likely from the 2ππ* state which has predominantly Feshbach resonance character, to the 1ππ* state.

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.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Interplay between Locally Excited and Charge Transfer States Governs the Photoswitching Mechanism in the Fluorescent Protein Dreiklang

We present the results of high-level electronic structure and dynamics simulations of the photoactive protein Dreiklang. With the goal of understanding the details of the Dreiklang photocycle, we carefully characterize the excited states of the ON- and OFF-forms of Dreiklang. The key finding of our study is the existence of a low-lying excited state of a charge-transfer character in the neutral ON form and that population of this state, which is nearly isoenergetic with the locally excited bright state, initiates a series of steps that ultimately lead to the formation of the hydrated dark chromophore (OFF state). These results allow us to refine the mechanistic picture of Dreiklang's photocycle and photoactivation.

Shining light on the electronic structure and relaxation dynamics of the isolated oxyluciferin anion

Firefly bioluminescence is exploited widely in imaging in the biochemical and biomedical sciences; however, our fundamental understanding of the electronic structure and relaxation processes of the oxyluciferin that emits the light is still rudimentary. Here, we employ photoelectron spectroscopy and quantum chemistry calculations to investigate the electronic structure and relaxation of a series of model oxyluciferin anions. We find that changing the deprotonation site has a dramatic influence on the relaxation pathway following photoexcitation of higher lying electronically excited states. The keto form of the oxyluciferin anion is found to undergo internal conversion to the fluorescent S₁ state, whereas we find evidence to suggest that the enol and enolate forms undergo internal conversion to a dipole bound state, possibly via the fluorescent S₁ state. Partially resolved vibrational structure points towards the involvement of out-of-plane torsional motions in internal conversion to the dipole bound state, emphasising the combined electronic and structural role that the microenvironment plays in controlling the electronic relaxation pathway in the enzyme.

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.

Dual-Self-Restricted GFP Chromophore Analogues with Significantly Enhanced Emission

The tremendous gap of fluorescence emission of synthetic green fluorescent protein (GFP) chromophore to the protein itself makes it impossible to use for applications in molecular and cellular imaging. Here, we developed an efficient methodology to enhance the photoluminescence response of synthetic GFP chromophore analogues by constructing dual-self-restricted chromophores. Single self-restricted chromophores were first generated with 2,5-dimethoxy substitution on the aromatic ring, which were further modified with phenyl or 2,5-dimethoxy phenyl to form dual-self-restricted chromophores. These two chromophores showed an obvious solvatofluorochromic color palette across blue to yellow with a maximum emission Stokes shift of 95 nm and dramatically enhanced fluorescence emission in various aprotic solvents, especially in hexane, where the QY reached around 0.6. Importantly, in acetonitrile and dimethyl sulfoxide, the fluorescence QYs of both chromophores were over 0.22, which were the highest reported so far in high polar organic solvents. Meanwhile, the fluorescence lifetimes also improved obviously with the maximum of around 4.5 ns. Theoretical calculations revealed a more favorable Mülliken atomic charge translocation over the double-bond bridge and illustrated much higher energy barriers for the Z/E photoisomerization together with larger bond orders compared with the GFP core chromophore, p-HBDI. Our work significantly improved the fluorescence emission of synthetic GFP chromophore analogues in polar solvents while reserved the multicolor emitting function, providing a solid molecular motif for engineering high-performance fluorescent probes.

Assessment of Functionals for TDDFT Calculations of One- and Two-Photon Absorption Properties of Neutral and Anionic Fluorescent Proteins Chromophores

Performance of DFT functionals with different percentages of exact Hartree-Fock exchange energy (EX) is assessed for recovery of the CC2 reference one- (OPA) and two-photon absorption (TPA) spectra of fluorescent proteins chromophores in vacuo. The investigated DFT functionals, together with their EX contributions are BLYP (0%), B3LYP (20%), B1LYP (25%), BHandHLYP (50%), and CAM-B3LYP (19% at short range and 65% at long range). Our test set consists of anionic and neutral chromophores as naturally occurring in the fluorescent proteins. For the first time, we compare TDDFT and CC2 methods for higher excited states than the S₁ state, exhibiting relatively large TPA intensity. Our TDDFT results for neutral chromophores reveal an increase in excitation energies as well as TPA and OPA intensities errors, compared to CC2-derived results, as the DFT functional contains less exact exchange. The long-range-corrected CAM-B3LYP functional performs the best, closely followed by BHandHLYP, while BLYP usually significantly underestimates all investigated spectral properties, hence being the worst in reproducing the reference CC2 results. The hybrid B3LYP and B1LYP functionals can be roughly placed in between. We propose that TDDFT may underestimate the TPA intensities for neutral chromophores of fluorescent proteins due to underestimated oscillator strengths between some excited states. In the case of anionic chromophores, we find that B3LYP and B1LYP functionals overcome others in terms of reproducing CC2 excitation energies. On the other hand, however, TPA intensity is usually significantly underestimated, and in this respect, CAM-B3LYP functional seems to be again superior. In contrast to the case of neutral chromophores, it seems that a large magnitude of excited-state dipole moments or changes in dipole moments upon excitation may be the driving force behind high TPA transition moments.

Rational designing of 8-hydroxyquinolin-imidazolinone-based fluorescent protein mutants with dramatically red shifted emission: A computational study

Engineering fluorescent proteins to be the customized in vivo labels for monitoring cellular dynamic events is critical in biochemical and biomedical studies. The design and development of novel red fluorescent proteins is one of the most important fronts in this field due to their potential of imaging the entire organism. A recent fluorescent protein mutant eqFP650-67-HqAla with the 8-hydroxyquinolin-imidazolinone (HQI) chromophore has the plausible bathochromic shift of ~30 nm in its emission spectrum wavelength comparing to the parent fluorescent protein eqFP650. However, molecular mechanism of this significant shift remains somewhat obscure. In this study, we carefully benchmarked our computational methods and performed extensive calculations to investigate various structural components' effect on the chromophore's emission energy and decipher the molecular origin of the spectral shift. The influences of conjugation size, substituent group, substituent site as well as the number of substituents have been examined by elaborately designed chromophore derivatives. Accordingly, we proposed several chromophore mutants with dramatic bathochromic shift of up to ~60 nm in their emission spectra. We further evaluated their structural stability in the protein using molecular dynamics simulations. Present theoretical study connects the structural feature of chromophore derivatives in red fluorescent proteins with their splendid performances in shifting the emission frequency and offer the molecular insight. The computational protocol and successive examination procedure to extract the structural effect utilized herein can also be widely applied to other fluorescent proteins in general. © 2018 Wiley Periodicals, Inc.

Revisiting the Dependence of Cu K-Edge X-ray Absorption Spectra on Oxidation State and Coordination Environment

X-ray absorption spectroscopy (XAS) at the Cu K-edge is an important tool for probing the properties of copper centers in transition-metal chemistry and catalysis. However, the interpretation of experimental XAS spectra requires a detailed understanding of the dependence of spectroscopic features on the local geometric and electronic structure, which can be established by theoretical X-ray spectroscopy. Here, we present a systematic computational study of the Cu K-edge XAS spectra of selected Cu complexes based on time-dependent density-functional theory in combination with a molecular orbital analysis of the relevant transitions. For a series of Cu ammine model complexes as well as a comprehensive test set of 12 Cu(I) and 5 Cu(II) complexes, we revisit the dependence of the pre-edge region in Cu K-edge XAS spectra on oxidation state and coordination geometry. While our calculations confirm earlier experimental assignments, we can also reveal additional signatures of the ligand orbitals and identify the underlying orbital interactions. The comprehensive picture revealed by this study will provide a reliable basis for the interpretation of in situ Cu K-edge XAS spectra of catalytic intermediates.

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.

Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution

The correct treatment of vibronic effects is vital for the modeling of absorption spectra of many solvated dyes. Vibronic spectra for small dyes in solution can be easily computed within the Franck-Condon approximation using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions on the electronic excited state. On the other hand, a straightforward way to account for solute-solvent interactions and temperature-dependent broadening is by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. Ensemble approaches usually do not account for vibronic transitions and thus often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining zero-temperature vibronic fine structure with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is treated classically through the sampling of configurations and quantum mechanical vibronic contributions are included as a zero-temperature correction to each vertical transition. In our calculation of the vertical excitations, significant regions of the solvent environment are treated fully quantum mechanically to account for solute-solvent polarization and charge-transfer. For the Franck-Condon calculations, a small amount of frozen explicit solvent is considered in order to capture solvent effects on the vibronic shape function. We test the proposed method by comparing calculated and experimental absorption spectra of Nile red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interactions, the combined approach yields significant improvements over the ensemble approach. For systems with weak to moderate solute-solvent interactions, both the high-energy vibronic tail and the width of the spectra are in excellent agreement with experiments.

Pyridinium Analogues of Green Fluorescent Protein Chromophore: Fluorogenic Dyes with Large Solvent-Dependent Stokes Shift

Novel fluorogenic dyes based on the GFP chromophore are developed. The compounds contain a pyridinium ring instead of phenolate and feature large Stokes shifts and solvent-dependent variations in the fluorescence quantum yield. Electronic structure calculations explain the trends in solvatochromic behavior in terms of the increase of the dipole moment upon excited-state relaxation in polar solvents associated with the changes in bonding pattern in the excited state. A unique combination of such optical characteristics and lipophilic properties enables using one of the new dyes for imaging the membrane structure of endoplasmic reticulum. An extremely high photostability (due to a dynamic exchange between the free and absorbed states) and selectivity make this compound a promising label for this type of cellular organelles.

Photoinduced Ultrafast Charge Transfer and Charge Migration in Small Gold Clusters Passivated by a Chromophoric Ligand

Because the development of attopulses, charge migration induced by short optical pulses has been extensively investigated. We report a computational purely electronic dynamical study of ultrafast few femtoseconds (fs) charge transfer and charge migration in realistic passivated stoichiometric Au₁₁ and Au₂₀ gold nanoclusters functionalized by a bipyridine ligand. We show that a net significant amount of electronic charge (0.1 to 0.4 |e| where |e| is the electron charge) is permanently transferred from the bipyridine chromophore to the gold cluster during the short 5-6 fs UV-vis strong pulse. This electron transfer to the metallic core is induced by the optical excitation of electronic states with a partial charge transfer character involving the chromophore before the onset of nuclei motion. In addition, the photoexcitation by the strong fs pulse builds a nonequilibrium electronic density that beats between the chromophore and the metallic core around the average of the transferred value. Modular systems made of a donor chromophore that can be photoexcited in the UV-vis range coupled to an efficient acceptor that could trap the charge are of interest for applications to nanodevices. Our study provides understanding on the very early, purely electronic dynamics built by the fs optical excitation and the initial charge separation step.

Combined quantum-mechanical molecular mechanics calculations with NWChem and AMBER: Excited state properties of green fluorescent protein chromophore analogue in aqueous solution

Combined quantum mechanical molecular mechanics (QM/MM) calculations have become a popular methodology for efficient and accurate description of large molecular systems. In this work we introduce our development of a QM/MM framework based on two well-known codes-NWChem and AMBER. As an initial application area we are focused on excited state properties of small molecules in an aqueous phase using an analogue of the green fluorescent protein (GFP) chromophore as a particular test case. Our approach incorporates high level coupled cluster theory for the analysis of excited states providing a reliable theoretical analysis of effects of an aqueous solvation environment on the photochemical properties of the GFP chromophore. Using a systematic approach, which involves comparison of gas phase and aqueous phase results for different protonation states and conformations, we resolve existing uncertainties regarding the theoretical interpretation of experimental data. We observe that the impact of aqueous environment on charged states generally results in blue shifts of the absorption spectra, but the magnitude of the effect is sensitive to both protonation state and conformation and can be rationalized based on charge movement into the area of higher/lower external electrostatic potentials. At neutral pH levels the experimentally observed absorption signal is most likely coming from the phenol protonated form. Our results also show that the high level electron correlated method is essential for a proper description of excited states of GFP. © 2017 Wiley Periodicals, Inc.

High-level Ab Initio Absorption Spectra Simulations of Neutral, Anionic and Neutral+ Chromophore of Green Fluorescence Protein Chromophore Models in Gas Phase and Solution

Semiclassical ab initio simulations of the absorption spectra of neutral and anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI), a model chromophore of green fluorescent protein (GFP) and of a positively charged neutral (N+)-HBDI chromophore model, were performed in gas phase with the resolution-of-identity algebraic diagrammatic construction through second-order (RI-ADC(2)) method. The calculated absorption spectra in gas phase are composed of one band centered at 3.51 eV (HBDI), 2.50 eV (HBDI^- ) and 3.02 eV ((N+)-HBDI) owing to the absorption of the first ¹ ππ* transition. Band maxima are redshifted by ~0.1 eV with respect to the corresponding vertical energies. The COSMO-RI-ADC(2) calculations of the first vertical excitation energy of HBDI, HBDI^- and (N+)-HBDI forms in polar solution including microsolvation simulate the observed solvent redshift for neutral HBDI and the solvent blueshift of the HBDI^- and (N+)-HBDI forms. The state-specific solvation approach applied to TDDFT calculations reproduced the experimental solvent shifts for the three HBDI forms, demonstrating a more accurate theoretical description as compared to the linear-response TDDFT approach.

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.

Photoelectron spectroscopy of isolated luciferin and infraluciferin anions in vacuo: competing photodetachment, photofragmentation and internal conversion

The electronic structure and dynamics of luciferin and infraluciferin have been investigated using photoelectron spectroscopy and quantum chemistry calculations.

Resolving Electronic Transitions in Synthetic Fluorescent Protein Chromophores by Magnetic Circular Dichroism

The detailed electronic structures of fluorescent chromophores are important for their use in imaging of living cells. A series of green fluorescent protein chromophore derivatives is examined by magnetic circular dichroism (MCD) spectroscopy, which allows the resolution of more bands than plain absorption and fluorescence. Observed spectral patterns are rationalized with the aid of time-dependent density functional theory (TDDFT) computations and the sum-over-state (SOS) formalism, which also reveals a significant dependence of MCD intensities on chromophore conformation. The combination of organic and theoretical chemistry with spectroscopic techniques also appears useful in the rational design of fluorescence labels and understanding of the chromophore's properties. For example, the absorption threshold can be heavily affected by substitution on the phenyl ring but not much on the five-member ring, and methoxy groups can be used to further tune the electronic levels.

Self-Restricted Green Fluorescent Protein Chromophore Analogues: Dramatic Emission Enhancement and Remarkable Solvatofluorochromism

The confinement effect of the β-barrel defines the emission profiles of the chromophores of the green fluorescent protein (GFP) family. Here, we describe the design strategy and mimicking of confinement effects via the chromophore itself, termed the self-restricted effect. By systematically tailoring the GFP core, a family of 2,5-dialkoxy-substituted GFP chromophore analogues is found to be highly emissive and show remarkable solvatofluorochromism in fluid solvents. Fluorescence quantum yield (QY) and lifetime measurements, in combination with theoretical calculations, illustrate the mechanism relying on inhibition of torsional rotation around the exocyclic CC bond. Meanwhile, theoretical calculations further reveal that the electrostatic interaction between the solvent and the imidazolinone oxygen can contribute to suppress the radiationless decay channel around the exocyclic C═C double bond. Our findings put forward a universal approach toward unlocked highly emissive GFPc analogues, potentially promoting the understanding of the photophysics and biochemical application of GFP chromophore analogues.

Extension of the Effective Fragment Potential Method to Macromolecules

The effective fragment potential (EFP) approach, which can be described as a nonempirical polarizable force field, affords an accurate first-principles treatment of noncovalent interactions in extended systems. EFP can also describe the effect of the environment on the electronic properties (e.g., electronic excitation energies and ionization and electron-attachment energies) of a subsystem via the QM/EFP (quantum mechanics/EFP) polarizable embedding scheme. The original formulation of the method assumes that the system can be separated, without breaking covalent bonds, into closed-shell fragments, such as solvent and solute molecules. Here, we present an extension of the EFP method to macromolecules (mEFP). Several schemes for breaking a large molecule into small fragments described by EFP are presented and benchmarked. We focus on the electronic properties of molecules embedded into a protein environment and consider ionization, electron-attachment, and excitation energies (single-point calculations only). The model systems include chromophores of green and red fluorescent proteins surrounded by several nearby amino acid residues and phenolate bound to the T4 lysozyme. All mEFP schemes show robust performance and accurately reproduce the reference full QM calculations. For further applications of mEFP, we recommend either the scheme in which the peptide is cut along the Cα-C bond, giving rise to one fragment per amino acid, or the scheme with two cuts per amino acid, along the Cα-C and Cα-N bonds. While using these fragmentation schemes, the errors in solvatochromic shifts in electronic energy differences (excitation, ionization, electron detachment, or electron-attachment) do not exceed 0.1 eV. The largest error of QM/mEFP against QM/EFP (no fragmentation of the EFP part) is 0.06 eV (in most cases, the errors are 0.01-0.02 eV). The errors in the QM/molecular mechanics calculations with standard point charges can be as large as 0.3 eV.

pH-Sensitive fluorophores from locked GFP chromophores by a non-alternant analogue of the photochemical meta effect

We report the synthesis and characterization of a pH-sensitive fluorescence switch based on a conformationally-locked green fluorescent protein (GFP) chromophore.

Emission shaping in fluorescent proteins: role of electrostatics and π-stacking

We obtained the fluorescence spectrum of the GFP with trajectory simulations, and revealed the role of the protein sidechains in emission shifts.

On the Nature of an Extended Stokes Shift in the mPlum Fluorescent Protein

Far-red fluorescent proteins (FPs) enable deep-tissue in vivo imaging. Combining FPs with large and small Stokes shifts enables single-excitation/dual-emission multicolor applications. Using a quantum mechanics/molecular mechanics (QM/MM) scheme, we carried out a series of simulations to identify the origin of an extended Stokes shift (0.2 eV) observed in mPlum, one of the most far-red-shifted FPs. We demonstrated that the red shift of emission is largely due to the excited-state relaxation of the chromophore itself. Rigid protein environment suppresses the relaxation; however, if the hydrogen-bond network around the chromophore is sufficiently flexible, it can rearrange upon electronic excitation, allowing the chromophore to relax. The reorganization of the hydrogen-bond network is driven by changes in bonding and charge distributions of the chromophore in the excited state. The ILE65 and GLU16 residues play the most important role. The MD simulations reveal two ground-state populations with the direct (Chro-ILE65···GLU16) and water-mediated (Chro-ILE65···Wat321···GLU16) hydrogen-bond patterns. In the excited state, both populations relax to a single emitting state with the water-mediated (Chro-ILE65···Wat321···GLU16) hydrogen-bond pattern, which provides a better match for the excited-state charge distribution (the acylimine's oxygen has a larger negative charge in S1 than in S0). The extended Stokes shift arises due to the conversion of the direct hydrogen-bond pattern to the water-mediated one accompanied by large structural relaxation of the electronically excited chromophore. This conclusion is supported by calculations for the GLU16LEU mutant, which has only one hydrogen-bond pattern. Consequently, no interconversion is possible, and the computed Stokes shift is small, in agreement with the experiment. Our theoretical findings provide support to a recent study of the Stokes shifts in mPlum and its mutants.

Long- and Short-Range Electrostatic Fields in GFP Mutants: Implications for Spectral Tuning

The majority of protein functions are governed by their internal local electrostatics. Quantitative information about these interactions can shed light on how proteins work and allow for improving/altering their performance. Green fluorescent protein (GFP) and its mutation variants provide unique optical windows for interrogation of internal electric fields, thanks to the intrinsic fluorophore group formed inside them. Here we use an all-optical method, based on the independent measurements of transition frequency and one- and two-photon absorption cross sections in a number of GFP mutants to evaluate these internal electric fields. Two physical models based on the quadratic Stark effect, either with or without taking into account structural (bond-length) changes of the chromophore in varying field, allow us to separately evaluate the long-range and the total effective (short- and long-range) fields. Both types of the field quantitatively agree with the results of independent molecular dynamic simulations, justifying our method of measurement.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

High-resolution X-ray absorption spectroscopy of iron carbonyl complexes

We apply high-energy-resolution fluorescence-detected (HERFD) X-ray absorption near-edge structure (XANES) spectroscopy to study iron carbonyl complexes.

Assessment of mode-mixing and Herzberg-Teller effects on two-photon absorption and resonance hyper-Raman spectra from a time-dependent approach

A time-dependent approach is presented to simulate the two-photon absorption (TPA) and resonance hyper-Raman scattering (RHRS) spectra including Duschinsky rotation (mode-mixing) and Herzberg-Teller (HT) vibronic coupling effects. The computational obstacles for the excited-state geometries, vibrational frequencies, and nuclear derivatives of transition dipole moments, which enter the expressions of TPA and RHRS cross sections, are further overcome by the recently developed analytical excited-state energy derivative approaches in the framework of time-dependent density functional theory. The excited-state potential curvatures are evaluated at different levels of approximation to inspect the effects of frequency differences, mode-mixing and HT on TPA and RHRS spectra. Two types of molecules, one with high symmetry (formaldehyde, p-difluorobenzene, and benzotrifluoride) and the other with non-centrosymmetry (cis-hydroxybenzylidene-2,3-dimethylimidazolinone in the deprotonated anion state (HDBI(-))), are used as test systems. The calculated results reveal that it is crucial to adopt the exact excited-state potential curvatures in the calculations of TPA and RHRS spectra even for the high-symmetric molecules, and that the vertical gradient approximation leads to a large deviation. Furthermore, it is found that the HT contribution is evident in the TPA and RHRS spectra of HDBI(-) although its one- and two-photon transitions are strongly allowed, and its effect results in an obvious blueshift of the TPA maximum with respect to the one-photon absorption maximum. With the HT and solvent effects getting involved, the simulated blueshift of 1291 cm(-1) agrees well with the experimental measurement.

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.

Competition between photodetachment and autodetachment of the 2(1)ππ* state of the green fluorescent protein chromophore anion

Using a combination of photoelectron spectroscopy measurements and quantum chemistry calculations, we have identified competing electron emission processes that contribute to the 350-315 nm photoelectron spectra of the deprotonated green fluorescent protein chromophore anion, p-hydroxybenzylidene-2,3-dimethylimidazolinone. As well as direct electron detachment from S0, we observe resonant excitation of the 2(1)ππ* state of the anion followed by autodetachment. The experimental photoelectron spectra are found to be significantly broader than photoelectron spectrum calculated using the Franck-Condon method and we attribute this to rapid (∼10 fs) vibrational decoherence, or intramolecular vibrational energy redistribution, within the neutral radical.

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.