Ion Hydration: Thermodynamic and Structural Analysis with an Integral Equation Theory of Liquids

Structure and thermodynamics of a 2D Lennard-Jones hexagonal fluid

The thermodynamic and structural properties of the 2D hexagonal soft-sites fluid are examined by integral equation theory benchmarked against extensive Monte Carlo simulations. Hexamers are built of six equal Lennard-Jones segments. Site-site integral equation theory is used to compute site-site correlation functions, excess internal energies and isotherms over a wide range of conditions and compared with results obtained from Monte Carlo simulations. Various approaches for computing the pressure are discussed as well. Satisfactory qualitative agreement between theory and simulations is found with details depending on the applied closure relation.

Renormalized site density functional theory for models of ion hydration

The development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile the co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. The renormalized formulation of SDFT (RSDFT), which we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system-a problem that historically has always been one of the most difficult cases for SDFT applications. Using a two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional three-dimensional reference interaction site model implementations and may prove useful in coarse grained simulations based on two-site solvent models.

Characterization of Specific Ion Effects on PI(4,5)P2 Clustering: Molecular Dynamics Simulations and Graph-Theoretic Analysis

Numerous cellular functions mediated by phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2; PIP2) involve clustering of the lipid as well as colocalization with other lipids. Although the cation-mediated electrostatic interaction is regarded as the primary clustering mechanism, the ion-specific nature of the intermolecular network formation makes it challenging to characterize the clusters. Here we use all-atom molecular dynamics (MD) simulations of PIP2 monolayers and graph-theoretic analysis to gain insight into the phenomenon. MD simulations reveal that the intermolecular interactions preferentially occur between specific cations and phosphate groups (P1, P4, and P5) of the inositol headgroup with better-matched kosmotropic/chaotropic characters consistent with the law of matching water affinities (LMWA). Ca2+ is strongly attracted to P4/P5, while K+ preferentially binds to P1; Na+ interacts with both P4/P5 and P1. These specific interactions lead to the characteristic clustering patterns. Specificially, the size distributions and structures of PIP2 clusters generated by kosmotropic cations Ca2+ and Na+ are bimodal, with a combination of small and large clusters, while there is little clustering in the presence of only chaotropic K+; the largest clusters are obtained in systems with all three cations. The small-world network (a model with both local and long-range connections) best characterizes the clusters, followed by the random and the scale-free networks. More generally, the present results interpreted within the LMWA are consistent with the relative eukaryotic intracellular concentrations Ca2+ ≪ Na+ < Mg2+ < K+; that is, concentrations of Ca2+ and Na+ must be low to prevent damaging aggregation of lipids, DNA, RNA and phosphate-containing proteins.

Ion hydration: linking self-diffusion and reorientational motion to water structure

A link between water dynamics and the “water structure” has been established through the combination of the extended jump model, transition state theory and the Kirkwood-Buff theory.

Small molecule hydration energy and entropy from 3D-RISM

Implicit solvent models offer an attractive way to estimate the effects of a solvent environment on the properties of small or large solutes without the complications of explicit simulations. One common test of accuracy is to compute the free energy of transfer from gas to liquid for a variety of small molecules, since many of these values have been measured. Studies of the temperature dependence of these values (i.e. solvation enthalpies and entropies) can provide additional insights into the performance of implicit solvent models. Here, we show how to compute temperature derivatives of hydration free energies for the 3D-RISM integral equation approach. We have computed hydration free energies of 1123 small drug-like molecules (both neutral and charged). Temperature derivatives were also used to calculate hydration energies and entropies of 74 of these molecules (both neutral and charged) for which experimental data is available. While direct results have rather poor agreement with experiment, we have found that several previously proposed linear hydration free energy correction schemes give good agreement with experiment. These corrections also provide good agreement for hydration energies and entropies though simple extensions are required in some cases.

A non-Bornian analysis of the Gibbs energy of hydration for organic ions

The hydration energy of organic ions can be well evaluated from the distribution of surface field strength, by using a simple semi-empirical equation.

Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions

Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.

Integral Equation Theory of Biomolecules and Electrolytes

The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.

Molecular dynamics investigation of dipeptide-transition metal salts in aqueous solutions

Molecular dynamics (MD) simulations of glycylglycine dipeptide with transition metal ions (Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), and Zn(2+)) in aqueous solutions have been carried out to get an insight into the solvation structure, intermolecular interactions, and salt effects in these systems. The solvation structure and hydrogen bonding were described in terms of radial distribution function (RDF) and spatial distribution function (SDF). The dynamical properties of the solvation structure were also analyzed in terms of diffusion and residence times. The simulation results show the presence of a well-defined first hydration shell around the dipeptide, with water molecules forming hydrogen bonds to the polar groups of the dipeptide. This shell is, however, affected by the strong electric field of divalent metal ions, which at higher ion concentrations lead to the shift in the dipeptide-water RDFs. Higher salt concentrations lead also to increased residence times and slower diffusion rates. In general, smaller ions (Cu(2+), Zn(2+)) demonstrate stronger binding to dipeptide than the larger ones (Fe(2+), Mn(2+)). Simulations do not show any stronger association of peptide molecules indicating their dissolution in water. The above results may be of potential interest to future researchers on these molecular interactions.

Receptacle model of salting-in by tetramethylammonium ions

Water is a poor solvent for nonpolar solutes. Water containing ions is an even poorer solvent. According to standard terminology, the tendency of salts to precipitate oils from water is called salting-out. However, interestingly, some salt ions, such as tetramethylammonium (TMA), cause instead the salting-in of hydrophobic solutes. Even more puzzling, there is a systematic dependence on solute size. TMA causes the salting-out of small hydrophobes and the salting-in of larger nonpolar solutes. We study these effects using NPT Monte Carlo simulations of the Mercedes-Benz (MB) + dipole model of water, which was previously shown to account for hydrophobic effects and ion solubilities in water. The present model gives a structural interpretation for the thermodynamics of salting-in. The TMA structure allows deep penetration by a first shell of waters, the dipoles of which interact electrostatically with the ion. This first water shell sets up a second water shell that is shaped to act as a receptacle that binds the nonpolar solute. In this way, a nonpolar solute can actually bind more tightly to the TMA ion than to another hydrophobe, leading to the increased solubility and salting-in. Such structuring may also explain why molecular ions do not follow the same charge density series as atomic ions do.

Computational analysis of water residence on ceramide and sphingomyelin bilayer membranes

Many physical chemical properties of lipid membranes, for example, the thickness, phase state, order parameter, and fluidity, can be understood straightforwardly. Water residence on a membrane is, however, an exception. To tackle this problem, we have performed molecular dynamics simulations of the distribution of water normal to the surface of several lipid membranes and from this deduced the associated water residence time. Our analysis of the results clearly indicates that lipid membranes have hydration shells on their surface, just as a solute in an aqueous solution does, and that the water residence time can be estimated from the potential for the mean force field derived from the distribution function of the water. We have done this atomic-scale analysis for ceramide bilayers and contrasted the calculation results with those for sphingomyelin bilayers, revealing that sphingomyelin bilayers can retain water molecules longer than ceramide bilayers and that the total number of water molecules retained on the membrane surface of sphingomyelin is larger than that for ceramide. In addition, we find that not only polar atoms of lipid molecules, such as oxygen, but also non-polar atoms, such as carbon, influence the motion of water on the membranes.

An integral equation theory for inhomogeneous molecular fluids: the reference interaction site model approach

An integral equation theory which is applicable to inhomogeneous molecular liquids is proposed. The "inhomogeneous reference interaction site model (RISM)" equation derived here is a natural extension of the RISM equation to inhomogeneous systems. This theory makes it possible to calculate the pair correlation function between two molecules which are located at different density regions. We also propose approximations concerning the closure relation and the intramolecular susceptibility of inhomogeneous molecular liquids. As a preliminary application of the theory, the hydration structure around an ion is investigated. Lithium, sodium, and potassium cations are chosen as the solute. Using the Percus trick, the local density of solvent around an ion is expressed in terms of the solute-solvent pair correlation function calculated from the RISM theory. We then analyze the hydration structure around an ion through the triplet correlation function which is defined with the inhomogeneous pair correlation function and the local density of the solvent. The results of the triplet correlation functions for cations indicate that the thermal fluctuation of the hydration shell is closely related to the size of the solute ion. The triplet correlation function from the present theory is also compared with that from the Kirkwood superposition approximation, which substitutes the inhomogeneous pair correlation by the homogeneous one. For the lithium ion, the behavior of the triplet correlation functions from the present theory shows marked differences from the one calculated within the Kirkwood approximation.

Partial molar volume of proteins studied by the three-dimensional reference interaction site model theory

The three-dimensional reference interaction site model (3D-RISM) theory is applied to the analysis of hydration effects on the partial molar volume of proteins. For the native structure of some proteins, the partial molar volume is decomposed into geometric and hydration contributions using the 3D-RISM theory combined with the geometric volume calculation. The hydration contributions are correlated with the surface properties of the protein. The thermal volume, which is the volume of voids around the protein induced by the thermal fluctuation of water molecules, is directly proportional to the accessible surface area of the protein. The interaction volume, which is the contribution of electrostatic interactions between the protein and water molecules, is apparently governed by the charged atomic groups on the protein surface. The polar atomic groups do not make any contribution to the interaction volume. The volume differences between low- and high-pressure structures of lysozyme are also analyzed by the present method.

Cellular magnesium acquisition: an anomaly in embryonic cation homeostasis

The intracellular dominance of magnesium ion makes clinical assessment difficult despite the critical role of Mg(++) in many key functions of cells and enzymes. There is general consensus that serum Mg(++) levels are not representative of the growing number of conditions for which magnesium is known to be important. There is no consensus method or sample source for testing for clinical purposes. High intracellular Mg(++) in vertebrate embryos results in part from interactions of cations which influence cell membrane transport systems. These are functionally competent from the earliest stages, at least transiently held over from the unfertilized ovum. Kinetic studies with radiotracer cations, osmolar variations, media lacking one or more of the four biological cations, Na(+), Mg(++), K(+), and Ca(++), and metabolic poison 0.05 mEq/L NaF, demonstrated that: (1) all four cations influence the behavior of the others, and (2) energy is required for uptake and efflux on different time scales, some against gradient. Na(+) uptake is energy dependent against an efflux gradient. The rate of K(+) loss is equal with or without fluoride, suggesting a lack of an energy requirement at these stages. Ca(++) efflux took twice as long in the presence of fluoride, likely due in part to intracellular binding. Mg(++) is anomalous in that early teleost vertebrate embryos have an intracellular content exceeding the surrounding sea water, an isolated unaffected yolk compartment, and a clear requirement for energy for both uptake and efflux. The physiological, pathological, and therapeutic roles of magnesium are poorly understood. This will change: (1) when (28)Mg is once again generally available at a reasonable cost for both basic research and clinical assessment, and (2) when serum or plasma levels are determined simultaneously with intracellular values, preferably as part of complete four cation profiles. Atomic absorption spectrophotometry, energy-dispersive x-ray analysis, and inductively coupled plasma emission spectroscopy on sublingual mucosal and peripheral blood samples are potential methods of value for coordinated assessments.

Partial molar volume of anionic polyelectrolytes in aqueous solution

In this work the partial molar volumes (V) of different anionic polyelectrolytes and hydrophobically modified polyelectrolytes (PHM) were measured. Polymers like polymaleic acid-co-styrene, polymaleic acid-co-1-olefin, polymaleic acid-co-vinyl-2-pyrrolidone, and polyacrylic acid (abbreviated as MAS-n, PA-n-K2, AMVP, and PAA, respectively) were employed. These materials were investigated by density measurements in highly dilute aqueous solutions. The molar volume results allow us to discuss the effect of the carboxylic groups and the contributions from the comonomeric principal chain. The PAA presents the smaller V, while the largest V value was for AMVP. The V of PHM shows a linear relationship with the number of methylene groups in the lateral chain. It is found that the magnitude of the contribution per methylene group decreases as the hydrophobic character of the environment increases.

Structure of electrolyte solutions sorbed in carbon nanospaces, studied by the replica RISM theory

The replica RISM theory is used to investigate the structure of electrolyte solutions confined in carbonized polyvinylidene chloride (PVDC) nanoporous material, compared to bulk electrolyte solution. Comparisons are made between the models of electrolyte solution sorbed in the carbonized PVDC material and a single carbon nanosphere in bulk electrolyte solution. Particular attention is paid to the chemical potential balance between the species of the sorbed electrolyte solution and the bulk solution in contact with the nanoporous material. As a result of the strong hydrophobicity of the carbonized PVDC material in the absence of activating chemical groups, the densities of water and ions sorbed in the material are remarkably low compared to those in the ambient bulk solution. The interaction between water molecules and cations becomes strong in nanospaces. It turns out that, in carbon nanopores, a cation adsorbed at the carbon surface is fully surrounded by the hydration shell of water molecules which separates the cation and the surface. Distinctively, an anion is adsorbed in direct contact with the carbon surface, which squeezes a part of its hydration shell out. The tendency increases toward smaller cations, which are characterized as "positive hydration" ions. In the bulk, cations are not hydrated so strongly and behave similarly to anions. The results suggest that the specific capacitance of an electric double-layer supercapacitor with nanoporous electrodes is intimately related to the solvation structure of electrolyte solution sorbed in nanopores, which is affected by the microscopic structure of the nanoporous electrode.

Modeling Water, the Hydrophobic Effect, and Ion Solvation

Water plays a central role in the structures and properties of biomolecules--proteins, nucleic acids, and membranes--and in their interactions with ligands and drugs. Over the past half century, our understanding of water has been advanced significantly owing to theoretical and computational modeling. However, like the blind men and the elephant, different models describe different aspects of water's behavior. The trend in water modeling has been toward finer-scale properties and increasing structural detail, at increasing computational expense. Recently, our labs and others have moved in the opposite direction, toward simpler physical models, focusing on more global properties-water's thermodynamics, phase diagram, and solvation properties, for example-and toward less computational expense. Simplified models can guide a better understanding of water in ways that complement what we learn from more complex models. One ultimate goal is more tractable models for computer simulations of biomolecules. This review gives a perspective from simple models on how the physical properties of water-as a pure liquid and as a solvent-derive from the geometric and hydrogen bonding properties of water.

Hydrophobic effects on partial molar volume

The hydrophobic effects on partial molar volume (PMV) are investigated as a PMV change in the transfer of a benzenelike nonpolar solute from the nonpolar solvent to water, using an integral equation theory of liquids. The volume change is divided into two effects. One is the "packing" effect in the transfer from the nonpolar solvent to hypothetical "nonpolar water" without hydrogen bonding networks. The other is the "iceberg" effect in the transfer from nonpolar water to water. The results indicate that the packing effect is negative and a half compensated by the positive iceberg effect. The packing effect is explained by the difference in the solvent compressibility. Further investigation shows that the sign and magnitude of the volume change depend on the solute size and the solvent compressibility. The finding gives a significant implication that the exposure of a hydrophobic residue caused by protein denaturation can either increase or decrease the PMV of protein depending on the size of the residue and the fluctuation of its surroundings.

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.