A molecular dynamics study of Beta-Glucosidase B upon small substrate binding

Paenibacillus polymyxa β-glucosidase B (BglB), belongs to a GH family 1, is a monomeric enzyme that acts as an exo-β-glucosidase hydrolysing cellobiose and cellodextrins of higher degree of polymerization using retaining mechanism. A molecular dynamics (MD) simulation was performed at 300 K under periodic boundary condition for 5 ns using the complexes structure obtained from previous docking study, namely BglB-Beta-d-glucose and BglB-Cellobiose. From the root-mean-square deviation analysis, both enzyme complexes were reported to deviate from the initial structure in the early part of the simulation but it was stable afterwards. The root-mean-square fluctuation analysis revealed that the most flexible regions comprised of the residues from 26 to 29, 43 to 53, 272 to 276, 306 to 325 and 364 to 367. The radius of gyration analysis had shown the structure of BglB without substrate became more compact towards the end of the simulation compare to other two complexes. The residues His122 and Trp410 were observed to form stable hydrogen bond with occupancy higher than 10%. In conclusion, the behaviour of BglB enzyme towards the substrate binding was successfully explored via MD simulation approaches.