RESUMEN
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.
Asunto(s)
Resistencia a Antineoplásicos/efectos de los fármacos , Metotrexato/farmacología , Tetrahidrofolato Deshidrogenasa/metabolismo , Área Bajo la Curva , Sitios de Unión , Dominio Catalítico , Antagonistas del Ácido Fólico/química , Antagonistas del Ácido Fólico/metabolismo , Antagonistas del Ácido Fólico/farmacología , Antagonistas del Ácido Fólico/uso terapéutico , Humanos , Enlace de Hidrógeno , Ligandos , Metotrexato/análogos & derivados , Metotrexato/metabolismo , Metotrexato/uso terapéutico , Simulación del Acoplamiento Molecular , Simulación de Dinámica Molecular , Mutagénesis Sitio-Dirigida , Neoplasias/tratamiento farmacológico , Neoplasias/patología , Curva ROC , Tetrahidrofolato Deshidrogenasa/química , Tetrahidrofolato Deshidrogenasa/genética , TermodinámicaRESUMEN
: Mechanistically, neurotoxic insults provoke Ca2+-mediated calpain activation, which cleaves the cytoplasmic region of membrane-embedded p35 and produces its truncated form p25. Upon physical interaction, cyclin-dependent kinase 5 (Cdk5) and p25 forms hyperactivated Cdk5/p25 complex and causes severe neuropathological aberrations including hyperphosphorylated tau-mediated neurofibrillary tangles formation, Alzheimer's symptoms, and neuronal death. Therefore, the inhibition of Cdk5/p25 complex may relieve p-tau-mediated Alzheimer's pathology. Herein, computational simulations have identified pyrrolidine-2,3-dione derivatives as novel inhibitors of Cdk5/p25 complex. A ligand-based pharmacophore was designed and employed as 3D query to retrieve drug-like molecules from chemical databases. By molecular docking, drug-like molecules obtaining dock score > 67.67 (Goldcore of the reference compound) were identified. Molecular dynamics simulation and binding free energy calculation retrieved four pyrrolidine-2,3-dione derivatives as novel candidate inhibitors of Cdk5/p25. The root means square deviation of Cdk5/p25 in complex with candidate inhibitors obtained an average value of ~2.15 Å during the 30 ns simulation period. Molecular interactions analysis suggested that each inhibitor occupied the ATP-binding site of Cdk5/p25 and formed stable interactions. Finally, the binding free energy estimation suggested that each inhibitor had lowest binding energy than the reference compound (-113.10 kJ/mol) to recapitulate their strong binding with Cdk5/p25. Overall, these inhibitors could mitigate tau-mediated Alzheimer's phenotype.
RESUMEN
Dihydrofolate reductase (DHFR) is an essential cellular enzyme and thereby catalyzes thereduction of dihydrofolate to tetrahydrofolate (THF). In cancer medication, inhibition of humanDHFR (hDHFR) remains a promising strategy, as it depletes THF and slows DNA synthesis and cellproliferation. In the current study, ligand-based pharmacophore modeling identified and evaluatedthe critical chemical features of hDHFR inhibitors. A pharmacophore model (Hypo1) was generatedfrom known inhibitors of DHFR with a correlation coefficient (0.94), root mean square (RMS)deviation (0.99), and total cost value (125.28). Hypo1 was comprised of four chemical features,including two hydrogen bond donors (HDB), one hydrogen bond acceptor (HBA), and onehydrophobic (HYP). Hypo1 was validated using Fischer's randomization, test set, and decoy setvalidations, employed as a 3D query in a virtual screening at Maybridge, Chembridge, Asinex,National Cancer Institute (NCI), and Zinc databases. Hypo1-retrieved compounds were filtered byan absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment test andLipinski's rule of five, where the drug-like hit compounds were identified. The hit compounds weredocked in the active site of hDHFR and compounds with Goldfitness score was greater than 44.67(docking score for the reference compound), clustering analysis, and hydrogen bond interactionswere identified. Furthermore, molecular dynamics (MD) simulation identified three compounds asthe best inhibitors of hDHFR with the lowest root mean square deviation (1.2 Å to 1.8 Å), hydrogenbond interactions with hDHFR, and low binding free energy (-127 kJ/mol to -178 kJ/mol). Finally,the toxicity prediction by computer (TOPKAT) affirmed the safety of the novel inhibitors of hDHFRin human body. Overall, we recommend novel hit compounds of hDHFR for cancer and rheumatoidarthritis chemotherapeutics.
RESUMEN
Prolyl oligopeptidase (POP) is a potential therapeutic target for treatment of several neurological disorders and α-synucleinopathies including Parkinson's disease. Most of the known POP inhibitors failed in the clinical trials due to poor pharmacokinetic properties and blood-brain impermeability. Therefore, a training set of 30 structurally diverse compounds with a wide range of inhibitory activity against POP was used to generate a quantitative pharmacophore model, Hypo 3, to identify potential POP inhibitors with desirable drug-like properties. Validations through test set, cost analysis, and Fisher's randomization methods proved that Hypo 3 accurately predicted the known inhibitors among inactive compounds. Hypo 3 was employed as 3D query for virtual screening on an in-house drug-like chemical database containing compounds with good brain permeability and ADMET parameters. Database screening with Hypo 3 resulted in 99 compounds that were narrowed down to 21 compounds through molecular docking. Among them, five compounds were identified in our earlier studies, while two compounds showed in vitro POP inhibition. The current study proposed new 16 virtually screened compounds as potential inhibitors against POP that possess Gold docking score in the range of 64.61-75.74 and Chemscore of -32.25 to -38.35. Furthermore, the top scoring four hit compounds were subjected to molecular dynamics simulations to reveal their appropriate binding modes and assessing binding free energies. The hit compounds interacted with POP effectively via hydrogen bonds with important active site residues along with hydrophobic interactions. Moreover, the hit compounds had key inter-molecular interactions and better binding free energies as compared to the reference inhibitor. A potential new hydrogen bond interaction was discovered between Hit 2 with the Arg252 residue of POP. To conclude, we propose four hit compounds with new structural scaffolds against POP for the lead development of POP-based therapeutics for neurological disorders.