Solvation dynamics in protein environments: comparison of fluorescence upconversion measurements of coumarin 153 in monomeric hemeproteins with molecular dynamics simulations

The complexes of the fluorescence probe coumarin 153 with apomyoglobin and apoleghemoglobin are used as model systems to study solvation dynamics in proteins. Time-resolved Stokes shift experiments are compared with molecular dynamics simulations, and very good agreement is obtained. The solvation of the coumarin probe is very rapid with approximately 60% occurring within 300 fs and is attributed to interactions with water (or possibly to the protein itself). Differences in the solvation relaxation (or correlation) function C(t) for the two proteins are attributed to differences in their hemepockets.     

Fluorescence probing of temperature-dependent dynamics and friction in ionic liquid local environments

The solvation dynamics and local orientational friction for a series of four ionic liquids have been probed using coumarin 153 (C153) as a function of temperature. These ionic liquids are comprised of nonaromatic organic cations paired with a common anion, bis(trifluoromethylsulfonyl)imide (NTf(2)-). The specific liquids are as follows: N-methyl-tri-N-butylammonium NTf(2)- (N(1444)+/NTf(2)-), N-hexyl-tri-N-butylammonium NTf(2)- (N(6444)+/NTf(2)-), N-methyl-N-butylpyrrolidinium NTf(2)- (Pyrr(14)+/NTf(2)-), and N-methyl-N-ethoxyethylpyrrolidinium NTf(2)- (Pyrr(1(2O2))+/NTf(2)-). The observed solvation dynamics and fluorescence depolarization dynamics occur over a broad range of time scales that can only be adequately fit by functions including three or more exponential components. Stretched exponential distributions cannot adequately fit our data. The solvation and reorientational dynamics of the C153 probe are studied over a range of temperatures from 278.2 to 353.2 K. For both the solvation dynamics and the probe reorientational dynamics, the observed temperature dependence is well fit by a Vogel-Tammann-Fulcher law. To correlate the observed microscopic dynamics with macroscopic physical properties, temperature-dependent viscosities are also measured. Differential scanning calorimetry is used to study the thermodynamics of the phase transitions from the liquid to supercooled liquid to glassy states. For the two tetraalkylammonium liquids, the observed melting transitions occur near 300 K, so we are able to study the dynamics in a clearly supercooled regime. Very long time scale orientational relaxation time constants dynamics on the order of 100 ns are observed in the C153 fluorescence anisotropy. These are interpreted to arise from long-lived local structures in the environment surrounding the C153 probe.     

Exploring the picosecond time domain of the solvation dynamics of coumarin 153 within beta-cyclodextrins

We report molecular dynamics simulation results of equilibrium and dynamical characteristics pertaining to the solvation of the dye coumarin 153 (C153) trapped within hydrophobic cavities of di- and trimethylated beta-cyclodextrins (CD) in aqueous solutions. We found that stable configurations of the encapsulated probe are characterized by a slanted docking, in which the plane of the C153 lies mostly parallel to one of the glucose units of the CD. "In and out" dynamical modes of the encapsulated probe present very small amplitudes. The rotational dynamics of the trapped coumarin can be cast in terms of a simple model that includes diffusive motions within a local restrictive environment coupled to the overall rotational motion of the CD. We have examined the early stages of the solvation response of the environment following a vertical excitation of the probe. Regardless of the degree of CD methylation, the water dynamical response seems to be completed within 2-3 ps and does not differ substantially from that observed for nonencapsulated probes. The CD response is characterized by a single, subpicosecond relaxation that involves intramolecular motions. We also explored dynamical modes that could account for the recently reported persistence of Stokes shifts in the nanosecond time domain. In all cases, the only sources of ultraslow dynamics that we detected were those associated with gauche-trans interconversions in primary hydroxyl chains of the CD, which do not seem to be directly connected to the electronic excitation of the probe.     

Computer simulations of the solvation dynamics of Coumarin 153 in dimethylsulfoxide

We report molecular dynamics (MD) simulations of the solvation dynamics of Coumarin 153 in liquid dimethylsulfoxide using two distinct sets of partial charges for the coumarin probe. The excited state dipole moment of the coumarin and the dynamic Stokes shift in solution depend significantly on the type of charge distributions used. Nevertheless, the overall characteristics of the solvation responses obtained from both sets of charges are very similar and show good agreement with time-dependent Stokes shift experiments. Microscopic details of the solvent reorganization around the probe are discussed in light of the charge transfer upon photoexcitation.     

Diffusion of organic dyes in a niosome immobilized on a glass surface using fluorescence correlation spectroscopy

Giant multilameller niosomes containing cholesterol and triton X-100 are studied using fluorescence correlation spectroscopy (FCS). Dynamic light scattering (DLS) data indicates formation of niosomes of broadly two different sizes (diameter)--~150 nm and ~1300 nm. This is confirmed by field emission scanning electron microscopy (FE-SEM) and confocal microscopy. The diffusion coefficient (D(t)) of three organic dyes in the niosome immobilized on a glass surface is studied using fluorescence correlation spectroscopy. On addition of the room temperature ionic liquids (RTIL) (1-methyl-3-pentylimidazolium bromide, [pmim][Br] and 1-methyl- 3-pentylimidazolium tetra-fluoroborate, [pmim][BF(4)]) the size of the niosome particles increases. The D(t) of all the organic dyes (DCM, C343 and C480) increases on addition of RTILs, indicating faster diffusion. The viscosity calculated from the D(t) of the three dyes exhibits weak probe dependence. Unlike lipid or catanionic vesicle, the D(t) values in a niosome exhibit very narrow distribution. This indicates that the niosomes are fairly homogeneous with small variation of viscosity.     

Polar solvation and solvation dynamics in supercritical CHF3: Results from experiment and simulation

Solvation dynamics of the probe trans-4-(dimethylamino)-4'-cyanostilbene (DCS) have been measured in supercritical fluoroform at 310 K (1.04 Tc) and solvent densities over the range 1.4-2.0 rho(c) using optical Kerr-gated emission spectroscopy. Steady-state measurements and computer simulations of this and the related system coumarin 153 (C153) in fluoroform are used to help interpret the observed dynamics. The solvent contribution to the Stokes shift of DCS is estimated to be 2300 +/- 400 cm(-1) and nearly density independent over the range (0.7-2.0)rho(c). Spectral response functions are bimodal and can be fit to biexponential functions having time constants of approximately 0.5 ps (85%) and 3-10 ps (15%) over the observable range ((1.4-2.0)rho(c)). Computer simulations based on a 2-site model of fluoroform and assuming an electrostatic solvation mechanism appear to properly account for the magnitude and weak density dependence of the Stokes shifts but predict much faster solvation than is observed. Possible reasons for the discrepancy are discussed.     

Photophysics of trans-4-(dimethylamino)-4'-cyanostilbene and its use as a solvation probe

Electronic structure calculations, steady-state electronic spectroscopy, and femtosecond time-resolved emission spectroscopy are used to examine the photophysics of trans-4-(dimethylamino)-4'-cyanostilbene (DCS) and its solvent dependence. Semiempirical AM1/CI calculations suggest that an anilino TICT state is a potential candidate for the emissive state of DCS in polar solvents. But observation of large and solvent-independent absorption and emission transition moments in a number of solvents (M(abs) = 6.7 +/- 0.4 D and M(em) = 7.6 +/- 0.8 D) rule out the involvement of any such state, which would have a vanishingly small transition moment. The absorption and steady-state emission spectra of DCS evolve in a systematic manner with solvent polarity, approximately as would be expected for a single, highly polar excited state. Attempts to fit the solvatochromism of DCS using standard dielectric continuum models are only partially successful when values of the solute dipole moments suggested by independent measurements are assumed. The shapes of the absorption and emission spectra of DCS change systematically with solvent polarity in a manner that is semiquantitatively reproduced using a coupled-state model of the spectroscopy. Kerr-gate emission measurements show that the emission dynamics of DCS down to subpicosecond times reflect only solvent relaxation, rather than any more complicated electronic state kinetics. The spectral response functions measured with DCS are well correlated to those previously reported for the solvation probe coumarin 153, indicating DCS to be a useful alternative probe of solvation dynamics.     

Temperature dependence of anisotropy decay and solvation dynamics of coumarin 153 in gamma-cyclodextrin aggregates

Effect of temperature on the fluorescence anisotropy decay and the ultraslow component of solvation dynamics of coumarin 153 (C153) in a gamma-cyclodextrin (gamma-CD) nanocavity are studied using a picosecond set up. The steady-state anisotropy (0.13 +/- 0.01) and residual anisotropy (0.14 +/- 0.01) in fluorescence anisotropy decay in an aqueous solution containing 7 microM C153 and 40 mM gamma-CD are found to be quite large. This indicates formation of large linear nanotube aggregates of gamma-CD linked by C153. It is estimated that >53 gamma-CD units are present in each aggregate. In these aggregates with rise in temperature, the average solvation time ((obs)) decreases markedly from 680 ps at 278 K to 160 ps at 318 K. The dynamic Stokes shift is found to decrease from 800 cm(-1) at 278 K to 250 cm(-1) at 318 K. The fraction of dynamic Stokes shift (f(d)) detected in a picosecond set up is calculated using the Fee-Maroncelli procedure. The corrected solvation time ((corr) = f(d)<(tau(s)>(obs)) displays an Arrhenius type temperature dependence. From the temperature variation, the activation energy and entropy of the solvation process are determined to be 12.5 kcal M(-1) and 28 cal M(-1) K(-1), respectively. The ultraslow component and its temperature dependence are ascribed to a dynamic exchange between bound and free water molecules.     

The complex of apomyoglobin with the fluorescent dye coumarin 153

Understanding a protein's dielectric response requires both a theoretical model and a well-defined experimental system. The former has already been proposed by Song (J. Chem. Phys. 116, 9359 [2002]). We suggest that the latter is provided by the complex of coumarin 153 (C153) with apomyoglobin (ApoMb). C153 has been exhaustively studied and has proven to be an excellent probe of the solvation dynamics of polar solvents. Myoglobin is one of the most thoroughly studied proteins. Myoglobins from a wide range of species have been subject to X-ray structural analysis and site-directed mutagenesis. Here, we demonstrate the existence of a robust C153-apomyglobin system by means of molecular dynamics simulations, equilibrium binding studies using a Job's plot and capillary electrophoresis, circular dichroism and time-resolved fluorescence. The reorganization energy of C153 bound to ApoMb is compared with that of C153 in bulk solvent using the method of Jordanides et al. (J. Phys. Chem. B 103, 7995 [1999]).     

Probing Deuterium Isotope Effect on Structure and Solvation Dynamics of Human Serum Albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up‐conversion, FRET, and single‐molecule spectroscopy. The CD spectra suggest that D₂O affects the structure of HSA, leading to a 20 % decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H₂O and D₂O, and hence, the location of the probe in the protein remains the same in the two solvents. The single‐molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96 %) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D₂O compared with H₂O. In H₂O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30 %) and a long component of 350 ps (70 %). In D₂O, we detected a short component of 4 ps (41 %) and a long component of 950 ps (59 %). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D₂O (950 ps) is 2.5‐fold slower than that in H₂O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D₂O.     

Femtosecond study of partially folded states of cytochrome C by solvation dynamics

Using femtosecond time-resolved fluorescence spectroscopy, it is shown that the solvation dynamics in the two partially folded states (IS' and IS' ') of a protein, cytochrome C, are very different. In the case of IS' (formed by the addition of 2 mM sodium dodecyl sulfate, SDS) almost the entire dynamic solvent shift of coumarin 153 (C153) is captured in a picosecond setup and the contribution of the ultrafast component (0.5 ps) is very small (5%). Solvation dynamics of IS' ' (formed by 2 mM SDS and 5 M urea) displays a major component (47%) of 1.3 ps. This indicates that the structure of IS' ' is much more open and exposed compared to that of IS'. The difference in the dynamics of IS' and IS' ' is attributed to differences in their structure, particularly near the heme region, and the presence of urea in IS' '.     

Studies on the Solvation Dynamics of Coumarin 153 in 1‐Ethyl‐3‐Methylimidazolium Alkylsulfate Ionic Liquids: Dependence on Alkyl Chain Length

Steady‐state and time‐resolved fluorescence behavior of coumarin 153 (C153) is investigated in a series of 1‐ethyl‐3‐methylimidazolium alkylsulfate ([C₂mim][C_nOSO₃]) ionic liquids differing only in the length of the linear alkyl chain (n=4, 6, and 8) in the anion. The aim of the present study is to understand the role of alkyl chain length in solute rotation and solvation dynamics of C153 in these ionic liquids. The blueshift observed in the steady‐state absorption and emission maxima of C153 on going from the C₄OSO₃ to the C₈OSO₃ system indicates increasing nonpolar character of the microenvironment of the solute with increasing length of the alkyl side chain of the anion of the ionic liquids. The average solvation time is also found to increase on changing the substituent from butyl to octyl, and this is attributed to the increase in the bulk viscosity of the ILs. A steady blueshift of the time‐zero maximum of the fluorescence spectrum with increasing alkyl chain length also indicates that the probe molecule experiences a less polar environment in the early part of the dynamics. Rotational dynamics of C153 are also analyzed by using the Stokes–Einstein–Debye (SED), Gierer–Wirtz (GW), and Dote–Kivelson–Schwartz (DKS) theories. Analyses of the results seem to suggest decoupling of the rotational motion of the probe from solvent viscosity.     

Effect of ionic liquid on diffusion in P123 gel: fluorescence correlation spectroscopy

The effect of two room-temperature ionic liquids (RTILs) on the diffusion of three fluorescent dyes in the gel phase of a triblock copolymer, (PEO)(20)-(PPO)(70)-(PEO)(20) [Pluronic P123; poly ethylene oxide (PEO), poly propylene oxide (PPO)], was studied by using fluorescence correlation spectroscopy (FCS). We used three dyes, 4-(dicyanomethylene)-2-methyl-6-(4-dimethylaminostyryl)-4H-pyran (DCM), coumarin 480 (C480), and coumarin 343 (C343). By field-emission scanning electron microscopy (FESEM), it was observed that the macroscopic structure of the P123 gel remained unaffected upon addition of RTIL. In the absence of RTIL, the diffusion coefficient (D(t)) of the hydrophobic dye DCM (1 μm(2) s(-1) at the core) is smaller than that of the other two hydrophilic dyes (7 μm(2) s(-1) for C480 and C343). On addition of RTIL, the D(t) values of all of the dyes increase, indicating a decrease in local viscosity (η(eff)). The η(eff) of the core of the RTIL-P123 gel estimated from the D(t) of DCM is lower than that of both the P123 gel (at the core η=90 cP) and RTIL (η=110 cP). It is shown that the RTIL affects the structure of the gel by modifying the size of the micellar aggregates and by penetrating the core.     

Solute rotation and solvation dynamics in an alcohol-functionalized room temperature ionic liquid

Steady-state and time-resolved fluorescence behaviors of two dipolar solutes, coumarin 153 and 4-aminophthalimide, have been studied in an alcohol-functionalized room-temperature ionic liquid, 1-(hydroxyethyl)-3-methylimidazolium bis(trifluoromethanesulfonyl)imide. The steady-state fluorescence parameters have been exploited for the estimation of the polarity of this ionic liquid and to obtain information on the hydrogen bonding interaction between the ionic liquid and the probe molecules. The time-resolved measurements have been focused on the dynamics of solvation by studying the dynamic Stokes shift in the ps-ns time scale and solute rotation by measuring the time dependence of the fluorescence anisotropy. The time-resolved anisotropy studies reveal a significant slow down of the rotational motion of one of the probe molecules. The time-dependent fluorescence Stokes shift measurements suggest that the time-resolvable part of the dynamics is biphasic in nature, highly dependent on the probe molecule and the ultrafast component is comparatively less than that in other ionic liquids. The influence of the hydrogen bonding interaction between the probe molecules and the ionic liquids on the solute rotation and the various components of the solvation dynamics is carefully analyzed in an attempt to obtain further insight into the mechanism of solvation in these novel media.     

Ultrafast proton transfer to solvent: molecularity and intermediates from solvation- and diffusion-controlled regimes

Photoinduced proton transfer (PT) from cations 6-hydroxyquinolinium (6HQc) and 6-hydroxy-1-methylquinolinium (6MQc) to water and alcohols, and solvation of the zwitterionic conjugate base 1-methylquinolinium-6-olate (6MQz) were studied with stationary and transient absorption spectroscopy and by quantum chemical calculations. Transient emission spectra from 6MQz in acetonitrile and protic solvents shift dynamically to the red without changing their shape and intensity. The shift matches the solvation correlation function C(t) either measured with known solvatochromic probes coumarin 343 and coumarin 153 or derived from infrared/dielectric-loss data on neat solvents. This indicates that 6MQz monitors the solvation dynamics and that no intramolecular electron transfer occurs on a subpicosecond or longer time scale. The PT dynamics S(t) from 6HQc and 6MQc closely follows C(t), being initially 2-3 times slower. This allows for the conclusion that PT is controlled by solvation, with a barrier of 2 kJ/mol. In water, a pre-condition of this ultrafast reaction seems to be hydrogen-bonding between the negatively charged oxygen and two water molecules, resulting in a complex 6HQc:H2O:H2O. The complex is stable due to a high (47 kJ/mol) bonding energy between 6HQc and a water molecule. In acetonitrile, the reaction equilibrium is strongly shifted to the cation. There an intermediate PT state was detected, which may be ascribed to the cationic form 6HQc:H2O due to residual water impurities. In water-acetonitrile mixtures, the ultrafast solvent-controlled PT is followed by a diffusion-controlled reaction; the measured rate kD approximately 1010 s-1 M-1 is characteristic for simple bimolecular diffusion. The dependence of the short-time PT signal on water concentration can be fitted with a Poisson distribution of water molecules around the cation. Altogether, the short-time and long-time behaviors provide strong evidence that diffusion of only one water molecule is sufficient to detach the proton. Subsequent solvent stabilization of the products completes the PT reaction.     

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.     

Dynamics of solvent and rotational relaxation of coumarin-153 in room-temperature ionic liquid 1-butyl-3-methyl imidazolium tetrafluoroborate confined in poly(oxyethylene glycol) ethers containing micelles

We have investigated solvent and rotational relaxation of coumarin 153 (C-153) in room-temperature ionic liquid (RTILs) 1-butyl-3-methyl-imidazolium tetrafluoroborate ([bmim][BF(4)]) and the ionic liquid confined in alkyl poly(oxyethylene glycol) ethers containing micelles. We have used octaethylene glycol monotetradecyl ether (C(14)E(8)) and octaethylene glycol monododecyl ether (C(12)E(8)) as surfactants. In the [bmim][BF(4)]-C(14)E(8) micelle, we have observed only a 22% increase in solvation time compared to neat [bmim][BF(4)], whereas in the [bmim][BF(4)]-C(12)E(8) system, we have observed approximately 57% increase in average solvation time due to micelle formation. However, the slowing down in solvation time on going from neat RTIL to RTIL-confined micelles is much smaller compared to that on going from water to water confined micellar aggregates. The 22-57% increase in solvation time is attributed to the slowing down of collective motions of cations and anions in micelles. The rotational relaxation times become faster in both the micelles compare to neat [bmim][BF(4)].     

Solvation dynamics in triton-X-100 and triton-X-165 micelles: effect of micellar size and hydration

Dynamic Stokes' shift measurements using coumarin 153 as the fluorescence probe have been carried out to study solvation dynamics in two nonionic micelles, viz., triton-X-100 (TX-100) and triton-X-165 (TX-165). In both the micelles, the solvent relaxation dynamics is biexponential in nature. While the fast solvation time tau(s1) is seen to be almost similar for both the micelles, the slow solvation time tau(s2) is found to be appreciably smaller in TX-165 than in TX-100 micelle. Dynamic light scattering measurements indicate that the TX-165 micelles are substantially smaller in size than that of TX-100. Assuming similar core size for both the micelles, as expected from the similar chemical structures of the nonpolar ends for both the surfactants, the Palisade layer is also indicated to be substantially thinner for TX-165 micelles than that of TX-100. The aggregation number of TX-165 micelles is also found to be substantially smaller than that of TX-100 micelles. Fluorescence spectral studies of C153 dye in the two micelles indicate that the Palisade layer of TX-165 micelles is more polar than that of TX-100 micelles. Fluorescence anisotropy measurements indicate that the microviscosity in the Palisade layer of TX-165 micelles is also lower than that of TX-100 micelles. Based on these results it is inferred that the structure of the Palisade layer of TX-165 micelles is quite loose and have higher degree hydration in comparison to that of TX-100 micelles. Due to these structural differences in the Palisade layers of TX-165 and TX-100 micelles the solvation dynamics is faster in the former micelles than in the latter. It has been further inferred that in the present systems the collective response of the water molecules at somewhat away from the probes is responsible for the faster component of the solvation time, which does not reflect much of the structural changes of the micellar Palisade layer. On the contrary, the slower solvation time component, which is mainly due to the single particle response arising from water molecules adjacent to the probe in the micellar Palisade layer, is largely affected by the structural changes in the micellar Palisade layer.