We have developed a method for calculating the association energy of quaternary complexes starting from their atomic coordinates. The association energy is described as the sum of two solvation terms and an energy term to account for the loss of translational and rotational entropy. The calculated solvation energy, using atomic solvation parameters and the solvent accessible surface areas, has a correlation of 96% with experimentally determined values. We have applied this methodology to examine intermediates in viral assembly and to assess the contribution isomerization makes to the association energy of molecular complexes. In addition, we have shown that the calculated association can be used as a predictive tool for analyzing modeled molecular complexes.Keywords: hydrophobicity; protein structure; solvation Specific interactions between macromolecules are responsible for the assembly of complex biological structures and are essential to the regulation of events within a cell or organism. The association of molecules to form higher ordered oligomers is in many respects analogous to the block condensation model for protein folding where prefolded units associate to form higher order structures (Richmond & Richards, 1978). To understand the structural basis of recognition we must be able to relate solution measurements of the association process to the structure of the macromolecular complex. This paper considers the problem of calculating the free energies of forming protein complexes from preformed subunits as derived from crystallographic data and relating these values to experimentally obtained association constants. Kauzmann (1959) suggested that a major factor in the stable formation of protein complexes is a consequence of the hydrophobic effect. Using an empirical correlation between the accessible surface area and free energies of transfer of amino acids from water to octanol, Chothia and Janin (1975) found that the free energy required to form a stable complex was directly related to the amount of surface area buried in the interface. Eisenberg and McLachlan (1986) recognized that it is an oversimplificaReprint requests to: Mitchell Lewis, The Johnson Research Foundation, Department of Biochemistry and Biophysics, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104. tion to base the energy of association on surface area alone; polarity and charge must also be considered. By introducing five atomic solvation parameters, to account for the polar or apolar character of each atom type most frequently found in proteins, they could more accurately relate surface area to the free energy of transfer. Moreover, the solvation energy was shown to be useful for assessing protein stability. We assume that the forces that govern the association between two molecules are the same as the forces that are responsible for the folding of a protein in water. As such, the solvation energy should be equally useful as a gauge for evaluating the association energies of quaternary structu...