In the current paper we present a receptor-independent 4D-QSAR method based on self-organizing mapping (SOM-4D-QSAR) and in particular focus on its pharmacophore mapping ability. We use a novel stochastic procedure to verify the predictive ability of the method for a large population of 4D-QSAR models generated. This systematic study was conducted on a series of benzoic acids, azo dyes, and steroids that bind aromatase. We show that the 4D-QSAR method coupled with IVE-PLS provides a very stable and predictive modeling technique. The method enables us to identify the molecular motifs contributing the most to the fiber-dye affinity and the aromatase enzyme binding activity of the steroid. However, the method appeared much less effective for the benzoic acid series, in which the efficacy was limited by electronic effects strictly correlated to a single conformer.