Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulation

The solvation dynamics of ionic liquids have been the subject of intense experimental study but remain poorly understood. We present the results of molecular dynamics simulations of the solvation dynamics of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate in response to photoexcitation of the fluorescent dye coumarin-153. We reproduce the time-resolved fluorescence Stokes shift using linear response theory, then use novel statistical techniques to analyze cation and anion contributions to the signal. We find that the solvation dynamics are dominated by collective ionic motion and characterize the time scale for various features of the collective response. Further, we use the Steele analysis [Mol. Phys. 61, 1031 (1987)] to characterize the contributions to the observed Stokes shift made by translational and rovibrational degrees of freedom. Our results indicate that in contrast to molecular liquids, the rovibrational response is trivial and the observed fluorescence response arises almost entirely from ionic translation. Our results resolve previously open questions in the literature about the nature of the rapid dynamics in room-temperature ionic liquids and offer insight into the physical principles governing ionic liquid behavior on longer time scales.     

Polar solvation dynamics of lysozyme from molecular dynamics studies

The solvation dynamics of a protein are believed to be sensitive to its secondary structures. We have explored such sensitivity in this article by performing room temperature molecular dynamics simulation of an aqueous solution of lysozyme. Nonuniform long-time relaxation patterns of the solvation time correlation function for different segments of the protein have been observed. It is found that relatively slower long-time solvation components of the α-helices and β-sheets of the protein are correlated with lower exposure of their polar probe residues to bulk solvent and hence stronger interactions with the dynamically restricted surface water molecules. These findings can be verified by appropriate experimental studies.     

Intramolecular charge-transfer dynamics in covalently linked perylene-dimethylaniline and cyanoperylene-dimethylaniline

The excited-state dynamics of covalently linked electron donor-acceptor systems consisting of N, N-dimethylaniline (DMA) as electron donor and either perylene (Pe) or cyanoperylene (CNPe) as acceptor has been investigated in a large variety of solvents, including a room-temperature ionic liquid, by using femtosecond time-resolved fluorescence and absorption spectroscopy. The negligibly small solvent dependence of the absorption spectrum of both compounds and the strong solvatochromism of the fluorescence are interpreted by a model where optical excitation results in the population of a locally excited state (LES) and emission takes place from a charge-separated state (CSS). This interpretation is supported by the fluorescence up-conversion and the transient absorption measurements that reveal substantial spectral dynamics in polar solvents only, occurring on time scales going from a few hundreds of femtoseconds in acetonitrile to several tens of picoseconds in the ionic liquid. The early transient absorption spectra are similar to those found in nonpolar solvents and are ascribed to the LES absorption. The late spectra due to CSS absorption show bands that are red-shifted relative to those of the radical anion of the acceptor moiety by an amount that depends on solvent polarity, pointing to partial charge separation. Global analysis of the time-resolved data indicates that the charge separation dynamics in PeDMA is essentially solvent controlled, whereas that in CNPeDMA is faster than diffusive solvation, this difference being accounted for by a larger driving force for charge separation in the latter. On the other hand, the CSS lifetime of PeDMA is of the order of a few nanoseconds independently of the solvent, whereas that of CNPeDMA decreases with increasing solvent polarity from a few nanoseconds to a few hundreds of picoseconds. Comparison of these results with previously published data on the fluorescence quenching of Pe and CNPe in pure DMA shows that the charge separation and the ensuing charge recombination occur on similar time scales independently of whether these processes are intra- or intermolecular.     

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.     

On the microscopic theory of polar solvation dynamics

A microscopic model of the time-resolved Stokes shift is developed. The model calculates the solvation dynamics by combining the atomic resolution of the solute structure with dipolar dynamics from the polarization structure factors of the homogeneous solvent. Calculations are made for coumarin 153 and quinoxaline optical dyes with atomic geometries and charge distributions taken from quantum calculations. Stokes shift dynamics is calculated and compared to experiment in high-temperature acetonitrile and methanol and in low-temperature 2-methyl-tetrahydrofurane using dielectric relaxation data from experiment.     

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.     

Theory of solvation in polar nematics

We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phases. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates the experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra and makes several experimentally testable predictions. The equilibrium free energy of solvation, anisotropic in the nematic phase, is given by a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for the equilibrium Stokes shift and the Stokes shift time correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl and 4,4-n-heptyl-cyanopiphenyl solvents as a function of temperature in both the nematic and isotropic phases.     

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.     

Solvation dynamics of the electron produced by two-photon ionization of liquid polyols. 1. Ethylene glycol

Solvated electrons have been produced in ethylene glycol by two-photon ionization of the solvent with 263 nm femtosecond laser pulses. The two-photon absorption coefficient of ethylene glycol at 263 nm is determined to be beta = (2.1 +/- 0.2) x 10(-11) m W(-1). The dynamics of electron solvation in ethylene glycol has been studied by pump-probe transient absorption spectroscopy. So, time-resolved absorption spectra ranging from 430 to 710 nm have been measured. A blue shift of the spectra is observed for the first tens of picoseconds. Using the Bayesian data analysis method, the observed solvation dynamics are reconstructed with different models: stepwise mechanisms, continuous relaxation models, or combinations of stepwise and continuous relaxation. Comparison between models is in favor of continuous relaxation, which is mainly governed by solvent molecular motions.     

Solvation and rotational dynamics of coumarin 153 in ionic liquids: comparisons to conventional solvents

Steady-state and time-resolved emission spectroscopy with 25 ps resolution are used to measure equilibrium and dynamic aspects of the solvation of coumarin 153 (C153) in a diverse collection of 21 room-temperature ionic liquids. The ionic liquids studied here include several phosphonium and imidazolium liquids previously reported as well as 12 new ionic liquids that incorporate two homologous series of ammonium and pyrrolidinium cations. Steady-state absorption and emission spectra are used to extract solvation free energies and reorganization energies associated with the S0 <--> S1 transition of C153. These quantities, especially the solvation free energy, vary relatively little in ionic liquids compared to conventional solvents. Some correlation is found between these quantities and the mean separation between ions (or molar volume). Time-resolved anisotropies are used to observe solute rotation. Rotation times measured in ionic liquids correlate with solvent viscosity in much the same way that they do in conventional polar solvents. No special frictional coupling between the C153 and the ionic liquid solvents is indicated by these times. But, in contrast to what is observed in most low-viscosity conventional solvents, rotational correlation functions in ionic liquids are nonexponential. Time-resolved Stokes shift measurements are used to characterize solvation dynamics. The solvation response functions in ionic liquids are also nonexponential and can be reasonably represented by stretched-exponential functions of time. The solvation times observed are correlated with the solvent viscosity, and the much slower solvation in ionic liquids compared to dipolar solvents can be attributed to their much larger viscosities. Solvation times of the majority of ionic liquids studied appear to follow a single correlation with solvent viscosity. Only liquids incorporating the largest phosphonium cation appear to follow a distinctly different correlation.     

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.     

Site-specific solvation determined by intermolecular nuclear Overhauser effect--measurements and molecular dynamics

Site-specific solvation has been determined by intermolecular NOE measurements between solvent and solute. The experimental effect is shown on the four compounds 2-butanol, L-alanyl-L-tryptophan (Ala-Trp), adenosine and the disodium salt of adenosine 5'-monophosphate (5'-AMP) in the two solvents water and dimethyl sulfoxide (DMSO). The strength of NOE transfer correlates with the average distribution of solvent molecules around the corresponding solvation sites represented by the number of solvent molecules in a first solvation sphere, which can be obtained from molecular dynamics simulations in water. Saturation transfer between exchanging protons explains some deviations from this correlation. The NOE transfer measurements provide information on specific solute-solvent interactions and contribute to a better understanding of solvation phenomena. On the basis of a distinct relationship between steric solvation hindrance and the strength of NOE transfer, the application of such measurements for conformational analysis has been demonstrated for the first time.     

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.     

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.     

A surface science approach to ultrafast electron transfer and solvation dynamics at interfaces

Excess electrons in polar media, such as water or ice, are screened by reorientation of the surrounding molecular dipoles. This process of electron solvation is of vital importance for various fields of physical chemistry and biology as, for instance, in electrochemistry or photosynthesis. Generation of such excess electrons in bulk water involves either photoionization of solvent molecules or doping with e.g. alkali atoms, involving possibly perturbing interactions of the system with the parent-cation. Such effects are avoided when using a surface science approach to electron solvation: in the case of polar adsorbate layers on metal surfaces, the substrate acts as an electron source from where photoexcited carriers are injected into the adlayer. Besides the investigation of electron solvation at such interfaces, this approach allows for the investigation of heterogeneous electron transfer, as the excited solvated electron population continuously decays back to the metal substrate. In this manner, electron transfer and solvation processes are intimately connected at any polar adsorbate-metal interface. In this tutorial review, we discuss recent experiments on the ultrafast dynamics of photoinduced electron transfer and solvation processes at amorphous ice-metal interfaces. Femtosecond time-resolved two-photon photoelectron spectroscopy is employed as a direct probe of the electron dynamics, which enables the analysis of all elementary processes: the charge injection across the interface, the subsequent electron localization and solvation, and the dynamics of electron transfer back to the substrate. Using surface science techniques to grow and characterize various well-defined ice structures, we gain detailed insight into the correlation between adsorbate structure and electron solvation dynamics, the location (bulk versus surface) of the solvation site, and the role of the electronic structure of the underlying metal substrate on the electron transfer rate.     

Effect of nonpolar solvents on the solute rotation and solvation dynamics in an imidazolium ionic liquid

Recognizing the potential of the mixed solvent systems comprising ionic liquid as one of the constituents in real applications, the steady-state and time-resolved fluorescence behavior of C153 has been studied in neat 1-butyl-3-methylimidazolium hexafluorophosphate and its mixtures with nonpolar solvents, namely, toluene and 1,4-dioxane. No significant effect of the cosolvent on the steady-state absorption or fluorescence spectra of C153 in ionic liquid has been observed. Time-resolved fluorescence anisotropy measurements show a decrease of the rotational correlation time of C153 with gradual addition of the cosolvent. Solvation dynamics in ionic liquid-cosolvent mixtures is found to be biphasic, and a decrease of the average solvation time is observed with increasing amount of the cosolvent in solution. The time-zero spectrum of C153 is found to shift toward higher energy with gradual addition of the nonpolar solvent, suggesting that the probe molecule experiences a more nonpolar environment at the early stage of the dynamics in mixed solvents. The blue shift of the time-zero spectrum caused by the addition of the nonpolar solvent results in a larger Stokes shift of the time-dependent spectra due to solvent relaxation in mixed solvents. A comparison of the time-dependent spectral data of the ionic liquid-toluene and ionic liquid-dioxane systems shows that, while a small amount of toluene can significantly affect the dynamics, comparatively, a larger amount of dioxane is required to bring about the same effect. This is explained in terms of favorable interactions between toluene and the imidazolium ring system leading to a more effective solubilization of toluene in the cybotactic region of the probe.     

Van der Waals interactions dominate ligand-protein association in a protein binding site occluded from solvent water

In the present study we examine the enthalpy of binding of 2-methoxy-3-isobutylpyrazine (IBMP) to the mouse major urinary protein (MUP), using a combination of isothermal titration calorimetry (ITC), NMR, X-ray crystallography, all-atom molecular dynamics simulations, and site-directed mutagenesis. Global thermodynamics data derived from ITC indicate that binding is driven by favorable enthalpic contributions, rather than a classical entropy-driven signature that might be expected given that the binding pocket of MUP-1 is very hydrophobic. The only ligand-protein hydrogen bond is formed between the side-chain hydroxyl of Tyr120 and the ring nitrogen of the ligand in the wild-type protein. ITC measurements on the binding of IBMP to the Y120F mutant demonstrate a reduced enthalpy of binding, but nonetheless binding is still enthalpy dominated. A combination of solvent isotopic substitution ITC measurements and all-atom molecular dynamics simulations with explicit inclusion of solvent water suggests that solvation is not a major contributor to the overall binding enthalpy. Moreover, hydrogen/deuterium exchange measurements suggest that there is no significant contribution to the enthalpy of binding derived from "tightening" of the protein structure. Data are consistent with binding thermodynamics dominated by favorable dispersion interactions, arising from the inequality of solvent-solute dispersion interactions before complexation versus solute-solute dispersion interactions after complexation, by virtue of poor solvation of the binding pocket.