In silico and in vitro mapping of specificity patterns of glycosaminoglycans towards cysteine cathepsins B, L, K, S and V

Human cysteine cathepsins are lysosomal proteases, which are involved in different biological processes. Their enzymatic activity can be regulated by glycosaminoglycans (GAGs): long linear periodic negatively charged polysaccharides, which dimeric building blocks consist of uronic acid and hexosamine monosaccharide units. In this study, molecular docking simulations of chondroitin 4-sulfate, chondroitin 6-sulfate, heparin, heparan sulfate, dermatan sulfate and hyaluronic acid of various chain lengths were performed with cathepsins B, L, K, S and V and followed by molecular dynamics-based refinement and binding free energy analysis. We concluded that electrostatics might be a driving force for cathepsin-GAG interactions; indeed as in most of characterised systems, the increase of GAG chain length consequently leads to a more pronounced effect on the strength of cathepsin-GAG interactions. Results also suggest that binding of GAGs at different regions on cathepsins surface affect differently their enzymatic activity and could is dependent on cathepsin and GAG type. Present data contribute to systematic description of cathepsin-GAG interactions, which is helpful in understanding the subtle molecular mechanisms of protease regulation behind their biological functions.