In silico ligand binding studies of cyanogenic β-glucosidase, dhurrinase-2 from Sorghum bicolor

Dhurrinase, a cyanogenic β-glucosidase from Sorghum bicolor is the key enzyme responsible for the hydrolysis of dhurrin to produce toxic hydrogen cyanide, as a part of plant defence mechanism. Dhurrinase 1 (SbDhr1) and dhurrinase 2 (SbDhr2), two isozymes have been isolated and characterized from S. bicolor. However, there is no information in the literature about the three dimensional (3D) structure of SbDhr2 and molecular interactions involved between the protein and ligand. In this study, the three dimensional structure of SbDhr2 was built based on homology modeling by using the X-ray crystallographic structure of its close homologue SbDhr1 as the template. The generated 3D model was energy minimized and the quality was validated by Ramachndran plot, various bioinformatic tools and their relevant parameters. Stability, folding-unfolding and flexibility of the modeled SbDhr2 was evaluated on the basis of RMSD, radius of gyration (Rg) and RMSF values respectively, obtained through molecular dynamic (MD) simulation. Further, molecular docking was performed with its natural substrate dhurrin, one substrate analogue, three un-natural substrates, and one inhibitor. Analysis of molecular interactions in the SbDhr2-ligand complexes revealed the key amino acid residues responsible to stabilize the ligands within the binding pocket through non-bonded interactions and some of them were found to be conserved (Glu239, Tyr381, Trp426, Glu454, Trp511). Reasonably broader substrate specificity of SbDhr2 was explained through the wider entrance passage observed in comparison to SbDhr1.