Microtubules play crucial role in process of mitosis and cell proliferation, which have been considered as attractive drug targets for anticancer therapy. The aim of this study was to discover novel and chemically diverse tubulin inhibitors for treatment of cancer. In this investigation, the multilayer virtual screening methods, including common feature pharmacophore model, structure-based pharmacophore model and molecular docking, were developed to screen BioDiversity database with 30,000 compounds. A total of 102 compounds were obtained by the virtual screening, and further filtered by diverse chemical clusters with desired properties and PAINS analysis. Finally, 50 compounds were selected and submitted to the biological evaluation. Among these hits, hits 8 and 30 with novel scaffolds displayed stronger antiproliferative activity on four human tumor cells including Hela, A549, MCF-7, and HepG2. Moreover, the two hits were subsequently submitted to molecular dynamic simulations of 90 ns with the aim of exploring the stability of ligandprotein interactions into the binding pocket, and further probing the mechanism of the interaction between tubulin and hits. The molecular dynamic simulation results revealed there had stronger interactions between tubulin and hits in equilibrium state. Therefore, the hits 8 and 30 have been well characterized as lead compounds for developing new tubulin inhibitors with potential anticancer activity.