Molecular dynamics computer simulations of solvation dynamics at liquid/liquid interfaces

Chemical Reaction Dynamics at Liquid Interfaces: A Computational Approach

Recent advances in experimental and theoretical studies of liquid interfaces provide remarkable evidence for the unique properties of these systems. In this review we examine how these properties affect the thermodynamics and kinetics of chemical reactions which take place at the liquid/vapor interface and at the liquid/liquid interface. We demonstrate how the rapidly varying density and viscosity, the marked changes in polarity and the surface roughness manifest themselves in isomerization, electron transfer and photodissociation reactions.

Microelectrochemical Techniques for Probing Kinetics at Liquid/Liquid Interfaces

We provide an overview of recent advances in microelectrochemical approaches to investigate the kinetics of various physicochemical processes that occur at the interface between two immiscible electrolyte solutions (ITIES). To place the advances in context, background material on the structure of the ITIES, derived from both experimental studies and computer simulation, is also provided. The main focus of the article is micro-ITIES techniques, single droplet measurements, microelectrochemical measurements at expanding droplets (MEMED) and scanning electrochemical microscopy (SECM). Recent developments in a combined SECM-Langmuir trough technique for probing diffusion processes across Langmuir monolayers at the water/air (W/A) interface are also highlighted, by considering an organic monolayer at a water surface as a special case of a liquid/liquid interface.

Solvation dynamics and rotational relaxation of coumarin 153 in mixed micelles of Triton X-100 and cationic gemini surfactants: effect of composition and spacer chain length of gemini surfactants

To demonstrate simultaneously how the solvation dynamics and rotational relaxation in nonionic micelles change with the composition of a gemini surfactant and how this change depends on spacer chain length of gemini surfactants.

Calculation of the intrinsic solvation free energy profile of an ionic penetrant across a liquid-liquid interface with computer simulations

We introduce the novel concept of an intrinsic free energy profile, allowing one to remove the artificial smearing caused by thermal capillary waves, which renders difficulties for the calculation of free energy profiles across fluid interfaces in computer simulations. We apply this concept to the problem of a chloride ion crossing the interface between water and 1,2-dichloroethane and show that the present approach is able to reveal several important features of the free energy profile which are not detected with the usual, nonintrinsic calculations. Thus, in contrast to the nonintrinsic profile, a free energy barrier is found at the aqueous side of the (intrinsic) interface, which is attributed to the formation of a water "finger" the ion pulls with itself upon approaching the organic phase. Further, by the presence of a nonsampled region, the intrinsic free energy profile clearly indicates the coextraction of the first hydration shell water molecules of the ion when entering the organic phase.

The unique molecular behavior of water at the chloroform-water interface

The molecular bonding and orientation of water at the chloroform-water interface has been examined in this study using vibrational sum-frequency spectroscopy (VSFS). The results provide a key puzzle piece towards our understanding of the systematic changes in the interfacial bonding and orientation of water that occur with variations in the polarity of the organic phase, especially when compared with previous studies of different liquid-liquid interfacial systems. In these VSFS studies the OH spectral responses of interfacial water molecules are used to characterize the interactions between water and the organic phase. The spectral analysis, aided by isotopic dilution studies, shows that the moderate polarity of the chloroform phase results in a mixed interfacial region with stronger organic-water bonding and fewer bonding interactions between adjacent water molecules than was previously found for studies of non-polar organic liquid-water interfaces. Even with the more mixed interfacial region and stronger organic-water interactions, interfacial water retains a significant amount of orientational ordering. These results are compared with recent predictions from molecular dynamics simulations about how molecules behave at the chloroform-water interface.

Solvation of Transmembrane Proteins by Isotropic Membrane Mimetics: A Molecular Dynamics Study

Mixtures of organic solvents are often used as membrane mimetics in structure determination of transmembrane proteins by solution NMR; however, the mechanism through which these isotropic solvents mimic the anisotropic environment of cell membranes is not known. Here, we use molecular dynamics simulations to study the solvation thermodynamics of the c-subunit of Escherichia coli F1F0 ATP synthase in membrane mimetic mixtures of methanol, chloroform, and water with varying fractions of components as well as in lipid bilayers. We show that the protein induces a local phase separation of the solvent components into hydrophobic and hydrophilic layers, which provides the anisotropic solvation environment to stabilize the amphiphilic peptide. The extent of this effect varies with solvent composition and is most pronounced in the ternary methanol-chloroform-water mixtures. Analysis of the solvent structure, including the local mole fraction, density profiles, and pair distribution functions, reveals considerable variation among solvent mixtures in the solvation environment surrounding the hydrophobic transmembrane region of the protein. Hydrogen bond analysis indicates that this is primarily driven by the hydrogen-bonding propensity of the essential Asp(61) residue. The impact of the latter on the conformational stability of the solvated protein is discussed. Comparison with the simulations in explicit all-atom models of lipid bilayer indicates a higher flexibility and reduced structural integrity of the membrane mimetic solvated c-subunit. This was particularly true for the deprotonated form of the protein and found to be linked to solvent stabilization of the charged Asp(61).

Local environments of coumarin dyes within mesostructured silica-surfactant nanocomposites

The local environments surrounding dye molecules were studied with use of coumarin dyes in a mesostructured silica-surfactant nanocomposite, which was formed in a porous alumina membrane by a surfactant-templated method and has an average pore diameter of 3.4 nm. Coumarin dyes, such as coumarin 480 (C480), coumarin 343 (C343), and propylamide coumarin 343 (PAC343), were extracted into the silica-surfactant nanocomposite and time-resolved fluorescence spectra of these dyes were examined. C480 and C343 show slow dynamic Stokes shifts and the decay curve can be fitted by a biexponential function. The decay-time constants obtained from the fitting are almost identical for C480 and C343: 0.87 and 7.5 ns for C480, and 0.86 and 7.6 ns for C343. In contrast to these two coumarin dyes, short decay-time constants (0.50 and 4.8 ns) were obtained for PAC343 in the silica-surfactant nanocomposite. These results indicate that the local environments of C480 and C343 are almost identical but different from that of PAC343. By considering the origin of the dynamic Stokes shift and the mesostructure of the silica-surfactant nanocomposite, the location and microenvironment of coumarin dyes within the silica-surfactant nanocomposite are discussed.

Time-dependent stokes shifts of fluorescent dyes in the hydrophobic backbone region of a phospholipid bilayer: combination of fluorescence spectroscopy and ab initio calculations

We explored the time-dependent Stokes shifts of fluorescent dyes containing an anthroyloxy chromophore (2-AS, 9-AS, and 16-AP) in bilayers composed of palmitoyl-oleoyl-phosphatidylcholine. The obtained data revealed a nontrivial solvation response of these dyes, which are located in the backbone region of the bilayer with a gradually increasing depth. For comparison, steady-state emission spectra in the neat solvents of various polarities and viscosities were also recorded. The results indicate that on the short picosecond time scale the AS dyes undergo complex photophysics including formation of states with a charge-transfer character. This observation is supported by ab initio calculations of the excited states of 9-methylanthroate. The slower nanosecond part of the relaxation process can be attributed to the solvation response of the dyes. A slowdown in solvent relaxation is observed upon moving toward the center of the bilayer. A mechanism similar to preferential solvation present in the mixture of a polar and nonpolar solvent is considered to explain the obtained data.

Synergistic effect of dicarbollide anions in liquid-liquid extraction: a molecular dynamics study at the octanol-water interface

We report a molecular dynamics study of chlorinated cobalt bis(dicarbollide) anions [(B(9)C(2)H(8)Cl(3))(2)Co](-)"CCD(-)" in octanol and at the octanol-water interface, with the main aim to understand why these hydrophobic species act as strong synergists in assisted liquid-liquid cation extraction. Neat octanol is quite heterogeneous and is found to display dual solvation properties, allowing to well solubilize CCD(-), Cs(+) salts in the form of diluted pairs or oligomers, without displaying aggregation. At the aqueous interface, octanol behaves as an amphiphile, forming either monolayers or bilayers, depending on the initial state and confinement conditions. In biphasic octanol-water systems, CCD(-) anions are found to mainly partition to the organic phase, thus attracting Cs(+) or even more hydrophilic counterions like Eu(3+) into that phase. The remaining CCD(-) anions adsorb at the interface, but are less surface active than at the chloroform interface. Finally, we compare the interfacial behavior of the Eu(BTP)(3)(3+) complex in the absence and in the presence of CCD(-) anions and extractant molecules. It is found that when the CCD(-)'s are concentrated enough, the complex is extracted to the octanol phase. Otherwise, it is trapped at the interface, attracted by water. These results are compared to those obtained with chloroform as organic phase and discussed in the context of synergistic effect of CCD(-) in liquid-liquid extraction, pointing to the importance of dual solvation properties of octanol and of the hydrophobic character of CCD(-) for synergistic extraction of cations.

A polarizable continuum model for molecules at diffuse interfaces

In this work we illustrate an extension of the polarizable continuum model to describe solvation effects on molecules at the interface between two fluid phases (liquid/liquid, liquid/vapor). This extension goes beyond the naive picture of the interface as a plane dividing two distinct dielectrics, commonly employed in continuum models. The main feature of the model is the use of a diffuse interface with an electric permittivity depending on the position. This characteristic clearly allows the study of simple interfaces as well as more complex membrane or multilayer structures. Moreover the smooth variation of the permittivity in the diffuse interface, in contrast to the sharp boundary between two regions, overcomes the numerical divergences due to charges placed at the boundary. The implementation of the model relies on the integral equation formalism, which allows one to calculate the reaction field acting on a molecule immersed in a dielectric environment once the proper Green's function is known. In the present case, such a Green's function is obtained numerically, allowing a large flexibility in the choice of the dielectric permittivity profile. The applications have been selected with the aim of illustrating the capabilities of the model; its present limitations are also discussed.

A femtosecond study of excitation wavelength dependence of solvation dynamics in a PEO-PPO-PEO triblock copolymer micelle

Excitation wavelength (lambdaex) dependence of solvation dynamics of coumarin 480 (C480) in the micellar core of a water soluble triblock copolymer, PEO20-PPO70-PEO20 (Pluronic P123), is studied by femtosecond and picosecond time resolved emission spectroscopies. In the P123 micelle, the width of the emission spectrum of C480 is found to be much larger than that in bulk water. This suggests that the P123 micelle is more heterogeneous than bulk water. The steady state emission maximum of C480 in P123 micelle shows a significant red edge excitation shift by 25 nm from 453 nm at lambdaex=345 nm to 478 nm at lambdaex=435 nm. The solvation dynamics in the interior of the triblock copolymer micelle is found to depend strongly on the excitation wavelength. The excitation wavelength dependence is ascribed to a wide distribution of locations of C480 molecules in the P123 micelle with two extreme environments-a bulklike peripheral region with very fast solvent response and a very slow core region. With increase in lambdaex, contribution of the bulklike region having an ultrafast component (< or =2 ps) increases from 7% at lambdaex=375 nm to 78% at lambda(ex)=425 nm while the contribution of the ultraslow component (4500 ps) decreases from 79% to 17%.

Solvent relaxation in phospholipid bilayers: principles and recent applications

Although there exist a number of methods, such as NMR, X-ray, e.g., which explore the hydration of phospholipid bilayers, the solvent relaxation (SR) method has the advantage of simple instrumentation, easy data treatment and possibility of measuring fully hydrated samples. The main information gained from SR by the analysis of recorded "time-resolved emission spectra" (TRES) is micro-viscosity and micro-polarity of the dye microenvironment. Based on these parameters, one can draw conclusions about water structure in the bilayer. In this review, we focus on physical background of this method, on all the procedures that are needed in order to obtain relevant parameters, and on the requirements on the fluorescence dyes. Furthermore, a few recent applications (the effect of curvature, binding of antibacterial peptides and phase transition) illustrating the versatility of this method are mentioned. Moreover, limitations and potential problems are discussed.

Temperature dependence of solvation dynamics and anisotropy decay in a protein: ANS in bovine serum albumin

Temperature dependence of solvation dynamics and fluorescence anisotropy decay of 8-anilino-1-naphthalenesulfonate (ANS) bound to a protein, bovine serum albumin (BSA), are studied. Solvation dynamics of ANS bound to BSA displays a component (300 ps) which is independent of temperature in the range of 278-318 K and a long component which decreases from 5800 ps at 278 K to 3600 ps at 318 K. The temperature independent part is ascribed to a dynamic exchange of bound to free water with a low barrier. The temperature variation of the long component of solvation dynamics corresponds to an activation energy of 2.1 kcal mol(-1). The activation energy is ascribed to local segmental motion of the protein along with the associated water molecules and polar residues. The time scale of solvation dynamics is found to be very different from the time scale of anisotropy decay. The anisotropy decays are analyzed in terms of the wobbling motion of the probe (ANS) and the overall tumbling of the protein.

Fluorescence Anisotropy Decay and Solvation Dynamics in a Nanocavity: Coumarin 153 in Methyl β-Cyclodextrins

Fluorescence anisotropy decay and solvation dynamics of coumarin 153 (C153) are studied in dimethyl beta-cyclodextrin (DIMEB) and trimethyl beta-cyclodextrin (TRIMEB) nanocavity in water. C153 binds to DIMEB and TRIMEB to form both 1:1 and 1:2 (C153:cyclodextrin) complexes. The anisotropy decays of C153 in DIMEB and TRIMEB are found to be biexponential. The fast component of anisotropy decay (approximately 1000 ps) is attributed to the 1:1 complex and the slower one (approximately 2500 ps) to the 1:2 complex. From the components of the anisotropy decay, the length of the 1:1 and 1:2 complexes are estimated. Solvation dynamics of C153 in DIMEB exhibits a very fast (2.4 ps) component (41%) and two slower components of 50 ps (29%) and 1450 ps (30%). Solvation dynamics in TRIMEB is described by three slow components of 10.3 ps (24%), 240 ps (45%), and 2450 ps (31%). Possible origins of the ultraslow components are discussed.

Computer Simulation Investigation of the Water−Benzene Interface in a Broad Range of Thermodynamic States from Ambient to Supercritical Conditions

The dependence of the properties of the water-benzene system on the thermodynamic conditions in a broad range of temperatures and pressures has been investigated by computer simulation methods. For this purpose, Monte Carlo simulations have been performed at 23 different thermodynamic states, ranging from ambient to supercritical conditions. The density profiles of the water and benzene molecules have been determined at each of the thermodynamic states investigated. Information on the dependence of the mutual solubility of the two components in each other as well as of the width of the interface on the temperature and pressure has been extracted from these profiles. The width of the interface has been found to increase with increasing temperature up to a certain point, where it diverges. The temperature of this divergence corresponds to the mixing of the two phases. The determination of the critical mixing temperature at various pressures allowed us to estimate the upper critical curve, separating the two-phase and one-phase liquid systems, of the phase diagram of the simulated water-benzene system. In analyzing the preferential orientation of the interfacial molecules relative to the interface, it has been found that the main orientational preference of the benzene molecules is to lie parallel with the plane of the interface, and the water molecules penetrated deepest into the benzene phase prefer to stay perpendicular to the interface, pointing by one of their O-H bonds almost straight toward the benzene phase, whereas the waters located at the aqueous side of the interface are preferentially aligned parallel with the interfacial plane. Although the strength of the observed orientational preferences decreases rapidly with increasing temperature, the preferred orientations themselves are found to be independent of the thermodynamic conditions. Remains of the orientational preferences of the molecules are found to be present up to temperatures as high as 650 K. The analysis of the relative orientation of the neighboring water-benzene pairs has revealed that the radius of the first hydration shell of the benzene molecules is independent of the thermodynamic conditions, even if the system consists of one single phase. It has been found that the nearest water neighbors of the benzene molecules are preferentially located above and below the benzene ring, whereas more distant water neighbors, belonging still to the first hydration shell, prefer to stay within the plane of the benzene molecule. In the two-phase systems the dipole vector of the nearest waters has been found to be preferentially perpendicular to the vector pointing from the center of the benzene molecule to the water O atom.

Hydration Shell Exchange Dynamics during Ion Transfer Across the Liquid/Liquid Interface

We examine using molecular dynamics simulations the rate and mechanism of water molecules exchange around the Li(+) and Na(+) ions during ion transfer across the interface between water and nitrobenzene. As the ions are transferred from the water to the organic phase, they keep their first hydration shell and an incomplete second shell. The rate of water exchange between the first shell and the rest of the interfacial water molecule decreases during the transfer, which is consistent with an increase in the barrier along the ion-water potential of mean force. While in bulk water the exchange of water molecules around the Li(+) follows an associative (A) or associative interchange (I(a)) type mechanism, the fraction of exchange events of type A increases at the interface. In contrast, while in bulk water the exchange of water molecules around the six coordinated Na(+) hydrated species mainly follows a dissociative mechanism, the situation at the interface involves an equilibrium interchange between the four- and five-coordinated hydrated ion. Simulation of the reversed process, in which the hydrated Li(+) ion is transferred to the aqueous phase, shows the same general behavior as a function of location from the interface.

Nucleophilic Substitution Reactions at Liquid/Liquid Interfaces: Molecular Dynamics Simulation of a Model SN1 Dissociation Reaction at the Water/Carbon Tetrachloride Interface

The ionic dissociation step of the nucleophilic substitution reaction t-BuCl --> t-Bu(+) + Cl(-) is studied at the water/carbon tetrachloride interface using molecular dynamics computer simulations. The empirical valence bond approach is used to couple two diabatic states, covalent and ionic, in the electronically adiabatic limit. The umbrella sampling technique is used to calculate the potential of mean force along the reaction coordinate (defined as the t-Bu to Cl distance) at several interface regions of varying distances from the Gibbs dividing surface. We find a significant increase of the ionic dissociation barrier height and of the reaction free energy at the interface relative to bulk water. This is shown to be due to the reduced polarity of the interface which causes a destabilization of the pure ionic state. However, deformation to the neat interface structure in the form of water protrusions into the organic phase may provide partial stabilization of the ionic species. The importance of these structural effects is examined by repeating the calculations with an artificially smooth interface. The destabilization of the ionic state at the interface also manifests itself with a rapid (picosecond time scale) recombination dynamics of the ions to form the parent molecule followed by a slow vibrational relaxation.

Temperature Dependence of Anisotropy Decay and Solvation Dynamics of Coumarin 153 in γ-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 (<tau(s)>(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 (<tau(s)>(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.

Hydrogen Bond Dynamics at Water/Organic Liquid Interfaces

Hydrogen bond dynamics at the neat interface between water and a series of organic liquids are studied with molecular dynamics computer simulation. The organic liquids are nonpolar (carbon tetrachloride), weakly polar (1,2-dichloroethane), and polar (nitrobenzene). The effect of surface polarity and surface roughness is examined. The dynamics are expressed in terms of the hydrogen bond population autocorrelation functions and are found to be nonexponential and strongly dependent on the nature of the organic phase. In particular, at all interfaces, the dynamics are slower at the interface than in the bulk and sensitive to the location of the water molecules along the interface normal.

Slow Solvation Dynamics at the Active Site of an Enzyme: Implications for Catalysis

Solvation dynamics at the active site of an enzyme, glutaminyl-tRNA synthetase (GlnRS), was studied using a fluorescence probe, acrylodan, site-specifically attached at cysteine residue C229, near the active site. The picosecond time-dependent fluorescence Stokes shift indicates slow solvation dynamics at the active site of the enzyme, in the absence of any substrate. The solvation dynamics becomes still slower when the substrate (glutamine or tRNA(Gln)) binds to the enzyme. A mutant Y211H-GlnRS was constructed in which the glutamine binding site is disrupted. The mutant Y211H-GlnRS labeled at C229 with acrylodan exhibited significantly different solvent relaxation, thus demonstrating that the slow dynamics is indeed associated with the active site. Implications for catalysis and specificity have been discussed.

Proton Transfer in Nanoconfined Polar Solvents. 1. Free Energies and Solute Position

The reaction free energy curves for a model phenol-amine proton-transfer system in a confined CH3Cl solvent have been calculated by Monte Carlo simulations. The free energy curves, as a function of a collective solvent coordinate, have been obtained for several fixed reaction complex radial positions (based on the center-of-mass). A smooth, hydrophobic spherical cavity was used to confine the solvent, and radii of 10 and 15 A have been considered. Quantum effects associated with the transferring proton have been included by adding the proton zero-point energy to the classical free energy. The results indicate the reaction complex position can be an important component of the reaction coordinate for proton-transfer reactions in nanoconfined solvents.