A combined physicochemical and 3D technique is used to establish the QSAR of four series of non-benzodiazepines towards BzR. In physicochemical QSAR study, the semi-empirical PM3 based parameters like hardness, electronegativity, electrophilicity index, molar refractivity, heat of formation, solvent assessable surface area and logP were used as descriptors. The heat of formation (HF) and logP are recognized as the most important descriptors for binding affinity of such compounds to BzR. The 3D QSAR study reveals that beta-Carbolines of series "A" and Imidazo[1, 2-alpha]pyrimidines of series "C" show steric bulk interaction, the beta-Carbolines of series "B" and di-substituted purines of series "D" have dominance of electrostatic interaction with BzR. The CoMSIA also indicate the same trend in terms of steric and electrostatic interaction. The CoMSIA reveals that only in case of series "D", the hydrophobic field effect is important. Some possible candidates for all four series were designed and their probable activities were estimated by using different QSAR models. The activities of designed molecules are in better range and developed models might be helpful to design the potent ligands of BzR.