Structure and conformational behaviour of a putative β-1,4-glucosidase of glycoside hydrolase family 3 (
Ps
GH3) from
Pseudopedobacter saltans
was predicted by using
in-silico
tools.
Ps
GH3 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
Ps
GH3 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
Ps
GH3 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
Ps
GH3 structure for 80 ns suggested stable and compact structure. Molecular docking studies revealed deeper active site conformation of
Ps
GH3 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
Ps
GH3 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
Ps
GH3.