Abstract: Despite the efficacy of imatinib therapy in chronic myelogenous leukemia, the development of drug-resistant Abl mutants, especially the most difficult overcoming T315I mutant, makes the search for new Abl T315I inhibitors a very interesting challenge in medicinal chemistry. In this work, a multistep computational framework combining the three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, molecular dynamics (MD) simulation and binding free energy calculation, was performed to explore the structural requirements for the Abl T315I activities of benzimidazole/benzothiazole derivatives and the binding mechanism between the inhibitors and Abl T315I. The established 3D-QSAR models exhibited satisfactory internal and external predictability. Docking study elucidated the comformations of compounds and the key amino acid residues at the binding pocket, which were confirmed by MD simulation. The binding free energies correlated well with the experimental activities. The MM-GBSA energy decomposition revealed that the van der Waals interaction was the major driving force for the interaction between the ligands and Abl T315I. The hydrogen bond interactions between the inhibitors and Met318 also played an important role in stablizing the binding of compounds to Abl T315I. Finally, four new compounds with rather high Abl T315I activities were designed and presented to experimenters for reference.