Structural basis of oligosaccharide processing by glycosaminoglycan sulfotransferases

Heparan sulfate (HS) is a sulfated polysaccharide that plays a key role in morphogenesis, physiology and pathogenesis. The biosynthesis of HS takes place in the Golgi apparatus by a group of enzymes that polymerize, epimerize and sulfate the sugar chain. This biosynthetic process introduces varying degrees of sulfate substitution, which are tightly regulated and directly dictate binding specificity to different cytokines, morphogens and growth factors. Here, we report the use of molecular dynamics simulations to investigate the dynamics of substrate recognition of two glycosaminoglycan (GAG) sulfotransferases, N-deacetylase-N-sulfotransferase and 2-O-sulfotransferase to the HS chain during the biosynthetic process. We performed multiple simulations of the binding of the sulfotransferase domains to both the HS oligosaccharide substrate and sulfate donor, 3'-phosphoadenosine-5'-phosphosulfate. Analysis of extended simulations provide detailed and useful insights into the atomic interactions that are at work during oligosaccharide processing. The fast information matching method was used to detect the enzyme global dynamics and to predict the pairwise contact of residues responsible for GAG-enzyme binding and unbinding. The correlation between HS displacement and the location of the modified GAG chain were calculated, indicating a possible route for HS and heparin during sulfotransferase processing. Our data also show sulfotransferases contain a conserved interspaced positively charged amino acid residues that form a patch which controls the protein-GAG binding equilibrium. Together, our findings provide further understanding on the fine-tuned complex mechanism of GAG biosynthesis. Our findings can also be extrapolated to other systems for calculating rates of protein-GAG binding.