Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.     

Single molecule photon emission statistics in the slow modulation limit

A framework for calculating photon emission statistics for single chromophores perturbed by slow environmental fluctuations is introduced. When internal chromophore dynamics are significantly faster than time scales for environmental modulation it becomes possible to invoke a type of adiabatic approximation, allowing for straightforward calculation of photon counting moments including explicitly quantum effects. Unlike previous exact treatments, the present methodology involves calculation of dynamics reflecting only the modulation characteristics of the environment and quantum dynamics of an isolated chromophore separately, i.e., the complicated intermingling of chromophore quantum dynamics and the environmental modulation are suppressed via the adiabatic approximation. This leads to significant conceptual and computational simplifications. Within its regime of applicability, the present approximation reproduces exact calculations quantitatively. We demonstrate this accuracy explicitly for the case of a two-level chromophore modulated by a number of different stochastic models.     

Existence of a new emitting singlet state of proflavine: femtosecond dynamics of the excited state processes and quantum chemical studies in different solvents

Proflavine (3,6-diaminoacridine) shows fluorescence emission with lifetime, 4.6 ± 0.2 ns, in all the solvents irrespective of the solvent polarity. To understand this unusual photophysical property, investigations were carried out using steady state and time-resolved fluorescence spectroscopy in the pico- and femtosecond time domain. Molecular geometries in the ground and low-lying excited states of proflavine were examined by complete structural optimization using ab initio quantum chemical computations at HF/6-311++G** and CIS/6-311++G** levels. Time dependent density functional theory (TDDFT) calculations were performed to study the excitation energies in the low-lying excited states. The steady state absorption and emission spectral details of proflavine are found to be influenced by solvents. The femtosecond fluorescence decay of the proflavine in all the solvents follows triexponential function with two ultrafast decay components (τ(1) and τ(2)) in addition to the nanosecond component. The ultrafast decay component, τ(1), is attributed to the solvation dynamics of the particular solvent used. The second ultrafast decay component, τ(2), is found to vary from 50 to 215 ps depending upon the solvent. The amplitudes of the ultrafast decay components vary with the wavelength and show time dependent spectral shift in the emission maximum. The observation is interpreted that the time dependent spectral shift is not only due to solvation dynamics but also due to the existence of more than one emitting state of proflavine in the solvent used. Time resolved area normalized emission spectral (TRANES) analysis shows an isoemissive point, indicating the presence of two emitting states in homogeneous solution. Detailed femtosecond fluorescence decay analysis allows us to isolate the two independent emitting components of the close lying singlet states. The CIS and TDDFT calculations also support the existence of the close lying emitting states. The near constant lifetime observed for proflavine in different solvents is suggested to be due to the similar dipole moments of the ground and the evolved emitting singlet state of the dye from the Franck-Condon excited state.     

Ultrafast relaxation dynamics of a biologically relevant probe dansyl at the micellar surface

We report picosecond-resolved measurement of the fluorescence of a well-known biologically relevant probe, dansyl chromophore at the surface of a cationic micelle (cetyltrimethylammonium bromide, CTAB). The dansyl chromophore has environmentally sensitive fluorescence quantum yields and emission maxima, along with large Stokes shift. In order to study the solvation dynamics of the micellar environment, we measured the fluorescence of dansyl chromophore attached to the micellar surface. The fluorescence transients were observed to decay (with time constant approximately 350 ps) in the blue end and rise with similar timescale in the red end, indicative of solvation dynamics of the environment. The solvation correlation function is measured to decay with time constant 338 ps, which is much slower than that of ordinary bulk water. Time-resolved anisotropy of the dansyl chromophore shows a bi-exponential decay with time constants 413 ps (23%) and 1.3 ns (77%), which is considerably slower than that in free solvents revealing the rigidity of the dansyl-micelle complex. Time-resolved area-normalized emission spectroscopic (TRANES) analysis of the time dependent emission spectra of the dansyl chromophore in the micellar environment shows an isoemissive point at 21066 cm-1. This indicates the fluorescence of the chromophore contains emission from two kinds of excited states namely locally excited state (prior to charge transfer) and charge transfer state. The nature of the solvation dynamics in the micellar environments is therefore explored from the time-resolved anisotropy measurement coupled with the TRANES analysis of the fluorescence transients. The time scale of the solvation is important for the mechanism of molecular recognition.     

Solvation dynamics in reverse micelles: the role of headgroup-solute interactions

We present molecular dynamics simulation results for solvation dynamics in the water pool of anionic-surfactant reverse micelles (RMs) of varying water content, w(0). The model RMs are designed to represent water/aerosol-OT/oil systems, where aerosol-OT is the common name for sodium bis(2-ethylhexyl)sulfosuccinate. To determine the effects of chromophore-headgroup interactions on solvation dynamics, we compare the results for charge localization in model ionic diatomic chromophores that differ only in charge sign. Electronic excitation in both cases is modeled as charge localization on one of the solute sites. We find dramatic differences in the solvation responses for anionic and cationic chromophores. Solvation dynamics for the cationic chromophore are considerably slower and more strongly w(0)-dependent than those for the anionic chromophore. Further analysis indicates that the difference in the responses can be ascribed in part to the different initial locations of the two chromophores relative to the surfactant interface. In addition, slow motion of the cationic chromophore relative to the interface is the main contributor to the longer-time decay of the solvation response to charge localization in this case.     

Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform

We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.     

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.     

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.     

Early time hydrogen-bonding dynamics of photoexcited coumarin 102 in hydrogen-donating solvents: theoretical study

To study the early time hydrogen-bonding dynamics of chromophore in hydrogen-donating solvents upon photoexcitation, the infrared spectra of the hydrogen-bonded solute-solvent complexes in electronically excited states have been calculated using the time-dependent density functional theory (TDDFT) method. The hydrogen-bonding dynamics in electronically excited states can be widely monitored by the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds. In this study, we have demonstrated that the intermolecular hydrogen bonds between coumarin 102 (C102) and hydrogen-donating solvents are strengthened in the early time of photoexcitation to the electronically excited state by theoretically monitoring the stretching modes of C=O and H-O groups. This is significantly contrasted with the ultrafast hydrogen bond cleavage taking place within a 200-fs time scale upon electronic excitation, proposed in many femtosecond time-resolved vibrational spectroscopy experiments. The transient hydrogen bond strengthening behaviors in excited states of chromophores in hydrogen-donating solvents, which we have demonstrated here for the first time, may take place widely in many other systems in solution and are very important to explain the fluorescence-quenching phenomena associated with some radiationless deactivation processes, for example, the ultrafast solute-solvent intermolecular electron transfer and the internal conversion process from the fluorescent state to the ground state.     

Site-selective photoinduced electron transfer from alcoholic solvents to the chromophore facilitated by hydrogen bonding: a new fluorescence quenching mechanism

Solute-solvent intermolecular photoinduced electron transfer (ET) reaction was proposed to account for the drastic fluorescence quenching behaviors of oxazine 750 (OX750) chromophore in protic alcoholic solvents. According to our theoretical calculations for the hydrogen-bonded OX750-(alcohol)(n) complexes using the time-dependent density functional theory (TDDFT) method, we demonstrated that the ET reaction takes place from the alcoholic solvents to the chromophore and the intermolecular ET passing through the site-specific intermolecular hydrogen bonds exhibits an unambiguous site selectivity. In our motivated experiments of femtosecond time-resolved stimulated emission pumping fluorescence depletion spectroscopy (FS TR SEP FD), it could be noted that the ultrafast ET reaction takes place as fast as 200 fs. This ultrafast intermolecular photoinduced ET is much faster than the diffusive solvation process, and even significantly faster than the intramolecular vibrational redistribution (IVR) process of the OX750 chromophore. Therefore, the ultrafast intermolecular ET should be coupled with the hydrogen-bonding dynamics occurring in the sub-picosecond time domain. We theoretically demonstrated for the first time that the selected hydrogen bonds are transiently strengthened in the excited states for facilitating the ultrafast solute-solvent intermolecular ET reaction.     

Exploration of the secondary structure specific differential solvation dynamics between the native and molten globule states of the protein HP-36

Recent experiments have shown that the time dependence of fluorescence Stokes shift of a chromophore is substantially different when the chromophore is located in a molten globule (MG) state and in the native state of the same protein. To understand the origin of this difference, particularly the role of water in the differential solvation of the protein in the native and the MG states, we have carried out fully atomistic molecular dynamics simulations with explicit water of a partially unfolded MG state of the protein HP-36 and compared the results with the solvation dynamics of the protein in the folded native state. It is observed that the polar solvation dynamics of the three helical segments of the protein is influenced in a nonuniform heterogeneous manner in the MG state. While the equilibrium solvation time correlation function for helix-3 has been found to relax faster in the MG state as compared to that in the native state, the decay of the corresponding function for the other two helices slows down in the MG state. A careful analysis shows that the origin of such heterogeneous relative solvation behavior lies in the differential location of the polar probe residues and their exposure to bulk solvent. We find a significant negative cross-correlation between the contribution (to the solvation energy of a tagged amino acid residue) of water and the other groups of the protein, indicating a competing role in solvation. The sensitivity of solvation dynamics to the secondary structure and the immediate environment can be used to discriminate the partially unfolded and folded states. These results therefore should be useful in explaining recent solvation dynamics experiments on native and MG states of proteins.     

Theoretical Computer‐Aided Mutagenic Study on the Triple Green Fluorescent Protein Mutant S65T/H148D/Y145F

Green fluorescent protein (GFP) mutant S65T/H148D has been proposed to host a photocycle that involves an excited‐state proton transfer between the chromophore (Cro) and the Asp148 residue and takes place in less than 50 fs without a measurable kinetic isotope effect. It has been suggested that the interaction between the unsuspected Tyr145 residue and the chromophore is needed for the ultrafast sub‐50 fs rise in fluorescence. To verify this, we have performed a computer‐aided mutagenic study to introduce the additional mutation Y145F, which eliminates this interaction. By means of QM/MM molecular dynamics simulations and time‐dependent density functional theory studies, we have assessed the importance of the Cro–Tyr145 interaction and the solvation of Asp148 and shown that in the triple mutant S65T/H148D/Y145F a significant loss in the ultrafast rise of the Stokes‐shifted fluorescence should be expected.     

Simulating the formation of sodium:electron tight-contact pairs: watching the solvation of atoms in liquids one molecule at a time

The motions of solvent molecules during a chemical transformation often dictate both the dynamics and the outcome of solution-phase reactions. However, a microscopic picture of solvation dynamics is often obscured by the concerted motions of numerous solvent molecules that make up a condensed-phase environment. In this study, we use mixed quantum/classical molecular dynamics simulations to furnish the molecular details of the solvation dynamics that leads to the formation of a sodium cation-solvated electron contact pair, (Na(+), e(-)), in liquid tetrahydrofuran following electron photodetachment from sodide (Na(-)). Our simulations reveal that the dominant solvent response is comprised of a series of discrete solvent molecular events that work sequentially to build up a shell of coordinating THF oxygen sites around the sodium cation end of the contact pair. With the solvent response described in terms of the sequential motion of single molecules, we are then able to compare the calculated transient absorption spectroscopy of the sodium species to experiment, providing a clear microscopic interpretation of ultrafast pump-probe experiments on this system. Our findings suggest that for solute-solvent interactions similar to the ones present in our study, the solvation dynamics is best understood as a series of kinetic events consisting of reactions between chemically distinct local structures in which key solvent molecules must be considered to be part of the identity of the reacting species.     

Cis-trans photoisomerization of fluorescent-protein chromophores

Photochromic variants of fluorescent proteins are opening the way to a number of opportunities for high-sensitivity regioselective studies in the cellular environment and may even lead to applications in information and communication technology. Yet, the detailed photophysical processes at the basis of photoswitching have not been fully clarified. In this paper, we used synthetic FP chromophores to clarify the photophysical processes associated with the photochromic behavior. In particular, we investigated the spectral modification of synthetic chromophore analogues of wild-type green fluorescent protein (GFP), Y66F GFP (BFPF), and Y66W GFP (CFP) upon irradiation in solutions of different polarities. We found that the cis-trans photoisomerization mechanism can be induced in all the chromophores. The structural assignments were carried out both by NMR measurements and DFT calculations. Remarkably, we determined for the first time the spectra of neutral trans isomers in different solvents. Finally, we calculated the photoconversion quantum yields by absorption measurements under continuous illumination at different times and by a nanosecond laser-flash photolysis method. Our results indicate that cis-trans photoisomerization is a general mechanism of FP chromophores whose efficiency is modulated by the detailed mutant-specific protein environment.     

Photochemistry of visual pigment chromophore models by ab initio molecular dynamics

Ab initio excited-state molecular dynamics calculations have been performed to study the effect of methyl substitution and chromophore distortion on the photoreaction of different four-double-bond retinal model chromophores. Randomly distributed starting geometries were generated by zero-point energy sampling; after Franck-Condon excitation the reaction was followed on the S1 surface. For determining the photoproduct and its configuration, a simplified approach--torsion angle following--is discussed and applied. We find that chromophore distortion significantly affects the outcome of the photoreaction: with dihedral angles taken from the rhodopsin-embedded 11-cis-retinal chromophore, the reaction rate of the model chromophore is increased by a factor of 3 compared to that of the relaxed chromophore. Also, the reaction proceeds in a completely stereoselective manner involving only the cis double bond and with a minimum quantum yield of 72%. Bond torsion is more effective than methyl substitution for fast and selective photochemistry, which is in agreement with photophysical measurements on rhodopsin analogues. We conclude that apart from the geometric distortions caused by the protein pocket it is not necessary to postulate other specific interactions between the protein and the chromophore to effect the selective and ultrafast photoreaction in rhodopsin.     

Early Time Solvation Dynamics Probed by Spectrally Resolved Degenerate Pump-Probe Spectroscopy

The dynamic role of solvent in influencing the rates of physico-chemical processes (for example, polar solvation and electron transfer) has been extensively studied using time-resolved fluorescence spectroscopy. Here we study ultrafast excited state relaxation dynamics of three different fluorescent probes (DNTTCI, IR-140 and IR-144) in two polar solvents, ethanol and ethylene glycol, using spectrally resolved degenerate pump-probe spectroscopy. We discuss how time-resolved emission spectra can be directly used for constructing relaxation correlation function, obviating spectral reconstruction and estimation of time-zero spectrum in non-polar solvents. We show that depending on the specific probe used, the relaxation dynamics is governed either by intramolecular vibrational relaxation (for IR140) or by intermolecular solvation (for DNTTCI) or by both (for IR144). We further show (using DNTTCI as a probe) that major differences in solvation by ethanol and ethylene glycol is contributed by early time (<1 ps) dynamics.     

A molecular spring for vision

Light absorption by the visual pigment rhodopsin leads to vision via a complex signal transduction pathway that is initiated by the ultrafast and highly efficient photoreaction of its chromophore, the retinal protonated Schiff base (RPSB). Here, we investigate this reaction in real time by means of unrestrained molecular dynamics simulations of the protein in a membrane mimetic environment, treating the chromophore at the density functional theory level. We demonstrate that a highly strained all-trans RPSB is formed starting from the 11-cis configuration (dark state) within approximately 100 fs by a minor rearrangement of the nuclei under preservation of the saltbridge with Glu113 and virtually no deformation of the binding pocket. Hence, the initial step of vision can be understood as the compression of a molecular spring by a minor change of the nuclear coordinates. This spring can then release its strain by altering the protein environment.