Structure and dynamics analysis of multi-domain putative β-1,4-glucosidase of family 3 glycoside hydrolase (PsGH3) from Pseudopedobacter saltans

Structure and conformational behaviour of a putative β-1,4-glucosidase of glycoside hydrolase family 3 (PsGH3) from Pseudopedobacter saltans was predicted by using in-silico tools. PsGH3 modeled structure constructed using Phyre2 displayed multidomain architecture comprising an N-terminal (β/α)(8)-fold domain followed by (α/β)(6)-sandwich domain, PA14 domain, and a C-terminal domain resembling an immunoglobulin fold. Ramachandran plot displayed 99.3% of amino acids in the allowed region and 0.7% residues in the disallowed region. Multiple sequence alignment (MSA) and structure superposition of PsGH3 with other homologues from GH3 family revealed the conserved residues, Asp274 and Glu624 present in loops LA and LB, respectively originating from N-terminal domain act as catalytic residues. The volume and area calculated for PsGH3 displayed a deep active-site conformation comparable with its homologues, β-1,4-glucosidases (GH3) of Kluyveromyces marxianus and Streptomyces venezuelae. Molecular dynamic (MD) simulation of PsGH3 structure for 80 ns suggested stable and compact structure. Molecular docking studies revealed deeper active site conformation of PsGH3 that could house larger cellooligosaccharides up to 7° of polymerization (DP7). The amino acid residues, Ala86, Leu88, Cys275, Pro483, Phe493, Asn417, Asn491, Pro492, and Leu495 created a binding pocket near the catalytic cleft, crucial for ligand binding. MD simulation of PsGH3 in the presence of cellooligosaccharides, viz., cellobiose and celloheptaose showed stability in terms of RMSD, R(g), and SASA values till 80 ns. The calculation of average number of hydrogen bond (H-bond), interaction energy, and binding free energy confirmed the stronger binding affinity of the larger cellooligosaccharides such as celloheptaose in the binding cavity of PsGH3.