Fructose 1,6-bisphosphatase (FBPase) has been identified as a drug discovery target for lowering glucose in type 2 diabetes mellitus. In this study, a large series of 105 FBPase inhibitors were studied using a combinational method by 3D-QSAR, molecular docking and molecular dynamics simulations for a further improvement in potency. The optimal 3D models exhibit high statistical significance of the results, especially for the CoMFA results with r(ncv)(2), q(2) values of 0.986, 0.514 for internal validation, and r(pred)(2), r(m)(2) statistics of 0.902, 0.828 statistics for external validation. Graphic representation of the results, as contoured 3D coefficient plots, also provides a clue to the reasonable modification of molecules. (1) Substituents with a proper length and size at the C5 position of the thiazole core are required to enhance the potency; (2) A small and electron-withdrawing group at the C2 position linked to the thiazole core is likely to help increase the FBPase inhibition; (3) Substituent groups as hydrogen bond acceptors at the C2 position of the furan ring are favored. In addition, the agreement between 3D-QSAR, molecular docking and molecular dynamics simulation proves the rationality of the developed models. These results, we hope, may be helpful in designing novel and potential FBPase inhibitors.