Virtual screening, structure based pharmacophore mapping, and molecular simulation studies of pyrido[2,3-d]pyrimidines as selective thymidylate synthase inhibitors.

Human thymidylate synthase is the rate-limiting enzyme in the de novo synthesis of 2'-deoxythymidine-5'-monophosphate. dUMP (pyrimidine) and folate binding site hTS inhibitors showed resistance in colorectal cancer (CRC). In the present study, we have performed virtual screening of the pyrido[2,3-d]pyrimidine database, followed by binding free energy calculations, and pharmacophore mapping to design novel pyrido[2,3-d]pyrimidine derivatives to stabilize inactive confirmation of hTS. A library of 42 molecules was designed. Based on the molecular docking studies, four ligands (T36, T39, T40, and T13) were identified to have better interactions and docking scores with the catalytic sites [dUMP (pyrimidine) and folate binding sites] of hTS protein than standard drug, raltitrexed. To validate efficacy of the designed molecules, we performed molecular dynamics simulation studies at 1000 ns with principal component analysis and binding free energy calculations on the hTS protein, also drug likeness properties of all hits were in acceptable range. Compounds T36, T39, T40, and T13 interacted with the catalytic amino acid (Cys195), an essential amino acid for anticancer activity. The designed molecules stabilized the inactive conformation of hTS, resulting in the inhibition of hTS. The designed compounds will undergo synthesis and biological evaluation, which may yield selective, less toxic, and highly potent hTS inhibitors.Communicated by Ramaswamy H. Sarma.