To improve the thermostability of α-L-rhamnosidase (r-Rha1), an enzyme previously identified from Aspergillus niger JMU-TS528, multiple arginine (Arg) residues were introduced into the r-Rha1 sequence to replace several lysine (Lys) residues that located on the surface of the folded r-Rha1. Hinted by in silico analysis, five surface Lys residues (K134, K228, K406, K440, K573) were targeted to produce a list of 5 single-residue mutants and 4 multiple-residue mutants using site-directed mutagenesis. Among these mutants, a double Lys to Arg mutant, i.e. K406R/K573R, showed the best thermostability improvement. The half-life of this mutant's enzyme activity increased 3 h at 60 °C, 23 min at 65 °C, and 3.5 min at 70 °C, when compared with the wild type. The simulated protein structure based interaction analysis and molecular dynamics calculation indicate that the thermostability improvement of the mutant K406R-K573R was possibly due to the extra hydrogen bonds, the additional cation-π interactions, and the relatively compact conformation. With the enhanced thermostability, the α-L-rhamnosidase mutant, K406R-K573R, has potentially broadened the r-Rha1 applications in food processing industry.