In this research, molecular docking and 3D-QSAR studies were carried out on a series of 79 thiazoloquin(az)olin(on)es as CD38 inhibitors. Based on docking results, four interactions including hydrogen bonding with main chain of GLU-226 (H-M-GLU-226), Van der Waals interactions with side chain of TRP-125 (V-S-TRP-125), TRP-189 (V-S-TRP- 189), and THR-221 (V-S-THR-221) were considered as pharmacological interactions. Active conformation of each ligand was extracted from docking studies and was used for carrying out 3D-QSAR modeling. Comparative molecular field analysis (CoMFA) was performed on CD38 inhibitory activities of these compounds on human and mouse. We developed CoMFA models with five components as optimum models for both data-sets. For human data-set, a model with high predictive power was developed. R2, RMSE, and F-test values for training set of this model were .94, .24, and 179.58, respectively, and R2 and RMSE for its test set were .92 and .32, respectively. The q2 and RMSE values for leave-one-out cross validation test on training set were .78 and .46, respectively, that demonstrate created model is robust. Based on extracted steric and electrostatic contour maps for this model, three inhibitors with pIC50 larger than 8.85 were designed.