D-Allulose 3-epimerase (DAEase) is a key enzyme in D-allulose bioproduction. DAEase from Thermoclostridium caenicola suffers from poor thermostability, hampering its large-scale applications in industry. In this study, mutants A70P, G107P, F155Y, and D162T with increased melting point temperature (T m ) were obtained by targeted mutagenesis based on the calculation of the free energy of folding. The optimal single-point mutant G107P showed 11.08 h, 5, and 5.70 °C increases in the values of halflife (t 1/2 ) at 60 °C, the optimum temperature (T opt ), and T m , respectively. Beneficial mutations were combined by ordered recombination mutagenesis, and the combinational mutant Var3 (G107P/F155Y/D162T/A70P) was generated with ΔT opt of 10 °C and ΔT m of 12.25 °C. Its t 1/2 value at 65 °C was more than 140 times higher than that of the wild-type enzyme. Molecular dynamics simulations and homology modeling analysis indicated that the enhanced overall rigidity, increased hydrogen bonds between subunits, and redistributed surface electrostatic charges might be responsible for the improved thermostability of the mutant Var3.