Origin, nature, and fate of the fluorescent state of the green fluorescent protein chromophore at the CASPT2//CASSCF resolution

Delayed vibrational modulation of the solvated GFP chromophore into a conical intersection

Green fluorescent protein (GFP) has revolutionized bioimaging and life sciences. Its successes have inspired modification of the chromophore structure and environment to tune emission properties, but outside the protein cage, the chromophore is essentially non-fluorescent. In this study, we employ the tunable femtosecond stimulated Raman spectroscopy (FSRS) and transient absorption (TA) to map the energy dissipation pathways of GFP model chromophore (HBDI) in basic aqueous solution. Strategic tuning of the Raman pump to 550 nm exploits the stimulated emission band to enhance excited state vibrational motions as HBDI navigates the non-equilibrium potential energy landscape to pass through a conical intersection. The time-resolved FSRS uncovers prominent anharmonic couplings between a global out-of-plane bending mode of ∼227 cm-1 and two modes at ∼866 and 1572 cm-1 before HBDI reaches the twisted intramolecular charge transfer (TICT) state on the ∼3 ps time scale. Remarkably, the wavelet transform analysis reveals a ∼500 fs delayed onset of the coupling peaks, in correlation with the emergence of an intermediate charge-separated state en route to the TICT state. This mechanism is corroborated by the altered coupling matrix for the HBDI Raman modes in the 50% (v/v) water-glycerol mixture, and a notable lengthening of the picosecond time constant. The real-time molecular "movie" of the general rotor-like HBDI isomerization reaction following photoexcitation represents a significant advance in comprehending the photochemical reaction pathways of the solvated GFP chromophore, therefore providing a crucial foundation to enable rational design of diverse nanomachines from efficient molecular rotors to bright fluorescent probes.

Volume-conserving photoisomerization of a nonplanar GFP chromophore derivative: Nonadiabatic dynamics simulation

Nonadiabatic dynamics of a nonplanar green fluorescent protein (GFP) chromophore derivative was examined theoretically by trajectory surface hopping approach at the CASSCF level based on Zhu-Nakamura theory. The geometry optimizations show that there are four ground-state minima and four conical intersections between the ground (S₀) and first excited (S₁) states. Four S₁-state minima were found at a perpendicularly twisted conformation around the imidazolinone-bridged bond. Upon excitation to S₁ state, the main decay pathways of four isomers involve different S₀/S₁ potential energy surface conical intersections. The dominant excited-state relaxation mechanism of this GFP chromophore derivative is the twists of two bridging bonds and the methyl group in the bridge combined with the pyramidalization character of the central carbon atom. Further twists of two bridging bonds and the methyl group occur sequentially in the S₀ state. It is worth to mention that the special volume-conserving motion of this molecule is attributed to twists of two bridging bonds in the same direction during the whole photoisomerization processes. The theoretical investigation presented herein provides important insights into the volume-conserving photoisomerization mechanisms of a nonplanar GFP chromophore derivative. We believe that the present work can benefit the design of the photochromic molecule devices in confined media.

Subatomic resolution X-ray structures of green fluorescent protein

Green fluorescent protein (GFP) is a light-emitting protein that does not require a prosthetic group for its fluorescent activity. As such, GFP has become indispensable as a molecular tool in molecular biology. Nonetheless, there has been no subatomic elucidation of the GFP structure owing to the structural polymorphism around the chromophore. Here, subatomic resolution X-ray structures of GFP without the structural polymorphism are reported. The positions of H atoms, hydrogen-bonding network patterns and accurate geometric parameters were determined for the two protonated forms. Compared with previously determined crystal structures and theoretically optimized structures, the anionic chromophores of the structures represent the authentic resonance state of GFP. In addition, charge-density analysis based on atoms-in-molecules theory and noncovalent interaction analysis highlight weak but substantial interactions between the chromophore and the protein environment. Considered with the derived chemical indicators, the lone pair-π interactions between the chromophore and Thr62 should play a sufficient role in maintaining the electronic state of the chromophore. These results not only reveal the fine structural features that are critical to understanding the properties of GFP, but also highlight the limitations of current quantum-chemical calculations.

Does the wavelength dependent photoisomerization process of the p‑coumaric acid come out from the electronic state dependent pathways?

Similar to the anion photoactive yellow protein (PYP) chromophore, the neutral form of the PYP chromophore was also found to exhibit a the wavelength-dependent photoisomerization quantum yield. The isomerization quantum yield increases with the increasing excitation energy on the S₁ state, while decreases when being excited to the S₂ state. Does this wavelength dependent product yield come out from the specific reaction pathways of the S₁ and S₂ states? This would mean that, the relaxation pathway of the S₂ state is distinct from that of the S₁ state and does not involve twisting motion. Does it break Kasha's rule by exhibiting a direct transition from the S₂ state to the ground state? The underlying mechanism needs further in. In this article, we employed the on-the-fly dynamics simulations and static electronic structure calculations to reveal the deactivation mechanism of the neutral form of the PYP chromophore. Our results indicated that the CC twisting motion dominates the S₁ state decay process. In contrast, for the decay process of the S₂ state, an ultrafast transition from the S₂ to the S₁ state through a planar conical intersection is observed, and the excess energy activates a new reaction channel to the ground state characterized by a puckering distortion of the ring. This pathway competes with the photoisomerization channel. No direct transition from S₂ to S₀ is observed, hence Kasha's rule is valid for this process. Our calcualtions can provide a reasonable explanation of the wavelength-dependent isomerization quantum yield of neutral PYP chromophore, and we hope it can provide theoretical foundations for comparing the effect of protonation state on the dynamcal behaviors of PYP chromophore.

Mechanism of the Chemiluminescent Reaction between Nitric Oxide and Ozone

The gas phase reaction of nitric oxide with ozone to give chemiluminescence is used extensively for detection of nitrogen oxides. The molecular mechanism of chemiluminescence in this reaction is not known. So far, the only chemiluminescent systems studied in depth are certain cycloperoxides, which emit light following decomposition. Given our understanding of the mechanism of chemiluminescence in those molecules, one would expect by extension that in the NO + O₃ reaction the chemiluminescent species (NO₂ in this case) is formed in the excited state through a reaction pathway that diverges from the ground state pathway near the transition state. A systematic search for such a pathway leads us to conclude that such a mechanism is unlikely. Instead, our study suggests that chemiluminescence in the NO + O₃ reaction is due to emission from the NO₂ vibronic states associated with the ground (X̃ ²A₁) and first excited (à ²B₂) electronic states, which are populated in the nascent NO₂ produced in the reaction. The vibronic coupling between the X̃ ²A₁ and à ²B₂ states of NO₂ is due to a conical intersection (CI), which is geometrically and energetically close to the à ²B₂ minimum energy geometry and only 1.3 eV higher than ground state NO₂. Further, the CI is 1.2 eV lower than the energy of the NO + O₃ reactants and therefore thermodynamically accessible following the reaction. An analysis of the product energy distribution indicates that the major fraction of the reaction energy is channeled into the vibrational modes of NO₂, sufficient to populate the vibronic states of NO₂ around the X̃/à CI. These vibronic states show dipole-allowed emission in a frequency range that is consistent with the observed broad chemiluminescence spectrum.

Role of Gln222 in Photoswitching of Aequorea Fluorescent Proteins: A Twisting and H-Bonding Affair?

Reversibly photoswitchable fluorescent proteins (RSFPs) admirably combine the genetic encoding of fluorescence with the ability to repeatedly toggle between a bright and dark state, adding a new temporal dimension to the fluorescence signal. Accordingly, in recent years RSFPs have paved the way to novel applications in cell imaging that rely on their reversible photoswitching, including many super-resolution techniques such as F-PALM, RESOLFT, and SOFI that provide nanoscale pictures of the living matter. Yet many RSFPs have been engineered by a rational approach only to a limited extent, in the absence of clear structure-property relationships that in most cases make anecdotic the emergence of the photoswitching. We reported [ Bizzarri et al. J. Am Chem Soc. 2010 , 102 , 85 ] how the E222Q replacement is a single photoswitching mutation, since it restores the intrinsic cis-trans photoisomerization properties of the chromophore in otherwise nonswitchable Aequorea proteins of different color and mutation pattern (Q-RSFPs). We here investigate the subtle role of Q222 on the excited-state photophysics of the two simplest Q-RSFPs by a combined experimental and theoretical approach, using their nonswitchable anacestor EGFP as benchmark. Our findings link indissolubly photoswitching and Q222 presence, by a simple yet elegant scenario: largely twisted chromophore structures around the double bond (including hula-twist configurations) are uniquely stabilized by Q222 via H-bonds. Likely, these H-bonds subtly modulate the electronic properties of the chromophore, enabling the conical intersection that connects the excited cis to ground trans chromophore. Thus, Q222 belongs to a restricted family of single mutations that change dramatically the functional phenotype of a protein. The capability to distinguish quantitatively T65S/E222Q EGFP ("WildQ", wQ) from the spectrally identical EGFP by quantitative Optical Lock-In Detection (qOLID) witnesses the relevance of this mutation for cell imaging.

A direct-compatible formulation of the coupled perturbed complete active space self-consistent field equations on graphical processing units

We recently developed an algorithm to compute response properties for the state-averaged complete active space self-consistent field method (SA-CASSCF) that capitalized on sparsity in the atomic orbital basis. Our original algorithm was limited to treating small to moderate sized active spaces, but the recent development of graphical processing unit (GPU) based direct-configuration interaction algorithms provides an opportunity to extend this to large active spaces. We present here a direct-compatible version of the coupled perturbed equations, enabling us to compute response properties for systems treated with arbitrary active spaces (subject to available memory and computation time). This work demonstrates that the computationally demanding portions of the SA-CASSCF method can be formulated in terms of seven fundamental operations, including Coulomb and exchange matrix builds and their derivatives, as well as, generalized one- and two-particle density matrix and σ vector constructions. As in our previous work, this algorithm exhibits low computational scaling and is accelerated by the use of GPUs, making possible optimizations and nonadiabatic dynamics on systems with O(1000) basis functions and O(100) atoms, respectively.

Systematic Excited State Studies of Reversibly Switchable Fluorescent Proteins

The reversibly switchable fluorescent proteins Dronpa, rsFastLime, rsKame, Padron, and bsDronpa feature the same chromophore but display a 40 nm variation in absorption maxima and an only 18 nm variation in emission maxima. In the present contribution, we employ QM/MM models to investigate the mechanism of such remarkably different spectral variations, which are caused by just a few amino acid replacements. We show that the models, which are based on CASPT2//CASSCF level of QM theory, reproduce the observed trends in absorption maxima, with only a 3.5 kcal/mol blue-shift, and in emission maxima, with an even smaller 1.5 kcal/mol blue-shift with respect to the observed quantities. In order to explain the variations across the series, we look at the chromophore's electronic structure change during absorption and emission. Such analysis indicates that a change in charge-transfer character, which is more pronounced during absorption, triggers a cascade of hydrogen-bond-network rearrangements, suggesting preparation to an isomerization event. We also show how the contribution of Arg 89 and Arg 64 residues to the chromophore conformational changes correlate with the spectral variations in absorption and emission. Furthermore, we describe how the conical intersection stability is related to the protein's photophysical properties. While for the Dronpa, rsFastLime, and rsKame triad, the stability correlates with the photoswitching speed, this does not happen for bsDronpa and Padron, suggesting a less obvious photoisomerization mechanism.

Ab Initio Study of Decay Dynamics of 1-Nitronaphthalene Initiated from the S2(ππ* + nNOπ*) State

Irradiation of nitro-PAHs in solution at ambient conditions leads to formation of its lowest excited triplet, dissociation intermediates nitrogen oxide (NO^•) and aryloxy radical (Ar-O^•). Experimental and theoretical studies demonstrated that Franck-Condon excited singlet state S_FC(ππ*) to a receiver, higher-energy triplet state T_n(nπ*) controlled the ultrafast population of the triplet state and, hence, the slight fluorescence yield of nitronaphthalenes. However, the detailed information about the curve-crossings of potential energy surfaces and the major channels for forming T₁ species and Ar-O^• radical were unclear. Here, by using the CASSCF//CASPT2 method, an efficient decay channel is revealed: S_2-FC-1NN → S_2-MIN-1NN or S₂T_3-MIN-1NN → T_3-MIN-1NN or T₃T_2-MIN-1NN→ T_2-MIN-1NN or T₂T_1-MIN-1NN → T_1-MIN-1NN. This explains the high yield of T_1-1NN species and minor yield of Ar-O^• and NO^• radicals. The calculation results suggest the bifurcation processes take place predominantly after the internal conversion to the T_1-1NN state via T₂T_1-MIN-1NN, one leads to T_1-MIN-1NN, while the other to T_1-MIN-ISO to produce Ar-O^• and NO^• radicals.

Fluorescent Proteins Detect Host Structural Rearrangements via Electrostatic Mechanism

The rational design of genetically encoded fluorescent biosensors, which can detect rearrangements of target proteins via interdomain allostery, is hindered by the absence of mechanistic understanding of the underlying photophysics. Here, we focus on the modulation of fluorescence by mechanical perturbation in a popular biological probe: enhanced Green Fluorescent Protein (eGFP). Using a combination of molecular dynamics (MD) simulations and quantum chemistry, and a set of physically motivated assumptions, we construct a map of fluorescence quantum yield as a function of a 2D electric field imposed by the protein environment on the fluorophore. This map is transferable between Tsien's Class 2 GFP's, and it allows one to estimate the shifts in fluorescence intensity due to mechanical perturbations directly from MD simulations. We use it in combination with steered MD simulations to put forward a hypothesis for the mechanism of a genetically encoded voltage probe (ArcLight) whose mechanism is currently under debate.

A new twist in the photophysics of the GFP chromophore: a volume-conserving molecular torsion couple

Dynamics of a nonplanar GFP chromophore are studied experimentally and theoretically. Coupled torsional motion is responsible for the ultrafast decay.

The simple structure of the chromophore of the green fluorescent protein (GFP), a phenol and an imidazolone ring linked by a methyne bridge, supports an exceptionally diverse range of excited state phenomena. Here we describe experimentally and theoretically the photochemistry of a novel sterically crowded nonplanar derivative of the GFP chromophore. It undergoes an excited state isomerization reaction accompanied by an exceptionally fast (sub 100 fs) excited state decay. The decay dynamics are essentially independent of solvent polarity and viscosity. Excited state structural dynamics are probed by high level quantum chemical calculations revealing that the fast decay is due to a conical intersection characterized by a twist of the rings and pyramidalization of the methyne bridge carbon. The intersection can be accessed without a barrier from the pre-twisted Franck–Condon structure, and the lack of viscosity dependence is due to the fact that the rings twist in the same direction, giving rise to a volume-conserving decay coordinate. Moreover, the rotation of the phenyl, methyl and imidazolone groups is coupled in the sterically crowded structure, with the methyl group translating the rotation of one ring to the next. As a consequence, the excited state dynamics can be viewed as a torsional couple, where the absorbed photon energy leads to conversion of the out-of-plane orientation from one ring to the other in a volume conserving fashion. A similar modification of the range of methyne dyes may provide a new family of devices for molecular machines, specifically torsional couples.

Improving the Design of the Triple-Decker Motif in Red Fluorescent Proteins

We characterize computationally a red fluorescent protein (RFP) with the chromophore (Chro) sandwiched between two aromatic tyrosine rings in a triple-decker motif. According to the original proposal [ J. Phys. Chem. Lett. 2013 , 4 , 1743 ], such a tyrosine-chromophore-tyrosine π-stacked construct can be accommodated in the green fluorescent protein (GFP). A recent study [ ACS Chem. Biol. 2016 , 11 , 508 ] attempted to realize the triple-decker motif and obtained an RFP variant called mRojoA-VYGV with two tyrosine residues surrounding the chromophore. The crystal structure showed that only a tyrosine-chromophore pair was involved in π-stacking, whereas the second tyrosine was oriented perpendicularly, edge-to-face with respect to the chromophore. We propose a more promising variant of this RFP with a perfect triple-decker unit achieved by introducing additional mutations in mRojoA-VYGV. The structures and optical properties of model proteins based on the structures of mCherry and mRojoA are characterized computationally by QM(DFT)/MM. The electronic transitions in the protein-bound chromophores are computed by high-level quantum chemical methods. According to our calculations, the triple-decker chromophore unit in the new RFP variant is stable within the protein and its optical bands are red-shifted with respect to the parent mCherry and mRojoA species.

DNA mimics of red fluorescent proteins (RFP) based on G-quadruplex-confined synthetic RFP chromophores

Red fluorescent proteins (RFPs) have emerged as valuable biological markers for biomolecule imaging in living systems. Developing artificial fluorogenic systems that mimic RFPs remains an unmet challenge. Here, we describe the design and synthesis of six new chromophores analogous to the chromophores in RFPs. We demonstrate, for the first time, that encapsulating RFP chromophore analogues in canonical DNA G-quadruplexes (G4) can activate bright fluorescence spanning red and far-red spectral regions (Em = 583-668 nm) that nearly match the entire RFP palette. Theoretical calculations and molecular dynamics simulations reveal that DNA G4 greatly restricts radiationless deactivation of chromophores induced by a twisted intramolecular charge transfer (TICT). These DNA mimics of RFP exhibit attractive photophysical properties comparable or superior to natural RFPs, including high quantum yield, large Stokes shifts, excellent anti-photobleaching properties, and two-photon fluorescence. Moreover, these RFP chromophore analogues are a novel and distinctive type of topology-selective G4 probe specific to parallel G4 conformation. The DNA mimics of RFP have been further exploited for imaging of target proteins. Using cancer-specific cell membrane biomarkers as targets, long-term real-time monitoring in single live cell and two-photon fluorescence imaging in tissue sections have been achieved without the need for genetic coding.

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.

Origin of the Intrinsic Fluorescence of the Green Fluorescent Protein

Green fluorescent protein, GFP, has revolutionized biology, due to its use in bioimaging. It is widely accepted that the protein environment makes its chromophore fluoresce, whereas the fluorescence is completely lost when the native chromophore is taken out of GFP. By the use of a new femtosecond pump-probe scheme, based on time-resolved action spectroscopy, we demonstrate that the isolated deprotonated GFP chromophore can be trapped in the first excited state when cooled to 100 K. The trapping is shown to last for 1.2 ns, which is long enough to establish conditions for fluorescence and consistent with calculated trapping barriers in the electronically excited state. Thus, GFP fluorescence is traced back to an intrinsic chromophore property, and by improving excited-state trapping, protein interactions enhance the molecular fluorescence.

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.

Mechanism and bottlenecks in strand photodissociation of split green fluorescent proteins (GFPs)

Split GFPs have been widely applied for monitoring protein-protein interactions by expressing GFPs as two or more constituent parts linked to separate proteins that only fluoresce on complementing with one another. Although this complementation is typically irreversible, it has been shown previously that light accelerates dissociation of a noncovalently attached β-strand from a circularly permuted split GFP, allowing the interaction to be reversible. Reversible complementation is desirable, but photodissociation has too low of an efficiency (quantum yield <1%) to be useful as an optogenetic tool. Understanding the physical origins of this low efficiency can provide strategies to improve it. We elucidated the mechanism of strand photodissociation by measuring the dependence of its rate on light intensity and point mutations. The results show that strand photodissociation is a two-step process involving light-activated cis-trans isomerization of the chromophore followed by light-independent strand dissociation. The dependence of the rate on temperature was then used to establish a potential energy surface (PES) diagram along the photodissociation reaction coordinate. The resulting energetics-function model reveals the rate-limiting process to be the transition from the electronic excited-state to the ground-state PES accompanying cis-trans isomerization. Comparisons between split GFPs and other photosensory proteins, like photoactive yellow protein and rhodopsin, provide potential strategies for improving the photodissociation quantum yield.

Deciphering the molecular mechanism responsible for GCaMP6m's Ca2+-dependent change in fluorescence

The goal of this work is to determine how GCaMP6m's fluorescence is altered in response to Ca2+-binding. Our detailed spectroscopic study reveals the simplest explanation for how GCaMP6m changes fluorescence in response to Ca2+ is with a four-state model, in which a Ca2+-dependent change of the chromophore protonation state, due to a shift in pKa, is the predominant factor. The pKa shift is quantitatively explained by a change in electrostatic potential around the chromophore due to the conformational changes that occur in the protein when calmodulin binds Ca2+ and interacts with the M13 peptide. The absolute pKa values for the Ca2+-free and Ca2+-saturated states of GCaMP6m are critical to its high signal-to-noise ratio. This mechanism has important implications for further improvements to GCaMP6m and potentially for other similarly designed biosensors.

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.

Fluorescent and “breathable” CO2 responsive vesicles inspired from green fluorescent protein

CO₂ responsive fluorescent vesicles from a GFP chromophore labeled block-copolymer could change their size and fluorescence to mimic jellyfish breathing.

Thermal Isomerization Mechanism in Dronpa and Its Mutants

The photoswitching speed of the reversibly switchable fluorescent proteins (RSFPs) from the family of green fluorescent proteins (GFPs) changes upon mutation which is of direct importance for various high-resolution techniques. Dronpa is one of the most used RSFPs. Its point mutants rsFastLime (Dronpa V157G) and rsKame (Dronpa V157L) exhibit a striking difference in their photoswitching speed. Here the QM/MM on-the-fly string method is used in order to explore the details of the thermal isomerization mechanism. The four principal ways in which isomerization may occur have been scrutinized for each of the three proteins. It has been shown that thermal isomerization occurs via a one-bond-flip mechanism in all three proteins, although, in rsKame, where the chromophore is constrained more, the activation free energy difference between hula-twist and one-bond-flip is significantly smaller. Functional mode analysis has been applied to examine the motions of the amino acids during the isomerization. It clearly identifies the importance of Val/Leu 157 as well as the amino acids in the α-helix during the isomerization.

Photoinduced Chemistry in Fluorescent Proteins: Curse or Blessing?

Photoinduced reactions play an important role in the photocycle of fluorescent proteins from the green fluorescent protein (GFP) family. Among such processes are photoisomerization, photooxidation/photoreduction, breaking and making of covalent bonds, and excited-state proton transfer (ESPT). Many of these transformations are initiated by electron transfer (ET). The quantum yields of these processes vary significantly, from nearly 1 for ESPT to 10^-4-10^-6 for ET. Importantly, even when quantum yields are relatively small, at the conditions of repeated illumination the overall effect is significant. Depending on the task at hand, fluorescent protein photochemistry is regarded either as an asset facilitating new applications or as a nuisance leading to the loss of optical output. The phenomena arising due to phototransformations include (i) large Stokes shifts, (ii) photoconversions, photoactivation, and photoswitching, (iii) phototoxicity, (iv) blinking, (v) permanent bleaching, and (vi) formation of long-lived intermediates. The focus of this review is on the most recent experimental and theoretical work on photoinduced transformations in fluorescent proteins. We also provide an overview of the photophysics of fluorescent proteins, highlighting the interplay between photochemistry and other channels (fluorescence, radiationless relaxation, and intersystem crossing). The similarities and differences with photochemical processes in other biological systems and in dyes are also discussed.

Multiscale Simulations on Spectral Tuning and the Photoisomerization Mechanism in Fluorescent RNA Spinach

Fluorescent RNA aptamer Spinach can bind and activate a green fluorescent protein (GFP)-like chromophore (an anionic DFHBDI chromophore) displaying green fluorescence. Spectroscopic properties, spectral tuning, and the photoisomerization mechanism in the Spinach-DFHBDI complex have been investigated by high-level QM and hybrid ONIOM(QM:AMBER) methods (QM method: (TD)DFT, SF-BHHLYP, SAC-CI, LT-DF-LCC2, CASSCF, or MS-CASPT2), as well as classical molecular dynamics (MD) simulations. First, our benchmark calculations have shown that TD-DFT and spin-flip (SF) TD-DFT (SF-BHHLYP) failed to give a satisfactory description of absorption and emission of the anionic DFHBDI chromophore. Comparatively, SAC-CI, LT-DF-LCC2, and MS-CASPT2 can give more reliable transition energies and are mainly used to further study the spectra of the anionic DFHBDI chromophore in Spinach. The RNA environmental effects on the spectral tuning and the photoisomerization mechanism have been elucidated. Our simulations show that interactions of the anionic cis-DFHBDI chromophore with two G-quadruplexes as well as a UAU base triple suppress photoisomerization of DFHBDI. In addition, strong hydrogen bonds between the anionic cis-DFHBDI chromophore and nearby nucleotides facilitate its binding to Spinach and further inhibit the cis-trans photoisomerization of DFHBDI. Solvent molecules, ions, and loss of key hydrogen bonds with nearby nucleotides could induce dissociation of the anionic trans-DFHBDI chromophore from the binding site. These results provide new insights into fluorescent RNA Spinach and may help rational design of other fluorescent RNAs.

Electric Field Keeps Chromophore Planar and Produces High Yield Fluorescence in Green Fluorescent Protein

The green fluorescent protein and its designed variants fluoresce efficiently. Because the isolated chromophore is not fluorescent in a practical sense, it is apparent that the protein environment plays a crucial role in its efficiency. Because of various obstacles in studying excited state dynamics of complex systems, however, the detailed mechanism of this efficiency enhancement is not yet clearly elucidated. Here, by adopting excited state nonadiabatic molecular dynamics simulations together with an interpolated quantum chemical potential model of the chromophore, we find that the strong electric field from the protein matrix contributes dominantly to the motional restriction of the chromophore. The delay in twisting motion subsequently obstructs the nonradiative decay that competes with fluorescence, leading naturally to an enhancement in light-emitting efficiency. Surprisingly, steric constraints make only a minor contribution to these aspects. Through residue specific analyses, we identify a group of key residues that control the excited state behavior. Testing a series of mutant GFPs with different brightnesses also supports the view regarding the importance of protein electrostatics. Our findings may provide a useful guide toward designing new fluorescent chemical systems in the future.

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.

Theoretical Studies on the Photochemistry of Pentose Aminooxazoline, a Hypothetical Intermediate Product in the Prebiotic Synthetic Scenario of RNA Nucleotides

2-Aminooxazole is generally considered a prebiotic precursor of ribonucleotides on the early earth. Its pentose compound, pentose aminooxazoline, has been suggested to be a key intermediate in the prebiotic synthetic scenario. In this article, detailed mechanism of the photochemistry of pentose aminooxazoline has been studied by performing density functional theory and multireference complete active space self-consistent field calculations. Parallel to the "ring-puckering" process, which leads to ultrafast nonradiative deactivation, several other photodissociation channels are explored in detail. In addition, the influences of the pentose structure and solvation effects with both implicit and explicit water models have been uncovered for both neutral and protonated forms. The current theoretical results provide very important information not only for the photostability of RNA nucleotides but also for an in-depth understanding of the synthesis of other prebiotic nucleotides.