Benzodiazepines and their derivatives belong to a category of new psychoactive substances that have been introduced into the continually expanding illicit market. However, there is a notable absence of available pharmacological data for these substances. To gain a deeper understanding of their pharmacology, we employed the Monte Carlo optimization conformation-independent method as a tool for developing QSAR models. These models were built using optimal molecular descriptors derived from both SMILES notation and molecular graph representations. The resulting QSAR model demonstrated robustness and a high degree of predictability, proving to be very reliable. Moreover, we were able to identify specific molecular fragments that exerted both positive and negative effects on binding activity. This discovery paves the way for the swift prediction of binding activity for emerging benzodiazepines, offering a faster and more cost-effective alternative to traditional in vitro/in vivo analyses.