We use ab initio results to reparameterize a multi-reference semiempirical method to reproduce the ground and excited state potential energy surfaces (PESs) for the chromophore of Green Fluorescent Protein (GFP). The validity of the new parameter set is tested, and the new method is combined with a quantum mechanical/molecular mechanical (QM/MM) treatment so that it can be applied in the solution phase. Solvent effects on the energetics of the relevant conical intersections are explored. We then combine this representation of the ground and excited state PESs with the full multiple spawning (FMS) nonadiabatic wavepacket dynamics method to simulate the photodynamics of the neutral GFP chromophore in both gas and solution phases. In these calculations, the PESs and their nonadiabatic couplings are evaluated simultaneously with the nuclear dynamics, ie. "on-the-fly". The effect of solvation is seen to be quite dramatic, resulting in an order of magnitude decrease in the excited state lifetime. We observe a correlated torsion about a double bond and its adjacent single bond in both gas and solution phases. This is discussed in the context of previous proposals about minimal volume isomerization mechanisms in protein environments.

The chromophores of the most widely known fluorescent proteins (FPs) are derivatives of a core p-hydroxybenzylidene-imidazolinon-5-one (HBI) motif, which usually occurs as a phenolate anion. Double bond photoisomerization of the exocyclic bridge of HBI is widely held to be an important internal conversion mechanism for FP chromophores. Herein we describe the ground and excited-state electronic structures and potential energy surfaces of two model chromophores: 4- p-hydroxybenzylidiene-1,2-dimethyl-imidazolin-5-one anion (HBDI), representing green FPs (GFPs), and 2-acetyl-4-hydroxybenylidene-1-methyl-imidazolin-5-one anion (AHBMI), representing kindling FPs (KFPs). These chromophores differ by a single substitution, but we observe qualitative differences in the potential energy surfaces which indicate inversion of bond selection in the photoisomerization reaction. Bond selection is also modulated by whether the reaction proceeds from a Z or an E conformation. These configurations correspond to fluorescent and nonfluorescent states of structurally characterized FPs, including some which can be reversibly switched by specific illumination regimes. We explain the difference in bond selectivity via substituent stabilization effects on a common set of charge-localized chemical structures. Different combinations of these structures give rise to both optically active (planar) and twisted intramolecular charge-transfer (TICT) states of the molecules. We offer a prediction of the gas-phase absorption of AHBMI, which has not yet been measured. We offer a hypothesis to explain the unusual fluorescence of AHBMI in DMF solution, as well as an experimental proposal to test our hypothesis.

We use CASSCF and MRPT2 calculations to characterize the bridge photoisomerization pathways of a model red fluorescent protein (RFP) chromophore model. RFPs are homologues of the green fluorescent protein (GFP). The RFP chromophore differs from the GFP chromophore via the addition of an N-acylimine substitution to a common hydroxybenzylidene-imidazolinone (HBI) motif. We examine the substituent effects on the manifold of twisted intramolecular charge-transfer (TICT) states which mediates radiationless decay via bridge isomerization in fluorescent protein chromophore anions. We find that the substitution destabilizes states associated with isomerization about the imidazolinone-bridge bond and stabilizes states associated with phenoxy-bridge bond isomerization. We discuss the results in the context of chromophore conformation and quantum yield trends in the RFP subfamily, as well as recent studies on synthetic models where the acylimine has been replaced with an olefin.

Reversible photoswitching fluorescent proteins can be photoswitched between fluorescent and nonfluorescent states by different irradiation regimes. Accumulating spectroscopic and crystallographic evidence suggest a correlated change in protonation state and methine bridge isomerism of the chromophore. The anion can decay by photoisomerization of either of the methine bonds, but only one channel can act as a switch. Using ab initio multiple spawning dynamics simulations, we show that protonation is sufficient to change the photoisomerization channel in the chromophore. We propose that this behavior can underlie a switch given certain other conditions. We also propose a basis for coupling between excited-state basicity changes and selection of the photoisomerization channel based on the polarity of twisted charge-transfer states for neutral and anionic forms of the chromophore.

No abstract