Color hues in red fluorescent proteins are due to internal quadratic Stark effect

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.

BF2-Functionalized Benzothiazole Amyloid Markers: Effect of Donor Substituents on One- and Two-Photon Properties

Investigation of amyloids with the aid of fluorescence microscopy provides crucial insights into the development of numerous diseases associated with the formation of aggregates. Here, we present a series of BF2-functionalized benzothiazoles with electron-donating methoxy group(s), which are tested as amyloid fluorescent markers. We evaluate how the position of donor functional group(s) influences optical properties (fluorescence lifetime (τ) and fluorescence quantum yield (FQY)) in a solution and upon binding to amyloids. We elucidate the importance of surrounding environmental factors (hydrogen-bonding network, polarity, and viscosity) on the observed changes in FQY and evaluate how the localization of a donor influences radiative and nonradiative decay pathways. We conclude that a donor attached to the benzothiazole ring contributes to the increment of radiative decay pathways upon binding to amyloids (kr), while the donor attached to the flexible part of a molecule (with rotational freedom) contributes to a decrease in nonradiative decay pathways (knr). We find that the donor-acceptor-donor architecture allows us to obtain 58 times higher FQY of the dye upon binding to bovine insulin amyloids. Finally, we measure two-photon absorption (2PA) cross sections (σ2) of the dyes and their change upon binding by the two-photon excited fluorescence (2PEF) technique. Measurements reveal that dyes that exhibit the increase/decrease of σ2 values when transferred from highly polar solvents to CHCl3 present a similar behavior upon amyloid binding. Our 2PA experimental values are supported by quantum mechanics/molecular mechanics (QM/MM) simulations. Despite this trend, the values of σ2 are not the same, which points out the importance of two-photon absorption measurements of amyloid-dye complexes in order to understand the performance of 2P probes upon binding.

Impact of the Protein Environment on Two-Photon Absorption Cross-Sections of the GFP Chromophore Anion Resolved at the XMCQDPT2 Level of Theory

The search for fluorescent proteins with large two-photon absorption (TPA) cross-sections and improved brightness is required for their efficient use in bioimaging. Here, we explored the impact of a single-point mutation close to the anionic form of the GFP chromophore on its TPA activity. We considered the lowest-energy transition of EGFP and its modification EGFP T203I. We focused on a methodology for obtaining reliable TPA cross-sections for mutated proteins, based on conformational sampling using molecular dynamics simulations and a high-level XMCQDPT2-based QM/MM approach. We also studied the numerical convergence of the sum-over-states formalism and provide direct evidence for the applicability of the two-level model for calculating TPA cross-sections in EGFP. The calculated values were found to be very sensitive to changes in the permanent dipole moments between the ground and excited states and highly tunable by internal electric field of the protein environment. In the case of the GFP chromophore anion, even a single hydrogen bond was shown to be capable of drastically increasing the TPA cross-section. Such high tunability of the nonlinear photophysical properties of the chromophore anions can be used for the rational design of brighter fluorescent proteins for bioimaging using two-photon laser scanning microscopy.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force  field, acting on the CRO. We  find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO-environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism.  We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO-environment polarization in quantitative calculations of absorption spectra in FPs.

The role of hydrogen bonds and electrostatic interactions in enhancing two-photon absorption in green and yellow fluorescent proteins

The spectral properties of fluorescent proteins (FPs) depend on the protein environment of the chromophore (CRO). A deeper understanding of the CRO - environment interactions in terms of FPs spectral characteristics will allow for a rational design of novel markers with desired properties. Here, we are taking a step towards achieving this important goal. With the time-dependent density functional theory (TDDFT), we calculate one- and two-photon absorption (OPA and TPA) spectra for 5 green FPs (GFPs) and 3 yellow FPs (YFPs) differing in amino acid sequence. The goal is to reveal a role of: (i) electrostatic interactions, (ii) hydrogen-bonds (h-bonds), and (iii) h-bonds together with distant electrostatic field in absorption spectra tuning. Our results point to design hypothesis towards FPs optimised for TPA-based applications. Both h-bonds and electrostatic interactions co-operate in enhancing TPA cross-section (σ TPA ) for the S 0 ->S 1   transition in GFPs. Furthermore, it seems that details of h-bonds network in the CRO's vicinity influences σ TPA   response to CRO - environment electrostatic interactions in YFPs. We postulate that engineering FPs with more hydrophilic CRO's environment can lead to greater σ TPA . We also find that removing h-bonds formed with the CRO's phenolate leads to TPA enhancement for transition to higher excited states than S 1 . Particularly Y145 and T203 residues are important in this regard.

Internal Stark effect of single-molecule fluorescence

The optical properties of chromophores can be efficiently tuned by electrostatic fields generated in their close environment, a phenomenon that plays a central role for the optimization of complex functions within living organisms where it is known as internal Stark effect (ISE). Here, we realised an ISE experiment at the lowest possible scale, by monitoring the Stark shift generated by charges confined within a single chromophore on its emission energy. To this end, a scanning tunneling microscope (STM) functioning at cryogenic temperatures is used to sequentially remove the two central protons of a free-base phthalocyanine chromophore deposited on a NaCl-covered Ag(111) surface. STM-induced fluorescence measurements reveal spectral shifts that are associated to the electrostatic field generated by the internal charges remaining in the chromophores upon deprotonation.

The internal Stark effect, a shift of the spectral lines of a chromophore induced by electrostatic fields in its close environment, plays an important role in nature. Here the authors observe a Stark shift in the fluorescence spectrum of a phthalocyanine molecule upon charge modifications within the molecule itself, achieved by sequential removal of the central protons with a STM tip.

Modeling Spectral Tuning in Red Fluorescent Proteins Using the Dipole Moment Variation upon Excitation

We describe a model for spectral tuning in red fluorescent proteins (RFPs) based on the relation between an electronic structure descriptor, the dipole moment variation upon excitation (DMV), and the excitation energy of a protein. This approach aims to overcome the problem of accurate prediction of excitation energies in RFPs, which span a very narrow window of band maxima. The latter roughly corresponds to the energy range of 0.1 eV, which is comparable with typical errors in calculations of the excitation energy by conventional quantum chemistry methods. In this work, we demonstrate a strong quantitative correlation between DMV values, obtained computationally with modest efforts, and excitation energies ΔE_ex at the experimental excitation band maxima for a series of RFPs with bands between 570 and 605 nm. Protein models are constructed by motifs of the relevant crystal structures, and atomic coordinates are optimized in quantum mechanics/molecular mechanics (QM/MM) calculations with QM-subsystems composed of large chromophore-containing regions. DMV values are evaluated with the electron density computed at the time-dependent density functional theory (TDDFT) level using several functionals and basis sets. We show that the results obtained with the CAM-B3LYP, BHHLYP, and M06-2X functionals demonstrate favorable correlations between DMV and ΔE_ex with the mean absolute error less than 0.01 eV. Taking into account the solid theoretical grounds of the relation between the DMV and the excitation energy in fluorescent proteins, the described modeling strategy presents a rational tool for spectral tuning in these efficient markers for in vivo imaging.

Structural Origins of Altered Spectroscopic Properties upon Ligand Binding in Proteins Containing a Fluorescent Noncanonical Amino Acid

Fluorescent noncanonical amino acids (fNCAAs) could serve as starting points for the rational design of protein-based fluorescent sensors of biological activity. However, efforts toward this goal are likely hampered by a lack of atomic-level characterization of fNCAAs within proteins. Here, we describe the spectroscopic and structural characterization of five streptavidin mutants that contain the fNCAA l-(7-hydroxycoumarin-4-yl)ethylglycine (7-HCAA) at sites proximal to the binding site of its substrate, biotin. Many of the mutants exhibited altered fluorescence spectra in response to biotin binding, which included both increases and decreases in fluorescence intensity as well as red- or blue-shifted emission maxima. Structural data were also obtained for three of the five mutants. The crystal structures shed light on interactions between 7-HCAA and functional groups, contributed either by the protein or by the substrate, that may be responsible for the observed changes in the 7-HCAA spectra. These data could be used in future studies aimed at the rational design of fluorescent, protein-based sensors of small molecule binding or dissociation.

Local Electric Field Controls Fluorescence Quantum Yield of Red and Far-Red Fluorescent Proteins

Genetically encoded probes with red-shifted absorption and fluorescence are highly desirable for imaging applications because they can report from deeper tissue layers with lower background and because they provide additional colors for multicolor imaging. Unfortunately, red and especially far-red fluorescent proteins have very low quantum yields, which undermines their other advantages. Elucidating the mechanism of nonradiative relaxation in red fluorescent proteins (RFPs) could help developing ones with higher quantum yields. Here we consider two possible mechanisms of fast nonradiative relaxation of electronic excitation in RFPs. The first, known as the energy gap law, predicts a steep exponential drop of fluorescence quantum yield with a systematic red shift of fluorescence frequency. In this case the relaxation of excitation occurs in the chromophore without any significant changes of its geometry. The second mechanism is related to a twisted intramolecular charge transfer in the excited state, followed by an ultrafast internal conversion. The chromophore twisting can strongly depend on the local electric field because the field can affect the activation energy. We present a spectroscopic method of evaluating local electric fields experienced by the chromophore in the protein environment. The method is based on linear and two-photon absorption spectroscopy, as well as on quantum-mechanically calculated parameters of the isolated chromophore. Using this method, which is substantiated by our molecular dynamics simulations, we obtain the components of electric field in the chromophore plane for seven different RFPs with the same chromophore structure. We find that in five of these RFPs, the nonradiative relaxation rate increases with the strength of the field along the chromophore axis directed from the center of imidazolinone ring to the center of phenolate ring. Furthermore, this rate depends on the corresponding electrostatic energy change (calculated from the known fields and charge displacements), in quantitative agreement with the Marcus theory of charge transfer. This result supports the dominant role of the twisted intramolecular charge transfer mechanism over the energy gap law for most of the studied RFPs. It provides important guidelines of how to shift the absorption wavelength of an RFP to the red, while keeping its brightness reasonably high.

Conformational analysis by UV spectroscopy: the decisive contribution of environment-induced electronic Stark effects

UV chromophores are frequently used as probes of the molecular structure. In particular, they are sensitive to the electric field generated by the molecular environment, resulting in the observation of Stark effects on UV spectra. While these environment-induced electronic Stark effects (EI-ESE) are already used for conformational analysis in the condensed phase, this work explores the potential of such an approach when performed at much higher conformational resolution in the gas phase. By investigating model alkali benzylacetate and 4-phenylbutyrate ion pairs, where the electric field applied to the phenyl ring is chemically tuned by changing the nature of the alkali cation, this work demonstrates that precise conformational assignments can be proposed based on the correlation between the conformation-dependent calculated electric fields and the frequency of the electronic transitions observed in the experimental UV spectra. Remarkably, the sole analysis of Stark effects and fragmentation patterns in mass-selected UV spectra provided an accurate and complete conformational analysis, where spectral differences as small as a few cm-1 between electronic transitions were rationalized. This case study illustrates that the identification of EI-ESE together with their interpretation at the modest cost of a ground state electric field calculation qualify UV spectroscopy as a powerful tool for conformational analysis.

Dipole Moment Variation Clears Up Electronic Excitations in the π-Stacked Complexes of Fluorescent Protein Chromophores

We propose a quantitative structure-property relationship (QSPR) model for prediction of spectral tuning in cyan, green, orange, and red fluorescent proteins, which are engineered by motifs of the green fluorescent protein. Protein variants, in which their chromophores are involved in the π-stacking interaction with amino acid residues tyrosine, phenylalanine, and histidine, are prospective markers useful in bioimaging and super-resolution microscopy. In this work, we constructed training sets of the π-stacked complexes of four fluorescent protein chromophores (of the green, orange, red, and cyan series) with various substituted benzenes and imidazoles and tested the use of dipole moment variation upon excitation (DMV) as a descriptor to evaluate the vertical excitation energies in these systems. To validate this approach, we computed and analyzed electron density distributions of the π-stacked complexes and correlated the QSPR predictions with the reference values of the transition energies obtained using the high-level ab initio quantum chemistry methods. According to our results, the use of the DMV descriptor allows one to predict excitation energies in the π-stacked complexes with errors not exceeding 0.1 eV, which makes this model a practically useful tool in the development of efficient fluorescent markers for in vivo imaging.

What is the Optimal Size of the Quantum Region in Embedding Calculations of Two-Photon Absorption Spectra of Fluorescent Proteins?

We systematically investigate an impact of the size and content of a quantum (QM) region, treated at the density functional theory level, in embedding calculations on one- (OPA) and two-photon absorption (TPA) spectra of the following fluorescent proteins (FPs) models: Aequorea victoria green FP (avGFP) with neutral (avGFP-n) and anionic (avGFP-a) chromophore as well as Citrine FP. We find that amino acid (a.a.) residues as well as water molecules hydrogen-bonded (h-bonded) to the chromophore usually boost both OPA and TPA processes intensity. The presence of hydrophobic a.a. residues in the quantum region also non-negligibly affects both absorption spectra but decreases absorption intensity. We conclude that to reach a quantitative description of OPA and TPA spectra in multiscale modeling of FPs, the quantum region should consist of a chromophore and most of a.a. residues and water molecules in a radius of 0.30–0.35 nm (ca. 200–230 atoms) when the remaining part of the system is approximated by the electrostatic point-charges. The optimal size of the QM region can be reduced to 80–100 atoms by utilizing a more advanced polarizable embedding model. We also find components of the QM region that are specific to a FP under study. We propose that the F165 a.a. residue is important in tuning the TPA spectrum of avGFP-n but not other investigated FPs. In the case of Citrine, Y203 and M69 a.a. residues must definitely be part of the QM subsystem. Furthermore, we find that long-range electrostatic interactions between the QM region and the rest of the protein cannot be neglected even for the most extensive QM regions (ca. 350 atoms).

Electronic Stark Effect in Isolated Ion Pairs

Stark spectral shifts of a molecular probe are commonly used to estimate the local electric field in condensed media. The very large fields reported, typically in the 0.1-10 GV m^-1 range, are, however, difficult to reproduce in a controlled manner, limiting the calibration of these molecular probes to ranges below 0.1 GV m^-1. In this context, we investigated gas-phase, isolated, molecular ion pairs, where a phenyl ring is immersed in the electric field produced by the nearby ionic groups. The intensity of the electric field is chemically tuned in the 1 GV m^-1 range by changing the nature of the cations, and the phenyl ring response is monitored by UV spectroscopy. A quadratic Stark effect is observed, demonstrating the possibility to characterize molecular probes in a solvent-free environment and in the very large field range they typically meet in condensed media such as biological environments.

The key to the yellow-to-cyan tuning in the green fluorescent protein family is polarisation

Computational approaches have to date failed to fully capture the large (about 0.4 eV) excitation energy tuning displayed by the nearly identical anionic chromophore in different green fluorescent protein (GFP) variants. Here, we present a thorough comparative study of a set of proteins in this sub-family, including the most red- (phiYFP) and blue-shifted (mTFP0.7) ones. We employ a classical polarisable embedding through induced dipoles and combine it with time-dependent density functional theory and multireference perturbation theory in order to capture both state-specific induction contributions and the coupling of the polarisation of the protein to the chromophore transition density. The obtained results show that only upon inclusion of both these two effects generated by the mutual polarisation between the chromophore and the protein can the full spectral tuning be replicated. We finally discuss how this mutual polarisation affects the correlation between excitation energies, dipole moment variation, and molecular electrostatic field.

Quantifying Acute Fuel and Respiration Dependent pH Homeostasis in Live Cells Using the mCherryTYG Mutant as a Fluorescence Lifetime Sensor

Intracellular pH plays a key role in physiology, and its measurement in living specimens remains a crucial task in biology. Fluorescent protein-based pH sensors have gained widespread use, but there is limited spectral diversity for multicolor detection, and it remains a challenge to measure absolute pH values. Here we demonstrate that mCherryTYG is an excellent fluorescence lifetime pH sensor that significantly expands the modalities available for pH quantification in live cells. We first report the 1.09 Å X-ray crystal structure of mCherryTYG, exhibiting a fully matured chromophore. We next determine that it has an extraordinarily large dynamic range with a 2 ns lifetime change from pH 5.5 to 9.0. Critically, we find that the sensor maintains a p Ka of 6.8 independent of environment, whether as the purified protein in solution or expressed in live cells. Furthermore, the lifetime measurements are robustly independent of total fluorescence intensity and scatter. We demonstrate that mCherryTYG is a highly effective sensor using time-resolved fluorescence spectroscopy on live-cell suspensions, which has been previously overlooked as an easily accessible approach for quantifying intracellular pH. As a red fluorescent sensor, we also demonstrate that mCherryTYG is spectrally compatible with the ATeam sensor and EGFP for simultaneous dual-color measurements of intracellular pH, ATP, and extracellular pH. In a proof-of-concept, we quantify acute respiration-dependent pH homeostasis that exhibits a stoichiometric relationship with the ATP-generating capacity of the carbon fuel choice in E. coli. Broadly speaking, our work presents a previously unemployed methodology that will greatly facilitate continuous pH quantification.

Effects of Intrinsic Local Fields on Molecular Fluorescence and Energy Transfer: Dipole Mechanisms and Surface Potentials

A general theory is developed to identify the influence of local dipole fields on fluorescence and intermolecular electronic excitation transfer. The analysis, based on electrodynamical principles, identifies the fundamental quantum mechanisms and delivers full analytical results. The aim is to afford new physical insights, assisting the interpretation of measurements on the specific effects of local molecular dipoles on direct fluorescence and on fluorescence resonance energy transfer. Dipole field effects, which include those originating from intrinsically polar chromophores and surface field gradients, also prove to be manifest in electronic transitions of quadrupole symmetry character. The results have particular significance for fluorescence studies of cell membrane biophysics.

The whither of bacteriophytochrome-based near-infrared fluorescent proteins: Insights from two-photon absorption spectroscopy

We present one- and two-photon-absorption fluorescence spectroscopic analysis of biliverdin (BV) chromophore-based single-domain near-infrared fluorescent proteins (iRFPs). The results of these studies are used to estimate the internal electric fields acting on BV inside iRFPs and quantify the electric dipole properties of this chromophore, defining the red shift of excitation and emission spectra of BV-based iRFPs. The iRFP studied in this work is shown to fit well the global diagram of the red-shift tunability of currently available BV-based iRFPs as dictated by the quadratic Stark effect, suggesting the existence of the lower bound for the strongest red shifts attainable within this family of fluorescent proteins. The absolute value of the two-photon absorption (TPA) cross section of a fluorescent calcium sensor based on the studied iRFP is found to be significantly larger than the TPA cross sections of other widely used genetically encodable fluorescent calcium sensors.

Understanding the Fluorescence Change in Red Genetically Encoded Calcium Ion Indicators

For over 20 years, genetically encoded Ca²⁺ indicators have illuminated dynamic Ca²⁺ signaling activity in living cells and, more recently, whole organisms. We are just now beginning to understand how they work. Various fluorescence colors of these indicators have been developed, including red. Red ones are promising because longer wavelengths of light scatter less in tissue, making it possible to image deeper. They are engineered from a red fluorescent protein that is circularly permuted and fused to a Ca²⁺-sensing domain. When Ca²⁺ binds, a conformational change in the sensing domain causes a change in fluorescence. Three factors can contribute to this fluorescence change: 1) a shift in the protonation equilibrium of the chromophore, 2) a change in fluorescence quantum yield, and 3) a change in the extinction coefficient or the two-photon cross section, depending on if it is excited with one or two photons. Here, we conduct a systematic study of the photophysical properties of a range of red Ca²⁺ indicators to determine which factors are the most important. In total, we analyzed nine indicators, including jRGECO1a, K-GECO1, jRCaMP1a, R-GECO1, R-GECO1.2, CAR-GECO1, O-GECO1, REX-GECO1, and a new variant termed jREX-GECO1. We find that these could be separated into three classes that each rely on a particular set of factors. Furthermore, in some cases, the magnitude of the change in fluorescence was larger with two-photon excitation compared to one-photon because of a change in the two-photon cross section, by up to a factor of two.

Dendrimeric Nanoparticles for Two-Photon Photodynamic Therapy and Imaging: Synthesis, Photophysical Properties, Innocuousness in Daylight and Cytotoxicity under Two-Photon Irradiation in the NIR

The synthesis and the photophysical properties of a new class of fully organic monodisperse nanoparticles for combined two-photon imaging and photodynamic therapy are described. The design of such nanoparticles is based on the covalent immobilization of a dedicated quadrupolar dye that combines excellent two-photon absorbing (2PA) properties, fluorescence and singlet oxygen generation ability, in a phosphorous-based dendrimeric architecture. First, a bifunctional quadrupolar dye bearing two different grafting moieties, a phenol function and an aldehyde function, was synthesized. It was then covalently grafted through its phenol function to a phosphorus-based dendrimer scaffold of generation 1. The remaining aldehyde functions were then used to continue the dendrimer synthesis up to generation 2, introducing finally 24 water-solubilizing triethyleneglycol chains at its periphery. A dendrimer confining 12 photoactive quadrupolar units in its inner scaffold and showing water solubility was thus obtained. Interestingly, the G1 and G2 dendrimers retain some fluorescence as well as significant singlet oxygen production efficiencies while they were found to show very high 2PA cross-sections in a broad range of the NIR biological spectral window. Hydrophilic dendrimer G2 was tested in vitro on breast cancer cells, first in one- and two-photon microscopy, which allowed for visualization of their cell internalization, then in two-photon photodynamic therapy. While being nontoxic in the dark and, more importantly, under exposure to daylight, dendrimer G2 proved to be a very efficient cell-death inducer only under two-photon irradiation in the NIR.

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.

Two photon spectroscopy and microscopy of the fluorescent flavoprotein, iLOV

Homans et al. show that engineered LOV-domains are amenable to two photon activation both in vitro and in human cells.

Blue-Shifted Green Fluorescent Protein Homologues Are Brighter than Enhanced Green Fluorescent Protein under Two-Photon Excitation

Fluorescent proteins (FPs) are indispensable markers for two-photon imaging of live tissue, especially in the brains of small model organisms. The quantity of physiologically relevant data collected, however, is limited by heat-induced damage of the tissue due to the high intensities of the excitation laser. We seek to minimize this damage by developing FPs with improved brightness. Among FPs with the same chromophore structure, the spectral properties can vary widely due to differences in the local protein environment. Using a physical model that describes the spectra of FPs containing the anionic green FP (GFP) chromophore, we predict that those that are blue-shifted in one-photon absorption will have stronger peak two-photon absorption cross sections. Following this prediction, we present 12 blue-shifted GFP homologues and demonstrate that they are up to 2.5 times brighter than the commonly used enhanced GFP (EGFP).

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.

Bright blue-shifted fluorescent proteins with Cys in the GAF domain engineered from bacterial phytochromes: fluorescence mechanisms and excited-state dynamics

Near-infrared fluorescent proteins (NIR FPs) engineered from bacterial phytochromes (BphPs) are of great interest for in vivo imaging. They utilize biliverdin (BV) as a chromophore, which is a heme degradation product, and therefore they are straightforward to use in mammalian tissues. Here, we report on fluorescence properties of NIR FPs with key alterations in their BV binding sites. BphP1-FP, iRFP670 and iRFP682 have Cys residues in both PAS and GAF domains, rather than in the PAS domain alone as in wild-type BphPs. We found that NIR FP variants with Cys in the GAF or with Cys in both PAS and GAF show blue-shifted emission with long fluorescence lifetimes. In contrast, mutants with Cys in the PAS only or no Cys residues at all exhibit red-shifted emission with shorter lifetimes. Combining these results with previous biochemical and BphP1-FP structural data, we conclude that BV adducts bound to Cys in the GAF are the origin of bright blue-shifted fluorescence. We propose that the long fluorescence lifetime follows from (i) a sterically more constrained thioether linkage, leaving less mobility for ring A than in canonical BphPs, and (ii) that π-electron conjugation does not extend on ring A, making excited-state deactivation less sensitive to ring A mobility.

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.

Prediction of two-photon absorption enhancement in red fluorescent protein chromophores made from non-canonical amino acids

Two-photon spectroscopy of fluorescent proteins is a powerful bio-imaging tool known for deep tissue penetration and little cellular damage. Being less sensitive than the one-photon microscopy alternatives, a protein with a large two-photon absorption (TPA) cross-section is needed. Here, we use time-dependent density functional theory (TD-DFT) at the B3LYP and CAM-B3LYP/6-31+G(d,p) levels of theory to screen twenty-two possible chromophores that can be formed upon replacing the amino-acid Tyr66 that forms the red fluorescent protein (RFP) chromophore with a non-canonical amino acid. The two-level model for TPA was used to assess the properties (i.e., transition dipole moment, permanent dipole moment difference, and the angle between them) leading to the TPA cross-sections determined via response theory. Computing TPA cross-sections with B3LYP and CAM-B3LYP yields similar overall trends. Results using both functionals agree that the RFP-derived model of the Gold Fluorescent Protein chromophore (Model 20) has the largest intrinsic TPA cross-section at the optimized geometry. TPA was further computed for selected chromophores following conformational changes: variation of both the dihedral angle of the acylimine moiety and the tilt and twist angles between the rings of the chromophore. The TPA cross-section assumed an oscillatory trend with the rotation of the acylimine dihedral, and the TPA is maximized in the planar conformation for almost all models. Model 21 (a hydroxyquinoline derivative) is shown to be comparable to Model 20 in terms of TPA cross-section. The conformational study on Model 21 shows that the acylimine angle has a much stronger effect on the TPA than its tilt and twist angles. Having an intrinsic TPA ability that is more than 7 times that of the native RFP chromophore, Models 20 and 21 appear to be very promising for future experimental investigation.

Local electric fields and molecular properties in heterogeneous environments through polarizable embedding

In spectroscopies, the local field experienced by a molecule embedded in an environment will be different from the externally applied electromagnetic field, and this difference may significantly alter the response and transition properties of the molecule. The polarizable embedding (PE) model has previously been developed to model the local field contribution stemming from the direct molecule-environment coupling of the electromagnetic response properties of molecules in solution as well as in heterogeneous environments, such as proteins. Here we present an extension of this approach to address the additional effective external field effect, i.e., the manifestations of the environment polarization induced by the external field, which allows for the calculation of properties defined in terms of the external field. Within a response framework, we report calculations of the one- and two-photon absorption (1PA and 2PA, respectively) properties of PRODAN-methanol clusters as well as the fluorescent protein DsRed. Our results demonstrate the necessity of accounting for both the dynamical reaction field and effective external field contributions to the local field in order to reproduce full quantum chemical reference calculations. For the lowest π→π* transition in DsRed, inclusion of effective external field effects gives rise to a 1.9- and 3.5-fold reduction in the 1PA and 2PA cross-sections, respectively. The effective external field is, however, strongly influenced by the heterogeneity of the protein matrix, and the resulting effect can lead to either screening or enhancement depending on the nature of the transition under consideration.

Open-ended response theory with polarizable embedding: multiphoton absorption in biomolecular systems

We present the theory and implementation of an open-ended framework for electric response properties that includes effects from the molecular environment modeled by the polarizable embedding model.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

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.

Ultrafast excited-state dynamics and fluorescence deactivation of near-infrared fluorescent proteins engineered from bacteriophytochromes

Near-infrared fluorescent proteins, iRFPs, are recently developed genetically encoded fluorescent probes for deep-tissue in vivo imaging. Their functions depend on the corresponding fluorescence efficiencies and electronic excited state properties. Here we report the electronic excited state deactivation dynamics of the most red-shifted iRFPs: iRFP702, iRFP713 and iRFP720. Complementary measurements by ultrafast broadband fluorescence and absorption spectroscopy show that single exponential decays of the excited state with 600~700 ps dominate in all three iRFPs, while photoinduced isomerization was completely inhibited. Significant kinetic isotope effects (KIE) were observed with a factor of ~1.8 in D2O, and are interpreted in terms of an excited-state proton transfer (ESPT) process that deactivates the excited state in competition with fluorescence and chromophore mobility. On this basis, new approaches for rational molecular engineering may be applied to iRFPs to improve their fluorescence.

Two-photon directed evolution of green fluorescent proteins

Directed evolution has been used extensively to improve the properties of a variety of fluorescent proteins (FPs). Evolutionary strategies, however, have not yet been used to improve the two-photon absorption (2PA) properties of a fluorescent protein, properties that are important for two-photon imaging in living tissues, including the brain. Here we demonstrate a technique for quantitatively screening the two-photon excited fluorescence (2PEF) efficiency and 2PA cross section of tens of thousands of mutant FPs expressed in E. coli colonies. We use this procedure to move EGFP through three rounds of two-photon directed evolution leading to new variants showing up to a 50% enhancement in peak 2PA cross section and brightness within the near-IR tissue transparency wavelength range.

Exploring color tuning strategies in red fluorescent proteins

Red-emitting fluorescent proteins (RFPs) with fluorescence emission above 600 nm are advantageous for cell and tissue imaging applications for various reasons. Fluorescence from an RFP is well separated from cellular autofluorescence, which is in the green region of the spectrum, and red light is scattered less, which allows thicker specimens to be imaged. Moreover, the phototoxic response of cells is lower for red than blue or green light exposure. Further red-shifted FP variants can be obtained by genetic modifications causing an extension of the conjugated π-electron system of the chromophore, or by placing amino acids near the chromophore that stabilize its excited state or destabilize its ground state. We have selected the tetrameric RFP eqFP611 from Entacmaea quadricolor as a lead structure and discuss several rational design trials to generate RFP variants with improved photochemical properties.

Advances in fluorescence labeling strategies for dynamic cellular imaging

Synergistic advances in optical physics, probe design, molecular biology, labeling techniques and computational analysis have propelled fluorescence imaging into new realms of spatiotemporal resolution and sensitivity. This review aims to discuss advances in fluorescent probes and live-cell labeling strategies, two areas that remain pivotal for future advances in imaging technology. Fluorescent protein- and bio-orthogonal-based methods for protein and RNA imaging are discussed as well as emerging bioengineering techniques that enable their expression at specific genomic loci (for example, CRISPR and TALENs). Important attributes that contribute to the success of each technique are emphasized, providing a guideline for future advances in dynamic live-cell imaging.

Two-photon fluorescence spectroscopy and imaging of 4-dimethylaminonaphthalimide peptide and protein conjugates

We report detailed photophysical studies on the two-photon fluorescence processes of the solvatochromic fluorophore 4-DMN as a conjugate of the calmodulin (CaM) and the associated CaM-binding peptide M13. Strong two-photon fluorescence enhancement has been observed which is associated with calcium binding. It is found that the two-photon absorption cross-section is strongly dependent on the local environment surrounding the 4-DMN fluorophore in the CaM conjugates, providing sensitivity between sites of fluorophore attachment. Utilizing time-resolved measurements, the emission dynamics of 4-DMN under various environmental (solvent) conditions are analyzed. In addition, anisotropy measurements reveal that the 4-DMN-S38C-CaM system has restricted rotation in the calcium-bound calmodulin. To establish the utility for cellular imaging, two-photon fluorescence microscopy studies were also carried out with the 4-DMN-modified M13 peptide in cells. Together, these studies provide strong evidence that 4-DMN is a useful probe in two-photon imaging, with advantageous properties for cellular experiments.

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.

All-Optical Sensing of the Components of the Internal Local Electric Field in Proteins

Here, we present a new all-optical method of interrogation of the internal electric field vector inside proteins. The method is based on experimental evaluation of the permanent dipole moment change upon excitation and the pure electronic transition frequency of a fluorophore embedded in a protein matrix. The permanent dipole moment change can be obtained from two-photon absorption measurements. In addition, permanent dipole moment change, tensor of polarizability change, and transition frequency for the free chromophore should be calculated quantum-mechanically. This allows obtaining the components of the electric field by considering the second-order Stark shift. We use the fluorescent protein mCherry as an example to demonstrate the applicability of the method.

Primary role of the chromophore bond length alternation in reversible photoconversion of red fluorescence proteins

Rapid photobleaching of fluorescent proteins can limit their use in imaging applications. The underlying kinetics is multi-exponential and strongly depends on the local chromophore environment. The first, reversible, step may be attributed to a rotation around one of the two exocyclic C-C bonds bridging phenol and imidazolinone groups in the chromophore. However it is not clear how the protein environment controls this motion - either by steric hindrances or by modulating the electronic structure of the chromophore through electrostatic interactions. Here we study the first step of the photobleaching kinetics in 13 red fluorescent proteins (RFPs) with different chromophore environment and show that the associated rate strongly correlates with the bond length alternation (BLA) of the two bridge bonds. The sign of the BLA appears to determine which rotation is activated. Our results present experimental evidence for the dominance of electronic effects in the conformational dynamics of the RFP chromophore.

Cell-based and in vivo spectral analysis of fluorescent proteins for multiphoton microscopy

Multiphoton microscopy of cells and subcellular structures labeled with fluorescent proteins is the state-of-the-art technology for longitudinal imaging studies in tissues and living animals. Successful analysis of separate cell populations or signaling events by intravital microscopy requires optimal pairing of multiphoton excitation wavelengths with spectrally distinct fluorescent proteins. While prior studies have analyzed two photon absorption properties of isolated fluorescent proteins, there is limited information about two photon excitation and fluorescence emission profiles of fluorescent proteins expressed in living cells and intact tissues. Multiphoton microscopy was used to analyze fluorescence outputs of multiple blue, green, and red fluorescent proteins in cultured cells and orthotopic tumor xenografts of human breast cancer cells. It is shown that commonly used orange and red fluorescent proteins are excited efficiently by 750 to 760 nm laser light in living cells, enabling dual color imaging studies with blue or cyan proteins without changing excitation wavelength. It is also shown that small incremental changes in excitation wavelength significantly affect emission intensities from fluorescent proteins, which can be used to optimize multi-color imaging using a single laser wavelength. These data will direct optimal selection of fluorescent proteins for multispectral two photon microscopy.

Describing two-photon absorptivity of fluorescent proteins with a new vibronic coupling mechanism

Fluorescent proteins (FPs) are widely used in two-photon microscopy as genetically encoded probes. Understanding the physical basics of their two-photon absorption (2PA) properties is therefore crucial for creation of two-photon brighter mutants. On the other hand, it can give us better insight into molecular interactions of the FP chromophore with a complex protein environment. It is known that, compared to the one-photon absorption spectrum, where the pure electronic transition is the strongest, the 2PA spectrum of a number of FPs is dominated by a vibronic transition. The physical mechanism of such intensity redistribution is not understood. Here, we present a new physical model that explains this effect through the "Herzberg-Teller"-type vibronic coupling of the difference between the permanent dipole moments in the ground and excited states (Δμ) to the bond-length-alternating coordinate. This model also enables us to quantitatively describe a large variability of the 2PA peak intensity in a series of red FPs with the same chromophore through the interference between the "Herzberg-Teller" and Franck-Condon terms.