How Crystals Nucleate and Grow in Aqueous NaCl Solution

Large-scale molecular dynamics simulations (64 000 particles) are used to examine the microscopic mechanism of crystal nucleation and growth in a slightly supersaturated solution of NaCl in water at 300 K and 1 atm. Early-stage nucleation is observed, and the growth of a single crystal is followed for ∼140 ns. It is shown that the nucleation and growth process is better described by Ostwald's rule of stages than by classical nucleation theory. Crystal nucleation originates in a region where the local salt concentration exceeds that of the bulk solution. The early-stage nucleus is a loosely ordered arrangement of ions that retains a significant amount of water. The residual water is slowly removed as the crystal grows and evolves toward its stable anhydrous state.

A Critical Assessment of Two-Body and Three-Body Interactions in Water

The microscopic behavior of water under different conditions and in different environments remains the subject of intense debate. A great number of the controversies arise due to the contradictory predictions obtained within different theoretical models. Relative to conclusions derived from force fields or density functional theory, there is comparably less room to dispute highly correlated electronic structure calculations. Unfortunately, such ab initio calculations are severely limited by system size. In this study, a detailed analysis of the two- and three-body water interactions evaluated at the CCSD(T) level is carried out to quantitatively assess the accuracy of several force fields, DFT models, and ab initio based interaction potentials that are commonly used in molecular simulations. On the basis of this analysis, a new model, HBB2-pol, is introduced which is capable of accurately mapping CCSD(T) results for water dimers and trimers into an efficient analytical function. The accuracy of HBB2-pol is further established through comparison with the experimentally determined second and third virial coefficients.

Toward a Universal Water Model: First Principles Simulations from the Dimer to the Liquid Phase

A full-dimensional model of water, HBB2-pol, derived entirely from "first-principles", is introduced and employed in computer simulations ranging from the dimer to the liquid. HBB2-pol provides excellent agreement with the measured second and third virial coefficients and, by construction, reproduces the dimer vibration-rotation-tunneling spectrum. The model also predicts the relative energy differences between isomers of small water clusters within the accuracy of highly correlated electronic structure methods. Importantly, when combined with simulation methods that explicitly include zero-point energy and quantum thermal motion, HBB2-pol accurately describes both structural and dynamical properties of the liquid phase.

On the Role of Dewetting Transitions in Host-Guest Binding Free Energy Calculations

We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host–guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states. Estimating free energies by TI can be especially challenging when there is a considerable difference in hydration between the two states of interest. We investigate these aspects using the popular TIP3P and TIP4P water models. TI free energy estimates through MD largely depend on water-related interactions, and water dynamics significantly affect the convergence of binding free energy calculations. Our results indicate that wetting/dewetting transitions play a major role in slowing the convergence of free energy estimation. We employ two alternative approaches—one analytical and the other empirically based on actual MD sampling—to correct for the standard state free energy. This correction is sizable (up to 4 kcal/mol), and the two approaches provide corrections that differ by about 1 kcal/mol. For the system considered here, the TIP4P water model combined with an analytical correction for the standard state free energy provides higher overall accuracy. This observation might be transferable to other systems in which water-related contributions dominate the binding process.

A transferable ab initio based force field for aqueous ions

We present a new polarizable force field for aqueous ions (Li(+), Na(+), K(+), Rb(+), Cs(+), Mg(2 +), Ca(2 +), Sr(2 +), and Cl(-)) derived from condensed phase ab initio calculations. We use maximally localized Wannier functions together with a generalized force and dipole-matching procedure to determine the whole set of parameters. Experimental data are then used only for validation purposes and a good agreement is obtained for structural, dynamic, and thermodynamic properties. The same procedure applied to crystalline phases allows to parametrize the interaction between cations and the chloride anion. Finally, we illustrate the good transferability of the force field to other thermodynamic conditions by investigating concentrated solutions.

Diffusion coefficient and shear viscosity of rigid water models

We report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.

A flexible model for water based on TIP4P/2005

A new flexible water model, TIP4P/2005f, is developed. The idea was to add intramolecular degrees of freedom to the successful rigid model TIP4P/2005 in order to try to improve the predictions for some properties, and to enable the calculation of new ones. The new model incorporates flexibility by means of a Morse potential for the bond stretching and a harmonic term for the angle bending. The parameters have been fitted to account for the peaks of the infrared spectrum of liquid water and to produce an averaged geometry close to that of TIP4P/2005. As for the intermolecular interactions, only a small change in the σ parameter of the Lennard-Jones potential has been introduced. The overall predictions are very close to those of TIP4P/2005. This ensures that the new model may be used with the same confidence as its predecessor in studies where a flexible model is advisable.

The role of conserved waters in conformational transitions of Q61H K-ras

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 µs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the C_α atoms of G60 at switch 2 and G10 at the P-loop and the N-C_α-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.

K-ras belongs to the Ras family of G-proteins that regulate cell proliferation and development. To execute its function, K-ras adopts different conformational states when it is active and inactive. In addition to these two states, it samples many transient intermediate conformations as it makes the transition from one state to the other. Mutations that affect the population of these states can cause cancer or developmental disorder. Using simulation approaches, here we show that a number of water molecules buried within the structure of an oncogenic K-ras protein modulate the distribution of its conformational states. Moreover, a detailed analysis based on two novel structural parameters revealed the existence of long-range water-mediated interactions that facilitate a dynamic coupling between the two lobes of the protein. These findings pave the way for a dynamics-guided strategy to inhibit abnormal Ras signaling.