Higher thermostability is required for 1,3-1,4-β-glucanase to maintain high activity under harsh conditions in the brewing and animal feed industries. In this study, a comprehensive and comparative analysis of thermostability in bacterial β-glucanases was conducted through a method named spatial compartmentalization of mutational hotspots (SCMH), which combined alignment of homologous protein sequences, spatial compartmentalization, and molecular dynamic (MD) simulation. The overall/local flexibility of six homologous β-glucanases was calculated by MD simulation and linearly fitted with enzyme optimal enzymatic temperatures. The calcium region was predicted to be the crucial region for thermostability of bacterial 1,3-1,4-β-glucanases, and optimization of four residue sites in this region by iterative saturation mutagenesis greatly increased the thermostability of a mesophilic β-glucanase (BglT) from Bacillus terquilensis. The E46P/S43E/H205P/S40E mutant showed a 20 °C increase in optimal enzymatic temperature and a 13.8 °C rise in protein melting temperature (T m) compared to wild-type BglT. Its half-life values at 60 and 70 °C were 3.86-fold and 7.13-fold higher than those of wild-type BglT. The specific activity of E46P/S43E/H205P/S40E mutant was increased by 64.4 %, while its stability under acidic environment was improved. The rational design strategy used in this study might be applied to improve the thermostability of other industrial enzymes.
Keywords: 1,3-1,4-β-Glucanase; Iterative saturation mutagenesis; Molecular dynamics simulation; Spatial compartmentalization; Thermostability.