β-Glycosidase from Sulfolobus solfataricus (SS-BGL) is a highly effective biocatalyst for the synthesis of compound K (CK) from glycosylated protopanaxadiol ginsenosides. In order to improve the thermal stability of SS-BGL, molecular dynamics simulations were used to determine the residue-level binding energetics of ginsenoside Rd in the SS-BGL-Rd docked complex and to identify the top ten critical contributors. Target sites for mutations were determined using dynamic cross-correlation mapping of residues via the Ohm server to identify networks of distal residues that interact with the key binding residues. Target mutations were determined rationally based on site characteristics. Single mutants and then recombination of top hits led to the two most promising variants SS-BGL-Q96E/N97D/N302D and SS-BGL-Q96E/N97D/N128D/N302D with 2.5-fold and 3.3-fold increased half-lives at 95 °C, respectively. The enzyme activities relative to those of wild-type for ginsenoside conversion were 161 and 116%, respectively..
Keywords: MD simulation; bioconversion; compound K; glycosidase; thermal stability.