Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for 67 molecules of 2-amino-benzothiazole-6-anilide derivatives against lymphocyte-specific protein tyrosine kinase (P56 LCK). The molecular field analysis (MFA) and receptor surface analysis (RSA) were employed for QSAR studies and the predictive ability of the model was validated by 15 test set molecules. Structure-based investigations using molecular docking simulation were performed with the crystal structure of P56 LCK. Good correlation between predicted fitness scores versus observed activities was demonstrated. The results suggested that the nature of substitutions at the 2-amino and 6-anilide positions were crucial in enhancing the activity, thereby providing new guidelines for the design of novel P56 LCK inhibitors.