A molecular dynamics study of the dielectric properties of aqueous solutions of alanine and alanine dipeptide

Molecular dynamics simulations were used to compute the frequency-dependent dielectric susceptibility of aqueous solutions of alanine and alanine dipeptide. We studied four alanine solutions, ranging in concentration from 0.13-0.55 mol/liter, and two solutions of alanine dipeptide (0.13 and 0.27 mol/liter). In accord with experiment we find a strong dielectric increment for both solutes, whose molecular origin is shown to be the zwitterionic nature of the solutes. The dynamic properties were analyzed based on a dielectric component analysis into solute, a first hydration shell, and all remaining (bulk) waters. The results of this three component decomposition were interpreted directly, as well as by uniting the solute and hydration shell component to a "suprasolute" component. In both approaches three contributions to the frequency-dependent dielectric properties can be discerned. The quantitatively largest and fastest component arises from bulk water [i.e., water not influenced by the solute(s)]. The interaction between waters surrounding the solute(s) (the hydration shell) and bulk water molecules leads to a relaxation process occurring on an intermediate time scale. The slowest relaxation process originates from the solute(s) and the interaction of the solute(s) with the first hydration shell and bulk water. The primary importance of the hydration shell is the exchange of shell and bulk waters; the self-contribution from bound water molecules is comparatively small. While in the alanine solutions the solute-water cross-terms are more important than the solute self-term, the solute contribution is larger in the dipeptide solutions. In the latter systems a much clearer separation of time scales between water and alanine dipeptide related properties is observed. The similarities and differences of the dielectric properties of the amino acid/peptide solutions studied in this work and of solutions of mono- and disaccharides and of the protein ubiquitin are discussed.     

Dielectric relaxation of aqueous solutions of hydrophilic versus amphiphilic peptides

We report on molecular dynamics simulations of the frequency-dependent dielectric relaxation spectra at room temperature for aqueous solutions of a hydrophilic peptide and an amphiphilic peptide at two concentrations. We find that only the high-concentration amphiphilic peptide solution exhibits an anomalous dielectric increment over that of pure water, while the hydrophilic peptide exhibits a significant dielectric decrement. The dielectric component analysis carried out by decomposing these peptide solutions into peptide, hydration layer, and outer layer(s) of water clearly shows the presence of a unique dipolar component with a relaxation time scale on the order of approximately 25 ps (compared to the bulk water time scale of approximately 11 ps) that originates from the interaction between the hydration layer water and the outer layer(s) of water. Results obtained from the dielectric component analysis further show the emergence of a distinct and much lower frequency relaxation process for the high-concentration amphiphilic peptide compared to the hydrophilic peptide due to strong peptide dipolar couplings to all constituents, accompanied by a slowing of the structural relaxation in all water layers, giving rise to time scales close to approximately 1 ns. We suggest that the molecular origin of the dielectric relaxation anomalies is due to frustration in the water network arising from the amphiphilic chemistry of the peptide that does not allow it to reorient on the picosecond time scale of bulk water motions. This explanation is consistent with the idea of the "slaving" of residue side chain motions to protein surface water, and furthermore offers the possibility that the anomalous dynamics observed from a number of spectroscopies arises at the interface of hydrophobic and hydrophilic domains on the protein surface.     

Beyond the first hydration shell: dielectric study of LiCl aqueous solution

The dielectric constant for lithium chloride (LiCl) in aqueous solution with the entire concentration has been determined in the frequency range 0.5 GHz–50GHz at 298 K by dielectric relaxation spectroscopy (DRS). The system behaviour is described according to the Cole–Cole and two Debye-type relaxation functions whose evolution with composition is analysed. Combining the Cole–Cole and two Debye-type relaxation functions, the results of the dielectric properties of aqueous LiCl solution are presented and discussed. The maximum number of water molecules perturbed by ions in the hydration shell decays with its concentration. Using the extended Froehlich theory, it is concluded that the water structure is perturbed by ions beyond the first hydration shell in LiCl aqueous solution system.     

Hydration and dynamic behavior of cyclodextrins in aqueous solution

The hydration state and dynamics of plain and chemically modified cyclodextrins (CDs) in aqueous solution were investigated by using dielectric relaxation measurements at 25 degrees C over a wide frequency range up to 20 GHz, which is the relaxation frequency of pure liquid water molecules. The obtained dielectric relaxation spectra were decomposed into two major and one minor relaxation modes with relaxation times of approximately 8.3, 20-25, and 1000-2500 ps, respectively, depending on the CD species. The two major modes, fast and medium, were attributed to a rotational relaxation process of water molecules belonging to the bulk (free) state and an exchange of water molecules hydrated to CDs owing to hydrogen bond formation. The hydration numbers of the CDs strongly depend on the number of hydroxy (OH) groups controlled by chemical modification such as methylation. Increasing the number of methoxy or 2-hydroxypropoxy groups increases the hydration number of CD molecules, and results in higher solubilities of the chemically modified CDs than those of the plain CDs. The minor, slow mode was assigned to overall rotational relaxation for CDs with finite permanent dipole moments, which also depends on the number of OH groups.     

Solute rotational dynamics at the water liquid/vapor interface

The rotational dynamics of a number of diatomic molecules adsorbed at different locations at the interface between water and its own vapors are studied using classical molecular dynamics computer simulations. Both equilibrium orientational and energy correlations and nonequilibrium orientational and energy relaxation correlations are calculated. By varying the dipole moment of the molecule and its location, and by comparing the results with those in bulk water, the effects of dielectric and mechanical frictions on reorientation dynamics and on rotational energy relaxation can be studied. It is shown that for nonpolar and weekly polar solutes, the equilibrium orientational relaxation is much slower in the bulk than at the interface. As the solute becomes more polar, the rotation slows down and the surface and bulk dynamics become similar. The energy relaxation (both equilibrium and nonequilibrium) has the opposite trend with the solute dipole (larger dipoles relax faster), but here again the bulk and surface results converge as the solute dipole is increased. It is shown that these behaviors correlate with the peak value of the solvent-solute radial distribution function, which demonstrates the importance of the first hydration shell structure in determining the rotational dynamics and dependence of these dynamics on the solute dipole and location.     

Hydration Dynamics of Water near an Amphiphilic Model Peptide at Low Hydration Levels: A Dielectric Relaxation Study

A dielectric relaxation study of aqueous solutions of the amphiphilic model peptide N‐acetyl‐leucine amide (NALA) at 298 K over a wide range of hydration levels is presented. The experiments range from states where water builds up several hydration layers to states where single water molecules or small water clusters are shared by several NALA molecules. The dielectric spectra reveal two modes on the 10 and 100 ps timescales. These are largely broadened with regard to the Lorentzian shape caused by simple Debye‐type relaxation, and are well described by the Kohlrausch–Williams–Watts stretched exponential function. The fast mode is assigned to water reorientation comprising bulk water as well as hydration water. Even when all water molecules are in contact with the solute, this fast component is dominant, and its mean relaxation time is retarded by less than a factor of two relative to neat water. The amplitude of the slow process is far higher than expected for the dipolar reorientation of the solute. The observations are consistent with results from molecular dynamics simulations for a similar model peptide reported in the literature. They suggest that the slow relaxation mode is mainly founded in peptide–water dipolar couplings, with some additional contribution from slowly reorienting hydration water molecules. The results are discussed with regard to the hydration dynamics of proteins and the interpretation of dielectric spectra of protein solutions.     

Hydration of poly(ethylene oxide)s in aqueous solution as studied by dielectric relaxation measurements

The hydration state of poly(ethylene oxide)s (PEOs) in aqueous solutions was investigated using dielectric relaxation measurements at 25 degrees C over a frequency range up to 20 GHz, which is the relaxation frequency of water molecules in a bulk state. The dielectric relaxation spectra obtained indicated decomposition into two major and one minor relaxation modes with relaxation times of 8.3, 22, and 250 ps, respectively. The two major modes were attributed to rotational relaxation of water molecules belonging to the bulk state and water molecules hydrogen bonded to ethylene oxide (EO) monomer units. The number of hydration water molecules per EO unit depended on the molar mass of PEO (M) and reached a constant value of 3.7 at M > 1500, which agrees with the value obtained by other experiments.     

Hydrophobic hydration in an orientational lattice model

To shed light on the microscopic mechanism of hydrophobic hydration, we study a simplified lattice model for water solutions in which the orientational nature of hydrogen bonding as well as the degeneracy related to proton distribution are taken into account. Miscibility properties of the model are looked at for both polar (hydrogen bonding) and nonpolar (non-hydrogen bonding) solutes. A quasichemical solution for the pure system is reviewed and extended to include the different kinds of solute. A Monte Carlo study of our model yields a novel feature for the local structure of the hydration layer: energy correlation relaxation times for solvation water are larger than the corresponding relaxation times for bulk water. This result suggests the presence of ordering of water particles in the first hydration shell. A nonassociating model solvent, represented by a lattice gas, presents opposite behavior, indicating that this effect is a result of the directionality of the interaction. In presence of polar solutes, we find an ordered mixed pseudophase at low temperatures, indicating the possibility of closed loops of immiscibility.     

Amino acid side chain interactions in the presence of salts

The effects of salt on the intermolecular interactions between polar/charged amino acids are investigated through molecular dynamics simulations. The mean forces and associated potentials are calculated for NaCl salt in the 0-2 M concentration range at 298 K. It is found that the addition of salt may stabilize or destabilize the interactions, depending on the nature of the interacting molecules. The degree of (de)stabilization is quantified, and the origin of the salt-dependent modulation is discussed based upon an analysis of solvent density profiles. To gain insight into the molecular origin of the salt modulation, spatial distribution functions (sdf's) are calculated, revealing a high degree of solvent structuredness in all cases. The peaks in the sdf's are consistent with long-range hydrogen-bonding networks connecting the solute hydrophilic groups, and that contribute to their intermolecular solvent-induced forces. The restructuring of water around the solutes as they dissociate from close contact is analyzed. This analysis offers clues on how the solvent structure modulates the effective intermolecular interactions in complex solutes. This modulation results from a critical balance between bulk electrostatic forces and those exerted by (i) the water molecules in the structured region between the monomers, which is disrupted by ions that transiently enter the hydration shells, and (ii) the ions in the hydration shells in direct interactions with the solutes. The implications of these findings in protein/ligand (noncovalent) association/dissociation mechanisms are briefly discussed.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

Increased interfacial thickness of the NaF, NaCl and NaBr salt aqueous solutions probed with non-resonant surface second harmonic generation (SHG)

Specific ion effects on the nonlinear optical response from the water molecules at the air/sodium halide solution interfaces are measured using non-resonant surface second harmonic generation (SHG). Procedures have been developed to monitor and remove the impurities in the salt solution samples to ensure measurement of small changes in the SHG signal. Quantitative polarization analysis of the measured SHG data indicated that the average orientation of the interfacial water molecules changed only slightly around 40 degrees with the increase of the bulk concentration of the three sodium halides, namely NaF, NaCl and NaBr, from that of the neat air/water interface. The observed significant SHG signal increase with the bulk salt concentration is attributed to the overall increase of the thickness of the interfacial water molecular layer, following the order of NaBr > NaCl approximately NaF. The absence of the electric-field-induced SHG (EFISHG) effect indicated that the electric double layer at the salt aqueous solution interface is much weaker than that predicted from the molecular dynamics (MD) simulations. These results provided quantitative data to the specific anion effects on the interfacial water molecules of the electrolyte aqueous solution, not only for the larger and more polarizable Br(-) anion, but also for the smaller and less polarizable F(-) and Cl(-) anions.     

Influence of hydration on protein dynamics: combining dielectric and neutron scattering spectroscopy data

Combining dielectric spectroscopy and neutron scattering data for hydrated lysozyme powders, we were able to identify several relaxation processes and follow protein dynamics at different hydration levels over a broad frequency and temperature range. We ascribe the main dielectric process to protein's structural relaxation coupled to hydration water and the slowest dielectric process to a larger scale protein's motions. Both relaxations exhibit a smooth, slightly super-Arrhenius temperature dependence between 300 and 180 K. The temperature dependence of the slowest process follows the main dielectric relaxation, emphasizing that the same friction mechanism might control both processes. No signs of a proposed sharp fragile-to-strong crossover at T approximately 220 K are observed in temperature dependences of these processes. Both processes show strong dependence on hydration: the main dielectric process slows down by six orders with a decrease in hydration from h approximately 0.37 (grams of water per grams of protein) to h approximately 0.05. The slowest process shows even stronger dependence on hydration. The third (fastest) dielectric relaxation process has been detected only in samples with high hydration ( h approximately 0.3 and higher). We ascribe it to a secondary relaxation of hydration water. The mechanism of the protein dynamic transition and a general picture of the protein dynamics are discussed.     

Vibrational and orientational dynamics of water in aqueous hydroxide solutions

We report the vibrational and orientational dynamics of water molecules in isotopically diluted NaOH and NaOD solutions using polarization-resolved femtosecond vibrational spectroscopy and terahertz time-domain dielectric relaxation measurements. We observe a speed-up of the vibrational relaxation of the O-D stretching vibration of HDO molecules outside the first hydration shell of OH(-) from 1.7 ± 0.2 ps for neat water to 1.0 ± 0.2 ps for a solution of 5 M NaOH in HDO:H(2)O. For the O-H vibration of HDO molecules outside the first hydration shell of OD(-), we observe a similar speed-up from 750 ± 50 fs to 600 ± 50 fs for a solution of 6 M NaOD in HDO:D(2)O. The acceleration of the decay is assigned to fluctuations in the energy levels of the HDO molecules due to charge transfer events and charge fluctuations. The reorientation dynamics of water molecules outside the first hydration shell are observed to show the same time constant of 2.5 ± 0.2 ps as in bulk liquid water, indicating that there is no long range effect of the hydroxide ion on the hydrogen-bond structure of liquid water. The terahertz dielectric relaxation experiments show that the transfer of the hydroxide ion through liquid water involves the simultaneous motion of ~7 surrounding water molecules, considerably less than previously reported for the proton.     

Thermodynamic origin of hofmeister ion effects

Quantitative interpretation and prediction of Hofmeister ion effects on protein processes, including folding and crystallization, have been elusive goals of a century of research. Here, a quantitative thermodynamic analysis, developed to treat noncoulombic interactions of solutes with biopolymer surface and recently extended to analyze the effects of Hofmeister salts on the surface tension of water, is applied to literature solubility data for small hydrocarbons and model peptides. This analysis allows us to obtain a minimum estimate of the hydration b1 (H2O A(-2)), of hydrocarbon surface and partition coefficients Kp, characterizing the distribution of salts and salt ions between this hydration water and bulk water. Assuming that Na+ and SO4(2-) ions of Na2SO4 (the salt giving the largest reduction in hydrocarbon solubility as well as the largest increase in surface tension) are fully excluded from the hydration water at hydrocarbon surface, we obtain the same b1 as for air-water surface (approximately 0.18 H2O A(-2)). Rank orders of cation and anion partition coefficients for nonpolar surface follow the Hofmeister series for protein processes, but are strongly offset for cations in the direction of exclusion (preferential hydration). By applying a coarse-grained decomposition of water accessible surface area (ASA) into nonpolar, polar amide, and other polar surface and the same hydration b1 to interpret peptide solubility increments, we determine salt partition coefficients for amide surface. These partition coefficients are separated into single-ion contributions based on the observation that both Cl- and Na+ (also K+) occupy neutral positions in the middle of the anion and cation Hofmeister series for protein folding. Independent of this assignment, we find that all cations investigated are strongly accumulated at amide surface while most anions are excluded. Cation and anion effects are independent and additive, allowing successful prediction of Hofmeister salt effects on micelle formation and other processes from structural information (ASA).     

How ions affect the structure of water

We model ion solvation in water. We use the MB model of water, a simple two-dimensional statistical mechanical model in which waters are represented as Lennard-Jones disks having Gaussian hydrogen-bonding arms. We introduce a charge dipole into MB waters. We perform (NPT) Monte Carlo simulations to explore how water molecules are organized around ions and around nonpolar solutes in salt solutions. The model gives good qualitative agreement with experiments, including Jones-Dole viscosity B coefficients, Samoilov and Hirata ion hydration activation energies, ion solvation thermodynamics, and Setschenow coefficients for Hofmeister series ions, which describe the salt concentration dependence of the solubilities of hydrophobic solutes. The two main ideas captured here are (1) that charge densities govern the interactions of ions with water, and (2) that a balance of forces determines water structure: electrostatics (water's dipole interacting with ions) and hydrogen bonding (water interacting with neighboring waters). Small ions (kosmotropes) have high charge densities so they cause strong electrostatic ordering of nearby waters, breaking hydrogen bonds. In contrast, large ions (chaotropes) have low charge densities, and surrounding water molecules are largely hydrogen bonded.     

Influence of concentration and temperature on the dynamics of water in the hydrophobic hydration shell of tetramethylurea

We study the influence of the amphipilic compound tetramethylurea (TMU) on the dynamical properties of water, using dielectric relaxation spectroscopy in the regime between 0.2 GHz and 2 THz. This technique is capable of resolving different water species, their relative fractions, and their corresponding reorientation dynamics. We find that the reorientation dynamics of water molecules in the hydration shell of the hydrophobic groups of TMU is between 3 (at low concentrations) and 10 (at higher concentrations) times slower than the dynamics of bulk water. The data indicate that the effect of hydrophobic groups on water is strong but relatively short-ranged. With increasing temperature, the fraction of water contained in the hydrophobic hydration shell decreases, which implies that the overall effect of hydrophobic groups on water becomes smaller.     

Impact of Hydrophilic Surfaces on Interfacial Water Dynamics Probed with NMR Spectroscopy

In suspensions of Nafion beads and of cationic gel beads, NMR spectroscopy showed two water-proton resonances, one representing intimate water layers next to the polymer surface, the other corresponding to water lying beyond. Both resonances show notably shorter spin-lattice relaxation times (T1) and smaller self-diffusion coefficients (D) indicating slower dynamics than bulk water. These findings confirm the existence of highly restricted water layers adsorbed onto hydrophilic surfaces and dynamically stable water beyond the first hydration layers. Thus, aqueous regions on the order of micrometers are dynamically different from bulk water.     

Thermodynamic linkage of ion binding and hydration in myofibrillar protein solutions determined from multinuclear spin relaxation studies

The mechanisms for the anionic and cationic interactions with myofibrillar proteins in aqueous solutions were investigated by nuclear magnetic resonance over a wide range of salt concentration. Markedly nonlinear dependeces of the ¹⁷O and ²³Na NMR transverse relaxation rates on salt concentration were analyzed with a thermodynamic linkage model of salt-dependent solubility and hydration (ligand-induced association model), according to Wyman's theory of linked functions. Nonlinear regression analysis of both ¹⁷O and ²³Na NMR data suggested cooperative, reversible binding of hydrated ions to myofibrillar proteins. Both ions and water were found to exchange fast, on the NMR timescale, between the binding sites of the myofibrillar proteins and the aqueous solution. At sodium chloride concentrations higher than about 0.1 grams salt/gram water, ion activities have marked effects upon the NMR relaxation rates of both ions and water. A salt activity model allowed quantitative fitting of the NMR data at high salt concentrations. The effect of neglecting the ion activity in solutions of myofibrillar proteins was also estimated and compared with the ligand-induced, cooperative association model for myofibrillar proteins. The comparison between the ¹⁷O and ²³Na results strongly suggests that water is exchanged as the hydrated ion species between the myofibrillar protein binding sites and the bulk, aqueous solution.     

Dielectric relaxation behavior of colloidal suspensions of palladium nanoparticle chains dispersed in PVP/EG solution

Dielectric measurements were carried out on colloidal suspensions of palladium nanoparticle chains dispersed in poly(vinyl pyrrolidone)/ethylene glycol (PVP/EG) solution with different particle volume fractions, and dielectric relaxation with relaxation time distribution and small relaxation amplitude was observed in the frequency range from 10(5) to 10(7) Hz. By means of the method based on logarithmic derivative of the dielectric constant and a numerical Kramers-Kronig transform method, two dielectric relaxations were confirmed and dielectric parameters were determined from the dielectric spectra. The dielectric parameters showed a strong dependence on the volume fraction of palladium nanoparticle chain. Through analyzing limiting conductivity at low frequency, the authors found the conductance percolation phenomenon of the suspensions, and the threshold volume fraction is about 0.18. It was concluded from analyzing the dielectric parameters that the high frequency dielectric relaxation results from interfacial polarization and the low frequency dielectric relaxation is a consequence of counterion polarization. They also found that the dispersion state of the palladium nanoparticle chain in PVP/EG solution is dependent on the particle volume fraction, and this may shed some light on a better application of this kind of materials.     

Cation electric field is related to hydration energy

Electronic structure calculations on ions that use dielectric continuum theory to mimic solvent around the bare ionic solutes are often prone to make large errors in the hydration energies. It is found for cations that much of the error can be accounted for by a simple linear correlation with the maximum value found anywhere on the dielectric cavity surface of the solute potential or, even better, the outgoing normal electric field, thus mirroring analogous results previously obtained for anions. This correlation allows for significantly improved estimates of cation hydration energies while still retaining cavities of physically reasonable size in determination of the bulk dielectric contributions.