An analysis of molecular packing and chemical association in liquid water using quasichemical theory

Unraveling the Influence of Osmolytes on Water Hydrogen-Bond Network: From Local Structure to Graph Theory Analysis

Water structure in aqueous osmolyte solutions, deduced from the slight alteration in the water-water radial distribution function, the decrease in water-water hydrogen bonding, and tetrahedral ordering based only on the orientation of nearest water molecules derived from the molecular dynamics simulations, appears to have been perturbed. A careful analysis, however, reveals that the hydrogen bonding and the tetrahedral ordering around a water molecule in binary solutions remain intact as in neat water when the contribution of osmolyte-water interactions is appropriately incorporated. Furthermore, the distribution of the water binding energies and the water excess chemical potential of solvation in solutions are also pretty much the same as in neat water. Osmolytes are, therefore, well integrated into the hydrogen-bond network of water. Indeed, osmolytes tend to preferentially hydrogen bond with water molecules and their interaction energies are strongly correlated to their hydrogen-bonding capability. The graph network analysis, further, illustrates that osmolytes act as hubs in the percolated hydrogen-bond network of solutions. The degree of hydrogen bonding of osmolytes predominantly determines all of the network properties. Osmolytes like ethanol that form fewer hydrogen bonds than a water molecule disrupt the water hydrogen-bond network, while other osmolytes that form more hydrogen bonds effectively increase the connectivity among water molecules. Our observation of minimal variation in the local structure and the vitality of osmolyte-water hydrogen bonds on the solution network properties clearly imply that the direct interaction between protein and osmolytes is solely responsible for the protein stability. Further, the relevance of hydrogen bonds on solution properties suggests that the hydrogen-bonding interaction among protein, water, and osmolyte could be the key determinant of the protein conformation in solutions.

Thermodynamics of Hydration from the Perspective of the Molecular Quasichemical Theory of Solutions

The quasichemical organization of the potential distribution theorem, molecular quasichemical theory (QCT), enables practical calculations and also provides a conceptual framework for molecular hydration phenomena. QCT can be viewed from multiple perspectives: (a) as a way to regularize an ill-conditioned statistical thermodynamic problem; (b) as an introduction of and emphasis on the neighborship characteristics of a solute of interest; or (c) as a way to include accurate electronic structure descriptions of near-neighbor interactions in defensible statistical thermodynamics by clearly defining neighborship clusters. The theory has been applied to solutes of a wide range of chemical complexity, ranging from ions that interact with water with both long-ranged and chemically intricate short-ranged interactions, to solutes that interact with water solely through traditional van der Waals interations, and including water itself. The solutes range in variety from monatomic ions to chemically heterogeneous macromolecules. A notable feature of QCT is that, in applying the theory to this range of solutes, the theory itself provides guidance on the necessary approximations and simplifications that can facilitate the calculations. In this Perspective, we develop these ideas and document them with examples that reveal the insights that can be extracted using the QCT formulation.

How Alcoholic Disinfectants Affect Coronavirus Model Membranes: Membrane Fluidity, Permeability, and Disintegration

Atomistic molecular dynamics simulations have been carried out with a view to investigating the stability of the SARS-CoV-2 exterior membrane with respect to two common disinfectants, namely, aqueous solutions of ethanol and n-propanol. We used dipalmitoylphosphatidylcholine (DPPC) as a model membrane material and did simulations on both gel and liquid crystalline phases of membrane surrounded by aqueous solutions of varying alcohol concentrations (up to 17.5 mol %). While a moderate effect of alcohol on the gel phase of membrane is observed, its liquid crystalline phase is shown to be influenced dramatically by either alcohol. Our results show that aqueous solutions of only 5 and 10 mol % alcohol already have significant weakening effects on the membrane. The effects of n-propanol are always stronger than those of ethanol. The membrane changes its structure, when exposed to disinfectant solutions; uptake of alcohol causes it to swell laterally but to shrink vertically. At the same time, the orientational order of lipid tails decreases significantly. Metadynamics and grand-canonical ensemble simulations were done to calculate the free-energy profiles for permeation of alcohol and alcohol/water solubility in the DPPC. We found that the free-energy barrier to permeation of the DPPC liquid crystalline phase by all permeants is significantly lowered by alcohol uptake. At a disinfectant concentration of 10 mol %, it becomes insignificant enough to allow almost free passage of the disinfectant to the inside of the virus to cause damage there. It should be noted that the disinfectant also causes the barrier for water permeation to drop. Furthermore, the shrinking of the membrane thickness shortens the gap needed to be crossed by penetrants from outside the virus into its core. The lateral swelling also increases the average distance between head groups, which is a secondary barrier to membrane penetration, and hence further increases the penetration by disinfectants. At alcohol concentrations in the disinfectant solution above 15 mol %, we reliably observe disintegration of the DPPC membrane in its liquid crystalline phase.

Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F-) Anion Hydration

Accurate predictions of the hydration free energy for anions typically has been more challenging than that for cations. Hydrogen bond donation to the anion in hydrated clusters such as F(H₂O) _n ^- can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that even though anharmonic effects are interesting they need not be of troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed, and application of implicit solvation models to these structures can incur more serious errors than those for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues.

Mini-grand canonical ensemble: Chemical potential in the solvation shell

Quantifying the statistics of occupancy of solvent molecules in the vicinity of solutes is central to our understanding of solvation phenomena. Number fluctuations in small solvation shells around solutes cannot be described within the macroscopic grand canonical framework using a single chemical potential that represents the solvent bath. In this communication, we hypothesize that molecular-sized observation volumes such as solvation shells are best described by coupling the solvation shell with a mixture of particle baths each with its own chemical potential. We confirm our hypotheses by studying the enhanced fluctuations in the occupancy statistics of hard sphere solvent particles around a distinguished hard sphere solute particle. Connections with established theories of solvation are also discussed.

Perturbation theory for water with an associating reference fluid

The theoretical description of the thermodynamics of water is challenged by the structural transition towards tetrahedral symmetry at ambient conditions. As perturbation theories typically assume a spherically symmetric reference fluid, they are incapable of accurately describing the liquid properties of water at ambient conditions. In this paper we address this problem by introducing the concept of an associated reference perturbation theory (APT). In APT we treat the reference fluid as an associating hard sphere fluid which transitions to tetrahedral symmetry in the fully hydrogen bonded limit. We calculate this transition in a theoretically self-consistent manner without appealing to molecular simulations. This associated reference provides the reference fluid for a second order Barker-Henderson perturbative treatment of the long-range attractions. We demonstrate that this approach gives a significantly improved description of water as compared to standard perturbation theories.

Conditional solvation thermodynamics of isoleucine in model peptides and the limitations of the group-transfer model

The hydration thermodynamics of the amino acid X relative to the reference G (glycine) or the hydration thermodynamics of a small-molecule analog of the side chain of X is often used to model the contribution of X to protein stability and solution thermodynamics. We consider the reasons for successes and limitations of this approach by calculating and comparing the conditional excess free energy, enthalpy, and entropy of hydration of the isoleucine side chain in zwitterionic isoleucine, in extended penta-peptides, and in helical deca-peptides. Butane in gauche conformation serves as a small-molecule analog for the isoleucine side chain. Parsing the hydrophobic and hydrophilic contributions to hydration for the side chain shows that both of these aspects of hydration are context-sensitive. Furthermore, analyzing the solute–solvent interaction contribution to the conditional excess enthalpy of the side chain shows that what is nominally considered a property of the side chain includes entirely nonobvious contributions of the background. The context-sensitivity of hydrophobic and hydrophilic hydration and the conflation of background contributions with energetics attributed to the side chain limit the ability of a single scaling factor, such as the fractional solvent exposure of the group in the protein, to map the component energetic contributions of the model-compound data to their value in the protein. But ignoring the origin of cancellations in the underlying components the group-transfer model may appear to provide a reasonable estimate of the free energy for a given error tolerance.

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Ion-water clusters, bulk medium effects, and ion hydration

Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical approximation to the quasichemical theory provides a transparent framework to anchor such efforts. Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate the different roles of the bulk medium. We find that the bulk medium can stabilize configurations of the cluster that are usually not observed in the gas phase, while also simultaneously lowering the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the coordination number that minimizes the excess chemical potential of the ion is identified as the optimal or most probable coordination number, for such soft ions the optimum cluster size and the hydration thermodynamics obtained with and without account of the bulk medium on the ion-water clustering reaction can be different. The ideas presented in this work are expected to be relevant to experimental studies that translate thermochemistry of ion-water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches that combine high-level calculations on clusters with coarse-grained models of the medium.

A local entropic signature of specific ion hydration

Monovalent ion hydration entropies are analyzed via energetic partitioning of the potential distribution theorem free energy. Extensive molecular dynamics simulations and free energy calculations are performed over a range of temperatures to determine the electrostatic and van der Waals components of the entropy. The far-field electrostatic contribution is negative and small in magnitude, and it does not vary significantly as a function of ion size, consistent with dielectric models. The local electrostatic contribution, however, varies widely as a function of ion size; the sign yields a direct indication of the kosmotropic (strongly hydrated) or chaotropic (weakly hydrated) nature of the ion hydration. The results provide a thermodynamic signature for specific ion effects in hydration and are consistent with experiments that suggest minimal perturbations of water structure outside the first hydration shell. The hydration entropies are also examined in relation to the corresponding entropies for the isoelectronic rare gas pairs; an inverse correlation is observed, as expected from thermodynamic hydration data.

Probing the thermodynamics of competitive ion binding using minimum energy structures

Ion binding is known to affect the properties of biomolecules and is directly involved in many biochemical pathways. Because of the highly polar environments where ions are found, a quantum-mechanical treatment is preferable for understanding the energetics of competitive ion binding. Due to computational cost, a quantum mechanical treatment may involve several approximations, however, whose validity can be difficult to determine. Using thermodynamic cycles, we show how intuitive models for complicated ion binding reactions can be built up from simplified, isolated ion-ligand binding site geometries suitable for quantum mechanical treatment. First, the ion binding free energies of individual, minimum energy structures determine their intrinsic ion selectivities. Next, the relative propensity for each minimum energy structure is determined locally from the balance of ion-ligand and ligand-ligand interaction energies. Finally, the environment external to the binding site exerts its influence both through long-ranged dispersive and electrostatic interactions with the binding site as well as indirectly through shifting the binding site compositional and structural preferences. The resulting picture unifies field-strength, topological control, and phase activation viewpoints into a single theory that explicitly indicates the important role of solute coordination state on overall reaction energetics. As an example, we show that the Na(+) → K(+) selectivities can be recovered by correctly considering the conformational contribution to the selectivity. This can be done even when constraining configuration space to the neighborhood around a single, arbitrarily chosen, minimum energy structure. Structural regions around minima for K(+)- and Na(+)-water clusters are exhibited that display both rigid/mechanical and disordered/entropic selectivity mechanisms for both Na(+) and K(+). Thermodynamic consequences of the theory are discussed with an emphasis on the role of coordination structure in determining experimental properties of ions in complex biological environments.

Communication: Thermodynamics of water modeled using ab initio simulations

We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. Assuming classical statistical mechanical simulations at 350 K model the liquid at 298 K, the calculated free energy is found in fair agreement with experiments, but the excess internal energy and hence also the excess entropy are not. The utility of thermodynamic characterization in understanding the role of high temperatures to mimic nuclear quantum effects and in evaluating ab initio simulations is noted.

Ion selectivity from local configurations of ligands in solutions and ion channels

Probabilities of numbers of ligands proximal to an ion lead to simple, general formulae for the free energy of ion selectivity between different media. That free energy does not depend on the definition of an inner shell for ligand-counting, but other quantities of mechanistic interest do. If analysis is restricted to a specific coordination number, then two distinct probabilities are required to obtain the free energy in addition. The normalizations of those distributions produce partition function formulae for the free energy. Quasi-chemical theory introduces concepts of chemical equilibrium, then seeks the probability that is simplest to estimate, that of the most probable coordination number. Quasi-chemical theory establishes the utility of distributions of ligand-number, and sharpens our understanding of quasi-chemical calculations based on electronic structure methods. This development identifies contributions with clear physical interpretations, and shows that evaluation of those contributions can establish a mechanistic understanding of the selectivity in ion channels.

Thermodynamically dominant hydration structures of aqueous ions

The hydration free energy of an ion is separated into a chemical term, arising due to the interaction of the ion with water molecules within the defined coordination sphere (the inner shell), a packing contribution, accounting for forming an ion-free coordination sphere (the observation volume) in the solvent, and a long range correction, accounting for the interaction of the ion with the solvent outside the coordination sphere. The chemical term is recast as a sum over coordination states, with the nth term depending on the probability of observing n water molecules in the observation volume and the free energy of assembling the n water molecules around the ion in the presence of the outer-shell solvent. Each stepwise increment in the coordination number more fully accounts for the chemical contribution, and this molecular aufbau approach is used to interrogate the thermodynamic importance of various hydration structures X[H(2)O](n) of X(aq) (X = Na(+), K(+), F(-)) within a classical molecular mechanics framework. States with n less than (and at best equal to) the most probable coordination state ñ account for all of the chemical term and evince the role of the ion in drawing water molecules into the coordination sphere. For states with n > ñ, the influence of the ion is tempered and changes in coordination states due to density fluctuations in water also appear important. Thus the influence of the ion on the solvent matrix is local, and only a subset of water molecules (n < or = ñ) contributes dominantly to the hydration thermodynamics. The n = 4 state of Na(+) (ñ = 5) and K(+) (ñ = 7) and the n = 6 state of F(-) (ñ = 6) are thermodynamically dominant; adding a water molecule to the dominant state additionally contributes only about 2-3 k(B)T toward the chemical term, but removing a water molecule is very unfavorable.

Ab initio molecular dynamics calculations of ion hydration free energies

We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or "lambda-path" technique to compute the intrinsic hydration free energies of Li(+), Cl(-), and Ag(+) ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential (phi) contributions, we obtain absolute AIMD hydration free energies (DeltaG(hyd)) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model phi predictions. The sums of Li(+)/Cl(-) and Ag(+)/Cl(-) AIMD DeltaG(hyd), which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag(+)+Ni(+)-->Ag+Ni(2+) in water. The predictions for this reaction suggest that existing estimates of DeltaG(hyd) for unstable radiolysis intermediates such as Ni(+) may need to be extensively revised.

Dewetting and hydrophobic interaction in physical and biological systems

Hydrophobicity manifests itself differently on large and small length scales. This review focuses on large-length-scale hydrophobicity, particularly on dewetting at single hydrophobic surfaces and drying in regions bounded on two or more sides by hydrophobic surfaces. We review applicable theories, simulations, and experiments pertaining to large-scale hydrophobicity in physical and biomolecular systems and clarify some of the critical issues pertaining to this subject. Given space constraints, we cannot review all the significant and interesting work in this active field.

Unraveling water's entropic mysteries: a unified view of nonpolar, polar, and ionic hydration

[Figure: see text]. Most chemical processes on earth are intimately linked to the unique properties of water, relying on the versatility with which water interacts with molecules of varying sizes and polarities. These interactions determine everything from the structure and activity of proteins and living cells to the geological partitioning of water, oil, and minerals in the Earth's crust. The role of hydrophobic hydration in the formation of biological membranes and in protein folding, as well as the importance of electrostatic interactions in the hydration of polar and ionic species, are all well known. However, the underlying molecular mechanisms of hydration are often not as well understood. This Account summarizes and extends emerging understandings of these mechanisms to reveal a newly unified view of hydration and explain previously mystifying observations. For example, rare gas atoms (e.g., Ar) and alkali-halide ions (e.g., K+ and Cl-) have nearly identical experimental hydration entropies, despite the significant charge-induced reorganization of water molecules. Here, we explain how such previously mysterious observations may be understood as arising from Gibbs inequalities, which impose rigorous energetic upper and lower bounds on both hydration free energies and entropies. These fundamental Gibbs bounds depend only on the average interaction energy of a solute with water, thus providing a deep link between solute-water interaction energies and entropies. One of the surprising consequences of the emerging picture is the understanding that the hydration of an ion produces two large but nearly perfectly cancelling, entropic contributions: a negative ion-water interaction entropy and a positive water reorganization entropy. Recent work has also clarified the relationship between the strong cohesive energy of water and the free energy required to form an empty hole (cavity) in water. Here, we explain how linear response theory (whose roots may also be traced to Gibbs inequalities) can provide remarkably accurate descriptions of the process of filling aqueous cavities with nonpolar, polar, or charged molecules. The hydration of nonpolar molecules is well-described by first-order perturbation theory, which implies that turning on solute-water van der Waals interactions does not induce a significant change in water structure. The larger changes in water structure that are induced by polar and ionic solutes are well-described by second-order perturbation theory, which is equivalent to linear response theory. Comparisons of the free energies of nonpolar and polar or ionic solutes may be used to experimentally determine electrostatic contributions to water reorganization energies and entropies. The success of this approach implies that water's ability to respond to solutes of various polarities is far from saturated, as illustrated by simulations of acetonitrile (CH 3CN) in water, which reveal that even such a strongly dipolar solute only produces subtle changes in the structure of water.

Studies of the thermodynamic properties of hydrogen gas in bulk water

The thermodynamic properties of hydrogen gas in liquid water are investigated using Monte Carlo molecular simulation and the quasichemical theory of liquids. The free energy of hydrogen hydration obtained by Monte Carlo simulations agrees well with the experimental result, indicating that the classical force fields used in this work provide an adequate description of intermolecular interactions in the aqueous hydrogen system. Two estimates of the hydration free energy for hydrogen made within the framework of the quasichemical theory also agree reasonably well with experiment provided local anharmonic motions and distant interactions with explicit solvent are treated. Both quasichemical estimates indicate that the hydration free energy results from a balance between chemical association and molecular packing. Additionally, the results suggest that the molecular packing term is almost equally driven by unfavorable enthalpic and entropic components.

Non-van der Waals Treatment of the Hydrophobic Solubilities of CF4

A quasi-chemical theory implemented on the basis of molecular simulation is derived and tested for the hydrophobic hydration of CF4(aq). The theory formulated here subsumes a van der Waals treatment of solvation and identifies contributions to the hydration free energy of CF4(aq) that naturally arise from chemical contributions defined by quasi-chemical theory and fluctuation contributions analogous to Debye-Hückel or random phase approximations. The resulting Gaussian statistical thermodynamic model avoids consideration of hypothetical drying-then-rewetting problems and is physically reliable in these applications as judged by the size of the fluctuation contribution. The specific results here confirm that unfavorable tails of binding energy distributions reflect few-body close solute-solvent encounters. The solvent near-neighbors are pushed by the medium into unfavorable interactions with the solute, in contrast to the alternative view that a preformed interface is pulled by the solute-solvent attractive interactions into contact with the solute. The polyatomic model of CF4(aq) studied gives a satisfactory description of the experimental solubilities including the temperature dependence. The proximal distributions evaluated here for polyatomic solutes accurately reconstruct the observed distributions of water near these molecules which are nonspherical. These results suggest that drying is not an essential consideration for the hydrophobic solubilities of CF4, or of C(CH3)4 which is more soluble.