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

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Combined quantum-mechanical molecular mechanics calculations with NWChem and AMBER: Excited state properties of green fluorescent protein chromophore analogue in aqueous solution

Combined quantum mechanical molecular mechanics (QM/MM) calculations have become a popular methodology for efficient and accurate description of large molecular systems. In this work we introduce our development of a QM/MM framework based on two well-known codes-NWChem and AMBER. As an initial application area we are focused on excited state properties of small molecules in an aqueous phase using an analogue of the green fluorescent protein (GFP) chromophore as a particular test case. Our approach incorporates high level coupled cluster theory for the analysis of excited states providing a reliable theoretical analysis of effects of an aqueous solvation environment on the photochemical properties of the GFP chromophore. Using a systematic approach, which involves comparison of gas phase and aqueous phase results for different protonation states and conformations, we resolve existing uncertainties regarding the theoretical interpretation of experimental data. We observe that the impact of aqueous environment on charged states generally results in blue shifts of the absorption spectra, but the magnitude of the effect is sensitive to both protonation state and conformation and can be rationalized based on charge movement into the area of higher/lower external electrostatic potentials. At neutral pH levels the experimentally observed absorption signal is most likely coming from the phenol protonated form. Our results also show that the high level electron correlated method is essential for a proper description of excited states of GFP. © 2017 Wiley Periodicals, Inc.

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

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

Probing microhydration effect on the electronic structure of the GFP chromophore anion: Photoelectron spectroscopy and theoretical investigations

The photophysics of the Green Fluorescent Protein (GFP) chromophore is critically dependent on its local structure and on its environment. Despite extensive experimental and computational studies, there remain many open questions regarding the key fundamental variables that govern this process. One outstanding problem is the role of autoionization as a possible relaxation pathway of the excited state under different environmental conditions. This issue is considered in our work through combined experimental and theoretical studies of microsolvated clusters of the deprotonated p-hydroxybenzylidene-2,3-dimethylimidazolinone anion (HBDI(-)), an analog of the GFP chromophore. Through selective generation of microsolvated structures of predetermined size and subsequent analysis of experimental photoelectron spectra by high level ab initio methods, we are able to precisely identify the structure of the system, establish the accuracy of theoretical data, and provide reliable description of auto-ionization process as a function of hydrogen-bonding environment. Our study clearly illustrates the first few water molecules progressively stabilize the excited state of the chromophore anion against the autodetached neutral state, which should be an important trait for crystallographic water molecules in GFPs that has not been fully explored to date.

On the optical absorption of the anionic GFP chromophore in vacuum, solution, and protein

In spite of the large number of experimental and theoretical studies, the optical absorption trend of the green fluorescent protein (GFP) chromophore in several environments has not been fully understood. We calculated at the same level of time dependent density functional theory the vertical excitation energy of the anionic GFP chromophore in the protein and in ethanol, dioxane, methanol and water solutions. As result, we reproduced for the first time the experimental trend of the absorption peaks with 0.015 eV as the standard deviation of the accuracy. This systematic error allowed us to analyze with confidence the relative weight of several solvation effects on the vertical excitation energy. Experimental trends not correlated with the solvent polarity were therefore explained with a fine balance of different steric and electronic effects on the photophysics of the chromophore. As an indirect and remarkable result, the present analysis confirms that the optical absorption of the chromophore in the gas phase is close to the value of 2.84 eV extrapolated by Dong et al. (J. Am. Chem. Soc., 2006, 128, 12038), and, as a consequence, that the protein environment induces a red shift of 0.23 eV.

Effect of deuteration on hydrogen bonding: a comparative concentration dependent Raman and DFT study of pyridine in CH₃OH and CD₃OD and pyrimidine in H₂O and D₂O

The relative effect of hydrogen bonding of pyrimidine (Pyr) in H(2)O/D(2)O and pyridine (Py) in CH(3)OH/CD(3)OD has been analyzed using Raman Difference Spectroscopic (RDS) technique and DFT calculations. This study is focused on analyzing the concentration dependent variation of linewidth, peak position and intensity of ring breathing mode of Py and Pyr. The ring breathing mode of Pyr in H(2)O and D(2)O has three components; due to free Pyr, lighter complexes of mPyr+nH(2)O/D(2)O and heavier complexes of mPyr+nH(2)O/D(2)O. The pyridine molecules, however, show only two components in CH(3)OH and CD(3)OD. Of these two components, one corresponds to free Py and the other inhomogeneously broadened profile corresponds to all mPy+nCH(3)OH/CD(3)OD complexes. The variation of peak position and linewidth establishes the role of dipole moment of complexes and the diffusion in the mixture. In case of CD(3)OD solution splitting was observed in ∼1030 cm(-1) band of Py, where an additional band at ∼1034 cm(-1) appears at x(Py) ≤ 0.4. However, this band remains single at all concentrations in case of CH(3)OH solvent.

Intramolecular hydrogen bonding plays a crucial role in the photophysics and photochemistry of the GFP chromophore

In commonly studied GFP chromophore analogues such as 4-(4-hydroxybenzylidene)-1,2-dimethyl-1H-imidazol-5(4H)-one (PHBDI), the dominant photoinduced processes are cis-trans isomerization and subsequent S(1) → S(0) decay via a conical intersection characterized by a highly twisted double bond. The recently synthesized 2-hydroxy-substituted isomer (OHBDI) shows an entirely different photochemical behavior experimentally, since it mainly undergoes ultrafast intramolecular excited-state proton transfer, followed by S(1) → S(0) decay and ground-state reverse hydrogen transfer. We have chosen 4-(2-hydroxybenzylidene)-1H-imidazol-5(4H)-one (OHBI) to model the gas-phase photodynamics of such 2-hydroxy-substituted chromophores. We first use various electronic structure methods (DFT, TDDFT, CC2, DFT/MRCI, OM2/MRCI) to explore the S(0) and S(1) potential energy surfaces of OHBI and to locate the relevant minima, transition state, and minimum-energy conical intersection. These static calculations suggest the following decay mechanism: upon photoexcitation to the S(1) state, an ultrafast adiabatic charge-transfer induced excited-state intramolecular proton transfer (ESIPT) occurs that leads to the S(1) minimum-energy structure. Nearby, there is a S(1)/S(0) minimum-energy conical intersection that allows for an efficient nonadiabatic S(1) → S(0) internal conversion, which is followed by a fast ground-state reverse hydrogen transfer (GSHT). This mechanism is verified by semiempirical OM2/MRCI surface-hopping dynamics simulations, in which the successive ESIPT-GSTH processes are observed, but without cis-trans isomerization (which is a minor path experimentally with less than 5% yield). These gas-phase simulations of OHBI give an estimated first-order decay time of 476 fs for the S(1) state, which is larger but of the same order as the experimental values measured for OHBDI in solution: 270 fs in CH(3)CN and 230 fs in CH(2)Cl(2). The differences between the photoinduced processes of the 2- and 4-hydroxy-substituted chromophores are attributed to the presence or absence of intramolecular hydrogen bonding between the two rings.

ELECTRONIC EXCITATIONS OF GREEN FLUORESCENT PROTEINS: MODELING SOLVATOCHROMATIC SHIFTS OF CHROMOPHORE MODEL COMPOUNDS IN SOLUTIONS

Green fluorescent protein (GFP) is a spontaneously fluorescent protein due to its p-hydroxylbenzylideneimidazolidinone chromophore. In this work, we have investigated the electronic structures, liquid structures, and solvent shifts of the GFP chromophore model compounds 4′-hydroxylbenzylidene-2,3-dimethylimidazolin-5-one (HBDI) and 4′-hydroxylbenzylidene-2-methyl-imidazoin-5-one-3-acetate (HBMIA) in NaOH /aqueous solutions, in which both compounds are protonated at anionic state. The electronic structure calculations predict that both model compounds could adopt both cis and trans conformations in solutions. Moreover, liquid simulations elucidate an extensive well-defined hydrogen-bonding network between the solvent and the solute in the ground state. Furthermore, solvent shifts calculations indicate that contributions from the specific solute-solvent hydrogen-bonding interactions are negligible for the solvatochromatic shifts observed in the absorption spectrum of the model compounds in solutions; rather, the solvent shifts are dominated by the dipolar solvation in which both permanent charge–charge interactions and many-body polarizations contribute significantly. Self-Consistent Reaction-Field (SCRF) approach could be the efficient method for studying the unusual optical properties of the GFP chromophore in solutions and proteins.

PROTONATION STATES AND CONFORMATIONAL FLEXIBILITY OF THE RED FLUORESCENT PROTEIN CHROMOPHORE

Fluorescent proteins are important reporter molecules widely used in biotechnology and the biological sciences, generally. Their unusual spectrophotometric and fluorescence characteristics are controlled via protonation states of the chromophore, as for green fluorescent proteins (GFPs), and/or via cis–trans isomerization of the chromophore, as for red fluorescent proteins (RFPs). Here, we have performed both quantum mechanical calculations on substituted chromophores and structural comparisons of several RFPs (Rtms5, eqFP611, HcRed, and DsRed) and wild-type GFP. Our results indicate that the chromophore cis and trans isomers are comparably stable, and cis–trans isomerization has only a minor effect on electronic excitation. The chromophore is also found to exist in either the anionic or possibly zwitterionic protonation state. Structural comparisons of RFPs reveal that the conformation of the chromophore within a specific RFP is determined by only a few key residues, which could serve as mutation targets for engineering new fluorescent proteins for novel biotechnological applications.

Potential Energy Landscape of the Electronic States of the GFP Chromophore in Different Protonation Forms: Electronic Transition Energies and Conical Intersections

We present the results of quantum chemical calculations of the transition energies and conical intersection points for the two lowest singlet electronic states of the green fluorescent protein chromophore, 4'-hydroxybenzylidene-2,3-dimethylimidazolinone, in the vicinity of its cis conformation in the gas phase. Four protonation states of the chromophore, i.e., anionic, neutral, cationic, and zwitterionic, were considered. Energy differences were computed by the perturbatively corrected complete active space self-consistent field (CASSCF)-based approaches at the corresponding potential energy minima optimized by density functional theory and CASSCF (for the ground and excited states, respectively). We also report the EOM-CCSD and SOS-CIS(D) results for the excitation energies. The minimum energy S0/S1 conical intersection points were located using analytic state-specific CASSCF gradients. The results reproduce essential features of previous ab initio calculations of the anionic form of the chromophore and provide an extension for the neutral, cationic, and zwitterionic forms, which are important in the protein environment. The S1 PES of the anion is fairly flat, and the barrier separating the planar bright conformation from the dark twisted one as well as the conical intersection point with the S0 surface is very small (less than 2 kcal/mol). On the cationic surface, the barrier is considerably higher (∼13 kcal/mol). The PES of the S1 state of the zwitterionic form does not have a planar minimum in the Franck-Condon region. The S1 surface of the neutral form possesses a bright planar minimum; the energy barrier of about 9 kcal/mol separates it from the dark twisted conformation as well as from the conical intersection point leading to the cis-trans chromophore isomerization.

The 559-to-600 nm shift observed in red fluorescent protein eqFP611 is attributed to cis-trans isomerization of the chromophore in an anionic protein pocket

Fluorescent proteins are commonly used as molecular labels, noninvasive markers of gene expression, and reporters of environmental conditions in live cells. We investigate the structural and spectroscopic properties of the chromophore of a far-red fluorescent protein eqFP611. Both the cis and trans isomers of the chromophore are examined within the protein for which both anionic and neutral states of protonation are considered. Spectroscopic properties are examined using time-dependent density functional theory (TDDFT), employing the B3LYP, PBE and B3PW91 density functionals. Intermolecular and long-range contributions to the structure and spectroscopy were treated using the own n-layered integrated molecular orbital and molecular mechanics (ONIOM) approach. The results indicated that the chromophore before excitation is in an anionic, protonated state, with the long-range contributions inducing a blue shift in the absorption and fluorescence maxima of the chromophore. Moreover, the calculated changes of the lowest pi-pi* excitation energy upon isomerization match the observed shift from 559 to 600 nm in the absorption maximum of the system following prolonged irradiation. Furthermore, decomposition analysis of the electrostatic contributions from individual residues indicated that the interactions from four residues Arg92, Lys67, Glu145, and His197 to the chromophore play a key role in the absorption and fluorescence spectra of eqFP611, suggesting that mutations at these sites should provide very useful mechanistic information.

Time-dependent density functional theory study on hydrogen-bonded intramolecular charge-transfer excited state of 4-dimethylamino-benzonitrile in methanol

The time-dependent density functional theory (TDDFT) method was carried out to investigate the hydrogen-bonded intramolecular charge-transfer (ICT) excited state of 4-dimethylaminobenzonitrile (DMABN) in methanol (MeOH) solvent. We demonstrated that the intermolecular hydrogen bond C[triple bond]N...H-O formed between DMABN and MeOH can induce the C[triple bond]N stretching mode shift to the blue in both the ground state and the twisted intramolecular charge-transfer (TICT) state of DMABN. Therefore, the two components at 2091 and 2109 cm(-1) observed in the time-resolved infrared (TRIR) absorption spectra of DMABN in MeOH solvent were reassigned in this work. The hydrogen-bonded TICT state should correspond to the blue-side component at 2109 cm(-1), whereas not the red-side component at 2091 cm(-1) designated in the previous study. It was also demonstrated that the intermolecular hydrogen bond C[triple bond]N...H-O is significantly strengthened in the TICT state. The intermolecular hydrogen bond strengthening in the TICT state can facilitate the deactivation of the excited state via internal conversion (IC), and thus account for the fluorescence quenching of DMABN in protic solvents. Furthermore, the dynamic equilibrium of these electronically excited states is explained by the hydrogen bond strengthening in the TICT state.

Protonation and conformational dynamics of GFP mutants by two-photon excitation fluorescence correlation spectroscopy

GFP mutants are known to display fluorescence flickering, a process that occurs in a wide time range. Because serine 65, threonine 203, glutamate 222, and histidine 148 have been indicated as key residues in determining the GFP fluorescence photodynamics, we have focused here on the role of histidine 148 and glutamate 222 by studying the fluorescence dynamics of GFPmut2 (S65A, V68L, and S72A GFP) and its H148G (Mut2G) and E222Q (Mut2Q) mutants. Two relaxation components are found in the fluorescence autocorrelation functions of GFPmut2: a 10-100 micros pH-dependent component and a 100-500 micros laser-power-dependent component. The comparison of these three mutants shows that the mutation of histidine 148 to glycine induces a 3-fold increase in the protonation rate, thereby indicating that the protonation-deprotonation of the chromophore occurs via a proton exchange with the solution mediated by the histidine 148 residue. The power-dependent but pH-independent relaxation mode, which is not affected by the E222Q and H148G mutations, is due to an excited-state process that is probably related to conformational rearrangements of the chromophore after the photoexcitation, more than to the chromophore excited-state proton transfer.

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.

Novel infrared spectra for intermolecular dihydrogen bonding of the phenol-borane-trimethylamine complex in electronically excited state

The intermolecular dihydrogen bonding in the electronically excited states of the dihydrogen-bonded phenol-BTMA complex in gas phase was theoretically investigated using the time-dependent density functional theory method for the first time. It was theoretically demonstrated that the S(1) state of the dihydrogen-bonded phenol-BTMA complex is a locally excited state, in which only the phenol moiety is electronically excited. The infrared spectra of the dihydrogen-bonded phenol-BTMA complex in ground state and the S(1) state were calculated at both the O-H and B-H stretching vibrational regions. A novel infrared spectrum of the dihydrogen-bonded phenol-BTMA complex in the electronically excited state was found. The stretching vibrational absorption bands of the dihydrogen-bonded O-H and B-H groups are very strong in the ground state, while they are disappeared in the S(1) state. At the same time, a new strong absorption band appears at the C[Double Bond]O stretching region. From the calculated bond lengths, it was found that both the O-H and B-H bonds in the dihydrogen bond O-H...H-B are significantly lengthened in the S(1) state of the dihydrogen-bonded phenol-BTMA complex. However, the C-O bond in the phenol moiety is markedly shortened in the excited state, and then has the characteristics of C[Double Bond]O group. Furthermore, it was demonstrated that the intermolecular dihydrogen bonds in the electronically excited state of the dihydrogen-bonded phenol-BTMA complex are strengthened, since calculated H...H distance is drastically shortened in the S(1) state.

Ultrafast Hydrogen Bond Strengthening of the Photoexcited Fluorenone in Alcohols for Facilitating the Fluorescence Quenching

The time-dependent density functional theory (TDDFT) method was performed to investigate the excited-state hydrogen-bonding dynamics of fluorenone (FN) in hydrogen donating methanol (MeOH) solvent. The infrared spectra of the hydrogen-bonded FN-MeOH complex in both the ground state and the electronically excited states are calculated using the TDDFT method, since the ultrafast hydrogen-bonding dynamics can be investigated by monitoring the vibrational absorption spectra of some hydrogen-bonded groups in different electronic states. We demonstrated that the intermolecular hydrogen bond C=O...H-O between fluorenone and methanol molecules is significantly strengthened in the electronically excited-state upon photoexcitation of the hydrogen-bonded FM-MeOH complex. The hydrogen bond strengthening in electronically excited states can be used to explain well all the spectral features of fluorenone chromophore in alcoholic solvents. Furthermore, the radiationless deactivation via internal conversion (IC) can be facilitated by the hydrogen bond strengthening in the excited state. At the same time, quantum yields of the excited-state deactivation via fluorescence are correspondingly decreased. Therefore, the total fluorescence of fluorenone in polar protic solvents can be drastically quenched by hydrogen bonding.

Ground-State Structures and Vertical Excitations for the Kindling Fluorescent Protein asFP595

Geometry configurations of a large fraction of the kindling fluorescent protein asFP595 around the chromophore region were optimized by using the effective fragment potential quantum mechanical-molecular mechanical (QM/MM) method. The initial coordinates of heavy atoms were taken from the structure from the Protein Data Bank archive corresponding to the dark-adapted state of the Ala143 --> Gly mutant of asFP595. Optimization of geometry parameters was performed for all internal coordinates in the QM part composed of the chromophore unit and the side chains of His197, Glu215, and Arg92 as well as for positions of effective fragments constituting the MMpart. The structures corresponding to the anion trans, anion cis, and zwitterion trans moieties were considered among various alternatives for the chromophore unit inside the protein matrix. The QM/MM simulations show that the protein environment provides stabilization for the trans-zwitterion isomer compared to the gas-phase conditions. By using the multiconfigurational CASSCF and the time-dependent density functional theory calculations, we estimated positions of spectral bands corresponding to vertical S(0)-S(1) transitions. The results of simulations support the assumption that the dark state of asFP595 corresponds to the anionic or zwitterionic trans-conformation, while the kindled state corresponds to the anionic cis-conformation.