A new set of bond-level molecular descriptors (bond-based linear indices) are used here in QSAR (quantitative structure–activity relationship) studies of tyrosinase inhibitors, for finding functions that discriminate between the tyrosinase inhibitor compounds and inactive ones. A database of 246 compounds was collected for this study; all organic chemicals were reported as tyrosinase inhibitors; they had great structural diversity. This dataset can be considered as a helpful tool, not only for theoretical chemists but also for other researchers in this area. The set used as inactive has 412 drugs with other clinical uses. Twelve LDA-based QSAR models were obtained, the first six using the non-stochastic total and local bond-based linear indices as well as the last six ones, the stochastic molecular descriptors. The best two discriminant models computed using the non-stochastic and stochastic molecular descriptors (Eqs. 7 and 13, respectively) had globally good classifications of 98.95% and 89.75% in the training set, with high Matthews correlation coefficients ( C) of 0.98 and 0.78. The external prediction sets had accuracies of 98.89% and 89.44%, and ( C) values of 0.98 and 0.78, for models 7 and 13, respectively. A virtual screening of compounds reported in the literature with such activity was carried out, to prove the ability of present models to search for tyrosinase inhibitors, not included in the training or test set. At the end, the fitted discriminant functions were used in the selection/identification of new ethylsteroids isolated from herbal plants, looking for tyrosinase inhibitory activity. A good behavior is shown between the theoretical and experimental results on mushroom tyrosinase enzyme. It might be highlighted that all the compounds showed values under 10 μM and that ES2 (IC 50 = 1.25 μM) showed higher activity in the inhibition against the enzyme than reference compounds kojic acid (IC 50 = 16.67 μM) and l-mimosine (IC 50 = 3.68 μM). In addition, a comparison with other established methods was carried to prove the adequate discriminatory performance of the molecular descriptors used here. The present algorithm provided useful clues that can be used to speed up in the identification of new tyrosinase inhibitor compounds.