A mouse monoclonal anti-2,4,6-trinitrophenyl IgE (clone Lb4) was screened with a random set of over 2000 compounds, and several ligands were found to bind with affinities comparable to that of the immunizing hapten (KD in the microM range). An automated docking algorithm was used for the prediction of complex structures formed by 2,4-dinitrophenyl (DNP) and non-DNP ligands in the fragment variable region of IgE(Lb4). All ligands were found to dock in an L-shaped cavity of 15 x 16 x 10 A, surrounded by complementary-determining regions L1, L3, H2 and H3. The ligands were found to occupy the same binding site in different orientations. For rigid ligands the most stable orientation could be predicted with high probability, based on the calculated energy of binding and the occurrence frequencies of identical complexes obtained by repeated simulations. The localization of a flexible ligand (cycrimine-R) was more ambiguous, but it still docked in the same site. The results support a model for heteroligating antibody (Ab) binding sites, where different ligands utilize the total set of available contacts in different combinations. It is suggested that although pseudoenergies calculated by the docking algorithm do not correlate with experimentally measured binding energies, the screening-and-docking procedure can be useful for the mapping of Ab and other receptor binding sites ligating small molecules.