Molecular simulation study of water mobility in aerosol-OT reverse micelles

In this work, we present results from molecular dynamics simulations on the single-molecule relaxation of water within reverse micelles (RMs) of different sizes formed by the surfactant aerosol-OT (AOT, sodium bis(2-ethylhexyl)sulfosuccinate) in isooctane. Results are presented for RM water content w(0) = [H(2)O]/[AOT] in the range from 2.0 to 7.5. We show that translational diffusion of water within the RM can, to a good approximation, be decoupled from the translation of the RM through the isooctane solvent. Water translational mobility within the RM is restricted by the water pool dimensions, and thus, the water mean-squared displacements (MSDs) level off in time. Comparison with models of diffusion in confined geometries shows that a version of the Gaussian confinement model with a biexponential decay of correlations provides a good fit to the MSDs, while a model of free diffusion within a sphere agrees less well with simulation results. We find that the local diffusivity is considerably reduced in the interfacial region, especially as w(0) decreases. Molecular orientational relaxation is monitored by examining the behavior of OH and dipole vectors. For both vectors, orientational relaxation slows down close to the interface and as w(0) decreases. For the OH vector, reorientation is strongly affected by the presence of charged species at the RM interface and these effects are especially pronounced for water molecules hydrogen-bonded to surfactant sites that serve as hydrogen-bond acceptors. For the dipole vector, orientational relaxation near the interface slows down more than that for the OH vector due mainly to the influence of ion-dipole interactions with the sodium counterions. We investigate water OH and dipole reorientation mechanisms by studying the w(0) and interfacial shell dependence of orientational time correlations for different Legendre polynomial orders.     

Excited-state double proton transfer of 7-azaindole in water nanopools

The excited-state proton-transfer dynamics of 7-azaindole occurring in the water nanopools of reverse micelles has been investigated by measuring time-resolved fluorescence spectra and kinetics, as well as static absorption and emission spectra, with varying water content and isotope. 7-Azaindole molecules are found to exist in the bound-water regions of reverse micelles. The rate constant and the kinetic isotope effect of proton transfer are smaller than those in bulk water although both increase with the size of the water nanopool. The retardation of proton transfer in the bound regions is attributed to the increased free energy of prerequisite solvation to form a cyclically H-bonded 1:1 7-azaindole/water complex.     

Validation and divergence of the activation energy barrier crossing transition at the AOT/lecithin reverse micellar interface

In this report, the validity and divergence of the activation energy barrier crossing model for the bound to free type water transition at the interface of the AOT/lecithin mixed reverse micelle (RM) has been investigated for the first time in a wide range of temperatures by time-resolved solvation of fluorophores. Here, picosecond-resolved solvation dynamics of two fluorescent probes, ANS (1-anilino-8-naphthalenesulfonic acid, ammonium salt) and Coumarin 500 (C-500), in the mixed RM have been carefully examined at 293, 313, 328, and 343 K. Using the dynamic light scattering (DLS) technique, the size of the mixed RMs at different temperatures was found to have an insignificant change. The solvation process at the reverse micellar interface has been found to be the activation energy barrier crossing type, in which interface-bound type water molecules get converted into free type water molecules. The activation energies, Ea, calculated for ANS and C-500 are 7.4 and 3.9 kcal mol(-1), respectively, which are in good agreement with that obtained by molecular dynamics simulation studies. However, deviation from the regular Arrhenius type behavior was observed for ANS around 343 K, which has been attributed to the spatial heterogeneity of the probe environments. Time-resolved fluorescence anisotropy decay of the probes has indicated the existence of the dyes in a range of locations in RM. With the increase in temperature, the overall anisotropy decay becomes faster revealing the lability of the microenvironment at elevated temperatures.     

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.     

Probing the behavior of confined water by proton-transfer reactions

The picosecond dynamics of a bifunctional and H-bonding molecule, 7-hydroxyquinoline (7HQ), has been studied in a reverse micelle with increasing water content. The fluorescence kinetics has a complex behavior as the water content is changed. All reactions are irreversible, and a two-step mechanism is invoked to explain the observations. H2O/D2O exchange and excitation energy effects show that the second step has a higher barrier and that the corresponding reaction occurs through tunneling. The results clearly indicate two regimes of water nanopool behavior switching at W0 approximately 5 (W0 = [water]/[surfactant]). Water collective dynamics explains these observations. The lower fluidity of confined water within the reverse micelle with respect to normal bulk water alters the related H-bond network dynamics and therefore is responsible for the slower proton-transfer processes.     

Water dynamics and proton transfer in nafion fuel cell membranes

The dynamics of water and its effect on proton transport kinetics in Nafion membranes are compared at several hydration levels. Nafion is the most widely used polyelectrolyte membrane in fuel cells. Ultrafast infrared spectroscopy of the O-D stretch of dilute HOD in H2O provides a probe of the local environment and hydrogen bond network dynamics of water confined in the hydrophilic regions of Nafion. The kinetics of proton transfer in Nafion are tracked by following the excited-state proton transfer and recombination kinetics of a molecular probe, pyranine (HPTS). The hydrophilic domains of Nafion grow with increased hydration, and the interfacial regions reorganize, leading to a changing local environment for water near the interface. Swelling is not uniform throughout the membrane, and heterogeneity is observed in the fluorescence anisotropy decays of the methoxy derivative of pyranine. Measurements of the time-dependent anisotropy of water in Nafion provide a direct probe of the hydrogen bond network dynamics. These dynamics, as well as the rate of proton transport over nanoscopic distances, are observed to slow significantly as the hydration level of the membrane decreases. The results provide insights into the influence of changes in the dynamics of water on the proton-transfer processes.     

Rapid distinction of intracellular and extracellular proteins using NMR diffusion measurements

In-cell NMR spectroscopy offers a unique opportunity to begin to investigate the structures, dynamics, and interactions of molecules within their functional environments. An essential aspect of this technique is to define whether observed signals are attributable to intracellular species rather than to components of the extracellular medium. We report here the results of NMR measurements of the diffusion behavior of proteins expressed within bacterial cells, and find that these experiments provide a rapid and nondestructive probe of localization within cells and can be used to determine the size of the confining compartment. We show that diffusion can also be exploited as an editing method to eliminate extracellular species from high-resolution multidimensional spectra, and should be applicable to a wide range of problems. This approach is demonstrated here for a number of protein systems, using both (15)N and (13)C (methyl-TROSY) based acquisition.     

Femtosecond dynamics of piroxicam structures in solutions

We report on studies of femtosecond dynamics of a nonsteroidal anti-inflammatory drug, piroxicam (1), in water at three different pHs and for comparison in two aprotic solvents. An ultrafast excited-state proton-transfer (ESIPT) process takes place in neutral and cationic enol-type structures. Femtosecond emission and transient absorption experiments show that this reaction is barrierless, and the proton-transferred keto tautomer is formed in less than 100 fs in both organic solvents and acidic water. In neutral and alkaline water, the ESIPT is not present because of the prevalence of the anion structures at the ground state. For the excited anions (pH = 7, 11) and formed keto cations (pH approximately 3), an intramolecular charge-transfer process takes place in times shorter than 0.3 ps. The formed structures have a fluorescence lifetime of approximately 2-5 ps, depending on pH. In contrast, the internal twisting motion in organic solvents is slower (approximately 0.5-1 ps) and gives rotamers with lifetimes of tens of picoseconds. These results clearly show strong interactions of 1 with water, significantly distinct from those present in organic aprotic solvents. We believe that the results are important for a better understanding on short time interactions of drugs with their environment.     

Study of molecular behavior in a water nanocluster: size and temperature effect

Temperature and size effects on the behavior of nanoscale water molecule clusters are investigated by molecular dynamics simulations. The flexible three-centered (F3C) water potential is used to model the inter- and intramolecular interactions of the water molecule. The differences between the structural properties for the surface region and those for the interior region of the cluster are also investigated. It is found that as the temperature rises, the average number of hydrogen bonds per water molecule decreases, but the ratio of surface water molecules increases. After comparing the water densities in interior regions and the average number of hydrogen bonds in those regions, we find there is no apparent size effect on water molecules in the interior region, whereas the size of the water cluster has a significant influence on the behavior of water molecules at the surface region.     

How does the urea dynamics differ from water dynamics inside the reverse micelle?

In this study, the urea dynamics inside AOT reverse micelle (RM) has been monitored without intervention of water using time-resolved fluorescence techniques from the picosecond to nanosecond time regime. It has been observed that urea dynamics inside the reverse micelle is severely retarded compared to water RM due to the formation of highly networked urea cluster inside the RM. Time-resolved fluorescence anisotropy study also confirms the existence of a confined environment around the dye at higher concentrations of urea inside the reverse micelle. The dynamics of urea-water mixtures inside AOT reverse micelle has also been monitored with increasing urea concentration to get insight about the effect of urea on the overall solvation dynamics feature. It has been observed that with the increase in urea concentration, the overall dynamics becomes slower, and it infers the presence of few water or urea molecules, those strongly associated with surrounding urea and (or) water by hydrogen bonds.     

Modulation of dynamics and reactivity of water in reverse micelles of mixed surfactants

In this contribution, we attempt to correlate the change in water dynamics in a reverse micellar (RM) core caused by the modification of the interface by mixing an anionic surfactant, sodium bis(2-ethylhexyl) sulfosuccinate (AOT), and a nonionic surfactant, tetraethylene glycol monododecyl ether (Brij-30), at different proportions, and its consequent effect on the reactivity of water, measured by monitoring the solvolysis reaction of benzoyl chloride (BzCl). The dimension of the RM droplets at different mixing ratios of AOT and Brij-30 (X(Brij-30)) has been measured using dynamic light scattering (DLS) technique. The physical properties of the RM water have been determined using Fourier transform infrared spectroscopy (FTIR) and compressibility studies, which show that with increasing X(Brij-30), the water properties tend toward that of bulk-like water. The solvation dynamics, probed by coumarin 500 dye, gets faster with X(Brij-30). The rotational anisotropy studies along with a wobbling-in-cone analysis show that the probe experiences less restriction at higher X(Brij-30). The kinetics of the water-mediated solvolysis also gets faster with X(Brij-30). The increased rate of solvolysis has been correlated with the accelerated solvation dynamics, which is another consequence of surfactant headgroup-water interaction.     

Solvation dynamics under a microscope: single giant lipid vesicle

Picosecond spectroscopy under a confocal microscope is employed to study solvation dynamics of coumarin 153 (C153) inside a single giant lipid vesicle (1,2-dilauroyl-sn-glycero-3-phosphocholine, DLPC) of diameter 20 μm. Fluorescence correlation spectroscopy (FCS) indicates that the diffusion coefficient (D(t)) of the probe (coumarin153, C153) in the immobilized vesicle displays a wide distribution from ~3 to 21 μm(2) s(-1). The distribution of D(t) suggests that the microenvironment of the probe (C153) is highly heterogeneous and the local friction is different for probe molecules in different regions. The values of D(t) is significantly smaller than that for the same dye in bulk water (550 μm(2) s(-1)). This suggests that the probe is located in the interface or membrane region rather than in the water pool of the vesicle. The solvation time of C153 in different regions of the lipid vesicle varies between 750 to 1200 ps. This result clearly shows that a confocal microscope is able to resolve the spatial heterogeneity in local friction (i.e., D(t)) and solvation dynamics within a lipid vesicle.     

Hydration layer of a cationic micelle, C(10)TAB: structure, rigidity, slow reorientation, hydrogen bond lifetime, and solvation dynamics

We report a theoretical study of the structure and dynamics of the water layer (the hydration layer) present at the surface of the cationic micelle decyltrimethylammonium bromide (DeTAB) by using atomistic molecular dynamics simulations. The simulated micelle consisted of 47 surfactant molecules (and an equal number of bromide ions), in good agreement with the pioneering light scattering experiments by Debye which found an aggregation number of 50. In this micelle, three partially positively charged methyl groups of each surfactant headgroup face the surrounding water. The nature of the cationic micellar surface is found to play an important role in determining the arrangement of water which is quite different from that in the bulk or on the surface of an anionic micelle, like cesium perfluorooctanoate. Water molecules present in the hydration layer are found to be preferentially distributed in the region between the three partially charged methyl headgroups. It is found that both the translational and rotational motions of water exhibit appreciably slower dynamics in the layer than those in the bulk. The solvation time correlation function (TCF) of bromide ions exhibits a long time component which is found to originate primarily from the interaction of the probe with the micellar headgroups. Thus, the decay of the solvation TCF is controlled largely by the residence time of the probe in the surface. The residence time distribution of the water molecules also exhibits a slow time component. We also calculate the collective number density fluctuation in the layer and find a prominent slow component compared to the similar quantity in the bulk. This slow component demonstrates that water structure in the hydration layer is more rigid than that in the bulk. These results demonstrate that the slow dynamics of hydration layer water is generic to macromolecular surfaces of either polarity.     

Microscopic Symmetry Imposed by Rotational Symmetry Boundary Conditions in Molecular Dynamics Simulation

A large number of viral capsids, as well as other macromolecular assemblies, have icosahedral structure or structures with other rotational symmetries. This symmetry can be exploited during molecular dynamics (MD) to model in effect the full viral capsid using only a subset of primary atoms plus copies of image atoms generated from rotational symmetry boundary conditions (RSBC). A pure rotational symmetry operation results in both primary and image atoms at short range, and within nonbonded interaction distance of each other, so that nonbonded interactions can not be specified by the minimum image convention and explicit treatment of image atoms is required. As such an unavoidable consequence of RSBC is that the enumeration of nonbonded interactions in regions surrounding certain rotational axes must include both a primary atom and its copied image atom, thereby imposing microscopic symmetry for some forces. We examined the possibility of artifacts arising from this imposed microscopic symmetry of RSBC using two simulation systems: a water shell and human rhinovirus 14 (HRV14) capsid with explicit water. The primary unit was a pentamer of the icosahedron, which has the advantage of direct comparison of icosahedrally equivalent spatial regions, for example regions near a 2-fold symmetry axis with imposed symmetry and a 2-fold axis without imposed symmetry. Analysis of structural and dynamic properties of water molecules and protein atoms found similar behavior near symmetry axes with imposed symmetry and where the minimum image convention fails compared with that in other regions in the simulation system, even though an excluded volume effect was detected for water molecules near the axes with imposed symmetry. These results validate the use of RSBC for icosahedral viral capsids or other rotationally symmetric systems.     

What can you learn from a molecular probe? New insights on the behavior of C343 in homogeneous solutions and AOT reverse micelles

The behavior of C343, a common molecular probe utilized in solvation dynamics experiments, was studied in homogeneous media and in aqueous and nonaqueous reverse micelles (RMs). In homogeneous media, the Kamlet and Taft solvatochromic comparison method quantified solute-solvent interactions from the absorption and emission bands showing that the solvatochromic behavior of the dye depends not only on the polarity of the medium but also on the hydrogen-bonding properties of the solvent. Specifically, in the ground state the molecule displays a bathochromic shift with the polarity polarizability (pi) and the H-bond acceptor (beta) ability of the solvents and a hypsochromic shift with the hydrogen donor ability (alpha) of the media. The carboxylic acid group causes C343 to display greater sensitivity to the beta than to the pi polarity parameter; this sensitivity increases in the excited state, while the dependence on alpha vanishes. This demonstrates that C343 forms a stable H-bond complex with solvents with high H-bond acceptor ability (high beta) and low H-bond donor character (low alpha). Spectroscopy in nonpolar solvents reveals J-aggregate formation. With information from the Kamlet-Taft analysis, C343 was used to explore RMs composed of water or polar solvents/sodium 1,4-bis-2-ethylhexylsulfosuccinate (AOT)/isooctane using absorption, emission, and time-resolved spectroscopies. Sequestered polar solvents included ethylene glycol (EG), formamide (FA), N,N-dimethylformamide (DMF), and N,N-dimethylacetamide (DMA). Dissolved in the AOT RM systems at low concentration, C343 exists as a monomer, and when introduced to the RM samples in its protonated form, C343 remains protonated driving it to reside in the interface rather than the water pool. The solvathochromic behavior of the dye depends the specific polar solvent encapsulated in the RMs, revealing different types of interactions between the solvents and the surfactant. EG and water H-bond with the AOT sulfonate group destroying their bulk H-bonded structures. While water remains well segregated from the nonpolar regions, EG appears to penetrate into the oil side of the interface. In aqueous AOT RMs, C343 interacts with neither the sulfonate group nor the water, perhaps because of intramolecular H-bonding in the dye. DMF and DMA interact primarily through dipole-dipole forces, and the strong interactions with AOT sodium counterions destroy their bulk structure. FA also interacts with the Na+ counterions but retains its H-bond network present in bulk solvent. Surprisingly, FA appears to be the only polar solvent other than water forming a "polar-solvent pool" with macroscopic properties similar to the bulk.     

A molecular dynamics study of a nafion polyelectrolyte membrane and the aqueous phase structure for proton transport

A molecular dynamics simulation study of hydrated Nafion at water contents ranging from 5 to 20 wt % was performed to examine the structure and dynamics of the hydrated polyelectrolyte system. The simulations show that the system forms segregated hydrophobic regions consisting primarily of the polymer backbone and hydrophilic regions with an inhomogeneous water distribution. We find that the water clustering strongly depends on the water content. At low water content, only isolated small water clusters are formed. As the water content increases, it becomes increasingly possible that a predominant majority of water molecules form a single cluster, suggesting that the hydrophilic regions become connected. We characterize the atomic structures formed within the system by various atomic pair correlation functions. The water structure factor shows a peak at q values corresponding to an intercluster distance about 2.5 nm and greater. With increasing water content, the distance moves to larger values, consistent with findings from scattering experiments. We find that the degree of solvation of hydronium ions by water molecules is a strong function of water content. At 5 wt %, a majority of the hydronium ions are hydrated by no more than two water molecules, prohibiting structural diffusion. As water content increases, the hydronium ions continue to become increasingly hydrated, resulting in structures capable of forming eigen ions, a necessary step in structural diffusion. Addressing the experimentally observed fact that conductivity in these membranes abruptly drops near 5 wt %, we find that both the local structure of the poorly hydrated hydronium ions and the disconnected nature of the global morphology of the water nanonetwork at low water content should contribute to poor conductivity.     

Proton-transfer dynamics in the activation of cytochrome P450eryF

Molecular dynamics simulations are combined with quantum chemistry calculations of instantaneous proton-transfer energy profiles to investigate proton-transfer events in the transient pathway of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1), from the oxyferrous species to the catalytically active ferryl oxygen species (compound I). This reaction is one of the most fundamental unresolved aspects in the mechanism of oxidation that is common to all cytochrome P450s. We find that this process involves an ultrafast proton transfer from the crystallographic water molecule W519 to the distal oxygen bound to the heme group, and a subsequent proton-transfer event from W564 to W519. Both proton-transfer events are found to be endothermic in the oxyferrous state, suggesting that the oxyferrous reduction is mechanistically linked to the proton-transfer dynamics. These findings indicate that the hydrogen bond network, proximate to the O(2)-binding cleft, plays a crucial functional role in the enzymatic activation of P450s. Our results are consistent with the effect of mutations on the enzymatic efficacy.     

Crystal growth simulations of methane hydrates in the presence of silica surfaces

We present a molecular dynamics simulation study of the crystal growth of methane hydrates in the presence of model silica (SiO(2)) surfaces. The crystal growth under apparent steady-state conditions shows a clear preference for bulk solution. We observe rather disordered water arrangements very close to the silica surface within about 5 ? in both liquid and crystalline regions of the system. These disordered structures have dynamic and structural properties intermediate between those exhibited by molecules in bulk liquid and crystalline phases. The presence of methane molecules appears to help stabilize these structures. We observe that under appropriate conditions, the hydroxylated silica surfaces can serve as a source of methane molecules which can help promote hydrate growth near the surfaces.     

Experimental and computational evidence for an inversion in guest capacity in high-generation triazine dendrimer hosts

The synthesis, characterization, and host-guest chemistry of high-generation triazine dendrimers are described. With pyrene and camptothecin as guests, experiments revealed that the guest capacity of odd-generation triazine dendrimers increased until generation 7 but decreased at generation 9. Molecular dynamics simulations conducted in explicit solvent showed a useful fingerprint for this behavior in radial distribution functions of water molecules penetrating the interior of the dendrimers. A linear relationship between the guest capacity of dendrimers measured experimentally and the number of water molecules within the interior determined computationally was observed.     

The phenomenon of water penetration into sodium octanoate micelles studied by molecular dynamics computer simulation

 In this publication we have studied the penetration process of water molecules into the hydrophobic core of a sodium octanoate micelle. The analysis of this phenomenon was based on a molecular dynamics computer simulation. We calculated the probability to find water molecules within a specific sphere which was adjusted in the center of the micelle. It turned out that the position of the micellar mass and geometry center was not too different, so that this reference point was well characterized. Water penetration was observed within the whole aggregate but if the radius is smaller than 300 pm, polar solvent molecules are very rarely observed. The results of our computer simulations suggest that significant water diffusion into the micelle occurs at larger distances from the micellar center with a lower threshold value of about 400 pm. In addition to these calculations, we used the Connolly algorithm in order to determine the solvent accessible surfaces of different micellar equilibrium structures. We observed large dynamical fluctuations with the formation of pores and channels. These structures are occasionally filled with water molecules. Received: 29 April 1998 Accepted: 27 May 1998