Free energy decoupling

From (2.70), it follows that the free energy cannot be divided simply into two terms, associated with the interactions of type a and type b. There are also coupling terms, which would vanish only if fluctuations in AUa and AUb were uncorrelated. One might expect that such a decoupling could be accomplished by carrying out the transformations that involve interactions of type a and type 6 separately. In Sect. 2,8.4, we have already discussed such a case for electrostatic and van der Waals interactions in the context of single-topology alchemical transformations. Even then, however, correlations between these two types of interactions are not... 

Hamelberg, D. McCammon, J. A., Standard free energy of releasing a localized water molecule from the binding pockets of proteins Double-decoupling method, J. Am. Chem. Soc. 2004,126, 7683-7689. 

Rao and Singh32 calculated relative solvation free energies for normal alkanes, tetra-alkylmethanes, amines and aromatic compounds using AMBER 3.1. Each system was solvated with 216 TIP3P water molecules. The atomic charges were uniformly scaled down by a factor of 0.87 to correct the overestimation of dipole moment by 6-31G basis set. During the perturbation runs, the periodic boundary conditions were applied only for solute-solvent and solvent-solvent interactions with a non-bonded interaction cutoff of 8.5 A. All solute-solute non-bonded interactions were included. Electrostatic decoupling was applied where electrostatic run was completed in 21 windows. Each window included 1 ps of equilibration and 1 ps of data... 

There is a discussion in the literature about the effect of undulation entropy on the equilibrium membrane tension [14,15], Formally, undulations are included in the surface tension, and thus we need not worry about this. However, if in some model the two are artificially decoupled, one may allow for a very small (positive) surface tension as the equilibrium structure. In other words, the entropy (per unit area) from undulations should compensate for the tension (excess free energy per unit area). 

Previous studies have shown that there is a correlation between the enthalpy of hydration of alkanes and their accessible surface area [30,31] or related magnitudes. Moreover, relationships between the hydration numbers calculated from discrete simulations for hydrocarbons and both the free energy and enthalpy of hydration of these molecules have also been reported [32] and have been often used to evaluate solvation enthalpies. Analysis of our results, illustrates the existence of a linear relationship between A//n eie and the surface of the van der Waals cavity,. SVw, defined in MST computations for the calculation of the non-electrostatic contributions (Figure 4-1). In contrast, no relationship was found for the electrostatic component of the hydration enthalpy (A//eie data not shown). Clearly, in a first approximation, one can assume that the electrostatic interactions between solute and solvent can be decoupled from the interactions formed between uncharged solutes and solvent molecules. 

---