The Relevance of Cavity Creation for Several Phenomena Occurring in Water

: The solvent-excluded volume effect is an under-appreciated general phenomenon occurring in liquids and playing a fundamental role in many cases. It is quantiﬁed and characterized by means of the theoretical concept of cavity creation and its Gibbs free energy cost. The magnitude of the reversible work of cavity creation proves to be particularly large in water, and this fact plays a key role for, among other things, the poor solubility of nonpolar species, the formation of host–guest complexes, and the folding of globular proteins. An analysis of some examples is provided in the present review


Introduction
The starting point of any theory able to describe processes occurring in liquids (i.e., pure liquids or solutions) is the recognition that a suitable void space, a cavity, has to be created to allow solute insertion . This is the simple consequence of the fact that liquids are a condensed state of matter and each molecule possess its own body. Cavity creation leads to a decrease in the number of configurations accessible to liquid molecules and thus leads to a solvent-excluded volume effect [22][23][24][25]. These words are right, but the matter has to be spelled out in more detail to reach a correct understanding (note that cavity creation is a theoretical concept and cannot be studied by performing experiments). Keeping a constant temperature and pressure, the creation of a cavity leads to an increase in liquid volume by the partial molar volume of the cavity itself. This volume increase does not cancel the solvent-excluded volume effect. If the cavity is to exist, the center of liquid molecules (assumed to be spherical) cannot go beyond the solvent-accessible surface area, SASA, and WASA in water [26], of the cavity itself. This means that the shell between the cavity van der Waals surface and the SASA is excluded to liquid molecules, causing a decrease in accessible configurations for basic geometric reasons. The latter constraint does affect the translational motion of all the liquid molecules, not solely the ones in the first solvation shell of the cavity (i.e., of the solute molecule to be hosted).
The solvent-excluded volume effect can be measured by calculating the reversible work of cavity creation, ∆G C , by means of analytical theories or computer simulations. Classic scaled particle theory, SPT [23][24][25][27][28][29], is a simple, geometry-based statistical mechanical model providing analytical formulas to calculate ∆G C for cavities of simple shape (i.e., a sphere, a prolate spherocylinder, and others) in liquids made up of hard particles. Its use in the case of water may appear strange, but it works well because the real liquid density is used as input in classic SPT calculations (i.e., density provides indirect information on the strength of the intermolecular attractions existing among liquid molecules; in addition, on the H-bonds between water molecules [30,31]). The classic SPT formulas to create a spherical cavity in a liquid (neglecting the pressure-volume term for its smallness at P = 1 atm) are: where R is the gas constant, αP is the isobaric thermal expansion coefficient of the liquid, ξ is the volume packing density of the liquid, which is defined as the ratio of the physical volume of a mole of liquid molecules over the liquid molar volume, v 1 (i.e., ξ = π · σ 1 3 · N Av /6 · v 1 ); x = σ C /σ 1 , and σ 1 is the hard sphere diameter of liquid molecules; σ C is the cavity diameter, defined as the diameter of the spherical region from which any part of liquid molecules is excluded. The ∆G C (SPT) magnitude depends upon the volume packing density of the liquid, ξ, that is, the fraction of the total liquid volume really occupied by liquid molecules, and the effective hard sphere diameter, σ 1 , of liquid molecules [25,31]. On increasing ξ, the void volume decreases and ∆G C increases; on decreasing σ 1 , the void volume is partitioned into smaller pieces and ∆G C increases (i.e., a significant fraction of the liquid volume is void, but most of these voids are too small to host an atom or a molecule). This implies that the effective diameter of liquid molecules is a fundamental length-scale for the liquid itself. The validity of these arguments has been verified and confirmed in several cases over the years [31].

Solvation of Noble Gases
A further test is provided here, by analyzing the solvation (i.e., the transfer from a fixed position in the gas phase to a fixed position in the liquid phase) of noble gases in water, carbon tetrachloride, CCl 4 , and benzene, C 6 H 6 . Experimental thermodynamic data at 25 • C and 1 atm [32][33][34], reported in Table 1, emphasize that: (1) noble gases are poorly soluble in water, with them being characterized by large positive ∆G · values, caused by large negative entropy changes; (2) Ar is characterized by positive ∆G · values also in CCl 4 and C 6 H 6 , while Xe is characterized by a negative ∆G · value in benzene. The ∆G C (SPT) values calculated for noble gases in the three liquids, whose molecules are assumed to be spherical, at 25 • C and 1 atm, are listed in the eighth column of Table 1. They prove to be largely positive in all cases. Actually, they are significantly larger in water with respect to the other two liquids (see the trends reported in Figure 1).

Figure 1.
Trend of ∆G C versus the cavity diameter for water, CCl 4 , and C 6 H 6 calculated by means of classic SPT, at 25 • C and 1 atm. The data necessary to perform the calculations are reported in the notes of Table 1.  [32][33][34], according to the Ben-Naim standard (i.e., the transfer from a fixed position in the gas phase to a fixed position in the liquid phase), at 25 • C and 1 atm, of noble gases in water (a), CCl 4 (b), and C 6 H 6 (c); the values of the hard sphere diameter and the Lennard-Jones energy parameter come from [35,36] with small modifications; the values of ∆G C are calculated by means of classic SPT analytical formulas [27,28]; those of E a are calculated by means of Pierotti's analytical formula [29]. Additional data used to perform the calculations [37]. Water: This holds because, even though the volume packing density of water is the smallest, ξ = 0.383 for water versus 0.503 for CCl 4 , and 0.513 for C 6 H 6 , the water molecules are the smallest, σ = 2.80 Å for water, 5.37 Å for CCl 4 , and 5.26 Å for C 6 H 6 [35][36][37]. In this respect, it is important to underscore that the effective hard sphere diameter assigned to water molecules is physically reliable because it corresponds to the location of the first peak in the oxygen-oxygen radial distribution function of water [38], the distance between two H-bonded water molecules. The size effect prevails because it is the molecular cause of the markedly larger number density of water: at 25 • C and 1 atm, ρ (in moles per liter) = 55.3 for water, 10.3 for CCl 4 , and 11.2 for benzene [37]. The magnitude of the solvent-excluded volume effect associated with cavity creation depends strongly upon the liquid number density: the entropy loss is larger with the greater the number of affected molecules. This is why the size is so important. The reliability of using classic SPT for water has been further confirmed recently by the agreement between the ∆G C (SPT) values and those obtained by computer simulations in detailed water models [39].
Moreover, a simple formula devised by Pierotti [29] allows for the calculation of the interaction energy E a between the noble gases and the three liquids. The Pierotti's formula is: where σ 12 = (σ 1 + σ 2 )/2 and ε 12 = (ε 1 ε 2 ) 1/2 , and where ε 1 and ε 2 are the Lennard-Jones parameters for the liquid and solute, respectively. The E a estimates, listed in the ninth column of Table 1, are negative and, when added to the ∆G C (SPT) values, produce numbers that are close to the experimental ∆G · ones for all of the three liquids. The success is mainly because the solvent-excluded volume effect associated with solute insertion in a liquid is correctly accounted for by calculating the reversible work of cavity creation.
Estimates of the enthalpy change associated with cavity creation, ∆H C (SPT), calculated by means of Equation (2) and listed in the third column of Table 2, are positive in all of the three liquids and are close to the values of the ∆H · − E a difference, listed in the fourth column of Table 2. This suggests that the structural reorganization of liquid molecules upon noble gas insertion is an endothermic process at 25 • C and 1 atm (i.e., in the case of water, there is no indication of iceberg formation [40][41][42][43]). Actually, the ∆H C (SPT) values of water are significantly smaller than those of the other two liquids [37]; this is a consequence of the smaller isobaric thermal expansion coefficient αP of water with respect to those of the other two liquids [37] (look at the values reported in the notes of Table 1). The latter quantity, present in the classic SPT formula of ∆H C [29], is a measure of the ensemble correlation between fluctuations in volume and fluctuations in enthalpy, and so it can account for the liquid structural reorganization upon cavity creation. The smallness of the α P of water is due to the strength of water-water H-bonds, in comparison to the weakness of van-der-Waals-type interactions occurring among benzene and carbon tetrachloride molecules [37]. Therefore, the ∆H C (SPT) values indicate that cavity creation does not cause the breakage of water-water H-bonds [30,31], but a significant breakage of van der Waals interactions occurs in the other two liquids [37]. Table 2. Enthalpy and entropy changes associated with cavity creation in water (a), CCl 4 (b), and C 6 H 6 (c), calculated by means of the classic SPT relationships at 25 • C and 1 atm, to be compared with the reorganization enthalpy change and the total solvation entropy change, respectively. Estimates of the entropy change associated with cavity creation, ∆S C (SPT), listed in the fifth column of Table 2, are close to the total solvation entropy changes, listed in the last column of Table 2, in all of the three considered liquids. This finding indicates that the process of cavity creation is the main process responsible of the negative solvation entropy change [31,37]. In water, the ∆S C (SPT) values are largely negative, increasing in magnitude with the solute diameter [37]. This entropy loss cannot be due to an increase in water structural order [44,45], because it comes from a hard sphere approach. It is due to the decrease in the number of accessible configurations for water molecules because of cavity creation (i.e., the solvent-excluded volume effect). Such a decrease in the number of accessible configurations also occurs in the other two liquids, but it is masked by a largely positive entropy change due to the structural reorganization upon cavity creation [25,37]. The latter structural reorganization, however, has a markedly different magnitude in water and the two organic liquids; it is also characterized by a complete enthalpy-entropy compensation in all liquids [31,46] and does not affect the ∆G C magnitude.

Formation of Host-Guest Complexes
It is interesting that noble gases are able to bind macrocyclic hosts in aqueous solutions. In particular, thanks to specialized NMR experiments, it has been possible to measure the binding constants of noble gases to cucurbit [5]uril, a rigid, synthetic, and water-soluble macrocyclic host [47]. Specifically, at 22 • C, K b (in M −1 ) = 87 for He, 72 for Ne, 360 for Ar, 2390 for Kr and 8700 for Xe. These numbers imply that the binding process is spontaneous, and thus, there is the need to identify the driving force of this host-guest recognition [47][48][49]. It is important to underscore that the inner part of cucurbit [5]uril proved not to be filled by water molecules, on the basis of both specialized NMR measurements and MD simulations [47] (note that the inner part of cucurbit [5]uril has a volume of 68 Å 3 and can host very few water molecules, considering that the van der Waals volume of a water molecule is 11.5 Å 3 ). Researchers calculated with great accuracy, at DFT level, the dispersion energetic attractions of noble gases in bulk water and in the inner part of cucurbit [5]uril. The unexpected result was that the magnitude of attractive dispersion interactions was larger in bulk water than in the inner part of the rigid macrocyclic host [47]. As a consequence, researchers turned their attention to the reversible work of cavity creation. The transfer of a noble gas atom from water to the inner part of cucurbit [5]uril implies the following steps: the switching-off of the energetic dispersion attractions with water, the closure of the cavity in water, the creation of the cavity in the macrocyclic host interior, and the switching-on of the energetic dispersion attractions with the host.
However, the reversible work to create a cavity in the inner part of cucurbit [5]uril is zero because this region does not contain water molecules (i.e., it is empty); in addition, as a first approximation, the magnitude of the energetic dispersion attractions of a noble gas atom in bulk water and in the host interior can be assumed to be equal. This implies ∆G(binding) ≈ −∆G C (H 2 O), and the driving force is given by the decrease in solventexcluded volume for cavity closure in water (i.e., leading to a gain in configurationaltranslational entropy of water molecules).
The experimental ∆G(binding) values of noble gases to cucurbit [5]uril are reported in Figure 2, together with the values of minus ∆G C (H 2 O), calculated via classic SPT formulas and listed in Table 1. One could say that the agreement between the two sets of numbers is better than expected, considering their totally different origin. In the original article, the authors calculated ∆G C (H 2 O) by means of computer simulations, they also considered the contribution of the difference in energetic dispersion attractions between the bulk water and the host interior, and obtained a good agreement with the experimental data [47]. This example demonstrates the pivotal role played by the reversible work of cavity creation in driving host-guest recognition phenomena in water [50][51][52].  [5]uril, measured via NMR at 22 • C (black filled squares) [47], contrasted with minus the ∆G C values for noble gases in water, calculated via classic SPT analytical formulas, and listed in column eight of Table 1 (red filled circles).

Conformational Stability of Globular Proteins
The geometric explanation for the solvent-excluded volume effect implies that the ∆G C magnitude has to increase if the cavity shape is changed, by keeping its van der Waals volume fixed, V vdW , and increasing its WASA. This is a fundamental point. Passing from a spherical cavity to several prolate spherocylinders with the same V vdW of the sphere, it is possible to test the rightness of the geometric arguments. Classic SPT analytical formulas allow for the calculation of ∆G C for both spherical and prolate spherocylindrical cavities. Therefore, the test can readily be completed and the results have confirmed that ∆G C increases with cavity WASA, even though V vdW is fixed [24,53]. This holds true also with ∆G C calculated by means of computer simulations [54,55]. The results of classic SPT calculations in water, at 25 • C and 1 atm, for two sets of cavities, the first starting with a sphere of 6 Å radius and the second starting with a sphere of 9 Å radius, are listed in Table 3. It is evident that on lengthening the prolate spherocylinder, WASA increases and ∆G C also increases. The plot of ∆G C versus WASA, constructed with the numbers of Table 3, is shown in Figure 3. The ∆G C values scale linearly with cavity WASA, but the line slope is not unique; the slope magnitude depends upon the V vdW of the cavity. In fact, the largest spherocylinder of the first set has a WASA larger than that of the smallest spherocylinder of the second set, but the order is reversed in the case of ∆G C values (see Table 3 and Figure 3). This means also that the cavity V vdW plays a role [53,[56][57][58]. Table 3. ∆G C estimates associated with the creation of prolate spherocylindrical cavities, at 25 • C and 1 atm, in a hard sphere fluid, with the experimental density of water and particle diameter σ = 2.8 Å. By keeping the cavity V vdW fixed at the volume of 6 Å and 9 Å radius spheres, respectively (i.e., 904.8 Å 3 and 3053.6 Å 3 , respectively), the ∆G C numbers have been calculated on varying the cylindrical length by means of the classic SPT analytic formulas. The first row of the A and B sections contains the numbers for the two spherical cavities. The geometric formulas for a prolate spherocylinder of radius a and cylindrical length l are: V vdW = (4/3)π·a 3 + π·l·a 2 and WASA = 4π(a + r w ) 2 + 2π·l·(a + r w ), where r w is the radius of water molecules, fixed at 1.4 Å; by setting l = 0, such formulas become right for a sphere of radius a.  Anyway, the trend of ∆G C versus cavity WASA is important to shed light on the driving force of protein folding and on the main factor responsible for the conformational stability of globular proteins. Experimental measurements have proved that the difference in molecular volume between the folded state and the unfolded state ensemble is negligibly small [59,60]. Thus, the folding process can be viewed as a collapse from a set of elongated conformations toward a compact, almost spherical one, keeping the volume occupied by the polypeptide chain fixed [24,25]. Such a collapse is characterized by a large WASA decrease; that means a large ∆G C decrease, which corresponds to a significant gain in the configurational-translational entropy of water molecules. The numbers listed in the last column of Table 3 indicate that a large negative Gibbs free energy change is associated, at 25 • C and 1 atm, with the collapse from the longest spherocylinder to the sphere. Polypeptide chains are flexible and can populate different conformations, producing markedly different solvent-excluded volume effects. Water molecules push these chains to assume compact conformations in order to gain configuration-translational entropy. This is the geometric-molecular basis of what is called the hydrophobic effect, considered to be the main determinant of the conformational stability of globular proteins [24,25].

Conclusions
In the present article, I have tried to show that the solvent-excluded volume effect associated with cavity creation in all liquids (that are a condensed state of the matter) allows one to devise a common and general theoretical approach to rationalize several disparate phenomena occurring in liquids. In particular, the ∆G C (SPT) values are able to rationalize the low solubility of noble gases in water and its entropic origin, the driving force of the recognition between noble gases and cucurbit [5]uril in water, and, last but not least, a reliable driving force for protein folding and stability.