An effective free energy potential, developed originally for binding free energy calculation, is compared to calorimetric data on protein unfolding, described by a linear combination of changes in polar and nonpolar surface areas. The potential consists of a molecular mechanics energy term calculated for a reference medium (vapor or nonpolar liquid), and empirical terms representing solvation and entropic effects. It is shown that, under suitable conditions, the free energy function agrees well with the calorimetric expression. An additional result of the comparison is an independent estimate of the side-chain entropy loss, which is shown to agree with a structure-based entropy scale. These findings confirm that simple functions can be used to estimate the free energy change in complex systems, and that a binding free energy evaluation model can describe the thermodynamics of protein unfolding correctly. Furthermore, it is shown that folding and binding leave the sum of solute-solute and solute-solvent van der Waals interactions nearly invariant and, due to this invariance, it may be advantageous to use a nonpolar liquid rather than vacuum as the reference medium.Keywords: binding free energy; free energy calculation; heat capacity; protein unfolding Calculating the free energy change in protein folding and association is a classical problem in biophysical chemistry. In principle, free energy differences can be obtained by molecular dynamics and Monte Carlo simulations that allow for similar molecules to be interconverted, and the relative free energies determined by perturbation or integration techniques (Mezei & Beveridge, 1986; Reynolds et al., 1992). However, simulation methods are far too expensive computationally for free energy calculation in conformational search, docking, and design (Wilson et al., 1991). The simplest remedy is to neglect solvation and entropic contributions in applications, but it is well known that energy-type target functions are frequently unable to distinguish between correct and incorrect proteins folds (Novotny et al., 1988), or correct and incorrect docked conformations (Shoichet & Kuntz, 1991 becoming the standard in computer-aided molecular design (Ajay & Murcko, 1995).We have developed a relatively complete empirical free energy function, and evaluated it against a range of structural and thermodynamic data (Vajda et al., Gulukota et al., 1996;King et al., 1996; Weng et ai., 1996). The free energy change, AG = G2 -G I , between two states is calculated according to the expression:where AE, Act,, and AS, represent the energy change, the desolvation free energy, and the change in conformational entropy, respectively. The last term, AC,,t/lrrr includes all other free energy changes associated with translational, rotational, vibrational, cratic, and protonation/deprotonation effects (Novotny et al., 1989; Vajda et al., 1994).The function has been developed originally for calculating receptor-ligand binding free energies, and we used a number of simplifying assumptions t...