banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
An improved coarse-grained model of solvation and the hydrophobic effect
Rent this article for


Image of FIG. 1.
FIG. 1.

Schematic showing the solute and the large length-scale density field on a grid.

Image of FIG. 2.
FIG. 2.

Solvation free energies G of hard spheres of increasing radii, as calculated from explicit SPC/E water simulations (Ref. 19) (solid blue), from the coarse-grained model (solid black), and from one of the most common variants of GBSA (Ref. 61) (arrow at bottom right). Different parameterizations of GBSA yield nonpolar solvation energies that may differ by as much as a factor of ten, reflecting the crudeness of the SA portion of the model (Ref. 61). When the coarse-grained model has no unbalancing potential (a = 0, dashed gray), the intermediate-size regime is only qualitatively reproduced. For large spheres, the ratio of G to surface area tends to the liquid-vapor surface tension γ (horizontal red dots). Inset: illustration of lattice artifacts. The spheres are centered at different offsets from the lattice: a generic position (0.98, 0.79, 1.89 Å) that breaks all rotational and mirror symmetries (black), a lattice cell corner (blue), and a lattice cell center (red). All three curves are identical for R ⩽ 0.35 nm.

Image of FIG. 3.
FIG. 3.

Solvation free energies G of hexagonal plates, as a function of plate size, as calculated by the coarse-grained model (solid lines), by explicit SPC/E water simulations (points), and by the same GBSA variant as in Fig. 2 (arrow on right). Three values of the attractive interaction strength η are shown: 0.0 (black), 0.5 (red), and 1.0 (blue). Solvent-accessible surface areas (SASAs) were calculated using VMD (Ref. 62) with a particle radius of 1.97 Å, a solvent radius of 1.4 Å, and 1 000 000 samples per atom. The bulk liquid-vapor surface tension of water (horizontal red dots) is shown. Inset: detail of the hexagonal plate. The solvent-excluded volume of each oily site is a sphere of radius R 0 = 3.37 Å.

Image of FIG. 4.
FIG. 4.

Water number distribution in a 12 × 12 × 12 Å3 cube, as obtained using explicit SPC/E water simulations (Ref. 33), the present model, the present model without the unbalancing potential (a = 0), and a model with an Ising lattice gas and no unbalancing potential.

Image of FIG. 5.
FIG. 5.

Water number distributions in a probe volume of size 24 × 24 × 3 Å3 immediately adjacent to a model plate solute (inset) of varying attractive strength η, in the coarse-grained model (solid lines) and in explicit SPC/E water (points). Defining z = 0 to be the plane passing through the plate center, points in the probe volume (green) satisfy 5 Å < z < 8 Å, so that a water molecule touching the plate is located at the edge of the probe volume (Ref. 33).

Image of FIG. 6.
FIG. 6.

Setup for examining water fluctuations under confinement (here, d = 8 Å). The model hydrophobic plates (Ref. 33) (gray particles) are placed d Å apart: taking into account the van der Waals radii of about 2 Å of the plates’ oily particles and the 3 Å thickness of each plate, the center of the first plate is placed at z = 0, the center of the second plate is placed at z = d + 7 Å. The van der Waals radius of water (red and white sticks) being about 1.5 Å, the 24 × 24 × (d − 3) Å3 probe volume (green) extends from z = 5 Å to z = d + 2 Å. The plates are not perfectly flat, so some waters fit between the plates and the probe volume.

Image of FIG. 7.
FIG. 7.

Phase diagram for the interplate region of the system depicted in Fig. 6. For the explicit SPC/E water simulations (blue), each symbol corresponds to an individual P V (N) distribution that we have calculated (filled: wet state stable; open: dry state stable). The phase boundary (blue dashes) is estimated from a linear interpolation of the relative stability of the wet and dry states. The relative stability is determined from the relative depths of the basins in −ln P V (N). The phase boundary for the present model (red solid line) was estimated from a dense sampling of P V (N) distributions, and is accurate to ±0.1 Å in d and ±0.1 in η.

Image of FIG. 8.
FIG. 8.

Water density distribution of confined water 1 Å from coexistence. These distributions are for the system depicted in Fig. 6 when η = 0.5. Coexistence lines are shown in Fig. 7. The explicit water simulation data (black) correspond to d = 11 Å, while the coarse-grained model (red) results correspond to d = 9 Å. The remaining P V (N) distributions are included in the Supplementary Data (Ref. 72).

Image of FIG. 9.
FIG. 9.

Constructing n(r) from {n i }. The binary field specifies whether the density at the center of each lattice cell should be that of the liquid or that of the vapor. Between cell centers, the density is interpolated using the basis function ψ(x) (whose form for water is shown in the lower left panel). The dashed lines delineate the domain of integration of the local free energy h i given by Eq. (10).

Image of FIG. 10.
FIG. 10.

Solvation free energies G of hard spheres as a function of sphere radius, where the term H +[{n i }] (Eq. (D1)) has been included and the renormalization constant K has been set to 1 (solid black), compared to the simpler model in Eq. (13) (circles). The averages of −〈H int[{n i }]〉 (red) and 〈H +[{n i }]〉 (blue) are nearly proportional to each other. Left inset: implied renormalization constant K, equal to 〈H int[{n i }] + H +[{n i }]〉/〈H int[{n i }]〉. Note that both the numerator and denominator take on essentially zero value for R ≲ 0.4 nm. Right inset: implied value of K for hexagonal plate solute (Fig. 3) with η = 1.0. The implied value of K is similar for different η.

Image of FIG. 11.
FIG. 11.

Solvation free energies G of spheres in the model of Ref. 14 (black), for cell sizes λ = 2.1 Å (solid) and λ = 2.3 Å (dashes). The use of the Ising Hamiltonian causes the average value of H large[{n i }] (red) to significantly exceed the solvation free energy, but also leads to large excess entropies (blue, TS = 〈H〉 − G). At λ = 2.1 Å, but not at λ = 2.3 Å, a fortuitous cancellation leads to correct solvation free energies.


Generic image for table
Table I.

Relative interfacial free energy h i for each distinct neighboring cell configuration (diagrams after Ref. 95). Highlighted corners denote cells j with n j = 1, whereas the others refer to cells with n j = 0; cell i is the lower-left corner in the back. To aid the eye, a schematic of the implied liquid-vapor interface of each configuration is shown in orange. The values of h i are inversion-symmetric: interchanging highlighted and unhighlighted corners yields the same interface and interfacial energy. Also shown are the values of h i that would reproduce the energetics of the standard Ising lattice gas, namely γλ2ij(n i n j )2.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: An improved coarse-grained model of solvation and the hydrophobic effect