Search

Owing to several mutations, the oncogene Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is activated in the majority of cancers, and targeting it has been pharmacologically challenging. In this study, using an in silico approach comprised of pharmacophore modeling, molecular docking, and molecular dynamics simulations, potential KRAS G12D inhibitors were investigated. A ligand-based common feature pharmacophore model was generated to identify the framework necessary for effective KRAS inhibition. The chemical features in the selected pharmacophore model comprised two hydrogen bond donors, one hydrogen bond acceptor, two aromatic rings and one hydrophobic feature. This model was used for screening in excess of 214,000 compounds from InterBioScreen (IBS) and ZINC databases. Eighteen compounds from the IBS and ten from the ZINC database mapped onto the pharmacophore model and were subjected to molecular docking. Molecular docking results highlighted a higher affinity of four hit compounds towards KRAS G12D in comparison to the reference inhibitor, BI-2852. Sequential molecular dynamics (MD) simulation studies revealed all four hit compounds them possess higher KRAS G12D binding free energy and demonstrate stable polar interaction with key residues. Further, Principal Component Analysis (PCA) analysis of the hit compounds in complex with KRAS G12D also indicated stability. Overall, the research undertaken provides strong support for further in vitro testing of these newly identified KRAS G12D inhibitors, particularly Hit1 and Hit2.

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.