In silico approach to elucidate factors associated with GH1 β-Glucosidase thermostability