Drug resistance is a core issue in cancer chemotherapy. A known folate antagonist, methotrexate (MTX) inhibits human dihydrofolate reductase (hDHFR), the enzyme responsible for the catalysis of 7,8-dihydrofolate reduction to 5,6,7,8-tetrahydrofolate, in biosynthesis and cell proliferation. Structural change in the DHFR enzyme is a significant cause of resistance and the subsequent loss of MTX. In the current study, wild type hDHFR and double mutant (engineered variant) F31R/Q35E (PDB ID: 3EIG) were subject to computational study. Structure-based pharmacophore modeling was carried out for wild type (WT) and mutant (MT) (variant F31R/Q35E) hDHFR structures by generating ten models for each. Two pharmacophore models, WT-pharma and MT-pharma, were selected for further computations, and showed excellent ROC curve quality. Additionally, the selected pharmacophore models were validated by the Guner-Henry decoy test method, which yielded high goodness of fit for WT-hDHFR and MT-hDHFR. Using a SMILES string of MTX in ZINC15 with the selections of 'clean', in vitro and in vivo options, 32 MTX-analogs were obtained. Eight analogs were filtered out due to their drug-like properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment tests and Lipinski's Rule of five. WT-pharma and MT-pharma were further employed as a 3D query in virtual screening with drug-like MTX analogs. Subsequently, seven screening hits along with a reference compound (MTX) were subjected to molecular docking in the active site of WT- and MT-hDHFR. Through a clustering analysis and examination of protein-ligand interactions, one compound was found with a ChemPLP fitness score greater than that of MTX (reference compound). Finally, a simulation of molecular dynamics (MD) identified an MTX analog which exhibited strong affinity for WT- and MT-hDHFR, with stable RMSD, hydrogen bonds (H-bonds) in the binding site and the lowest MM/PBSA binding free energy. In conclusion, we report on an MTX analog which is capable of inhibiting hDHFR in wild type form, as well as in cases where the enzyme acquires resistance to drugs during chemotherapy treatment.

Non-small cell lung cancer (NSCLC) is a lethal non-immunogenic malignancy and proto-oncogene ROS-1 tyrosine kinase is one of its clinically relevant oncogenic markers. The ROS-1 inhibitor, crizotinib, demonstrated resistance due to the Gly2032Arg mutation. To curtail this resistance, researchers developed lorlatinib against the mutated kinase. In the present study, a receptor-ligand pharmacophore model exploiting the key features of lorlatinib binding with ROS-1 was exploited to identify inhibitors against the wild-type (WT) and the mutant (MT) kinase domain. The developed model was utilized to virtually screen the TimTec flavonoids database and the retrieved drug-like hits were subjected for docking with the WT and MT ROS-1 kinase. A total of 10 flavonoids displayed higher docking scores than lorlatinib. Subsequent molecular dynamics simulations of the acquired flavonoids with WT and MT ROS-1 revealed no steric clashes with the Arg2032 (MT ROS-1). The binding free energy calculations computed via molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) demonstrated one flavonoid (Hit) with better energy than lorlatinib in binding with WT and MT ROS-1. The Hit compound was observed to bind in the ROS-1 selectivity pocket comprised of residues from the β-3 sheet and DFG-motif. The identified Hit from this investigation could act as a potent WT and MT ROS-1 inhibitor.