Modeled 3d-structures of proteobacterial transglycosylases from glycoside hydrolase family 17 give insight in ligand interactions explaining differences in transglycosylation products

The structures of glycoside hydrolase family 17 (GH17) catalytic modules from modular proteins in the ndvB loci in Pseudomonas aeruginosa (Glt1), P. putida (Glt3) and Bradyrhizobium diazoefficiens (previously B. japonicum) (Glt20) were modeled to shed light on reported differences between these homologous transglycosylases concerning substrate size, preferred cleavage site (from reducing end (Glt20: DP2 product) or non-reducing end (Glt1, Glt3: DP4 products)), branching (Glt20) and linkage formed (1,3-linkage in Glt1, Glt3 and 1,6-linkage in Glt20). Hybrid models were built and stability of the resulting TIM-barrel structures was supported by molecular dynamics simulations. Catalytic amino acids were identified by superimposition of GH17 structures, and function was verified by mutagenesis using Glt20 as template (i.e., E120 and E209). Ligand docking revealed six putative subsites (−4, −3, −2, −1, +1 and +2), and the conserved interacting residues suggest substrate binding in the same orientation in all three transglycosylases, despite release of the donor oligosaccharide product from either the reducing (Glt20) or non-reducing end (Glt1, Gl3). Subsites +1 and +2 are most conserved and the difference in release is likely due to changes in loop structures, leading to loss of hydrogen bonds in Glt20. Substrate docking in Glt20 indicate that presence of covalently bound donor in glycone subsites −4 to −1 creates space to accommodate acceptor oligosaccharide in alternative subsites in the catalytic cleft, promoting a branching point and formation of a 1,6-linkage. The minimum donor size of DP5, can be explained assuming preferred binding of DP4 substrates in subsite −4 to −1, preventing catalysis.