Abstract
The atomic environments of 170 phenylalanine-residue aromatic rings from 28 protein crystal structures are transformed into a common orientation and combined to calculate an average three-dimensional environment. The spatial distribution of atom types in this environment reveals a preferred interaction between oxygen atoms and the edge of the planar aromatic rings. From the difference in frequency of interaction of oxygen atoms with the edge and the top of the ring, an apparent net free energy difference of interaction favoring the edge of the ring is estimated to be about -1 kcal/mol (1 cal = 4.184 J). Ab initio quantum mechanical calculations, performed on a model consisting of benzene and formamide, indicate that the observed geometry is stabilized by a favorable enthalpic interaction. Although benzene rings are considered to be nonpolar, the electron distribution is a complex multipole with no net dipole moment. The observed interaction orientation frequencies demonstrate that these multipolar electron distributions, when occurring at the short distances encountered in densely packed protein molecules, are significant determinants of internal packing geometries.