Computational strategy for tuning spectral properties of red fluorescent proteins

Theoretical Insight into the Fluorescence Spectral Tuning Mechanism: A Case Study of Flavin-Dependent Bacterial Luciferase

Bioluminescence of bacteria is widely applied in biological imaging, environmental toxicant detection, and many other situations. Understanding the spectral tuning mechanism not only helps explain the diversity of colors observed in nature but also provides principles for bioengineering new color variants for practical applications. In this study, time-dependent density functional theory (TD-DFT) and quantum mechanics and molecular mechanics (QM/MM) calculations have been employed to understand the fluorescence spectral tuning mechanism of bacterial luciferase with a focus on the electrostatic effect. The spectrum can be tuned by both a homogeneous dielectric environment and oriented external electric fields (OEEFs). Increasing the solvent polarity leads to a redshift of the fluorescence emission maximum, λF, accompanied by a substantial increase in density. In contrast, applying an OEEF along the long axis of the isoalloxazine ring (X-axis) leads to a significant red- or blue-shift in λF, depending on the direction of the OEEF, yet with much smaller changes in intensity. The effect of polar solvents is directionless, and the red-shifts can be attributed to the larger dipole moment of the S1 state compared with that of the S0 state. However, the effect of OEEFs directly correlates with the difference dipole moment between the S1 and S0 states, which is directional and is determined by the charge redistribution upon deexcitation. Moreover, the electrostatic effect of bacterial luciferase is in line with the presence of an internal electric field (IEF) pointing in the negative X direction. Finally, the key residues that contribute to this IEF and strategies for modulating the spectrum through site-directed point mutations are discussed.

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).

FPredX: Interpretable models for the prediction of spectral maxima, brightness, and oligomeric states of fluorescent proteins

Fluorescent protein (FP) design is among the challenging protein design problems due to the tradeoffs among multiple properties to be optimized. Despite the accumulated efforts in design and characterization, progress has been slow in gaining a full understanding of sequence-property relationships to tackle the multiobjective design problem in FPs. In this study, we approach this problem by developing FPredX, a collection of gradient-boosted decision tree models, which mapped FP sequences to four major design targets of FPs, including excitation maximum, emission maximum, brightness, and oligomeric state. By training using one-hot encoded multiple aligned sequences with hyperparameters optimization in each model, FPredX models showed excellent prediction performance for all target properties compared with existing methods. We further interpreted the FPredX models by comparing the importance of positions along the aligned FP sequence to the predictive performance and suggested positions, which showed differential importance deemed by FPredX models to the prediction of each target property.

Influence of the Environment on Shaping the Absorption of Monomeric Infrared Fluorescent Proteins

Infrared fluorescent proteins (iRFPs) are potential candidates for deep-tissue in vivo imaging. Here, we provide molecular-level insights into the role of the protein environment in the structural stability of the chromophore within the protein binding pocket through the flexible hydrogen-bonding network using molecular dynamics simulation. Furthermore, we present systematic excited-state analysis to characterize the nature of the first two excited states and the role of the environment in shaping the nature of the chromophore's excited states within the hybrid quantum mechanics/molecular mechanics framework. Our results reveal that the environment red-shifts the absorption of the chromophore by about 0.32 eV compared to the isolated counterpart, and besides the structural stability, the protein environment does not alter the nature of the excited state of the chromophore significantly. Our study contributes to the fundamental understanding of the excited-state processes of iRFPs in a complex environment and provides a design principle for developing iRFPs with desired spectral properties.

Efficient switching of mCherry fluorescence using chemical caging

Fluorophores with dynamic or controllable fluorescence emission have become essential tools for advanced imaging, such as superresolution imaging. These applications have driven the continuing development of photoactivatable or photoconvertible labels, including genetically encoded fluorescent proteins. These new probes work well but require the introduction of new labels that may interfere with the proper functioning of existing constructs and therefore require extensive functional characterization. In this work we show that the widely used red fluorescent protein mCherry can be brought to a purely chemically induced blue-fluorescent state by incubation with β-mercaptoethanol (βME). The molecules can be recovered to the red fluorescent state by washing out the βME or through irradiation with violet light, with up to 80% total recovery. We show that this can be used to perform single-molecule localization microscopy (SMLM) on cells expressing mCherry, which renders this approach applicable to a very wide range of existing constructs. We performed a detailed investigation of the mechanism underlying these dynamics, using X-ray crystallography, NMR spectroscopy, and ab initio quantum-mechanical calculations. We find that the βME-induced fluorescence quenching of mCherry occurs both via the direct addition of βME to the chromophore and through βME-mediated reduction of the chromophore. These results not only offer a strategy to expand SMLM imaging to a broad range of available biological models, but also present unique insights into the chemistry and functioning of a highly important class of fluorophores.

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.

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.

Estimating orientation factors in the FRET theory of fluorescent proteins: the TagRFP-KFP pair and beyond

The orientation factor κ(2), one of the key parameters defining Förster resonance energy transfer efficiency, is determined by the transition dipole moment orientations of the donor and acceptor species. Using the results of quantum chemical and quantum mechanical/molecular mechanical calculations for the chromophore-containing pockets in selected colored proteins of the green fluorescent protein family, we derived transition dipole moments corresponding to the S0,min → S1 excitation for green fluorescent protein, red fluorescent protein (TagRFP), and kindling fluorescent protein, and the S1,min → S0 emission for TagRFP. These data allowed us to estimate κ(2) values for the TagRFP-linker-kindling fluorescent protein tetrameric complex required for constructing novel sensors.

Studies of pH-Sensitive Optical Properties of the deGFP1 Green Fluorescent Protein Using a Unique Polarizable Force Field

The aim of this study is to identify the responsible molecular forms for the pH dependent optical properties of the deGFP1 green fluorescent protein mutant. We have carried out static and dynamic type calculations for all four protonation states of the chromophore to unravel the contributions due to finite temperature and the flexible protein backbone on the pH dependent optical properties. In particular, we have used a combined molecular dynamics and density functional-molecular mechanics linear response approach by means of which the optical property calculations were carried out for the chromophore in the explicitly treated solvent and bioenvironment. Two different models were used to describe the environment-electronic embedding and polarizable electronic embedding-accounting for the polarization of the chromophore and the mutual polarization between the chromophore and the environment, respectively. For this purpose a polarizable force field was derived quantum mechanically for the protein environment by use of analytical response theory. While the gas-phase calculations for the chromophore predict that the induced red shift going from low to high pH is attributed to the change of molecular forms from neutral to zwitterionic, the two more advanced models that explicitly account for the protein backbone attribute the pH shift to a neutral to anionic conversion. Some ramifications of the results for the use of GFPs as pH sensors are discussed.

Multiphoton photochemistry of red fluorescent proteins in solution and live cells

Genetically encoded fluorescent proteins (FPs), and biosensors based on them, provide new insights into how living cells and tissues function. Ultimately, the goal of the bioimaging community is to use these probes deep in tissues and even in entire organisms, and this will require two-photon laser scanning microscopy (TPLSM), with its greater tissue penetration, lower autofluorescence background, and minimum photodamage in the out-of-focus volume. However, the extremely high instantaneous light intensities of femtosecond pulses in the focal volume dramatically increase the probability of further stepwise resonant photon absorption, leading to highly excited, ionizable and reactive states, often resulting in fast bleaching of fluorescent proteins in TPLSM. Here, we show that the femtosecond multiphoton excitation of red FPs (DsRed2 and mFruits), both in solution and live cells, results in a chain of consecutive, partially reversible reactions, with individual rates driven by a high-order (3-5 photon) absorption. The first step of this process corresponds to a three- (DsRed2) or four-photon (mFruits) induced fast isomerization of the chromophore, yielding intermediate fluorescent forms, which then subsequently transform into nonfluorescent products. Our experimental data and model calculations are consistent with a mechanism in which ultrafast electron transfer from the chromophore to a neighboring positively charged amino acid residue triggers the first step of multiphoton chromophore transformations in DsRed2 and mFruits, consisting of decarboxylation of a nearby deprotonated glutamic acid residue.

Theoretical study on molecular design and optical properties of organic sensitizers

A series of organic sensitizers based on different configurations such as D-π-A, D-[π]n-A, D-π-[A]n, [D]n-π-A, D-π-A-π-D, D-π-[A]n-π-D and D-[π-A]n-π-D (where n = 1-4) are designed using theoretical methods. The effects of repeating π-linker, donor-acceptor moieties and the substitution of additional donor-acceptor moieties on the optoelectronic properties have been addressed. Our results show that the strength of the acceptor units changes the mono band absorption into dual band absorption in all the considered strategies. We found that repeating π-linker/donor moieties in the D-π-A series enhances the intensity of the absorption and can tune their absorption spectra from UV-to-visible and visible to near IR regions by repeating acceptor units. Also, the present results indicate that the D-π-A-π-D configuration shows improved optical properties than the conventional D-π-A structure. This theoretical study explores the new configurations and design strategies of organic dyes for developing efficient light harvesting devices working in the whole visible and near IR regions.

Hydrogen bond flexibility correlates with Stokes shift in mPlum variants

Fluorescent proteins have revolutionized molecular biology research and provide a means of tracking subcellular processes with extraordinary spatial and temporal precision. Species with emission beyond 650 nm offer the potential for deeper tissue penetration and lengthened imaging times; however, the origin of their extended Stokes shift is not fully understood. We employed spectrally resolved transient grating spectroscopy and molecular dynamics simulations to investigate the relationship between the flexibility of the chromophore environment and Stokes shift in mPlum. We examined excited state solvation dynamics in a panel of strategic point mutants of residues E16 and I65 proposed to participate in a hydrogen-bonding interaction thought responsible for its red-shifted emission. We observed two characteristic relaxation constants of a few picoseconds and tens of picoseconds that were assigned to survival times of direct and water-mediated hydrogen bonds at the 16-65 position. Moreover, variants of the largest Stokes shift (mPlum, I65V) exhibited significant decay on both time scales, indicating the bathochromic shift correlates with a facile switching between a direct and water-mediated hydrogen bond. This dynamic model underscores the role of environmental flexibility in the mechanism of excited state solvation and provides a template for engineering next-generation red fluorescent proteins.

Electric field effect on the ground state proton transfer in the H-bonded HBDI complex: an implication of the green fluorescent protein

In this paper, first-principles calculations were performed regarding the electric field effect on the ground state proton transfer (GSPT) in the H-bonded p-hydroxybenzylideneimidazolidinone (HBDI) network that represents the active site of the green fluorescent protein (GFP).

Molecular-Level Insight into the Spectral Tuning Mechanism of the DsRed Chromophore

We present a detailed study of the protein environmental effects on the one- and two-photon absorption (1PA and 2PA, respectively) properties of the S0-S1 transition in the DsRed protein using the polarizable embedding density functional theory formalism. We find that steric factors and chromophore-protein interactions act in concert to enhance the 2PA activity inside the protein while adversely blue-shifting the 1PA maximum. A two-state model reveals that the 2PA intensity gain is primarily governed by the increased change in the permanent dipole moment between the ground and the excited states acquired inside the protein. Our results indicate that this mainly is attributable to counter-directional contributions stemming from Lys163 and the conserved Arg95 with the former additionally identified as a key residue in the color tuning mechanism. The results provide new insight into the tuning mechanism of DsRed and suggest a possible strategy for simultaneous improvement of its 1PA and 2PA properties.

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.