The three-dimensional quantitative structure-activity relationship comparative molecular field analysis (3D-QSAR CoMFA) approach was applied to some classes of melatonin (MLT) membrane receptor ligands, with the principal aim of exploring the correlation between their steric features and MT(2)-selective antagonism. Binding data obtained from cloned MT(1) and MT(2) receptor subtypes were used to develop 3D-QSAR models for agonists and for antagonists at the two receptor subtypes, looking for the structural requirements for receptor subtype selectivity. In particular, we superposed the compounds showing antagonist activity, or very low intrinsic activity at the GTPgammaS test, following the hypothesis that the occupation of an additional pocket positioned out of the plane of MLT is one of the major determinants for MT(2) selectivity; the statistical models obtained confirmed this hypothesis. Structure-intrinsic activity relationship studies, applied to a set of compounds homogeneously tested, allowed the identification of the structural features whose modulation shifts the behavior from that of the agonist to that of the antagonist. The pocket out of the plane of MLT was identified as one of the key features for obtaining selective MT(2) antagonists. The reliability of our statistical models was further confirmed by the correct prediction of the pharmacological behavior of some N-substituted melatonin derivatives, which were prepared and tested on cloned receptor subtypes.