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.     

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.     

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.     

Effect of microhydration on the electronic structure of the chromophores of the photoactive yellow and green fluorescent proteins

Electronic structure calculations of microhydrated model chromophores (in their deprotonated anionic forms) of the photoactive yellow and green fluorescent proteins (PYP and GFP) are reported. Electron-detachment and excitation energies as well as binding energies of mono- and dihydrated isomers are computed and analyzed. Microhydration has different effects on the excited and ionized states. In lower-energy planar isomers, the interaction with one water molecule blueshifts the excitation energies by 0.1-0.2 eV, whereas the detachment energies increase by 0.4-0.8 eV. The important consequence is that microhydration by just one water molecule converts the resonance (autoionizing) excited states of the bare chromophores into bound states. In the lower-energy microhydrated clusters, interactions with water have negligible effect on the chromophore geometry; however, we also identified higher-energy dihydrated clusters of PYP in which two water molecules form hydrogen-bonding network connecting the carboxylate and phenolate moieties and the chromophore is strongly distorted resulting in a significant shift of excitation energies (up to 0.6 eV).     

"标记"芘的激基缔合物荧光在水溶性高分子研究中的应用

首先介绍了芘的荧光发射光谱的特点及其激基缔合物荧光、然后给出了用芘标记高分子的各种方式和方法,讨论了水溶性高分子的”标记”芘形成激基缔合物的原理,评述了芘标记高分子箕缔合物荧光在水溶性分子疏水相互作用,静电相互作用、氢键、络合等研究中的应用,也介绍了芘标记激基缔的荧光在高分子凝胶体积相变、分子链运动等领域的研究进展。     

Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore

We report that the photoinduced dynamics of the phytochrome chromophore is strongly dependent on the protonation/deprotonation states of the pyrrole ring. The on-the-fly surface hopping dynamics simulations were performed to study the photoisomerization of different protonation/deprotonation phytochrome chromophore models. The simulation results indicate that the deprotonations at the pyrrole rings significantly modify the photoinduced nonadiabatic dynamics, leading to distinctive population decay dynamics and different reaction channels. Such feature can be well explained by the formation of the different hydrogen bond network patterns. Therefore, the proper understanding of the photoisomerization mechanism of phytochrome chromophore must take the hydrogen bond network into account. This work provides the new insights into the photobiological functions of phytochrome chromophore and suggests the possible ideas to control of its photoconversion processes for further rational engineering in optical applications.     

The thiourea group modulates the fluorescence emission decay of fluorescein-labeled molecules

We have investigated the spectral properties and emission characteristics of fluorescein-5-thiocarbamoyl-N,N'-caproate (FITC-ACA) to examine the origin of the complex emission decay often observed in fluorescein-labeled molecules. The covalent attachment of fluorescein to epsilon-amino-n-caproic acid does not perturb the prototropic transitions of the chromophore or the general fluorescence characteristics of the various prototropic forms. However, both the monoanion and dianion forms of FITC-ACA are quenched relative to free fluorescein and exhibit a complex emission decay that is described by two discrete lifetimes. The thiourea group that links the chromophore to the caproic acid is shown to modulate the emission properties of the FITC-ACA. We show that the emission decay can also be analyzed using the asymmetric distribution model of Alcala et al. In this analysis, the tauL and tauu parameters that represent the lower and upper lifetime limits of the distribution reflect the quenched (0 ns) and unquenched lifetimes, respectively. The beta parameter that describes the distribution of lifetimes between the two limiting states can be related to the quenching efficiency of the thiourea group and to the structure and dynamics of the FITC-ACA molecule.     

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.     

A molecular dynamics study of Lys-Trp-Lys: structure and dynamics in solution following photoexcitation

We report studies of the structure and dynamics of a tripeptide Lys-Trp-Lys (KWK) in aqueous solution following photoexcitation by molecular dynamics simulations. For ground-state KWK, we observe three stable conformations with free energy differences of less than 5.2 kJ/mol. Each conformer is stabilized by a pi-cation interaction between one of three protonated amino groups and the indole moiety. For the excited state of tryptophan in KWK, the simulated molecular dynamics of the three isomers are similar, all in good agreement with recent femtosecond experiments (J. Phys. Chem. B 2005, 109, 16901). Specifically, we observe: (1) the fluorescence anisotropy is dominated by a single-exponential component and decays in approximately 130 ps, (2) the total dynamic Stokes shift reaches approximately 2700 cm(-1), and (3) the excited state relaxation dynamics occurs on several time scales ranging from femtoseconds to tens of picoseconds. The relaxation dynamics involve rapid initial response of neighboring water, followed by local motions of flexible peptide chains. These processes drive global restructuring of the tripeptide on a rather flat energy surface, inducing slower dynamics evident in both the water and protein contributions to the stabilization energy of the photoexcited chromophore. The water and protein dynamics are strongly correlated. On a still longer time scale, we observe isomerization of two excited state conformers to the other most stable one, an analogue for evolution of trajectories along the funnel on the rugged free energy landscape to the final "native" state. Our studies suggest new experiments to detect this unique dynamics.     

Dependence of tryptophan emission wavelength on conformation in cyclic hexapeptides

The wavelength of maximum emission of tryptophan depends on the local electrostatic environment of the indole chromophore. The time-resolved emission spectra of seven rigid cyclic hexapeptides containing a single tryptophan residue were measured. The emission maxima of the three decay-associated spectra for the seven peptides ranged from 341 to 359 nm, suggesting that different tryptophan rotamers have different emission maxima even in the case of solvent-exposed tryptophans. This conclusion is supported by quantum mechanical/molecular dynamics simulations of the six canonical side chain rotamers of tryptophan in solvated hexapeptides. The calculated range of emission maxima for the tryptophan rotamers of the seven peptides is 344-365 nm. The precision of the wavelength calculations and the peptide, water, and charged side chain contributions to the spectral shifts are examined. The results indicate that the emission maxima of decay-associated spectra can aid in the assignment of fluorescence lifetimes to tryptophan rotamers.     

Hydrogen Bond Fluctuations Control Photochromism in a Reversibly Photo‐Switchable Fluorescent Protein

Reversibly switchable fluorescent proteins (RSFPs) are essential for high‐resolution microscopy of biological samples, but the reason why these proteins are photochromic is still poorly understood. To address this problem, we performed molecular dynamics simulations of the fast switching Met159Thr mutant of the RSFP Dronpa. Our simulations revealed a ground state structural heterogeneity in the chromophore pocket that consists of three populations with one, two, or three hydrogen bonds to the phenolate moiety of the chromophore. By means of non‐adiabatic quantum mechanics/molecular dynamics simulations, we demonstrated that the subpopulation with a single hydrogen bond is responsible for off‐switching through photo‐isomerization of the chromophore, whereas two or more hydrogen bonds inhibit the isomerization and promote fluorescence instead. While rational design of new RSFPs has so far focused on structure alone, our results suggest that structural heterogeneity must be considered as well.     

Water Diffusion In And Out Of The β-Barrel Of GFP and The Fast Maturing Fluorescent Protein, TurboGFP

The chromophore of fluorescent proteins is formed by an internal cyclization of the tripeptide 65SYG67 fragment and a subsequent oxidation. The oxidation is slow - the kinetics of this step is presumably improved in fast maturing GFPs. Water molecules can aid in the chromophore formation. We have used 50 ns molecular dynamics simulations of the mature and immature forms of avGFP and TurboGFP to examine the diffusion of water molecules in-and-out of the protein β-barrel. Most crystal structures of GFPs have well-structured waters within hydrogen-bonding distance of Glu222 and Arg96. It has been proposed that they have an important role in chromophore formation. Stable waters are found in similar positions in all simulations conducted. The simulations confirm the existence of a pore that leads to the chromophore in the rapidly maturing TurboGFP; decreased water diffusion upon chromophore formation; and increased water diffusion due to the pore formation.     

Nonadiabatic molecular dynamics simulations of correlated electrons in solution. 1. Full configuration interaction (CI) excited-state relaxation dynamics of hydrated dielectrons

The hydrated dielectron is composed of two excess electrons dissolved in liquid water that occupy a single cavity; in both its singlet and triplet spin states there is a significant exchange interaction so the two electrons cannot be considered to be independent. In this paper and the following paper,we present the results of mixed quantum/classical molecular dynamics simulations of the nonadiabatic relaxation dynamics of photoexcited hydrated dielectrons, where we use full configuration interaction (CI) to solve for the two-electron wave function at every simulation time step. To the best of our knowledge, this represents the first systematic treatment of excited-state solvation dynamics where the multiple-electron problem is solved exactly. The simulations show that the effects of exchange and correlation contribute significantly to the relaxation dynamics. For example, spin-singlet dielectrons relax to the ground state on a time scale similar to that of single electrons excited at the same energy, but spin-triplet dielectrons relax much faster. The difference in relaxation dynamics is caused by exchange and correlation: The Pauli exclusion principle imposes very different electronic structure when the electrons' spins are singlet paired than when they are triplet paired, altering the available nonadiabatic relaxation pathways. In addition, we monitor how electronic correlation changes dynamically during nonadiabatic relaxation and show that solvent dynamics cause electron correlation to evolve quite differently for singlet and triplet dielectrons. Despite such differences, our calculations show that both spin states are stable to excited-state dissociation, but that the excited-state stability has different origins for the two spin states. For singlet dielectrons, the stability depends on whether the solvent structure can rearrange to create a second cavity before the ground state is reached. For triplet dielectrons, in contrast, electronic correlation ensures that the two electrons do not dissociate, even if the dielectron is artificially kept from reaching the ground state. In addition, both singlet and triplet dielectrons change shape dramatically during relaxation, so that linear response fails to describe the solvation dynamics for either spin state. In the following paper (Larsen, R. E.; Schwartz, B. J. J. Phys. Chem. B 2006, 110, 9692), we use these simulations to calculate the pump-probe spectroscopic signal expected for photoexcited hydrated dielectrons and to predict an experiment to observe hydrated dielectrons directly.     

Color tuning in photofunctional proteins

Depending on protein environment, a single photofunctional chromophore shows a wide variation of photoabsorption/emission energies. This photobiological phenomenon, known as color tuning, is observed in human visual cone pigments, firefly luciferase, and red fluorescent protein. We investigate the origin of color tuning by quantum chemical calculations on the excited states: symmetry-adapted cluster-configuration interaction (SAC-CI) method for excited states and a combined quantum mechanical (QM)/molecular mechanical (MM) method for protein environments. This Minireview summarizes our theoretical studies on the above three systems and explains a common feature of their color-tuning mechanisms. It also discuss the possibility of artificial color tuning toward a rational design of photoabsorption/emission properties.     

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.     

Open science grid study of the coupling between conformation and water content in the interior of a protein

Computational grids are a promising resource for modeling complex biochemical processes such as protein folding, penetration of gases or water into proteins, or protein structural rearrangements coupled to ligand binding. We have enabled the molecular dynamics program CHARMM to run on the Open Science Grid. The implementation is general, flexible, easily modifiable for use with other molecular dynamics programs and other grids and automated in terms of job submission, monitoring, and resubmission. The usefulness of grid computing was demonstrated through the study of hydration of the Glu-66 side chain in the interior of protein staphylococcal nuclease. Multiple simulations started with and without two internal water molecules shown crystallographically to be associated with the side chain of Glu-66 yielded two distinct populations of rotameric states of Glu-66 that differed by as much as 20%. This illustrates how internal water molecules can bias protein conformations. Furthermore, there appeared to be a temporal correlation between dehydration of the side chain and conformational transitions of Glu-66. This example demonstrated how difficult it is to get convergence even in the relatively simple case of a side chain oscillating between two conformations. With grid computing, we also benchmarked the self-guided Langevin dynamics method against the Langevin dynamics method traditionally used for temperature control in molecular dynamics simulations and showed that the two methods yield comparable results.     

Hydration dynamics and time scales of coupled water-protein fluctuations

We report experimental and theoretical studies on water and protein dynamics following photoexcitation of apomyoglobin. Using site-directed mutation and with femtosecond resolution, we experimentally observed relaxation dynamics with a biphasic distribution of time scales, 5 and 87 ps, around the site Trp7. Theoretical studies using both linear response and direct nonequilibrium molecular dynamics (MD) calculations reproduced the biphasic behavior. Further constrained MD simulations with either frozen protein or frozen water revealed the molecular mechanism of slow hydration processes and elucidated the role of protein fluctuations. Observation of slow water dynamics in MD simulations requires protein flexibility, regardless of whether the slow Stokes shift component results from the water or protein contribution. The initial dynamics in a few picoseconds represents fast local motions such as reorientations and translations of hydrating water molecules, followed by slow relaxation involving strongly coupled water-protein motions. We observed a transition from one isomeric protein configuration to another after 10 ns during our 30 ns ground-state simulation. For one isomer, the surface hydration energy dominates the slow component of the total relaxation energy. For the other isomer, the slow component is dominated by protein interactions with the chromophore. In both cases, coupled water-protein motion is shown to be necessary for observation of the slow dynamics. Such biologically important water-protein motions occur on tens of picoseconds. One significant discrepancy exists between theory and experiment, the large inertial relaxation predicted by simulations but clearly absent in experiment. Further improvements required in the theoretical model are discussed.     

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.