Protein engineering to improve the thermostability of glucoamylase from Aspergillus awamori based on molecular dynamics simulations

Twelve mutations were constructed to improve the thermostability of glucoamylase from Aspergillus awamori based on the results of molecular dynamics simulations. The thermal unfolding of the catalytic domain followed a putative hierarchical behavior. In addition, the unfolding of the 13 alpha-helices obeyed the random ordered mechanism, in which the alpha-helices 8, 1 and 11 unfolded more rapidly than the others. The catalytic center was well protected by the (alpha/alpha)(6)-barrel at simulation temperatures up to 600 K, whereas the catalytic base, E400, migrated from its original interior pocket to the surface of the catalytic domain by surmounting the hydrophobic barrier provided by alpha-helices 12 and 13 at 800 K. The disulfide bonds engineered to 'lock' the alpha-helix 11 on the surface of the catalytic domain dramatically increased the thermostability. Substituting G396 and G407 with Ala residues slightly increased the thermostability, whereas their specific activity and catalytic efficiency were reduced. This indicates that the introduced residues with higher hydrophobicity were favorable in the loop between alpha-helices 12 and 13, whereas they partially destroyed the hydrogen bond and salt linkage network in the catalytic center. alpha-Helices 12 and 13 can be stabilized by introducing residues with higher hydrophobicity, except for the H391M mutation.