Heparanase (Hpse) is an endo-β-D-glucuronidase capable of cleaving heparan sulfate side chains. Its upregulated expression is implicated in tumor growth, metastasis and angiogenesis, thus making it an attractive target in cancer therapeutics. Currently, a few small molecule inhibitors have been reported to inhibit Hpse, with promising oral administration and pharmacokinetic (PK) properties. In the present study, a ligand-based pharmacophore model was generated from a dataset of well-known active small molecule Hpse inhibitors which were observed to display favorable PK properties. The compounds from the InterBioScreen database of natural (69,034) and synthetic (195,469) molecules were first filtered for their drug-likeness and the pharmacophore model was used to screen the drug-like database. The compounds acquired from screening were subjected to molecular docking with Heparanase, where two molecules used in pharmacophore generation were used as reference. From the docking analysis, 33 compounds displayed higher docking scores than the reference and favorable interactions with the catalytic residues. Complex interactions were further evaluated by molecular dynamics simulations to assess their stability over a period of 50 ns. Furthermore, the binding free energies of the 33 compounds revealed 2 natural and 2 synthetic compounds, with better binding affinities than reference molecules, and were, therefore, deemed as hits. The hit compounds presented from this in silico investigation could act as potent Heparanase inhibitors and further serve as lead scaffolds to develop compounds targeting Heparanase upregulation in cancer.
Keywords: Heparanase; MM/PBSA; binding free energy calculations; molecular docking; molecular dynamics simulations; pharmacophore modeling; virtual screening.