Discovery of anti-colon cancer agents targeting wild-type and mutant p53 using computer-aided drug design

Mutations in the p53 gene are common and occur in over 50% of all cancers, as it is involved in DNA damage repair, cell cycle regulation and apoptosis. Moreover, the p53 gene is mutated in 70% of colon cancers. Therefore, the development of drugs to combat this mutation requires urgent attention. With this in mind, in silico drug design approaches were applied on quinoline derivatives with anticancer activity. In 3D-QSAR study, steric, electrostatic, hydrophobic and H-bond acceptor fields (SEHA) play an important role in prediction and design of new colon cancer compounds. Indeed, the two best CoMSIA/SEHA models with (Q2 = 0.737, R2 = 0.914, R pred 2 = 0.720) and (Q2 = 0.738, R2 = 0.919, R pred 2 = 0.739) show good prediction of human colon carcinoma HCT 116 (p53+/+) and (p53-/-) activities, respectively. Furthermore, the predictive ability and robustness of these models were tested by several validation methods. Molecular docking analyses reveal crucial interactions with the active sites of the p53 protein in both wild type and mutant. Based on these theoretical studies, we designed 10 new compounds with good anticancer activity potential, which were evaluated using ADMET properties. Molecular dynamics simulations were performed to confirm the detailed binding mode of the docking results. Finally, the MM-GBSA based on molecular dynamics simulation confirmed that the designed compounds were able to form stable hydrogen bonding interactions with the crucial residues, which are essential to overcome the p53 mutation in colon cancer.Communicated by Ramaswamy H. Sarma.

Identification of novel discoidin domain receptor 1 (DDR1) inhibitors using E-pharmacophore modeling, structure-based virtual screening, molecular dynamics simulation and MM-GBSA approaches

Dysregulation of the discoidin domain receptor (DDR1), a collagen-activated receptor tyrosine kinase, has been linked to several human cancer diseases including non-small cell lung carcinoma (NSCLC), ovarian cancer, glioblastoma, and breast cancer, in addition to several inflammatory and neurological conditions. Although there are some selective DDR1 inhibitors that have been discovered during the last two decades, a combination of elevated cytotoxicity, kinome selectivity and/or poor DMPK profile has prevented more in-depth studies from being performed. As such, no DDR1 inhibitor has reached clinical investigation to date, forming an urgent need to develop specific DDR1 inhibitor(s) using various drug discovery means. However, the recent discovery of VU6015929, a potent and selective DDR1 kinase inhibitor, with enhanced physiochemical and DMPK properties in addition to its clean kinome profile marked a milestone in the development of DDR1 inhibitors. Herein, VU6015929 was used to construct a 3D e-pharmacophore model which was validated via calculating the difference of score between the active compounds and decoys. The validated e-pharmacophore model was then utilized to screen 20 million drug-like compounds obtained from the freely accessible Zinc database. The generated hits were ranked using high throughput virtual screening technique (HTVS), and the top 8 small molecules were subjected to a molecular docking study and MM-GBSA calculations. Protein-ligand complexes of compounds 1, 2, 3 and the standard compound (VU6015929) were performed for 100 ns and compared with the DDR1 unbound protein state and the DDR1 bound to a co-crystallized ligand. The molecular docking, MD and MM-GBSA outputs revealed compounds 1-3 as potential DDR1 inhibitors, with compound 2 displaying superior binding affinity, comparable binding stability and average binding free energy for the ligand-enzyme complex compared to VU6015929.