Investigating hydration dependence of dynamics of confined water: monolayer, hydration water and Maxwell-Wagner processes

The dynamics of water confined in silica matrices MCM-41 C10 and C18, with pore diameter of 21 and 36 A, respectively, is examined by broadband dielectric spectroscopy (10(-2)-10(9) Hz) and differential scanning calorimetry for a wide temperature interval (110-340 K). The dynamics from capillary condensed hydration water and surface monolayer of water are separated in the analysis. Contrary to previous reports, the rotational dynamics are shown to be virtually independent on the hydration level and pore size. Moreover, a third process, also reported for other systems, and exhibiting a saddlelike temperature dependence is investigated. We argue that this process is due to a Maxwell-Wagner process and not to strongly bound surface water as previously suggested in the literature. The dynamics of this process is strongly dependent on the amount of hydration water in the pores. The anomalous temperature dependence can then easily be explained by a loss of hydration water at high temperatures in contradiction to previous explanations.     

Hofmeister anionic effects on hydration electric fields around water and peptide

Specific ion effects on water dynamics and local solvation structure around a peptide are important in understanding the Hofmeister series of ions and their effects on protein stability in aqueous solution. Water dynamics is essentially governed by local hydrogen-bonding interactions with surrounding water molecules producing hydration electric field on each water molecule. Here, we show that the hydration electric field on the OD bond of HOD molecule in water can be directly estimated by measuring its OD stretch infrared (IR) radiation frequency shift upon increasing ion concentration. For a variety of electrolyte solutions containing Hofmeister anions, we measured the OD stretch IR bands and estimated the hydration electric field on the OD bond to be about a hundred MV∕cm with standard deviation of tens of MV∕cm. As anion concentration increases from 1 to 6 M, the hydration electric field on the OD bond decreases by about 10%, indicating that the local H-bond network is partially broken by dissolved ions. However, the measured hydration electric fields on the OD bond and its fluctuation amplitudes for varying anions are rather independent on whether the anion is a kosmotrope or a chaotrope. To further examine the Hofmeister effects on H-bond solvation structure around a peptide bond, we examined the amide I' and II' mode frequencies of N-methylacetamide in various electrolyte D(2)O solutions. It is found that the two amide vibrational frequencies are not affected by ions, indicating that the H-bond solvation structure in the vicinity of a peptide remains the same irrespective of the concentration and character of ions. The present experimental results suggest that the Hofmeister anionic effects are not caused by direct electrostatic interactions of ions with peptide bond or water molecules in its first solvation shell. Furthermore, even though the H-bond network of water is affected by ions, thus induced change of local hydration electric field on the OD bond of HOD is not in good correlation with the well-known Hofmeister series. We anticipate that the present experimental results provide an important clue about the Hofmeister effect on protein structure and present a discussion on possible alternative mechanisms.     

Molecular modeling and simulation of water near model micelles: diffusion, rotational relaxation and structure at the hydration interface

This paper reports on molecular dynamics simulations of two hydrated micelles composed of C12E6 and LDAO surfactants. The simulations results provide a quantitative picture of the dynamics of the hydration water at the water/micelle interface. Both the residence time of water near the micelle surface and its retardation with respect to the bulk have been estimated. It is found that the water dynamics is radically different for the two micellar systems and depends on the physical nature of the micelle surface in contact with water. For C12E6 this interface is thicker and presents a stronger hydrophilic character than that of LDAO. Thus, in C12E6, surface water dynamics is 1-2 orders of magnitude slower than that of bulk water, compared with only 18% for the LDAO system. The simulations have also revealed the nature of the rotational landscape experienced by water at the micellar surface: In the C12E6 micelle water rotation occurs in a highly anisotropic space due to confinement of waters at the interface; in LDAO the rotational landscape is instead isotropic. These findings clearly indicate that the slowdown of interfacial water relaxation near complex micelles depends, case by case, on the structural properties of the interface itself, such as the ratio between hydrophobic/hydrophilic exposed regions and on the interface thickness and topography.     

Molecular dynamic simulation of the hydration and diffusion of chloride ions from bulk water to polypyrrole matrix

A molecular dynamic simulation of wet polypyrrole film was carried out, in both oxidized and reduced state. The system was modeled by two layers of polypyrrole, water and chloride ions (as counterions required for charge balance in the oxidized state) in atomic detail to provide an insight into some dynamic and steady properties of the system. Our simulations pointed to a swelling of the polymer matrix after oxidation due to electrostatic repulsions between charged sites of the oxidized polypyrrole, followed by penetration of the polypyrrole by counterions to maintain the electroneutrality of the system. Associated with this penetration of counterions toward the core of the oxidized polypyrrole, dehydration of the counterions was observed. This dehydration was compensated (in part) by a strong coordination with the charged sites of the polymer. The remaining hydrophobicity inside the polymer also contributed to the dehydration of these counterions. The translational diffusion coefficient of chloride ions was also calculated at different positions of the polypyrrole/water interface, from bulk water to the inner polymer matrix. A value of 4.1 x 10(-5) cm(2) s(-1) was measured in the bulk water compared to 5 x 10(-7) cm(2) s(-1) inside the polymer, representing a diminution of two orders of magnitude for the translational diffusion coefficient from bulk water to the core of a oxidized polypyrrole matrix. These results were in good agreement with experimental data.     

Molecular interpretation of electrokinetic potentials

A discussion is given of the electrokinetic, or ζ-potential in terms of the slip process and the composition of the electric double layer. Electrokinetically, only the outer parts of double layers are active. The existence of a stagnant part is generally observed for aqueous solutions adjacent to solid surfaces. It is claimed that this stagnancy is caused by the spontaneous structuring of fluids near solid surfaces. Hence, it is a ubiquitous phenomenon and the thickness of the stagnant layer does not significantly depend on the wettability and the surface charge of the surfaces.     

Molecular dynamics simulations of peptide carboxylate hydration

Aqueous solvation of carboxylate groups, as present in the glycine zwitterion and the dipeptide aspartylalanine, is studied employing a force-field that includes distributed multipole electrostatics and induction contributions (Amoebapro: P. Ren and J. W. Ponder, J. Comput. Chem., 2002, 23, 1497; P. Ren and J. W. Ponder, J. Phys. Chem. B, 2003, 107, 5933; J. W. Ponder and D. A. Case, Adv. Protein Chem., 2003, 66, 27). Radial and orientation distribution functions, as well as hydration numbers, are calculated and compared with existing simulation data derived from Car-Parrinello molecular dynamics (CPMD), and also distributed-charge force-fields. Connections are also made with experimental data for solvation of carboxylates in water. Our findings show that Amoebapro yields carboxylate solvation properties in very good agreement with CPMD results, significantly closer agreement than can be obtained from traditional force-fields. We also demonstrate that the influence of solvation on the conformation of the dipeptide is markedly different using Amoebapro compared with the other force-fields.     

Ion-specific hydration effects: Extending the Poisson-Boltzmann theory

In aqueous solutions, dissolved ions interact strongly with the surrounding water and surfaces, thereby modifying solution properties in an ion-specific manner. These ion-hydration interactions can be accounted for theoretically on a mean-field level by including phenomenological terms in the free energy that correspond to the most dominant ion-specific interactions. Minimizing this free energy leads to modified Poisson-Boltzmann equations with appropriate boundary conditions. Here, we review how this strategy has been used to predict some of the ways ion-specific effects can modify the forces acting within and between charged interfaces immersed in salt solutions.     

Kinetically stable complexes in water: the role of hydration and hydrophobicity

We describe here the synthesis and characterization of a molecular receptor that forms kinetically and thermodynamically stable host-guest complexes in water. This cavitand-based host is composed of a preorganized aromatic pocket whose rim is decorated with four negatively charged carboxylate groups. (1)H NMR and isothermal titration calorimetry have been used to characterize the behavior of the resulting complexes in response to changes in pH, buffer identity, and salt concentration and in the presence of sodium dodecyl sulfate micelles.     

Aqueous solutions of divalent chlorides: ions hydration shell and water structure

By combining neutron diffraction and Monte Carlo simulations, we have determined the microscopic structure of the hydration ions shell in aqueous solutions of MgCl(2) and CaCl(2), along with the radial distribution functions of the solvent. In particular the hydration shell of the cations, show cation specific symmetry, due to the strong and directional interaction of ions and water oxygens. The ions and their hydration shells likely form molecular moieties and bring clear signatures in the water-water radial distribution functions. Apart from these signatures, the influence of divalent salts on the microscopic structure of water is similar to that of previously investigated monovalent solutes, and it is visible as a shift of the second peak of the oxygen-oxygen radial distribution function, caused by distortion of the hydrogen bond network of water.     

Ab initio rigid water: effect on water structure, ion hydration, and thermodynamics

We investigate the liquid structure, ion hydration, and some thermodynamic properties associated with the rigid geometry approximation to water by applying ab initio molecular dynamics simulations (AIMD) with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional at T = 320 K. We vary the rigid water geometry in order to locate a class of practical water models that yield reasonable liquid structure and dynamics, and to examine the progression of AIMD-predicted water behavior as the OH bond length varies. Water constrained at the optimal PBE gas phase geometry yields reasonable pair correlation functions. The predicted liquid phase pressure, however, is large ( approximately 8.0 kbar). Although the O-H bond in water should elongate when transferred from gas to the condensed phase, when it is constrained to 0.02, or even just 0.01 A longer than the optimal gas phase value, liquid water is predicted to be substantially overstructured compared to experiments. Zero temperature calculations of the thermodynamic properties of cubic ice underscore the sensitivity toward small variations in the O-H bond length. We examine the hydration structures of potassium, chloride, and formate ions in one rigid PBE water model. The results are in reasonable agreement with unconstrained AIMD simulations.     

Synthesis and mesomorphic properties of esters derived from alkanediols

New alkane-α,ω-diyl esters, propane-1,3-diyl and butane-1,4-diyl bis[4-(4-alkoxybenzoyloxy)-3-bromobenzoates], were synthesized and were found to form nematic mesophase. The effects of the lateral substituent and the length of the central and terminal hydrocarbon chains on the mesomorphic properties of butane-1,4-diyl bis[4-(4-alkoxybenzoyloxy)-3-bromobenzoates] were studied. Structural factors responsible for mesogenic properties of these compounds were determined by conformation analysis.     

Hierarchical supramolecular structuring and dynamical properties of water soluble polyethylene glycol-perylene self-assemblies

The structural and dynamical properties of dilute aqueous solutions of poly(ethylene glycol)-perylene diimides (PEG(n)-PDI) have been investigated by means of static and dynamic light scattering, TEM microscopy, and small-angle X-ray scattering experiments. The amphiphilic PEG(n)-PDI molecules first self-assemble into stable and compact primary stacks of a few units of planar PDI through hydrophobic and π-π interactions. These primary stacks subsequently arrange in large and globular aggregates of typically 100-250 nm via weak PEG chain interpenetration. Surprisingly, the scattered electric field autocorrelation function g((1))(q,t) measured by dynamic light scattering evolves over very long periods of times (several months) and up to a bimodal distribution. The fast relaxation mechanism is associated to the diffusion of free primary stacks, whereas the slower relaxation still indicates the presence of large self-assemblies. Kinetic experiments show that the large supramolecular aggregates slowly release the free primary stacks whose proportion increases with time. This dissociation depends on several parameters such as PEG side chain length, total concentration, and shaking.     

Molecular simulation of the hydration Gibbs energy of barbiturates

In the present work, molecular dynamics calculations of the Gibbs energy of hydration of 10 different substituted barbiturates in SPC/E water were performed using thermodynamic integration. Given that experimental determination of the Gibbs hydration energy for this class of compounds is currently unfeasible, computer simulations appear as the only alternative for the estimation of this important quantity. Several simulation parameters are discussed and optimized based on calculations for barbituric acid. It is concluded that accounting for electrostatic interactions with the Reaction-Field method can be up to two times faster than with Particle-Mesh-Ewald method, without loss of accuracy. Different number of solvent molecules and simulation lengths were also tested. Lennard-Jones and electrostatic contributions were scaled down to zero in an independent way. It is shown that the electrostatic contribution is dominant (representing approximately 90% of the total Gibbs energy of hydration) and that barbiturate intra-molecular interactions cannot be neglected. The importance of the electrostatic contribution is attributed to the formation of hydrogen bonds between the barbiturates and water, which play an important role in the solvation process. The influence of the different substituents and their contribution to the Gibbs energy of hydration was assessed. Finally, the Lennard-Jones contributions and the total hydration Gibbs energy can both be correlated against molecular weight or partition coefficient data for mono- and di-substituted barbiturates.     

Molecular dynamics study of the hydration of lanthanum(III) and europium(III) including many-body effects

Lanthanides complexes are widely used as contrast agents in magnetic resonance imaging (MRI) and are involved in many fields such as organic synthesis, catalysis, and nuclear waste management. The complexation of the ion by the solvent or an organic ligand and the resulting properties (for example the relaxivity in MRI) are mainly governed by the structure and dynamics of the coordination shells. All of the MD approaches already carried out for the lanthanide(III) hydration failed due to the lack of accurate representation of many-body effects. We present the first molecular dynamics simulation including these effects that accounts for the experimental results from a structural and dynamic (water exchange rate) point of view.     

Structuring effects and hydration phenomena in poly(ethylene glycol)/water mixtures investigated by brillouin scattering

Aqueous solutions of poly(ethylene glycol) (PEG) of mean molecular mass of 600 g/mol (PEG600) are investigated by Brillouin scattering technique. At high PEG content, a relaxation phenomenon is observed, which is related to a local rearrangement of the polymer structure where the interaction, via hydrogen bonding, with the solvent molecules plays a role. The obtained values of the relaxation times match the literature data very well for a fast relaxation time revealed by dielectric relaxation measurements in very similar mixtures. The calculated concentration behaviors of the excess adiabatic compressibility turns out in good agreement with the previous findings from ultrasonic measurements at 3 MHz. The observed minimum in the adiabatic compressibility is interpreted as the result of the interaction between water and the EO units of the PEG chain, which results in a structure tighter then that typical of bulk water and of pure PEG600. Such a hypothesis is supported by the observation that volume fraction value of about 0.3 coincides with the concentration value at which full hydration of EO units takes place. The observation that at the same concentration, the polymer coils start to overlap each other further supports the idea that the adiabatic compressibility behavior is monitoring the structural evolution of the mixture. However, similar results are obtained for largely different binary mixture which suggests caution in taking this conclusion too literally. In particular, the hypothesis that the occurrence of an extreme in the excess adiabatic compressibility could be simply originated by statistical effects and that further work is required for disentangling entropic contribution from effects of hetero-association and self-aggregation of one or both the components.     

Molecular dynamics simulations of hydration, dissolution and nucleation processes at the alpha-quartz (0001) surface in liquid water

Molecular dynamics simulations have been employed to investigate the hydration and dissolution of alpha-quartz (0001) surfaces in a liquid water environment. Our study indicates that the structure of the water layers near the surfaces is affected by the nature of the substrate surface and by temperature. Ordered mono-layers of interfacial water molecules form in the region of the substrate where the surface is highly charged and built up of Si-O-Si bridges. As the temperature is increased this ordered mono-layer structure is gradually lost. When the surface is terminated by silanol groups, the water retains liquid-like properties even at low temperature and the molecules are distributed in a random manner, without the formation of distinct ordered mono-layers of water molecules near the surface. Taking into account the entropy of the system, the calculated energies of stepwise dissolution of a silicon species from the surface suggest that on thermodynamic grounds the complete dissolution of silicon atoms from the quartz surfaces in a liquid water environment is an endothermic process, but that the formation of a -Si(OH)3 species at the surface would be possible. In addition, if the Si(OH)(4) species were to be dissolved, it would remain near the surface, and re-deposition at the defect-free surface is thermodynamically preferred, although there is an activation enthalpy to the first step in the process of nucleation of Si(OH)4 at the perfect surface.     

Hydrogen-bonded network of hydration water around model solutes

Clustering of water molecules in the hydration shells of spherical structureless solutes was studied in dependence on thermodynamic state, solute radius R(sp) and strength U(0) of water-solute interaction. Two qualitatively different clustering states of hydration water have been found: an "ordered" state with a hydrogen-bonded (H-bonded) network, which includes most of the hydration water, and a "disordered" state with small H-bonded clusters of hydration water. The transition from the ordered to disordered state occurs upon increasing temperature and decreasing pressure. This percolation transition is rounded due to the finite solute size and occurs in some temperature (pressure) interval. A finite-size scaling was applied to determine the transition temperature T(∞) in the limit R(sp)→∞. Strengthening of the water-solute interaction strongly enhances the stability of the ordered state: the transition temperature increases by about 35 °C, when U(0) decreases by 1 kcal mol(-1). At T > T(∞) and fixed U(0), the stability of the H-bonded water network increases upon decreasing solute size.     

A cellular automata model of water structuring by a chiral solute

The organization of water around a solute molecule with surface features of varying hydropathic states is studied. A stationary solute molecule and mobile water solvent molecules are modeled using cellular automata dynamics. It is shown that varying hydropathic states of solute molecule surface features influence the relative affinities of water for these features. In the case of a simulated chiral solute, a chiral pattern of associated water molecules binding to the surface is produced. This finding is in agreement with published simulations and circular dichroism measurements. A pattern of water molecules at locations beyond the surface of the solute molecules is detected, evidence of an emergent property in this solvent-solute system.