Electronic excitations of green fluorescent proteins: modeling solvatochromatic shifts of red fluorescent protein chromophore model compound in aqueous solutions

While green fluorescent proteins (GFPs) have been widely used as tools in biochemistry, cell biology, and molecular genetics, novel red fluorescent proteins (RFPs) with red fluorescence emission have also been identified, as complements to the existing GFP technology. The unusual spectrophotometric and fluorescence properties of GFPs and RFPs are controlled by the protonation states and possibly cis/trans isomerization within their chromophores. In this work, we have investigated the electronic structures, liquid structures, and solvent shifts of the possible neutral and anionic protonated states and the cis/trans isomerization of a RFP chromophore model compound HBMPI in aqueous solutions. The calculations reproduced the experimental absorption solvatochromatic shifts of dilute HBMPI in water under neutral and anionic conditions. Unlike the GFP chromophore, the RFP chromophore model compound HBMPI in basic solution can only adopt a conformation where the C=C bond between the bridge group and the imidazolinone ring and the C-C bond between the imidazolinone and ethylene groups exist in cis and trans conformations, respectively. Moreover, the solvent-solute hydrogen-bonding interactions are found to contribute significantly to the total solvent shifts of pi-pi* excitations of aqueous HBMPI solutions, signifying the importance of protein environment in the determination of the conformation of the chromophores in red fluorescent proteins.     

Insight into the common mechanism of the chromophore formation in the red fluorescent proteins: the elusive blue intermediate revealed

Understanding the chromophore maturation process in fluorescent proteins is important for the design of proteins with improved properties. Here, we present the results of electronic structure calculations identifying the nature of a blue intermediate, a key species in the process of the red chromophore formation in DsRed, TagRFP, fluorescent timers, and PAmCherry. The chromophore of the blue intermediate has a structure in which the π-system of the imidazole ring is extended by the acylimine bond, which can be represented by the model N-[(5-hydroxy-1H-imidazole-2yl)methylidene]acetamide (HIMA) compound. Ab initio and QM/MM calculations of the isolated model and protein-bound (mTagBFP) chromophores identify the anionic form of HIMA as the only structure that has absorption that is consistent with the experiment and is stable in the protein binding pocket. The anion and zwitterion are the only protonation forms of HIMA whose absorption (421 and 414 nm, or 2.95 and 3.00 eV) matches the experimental spectrum of the blue form in DsRed (the absorption maximum is 408 nm or 3.04 eV) and mTagBFP (400 nm or 3.10 eV). The QM/MM optimization of the protein-bound anionic form results in a structure that is close to the X-ray one, whereas the zwitterionic chromophore is unstable in the protein binding pocket and undergoes prompt proton transfer. The computed excitation energy of the protein-bound anionic form of the mTagBFP-like chromophore (3.04 eV) agrees with the experimental absorption spectrum of the protein. The DsRed-like chromophore formation in red fluorescent proteins is revisited on the basis of ab initio results and verified by directed mutagenesis revealing a key role of the amino acid residue 70, which is the second after the chromophore tripeptide, in the formation process.     

Competition between energy and proton transfer in ultrafast excited-state dynamics of an oligomeric fluorescent protein red Kaede

We investigated femtosecond and picosecond time-resolved fluorescence dynamics of a tetrameric fluorescent protein Kaede with a red chromophore (red Kaede) to examine a relationship between the excited-state dynamics and a quaternary structure of the fluorescent protein. Red Kaede was obtained by photoconversion from green Kaede that was cloned from a stony coral Trachyphyllia geoffroyi. In common with other typical fluorescent proteins, a chromophore of red Kaede has two protonation states, the neutral and the anionic forms in equilibrium. Time-resolved fluorescence measurements clarified that excitation of the neutral form gives the anionic excited state with a time constant of 13 ps at pH 7.5. This conversion process was attributed to fluorescence resonance energy transfer (FRET) from the photoexcited neutral form to the ground-state anionic form that is located in an adjacent subunit in the tetramer. The time-resolved fluorescence data measured at different pH revealed that excited-state proton transfer (ESPT) also occurs with a time constant of 300 ps and hence that the FRET and ESPT take place simultaneously in the fluorescent protein as competing processes. The ESPT rate in red Kaede was significantly slower than the rate in Aequorea GFP, which highly likely arises from the different hydrogen bond network around the chromophore.     

Electronic excitations of green fluorescent proteins: protonation states of chromophore model compound in solutions

Green fluorescent proteins (GFPs) are widely used as tools in biochemistry, cell biology, and molecular genetics due to their unusual optical spectroscopic characteristics. The spectrophotometric and fluorescence properties of GFPs are controlled by the protonation states and possibly cis-trans isomerization of the chromophore (p-hydroxybenzylideneimidazolinone). In this work, we have investigated electronic structures, liquid structures, and solvent shifts of the three possible protonated states (neutral, anionic, and zwitterionic) and their cis-trans isomerization of a model compound 4'-hydroxybenzylidene-2-methyl-imidazolin-5-one-3-acetate (HBMIA) in aqueous solutions. Our calculated results suggest that HBMIA could adopt both cis and trans conformations in a solution, and it exists in three different protonation states depending on the pH conditions. The absorption spectrum observed in neutral solution is thus assigned to the electronic excitations within the neutral form and the cis isomer of the zwitterionic form, while the absorption band at 425 nm in the basic solution is due to the excitations within the anionic form and the trans isomer of the zwitterionic form. Some technical problems related to the computation of electronic excitations within the HBMIA at the anionic state are also discussed.     

Solvation of p-coumaric acid in water

The photoactive yellow protein (PYP) is an important model protein for many (photoactive) signaling proteins. Key steps in the PYP photocycle are the isomerization and protonation of its chromophore, p-coumaric acid (pCA). In the ground state of the protein, this chromophore is in the trans configuration with its phenolic oxygen deprotonated. For this paper, we studied four different configurations of pCA solvated in water with ab initio molecular dynamics simulations as implemented in CP2K/Quickstep. We researched the influence of the protonation and isomerization state of pCA on its hydrogen-bonding properties and on the Mulliken charges of the atoms in the simulation. The chromophore isomerization state influenced the hydrogen-bonding less than its protonation state. In general, more charge yielded a higher hydrogen-bond coordination number. Where deprotonation increases both the coordination number and the residence time of the water molecules around the chromophore, protonation showed a somewhat lower coordination number on two of the three pCA oxygens but much higher residence times on all of them. This could be explained by the increased polarization of the OH groups of the molecule. The presence of the chromophore also influenced the charge and polarization of the water molecules around it. This effect was different in the four systems studied and mainly localized in the first solvation shell. We also performed a proton-transfer reaction from hydronium through various other water molecules to the chromophore. In this small charge-separated system, the protonation occurred within 6.5 ps. We identified the transition state for the final step in this protonation series.     

Modeling photoabsorption of the asFP595 chromophore

The fluorescent protein asFP595 is a promising photoswitchable biomarker for studying processes in living cells. We present the results of a high level theoretical study of photoabsorption properties of the model asFP595 chromophore molecule in biologically relevant protonation states: anionic, zwitterionic, and neutral. Ground state equilibrium geometry parameters are optimized in the PBE0/(aug)-cc-pVDZ density functional theory approximation. An augmented version of multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) following the state-averaged CASSCF/(aug)-cc-pVDZ calculations is used to estimate the vertical S0-S1 excitation energies for all chromophore species. An accuracy of this approach is validated by comparing the computed estimates of the S0-S1 absorption maximum of the closely related chromophore from the DsRed protein to the known experimental value in the gas phase. An influence of the CASSCF active space on the aug-MCQDPT2 excitation energies is analyzed. The zwitterionic form of the asFP595 chromophore is found to be the most sensitive to a particular choice and amount of active orbitals. This observation is explained by the charge-transfer type of the S0-S1 transition involving the entire conjugated pi-electron system for the zwitterionic protonation state. According to the calculation results, the anionic form in the trans conformation is found to possess the most red-shifted absorption band with the maximum located at 543 nm. The bands of the zwitterionic and neutral forms are considerably blue-shifted compared to those of the anionic form. These conclusions are at variance with the results obtained in the TDDFT approximation for the asFP595 chromophore. The absorption wavelengths computed in the aug-MCQDPT2/CASSCF theory are as follows: 543 (535), 470 (476), and 415 (417) nm for the anionic, zwitterionic, and neutral forms of the trans and cis (in parentheses) isomers of the asFP595 chromophore. These data can be used as a reference for further theoretical studies of the asFP595 chromophore in different media and for modeling photoabsorption properties of the asFP595 fluorescent protein.     

Conversion of red fluorescent protein into a bright blue probe

We used a red chromophore formation pathway, in which the anionic red chromophore is formed from the neutral blue intermediate, to suggest a rational design strategy to develop blue fluorescent proteins with a tyrosine-based chromophore. The strategy was applied to red fluorescent proteins of the different genetic backgrounds, such as TagRFP, mCherry, HcRed1, M355NA, and mKeima, which all were converted into blue probes. Further improvement of the blue variant of TagRFP by random mutagenesis resulted in an enhanced monomeric protein, mTagBFP, characterized by the substantially higher brightness, the faster chromophore maturation, and the higher pH stability than blue fluorescent proteins with a histidine in the chromophore. The detailed biochemical and photochemical analysis indicates that mTagBFP is the true monomeric protein tag for multicolor and lifetime imaging, as well as the outstanding donor for green fluorescent proteins in F?rster resonance energy transfer applications.     

Ultrafast photoisomerization of photoactive yellow protein chromophore analogues in solution: influence of the protonation state.

We investigate solvent viscosity and polarity effects on the photoisomerization of the protonated and deprotonated forms of two analogues of the photoactive yellow protein (PYP) chromophore. These are trans-p-hydroxybenzylidene acetone and trans-p-hydroxyphenyl cinnamate, studied in solutions of different polarity and viscosity at room temperature, by means of femtosecond fluorescence up-conversion. The fluorescence lifetimes of the protonated forms are found to be barely sensitive to solvent viscosity, and to increase with increasing solvent polarity. In contrast, the fluorescence decays of the deprotonated forms are significantly slowed down in viscous media and accelerated in polar solvents. These results elucidate the dramatic influence of the protonation state of the PYP chromophore analogues on their photoinduced dynamics. The viscosity and polarity effects are, respectively, interpreted in terms of different isomerization coordinates and charge redistribution in S(1). A trans-to-cis isomerization mechanism involving mainly the ethylenic double-bond torsion and/or solvation is proposed for the anionic forms, whereas "concerted" intramolecular motions are proposed for the neutral forms.     

Trans-cis isomerization is responsible for the red-shifted fluorescence in variants of the red fluorescent protein eqFP611

An important class of red fluorescent proteins (RFPs) feature a 2-iminomethyl-5-(4-hydroxybenzylidene)imidazolinone chromophore. Among these proteins, eqFP611 has the chromophore in a coplanar trans orientation, whereas the cis isomer is preferred by other RFPs such as DsRed and its variants. In the photoactivatable protein asFP595, the chromophore can even be switched from the nonfluorescent trans to the fluorescent cis state by light. By using X-ray crystallography, we have determined the structure of dimeric eqFP611 at high resolution (up to 1.1 A). In the far-red emitting eqFP611 variant d2RFP630, which carries an additional Asn143Ser mutation, the chromophore resides predominantly (approximately 80%) in the cis isomeric state, and in RFP639, which has Asn143Ser and Ser158Cys mutations, the chromophore is found completely in the cis form. The pronounced red shift of excitation and emission maxima of RFP639 can thus unambiguously be assigned to trans-cis isomerization of the chromophore. Among RFPs, eqFP611 is thus unique because its chromophore is highly fluorescent in both the cis and trans isomeric forms.     

Determination of chromophore charge states in the low pH color transition of the fluorescent protein Rtms5 via time-dependent DFT

The red fluorescent protein Rtms5^H146S displays a transition from blue (absorbance λ_max 590 nm) to yellow (absorbance λ_max 453 nm) upon titration to low pH. The pK_a of the reaction depends on the concentration of halide, offering promise for new expressible halide sensors. The protonation state involved in the low pH form of the chromophore remains, however, ambiguous. We report calculated excitation energies of different protonation states of an RFP chromophore model. These suggest that the relevant titration site is the phenoxy moiety of the chromophore, and the relevant base and conjugate acid are anionic and neutral chromophore species, respectively.     

Electronic excitations of the green fluorescent protein chromophore in its protonation states: SAC/SAC-CI study

Two ground-state protonation forms causing different absorption peaks of the green fluorescent protein chromophore were investigated by the quantum mechanical SAC/SAC-CI method with regard to the excitation energy, fluorescence energy, and ground-state stability. The environmental effect was taken into account by a continuum spherical cavity model. The first excited state, HOMO-LUMO excitation, has the largest transition moment and thus is thought to be the source of the absorption. The neutral and anionic forms were assigned to the protonation states for the experimental A- and B-forms, respectively. The present results support the previous experimental observations.     

Structural and Vibrational Characterization of the Chromophore Binding Site of Bacterial Phytochrome Agp1

Agp1 is a prototypical bacterial phytochrome from Agrobacterium fabrum harboring a biliverdin cofactor which reversibly photoconverts between a red‐light‐absorbing (Pr) and a far‐red‐light‐absorbing (Pfr) states. The reaction mechanism involves the isomerization of the bilin‐chromophore followed by large structural changes of the protein matrix that are coupled to protonation dynamics at the chromophore binding site. Histidines His250 and His280 participate in this process. Although the three‐dimensional structure of Agp1 has been solved at high resolution, the precise position of hydrogen atoms and protonation pattern in the chromophore binding pocket has not been investigated yet. Here, we present protonated structure models of Agp1 in the Pr state involving appropriately placed hydrogen atoms that were generated by hybrid quantum mechanics/molecular mechanics‐ and electrostatic calculations and validated against experimental structural‐ and spectroscopic data. Although the effect of histidine protonation on the vibrational spectra is weak, our results favor charge neutral H250 and H280 both protonated at Nε. However, a neutral H250 with a proton at Nε and a cationic H280 may also be possible. Furthermore, the present QM/MM calculations of IR and Raman spectra of Agp1 containing isotope‐labeled BV provide a detailed vibrational assignment of the biliverdin modes in the fingerprint region.     

Protonation of the chromophore in the photoactive yellow protein

The photoactive yellow protein (PYP) acts as a light sensor to its bacterial host: it responds to light by changing shape. After excitation by blue light, PYP undergoes several transformations, to partially unfold into its signaling state. One of the crucial steps in this photocycle is the protonation of p-coumaric acid after excitation and isomerization of this chromophore. Experimentalists still debate on the nature of the proton donor and on whether it donates the hydrogen directly or indirectly. To obtain better knowledge of the mechanism, we studied this proton transfer using Car-Parrinello molecular dynamics, classical molecular dynamics, and computer simulations combining these two methods (quantum mechanics/molecular mechanics, QMMM). The simulations reproduce the chromophore structure and hydrogen-bond network of the protein measured by X-ray crystallography and NMR. When the chromophore is protonated, it leaves the assumed proton donor, glutamic acid 46, with a negative charge in a hydrophobic environment. We show that the stabilization of this charge is a very important factor in the mechanism of protonation. Protonation frequently occurs in simplified ab initio simulations of the chromophore binding pocket in vacuum, where amino acids can easily hydrogen bond to Glu46. When the complete protein environment is incorporated in a QMMM simulation on the complete protein, no proton transfer is observed within 14 ps. The hydrogen-bond rearrangements in this time span are not sufficient to stabilize the new protonation state. Force field molecular dynamics simulations on a much longer time scale have shown which internal rearrangements of the protein are needed. Combining these simulations with more QMMM calculations enabled us to check the stability of protonation states and clarify the initial requirements for the proton transfer in PYP.     

Peptide bond trans-cis isomerization and acylimine formation in chromophore maturation of the red fluorescent proteins

In red fluorescent proteins such as DsRed, an acylimine is formed from the Phe65-Gln66 linkage in GFP-like immature form, while it shows a cis configuration in its mature state. To date, the relationship between acylimine formation and trans-cis isomerization is still unresolved. We have calculated bond rotation profiles for mature and immature chromophores within the protein using our own n-layered integrated molecular orbital and molecular mechanism (ONIOM) approach. The results suggested that the isomerization is barrierless in acylimine formed in the mature state, suggesting that the acylimine formation precedes the trans-cis isomerization in DsRed chromophores. Further decomposition analysis of electrostatic contributions from individual residues has identified several residues and a specific water molecule which could play key roles in controlling the rate of the trans-cis isomerization of peptide bond in immature GFP-like protein. The results also highlight the importance of Gln66-like of tripeptide motif (chromophore) in the maturation of red fluorescent proteins. In view of the considerable interest in developing red fluorescent proteins for numerous biotechnological applications, these results should be useful for design of novel fluorescent proteins.     

Isomerization mechanism of the HcRed fluorescent protein chromophore

To understand how the protein achieves fluorescence, the isomerization mechanism of the HcRed chromophore is studied both under vacuum and in the solvated red fluorescent protein. Quantum mechanical (QM) and quantum mechanical/molecular mechanical (QM/MM) methods are applied both for the ground and the first excited state. The photoinduced processes in the chromophore mainly involve torsions around the imidazolinone-bridge bond (τ) and the phenoxy-bridge bond (φ). Under vacuum, the isomerization of the cis-trans chromophore essentially proceeds by τ twisting, while the radiationless decay requires φ torsion. By contrast, the isomerization of the cis-trans chromophore in HcRed occurs via simultaneous τ and φ twisting. The protein environment significantly reduces the barrier of this hula twist motion compared with vacuum. The excited-state isomerization barrier via the φ rotation of the cis-coplanar conformer in HcRed is computed to be significantly higher than that of the trans-non-coplanar conformer. This is consistent with the experimental observation that the cis-coplanar-conformation of the chromophore is related to the fluorescent properties of HcRed, while the trans-non-planar conformation is weakly fluorescent or non-fluorescent. Our study shows how the protein modifies the isomerization mechanism, notably by interactions involving the nearby residue Ile197, which keeps the chromophore coplanar and blocks the twisting motion that leads to photoinduced radiationless decay.     

Photoconversion in the red fluorescent protein from the sea anemone Entacmaea quadricolor: is cis-trans isomerization involved?

Proteins from the family of the green fluorescent protein (GFP) are presently extensively used in molecular and cellular biology. Recent studies suggest that isomerization of the chromophore occurs upon excitation and is involved in nonradiative deactivation. Using Raman spectroscopy, we report on photoinduced cis-trans isomerization in the red fluorescent protein eqFP611 from the sea anemone Entacmaea quadricolor. The crystal structure of eqFP611 shows that the chemical structure of the chromophore, p-hydroxybenzylidene-imidazolinone with an extended -conjugated system, is nearly identical to the chromophore of other red fluorescent proteins such as DsRed and HcRed. However, the chromophore of eqFP611 has a trans configuration whereas the chromophore of DsRed has a cis configuration. Upon irradiation with 532-nm light, the absorption of eqFP611 peaking at 559 nm diminished, and concomitantly a drastic decrease in the quantum yield of fluorescence as well as more complex decay kinetics was observed. Upon irradiation, changes in the Raman spectrum of eqFP611 were observed, and the relative intensities and peak positions of the irradiated eqFP611 showed striking similarity with the peaks in the Raman spectrum of DsRed. These observations are tentatively interpreted as trans-to-cis isomerization of the chromophore taking place upon irradiation together with the opening of new, nonradiative pathways.     

Structural basis of fluorescence fluctuation dynamics of green fluorescent proteins in acidic environments

Green fluorescent proteins (GFPs) have become powerful markers for numerous biological studies due to their robust fluorescence properties, site-specific labeling, pH sensitivity, and mutations for multiple-site labeling. Fluorescence correlation spectroscopy (FCS) studies have indicated that fluorescence blinking of anionic GFP mutants takes place on a time scale of 45-300 ms, depending on pH, and have been attributed to external proton transfer. Here we present experimental evidence indicating that conformational change in the protein &beta-barrel is a determining step for the external protonation of GFP-S65T (at low pH) using time-resolved fluorescence and polarization anisotropy measurements. While the average anionic fluorescence lifetime of GFP-S65T is reduced by approximately 18% over a pH range of 3.6-10.0, the fluorescence polarization anisotropy decays mostly as a single exponential with a rotational time of phi = 17 +/- 1 ns, which indicates an intact beta-barrel with a hydrodynamic volume of 78 +/- 5 nm3. In contrast, the total fluorescence (525 +/- 50 nm) of the excited neutral state of S65T reveals a strong correlation between the fluorescence lifetime, structural conformation, and pH. The average fluorescence lifetime of the excited neutral state of S65T as a function of pH yields pKa approximately 5.9 in agreement with literature values using steady-state techniques. In contrast to the intact beta-barrel at high pH, the anisotropy of neutral S65T (at pH 相似文献   

The nature of transient dark states in a photoactivatable fluorescent protein

Fluorescent proteins (FPs) of the green fluorescent protein family blink and bleach like all fluorophores. However, contrary to organic dyes, the mechanisms by which transient losses of fluorescence occur in FPs have received little attention. Here, we focus on the photoactivatable IrisFP, for which a transient non-fluorescent chromophoric state with distorted geometry was recently reported (Adam, V.; et al. J. Am. Chem. Soc. 009, 131, 18063). We investigated the chemical nature of this blinked state by employing quantum chemical/molecular mechanical calculations. Our findings suggest two previously unidentified dark states that display similar distorted chromophores with a transiently ruptured π-electron system. Both are protonated at atom C(α) of the chromophore methylene bridge. Transient protonation may occur via proton transfer from the nearby Arg66 either in the triplet state T(1) after intersystem crossing or in an anionic radical (doublet) ground state. As Arg66 is conserved in green-to-red photoconvertible FPs, these dark states are predicted to be common to all these proteins. We also suggest that C(α) protonated dark states may accelerate photobleaching by favoring decarboxylation of the fully conserved Glu212.     

Kinetics of proton transfer in a green fluorescent protein: A laser-induced pH jump study

The sub-millisecond protonation dynamics of the chromophore in S65T mutant form of the green fluorescent protein (GFP) was tracked after a rapid pH jump following laser-induced proton release from the caged photolabile compoundo-nitrobenzaldehyde. Following a jump in pH from 8 to 5 (which is achieved within 2 μs), the fluorescence of S65T GFP decreased as a single exponential with a time constant of ∼90 μs. This decay is interpreted as the conversion of the deprotonated fluorescent GFP chromophore to a protonated non-fluorescent species. The protonation kinetics showed dependence on the bulk viscosity of the solvent, and therefore implicates bulk solvent-controlled protein dynamics in the protonation process. The protonation is proposed to be a sequential process involving two steps: (a) proton transfer from solvent to the chromophore, and (b) internal structural rearrangements to stabilize a protonated chromophore. The possible implications of these observations to protein dynamics in general is discussed