Molecular Dynamics Unveils Multiple-Site Binding of Inhibitors with Reduced Activity on the Surface of Dihydrofolate Reductase.

The sensitivity to protein inhibitors is altered by modifications or protein mutations, as represented by drug resistance. The mode of stable drug binding to the protein pocket has been experimentally clarified. However, the nature of the binding of inhibitors with reduced sensitivity remains unclear at the atomic level. In this study, we analyzed the thermodynamics and kinetics of inhibitor binding to the surface of wild-type and mutant dihydrofolate reductase (DHFR) using molecular dynamics simulations combined with Markov state modeling. A strong inhibitor of methotrexate (MTX) showed a preference for the active site of wild-type DHFR with minimal binding to unrelated (secondary) sites. Deletion of a side-chain fragment in MTX largely destabilized the active site-bound state, with clear evidence of binding to secondary sites. Similarly, the F31V mutation in DHFR diminished the specificity of MTX binding to the active site. These results reveal the presence of multiple-bound states whose stabilities are comparable to or higher than those of the unbound state, suggesting that a reduction in the binding affinity for the active site significantly elevates the fractions of these states. This study presents a theoretical model that more accurately interprets the altered drug sensitivity than the traditional two-state model.

Design, synthesis, in vitro and in silico studies of novel piperidine derived thiosemicarbazones as inhibitors of dihydrofolate reductase.

Dihydrofolate reductase (DHFR), an essential enzyme in folate metabolism, presents a promising target for drug development against various diseases, including cancer and tuberculosis. Herein, we present an integrated approach combining in vitro biochemical assays with in silico molecular docking analysis to evaluate the inhibitory potential of 4-piperidine-based thiosemicarbazones 5(a-s) against DHFR. In our in vitro study, a novel series of 4-piperidine-based thiosemicarbazones 5(a-s) were assessed for their inhibitory activity against DHFR enzyme. The synthesized compounds 5(a-s) exhibited potent inhibition with IC50 values in the range of 13.70 ± 0.25 µM to 47.30 ± 0.86 µM. Among all the derivatives 5p displayed highest inhibitory activity. Simultaneously, in silico analysis were performed and compared with standard drug (Methotrexate) to predict the binding affinity and interaction pattern of synthesized compounds with DHFR active site. SAR analysis was done to elucidate how structural modifications impact compound's biological activity, guiding the rational design of potent and selective drug candidates for targeted diseases. These findings may provide a comprehensive assessment of 4-piperdine-based thiosemicarbazones as DHFR inhibitors and contribute to the development of novel therapeutics targeting DHFR-associated diseases.

Identification of Innovative Folate Inhibitors Leveraging the Amino Dihydrotriazine Motif from Cycloguanil for Their Potential as Anti-Trypanosoma brucei Agents.

Folate enzymes, namely, dihydrofolate reductase (DHFR) and pteridine reductase (PTR1) are acknowledged targets for the development of antiparasitic agents against Trypanosomiasis and Leishmaniasis. Based on the amino dihydrotriazine motif of the drug Cycloguanil (Cyc), a known inhibitor of both folate enzymes, we have identified two novel series of inhibitors, the 2-amino triazino benzimidazoles (1) and 2-guanidino benzimidazoles (2), as their open ring analogues. Enzymatic screening was carried out against PTR1, DHFR, and thymidylate synthase (TS). The crystal structures of TbDHFR and TbPTR1 in complex with selected compounds experienced in both cases a substrate-like binding mode and allowed the rationalization of the main chemical features supporting the inhibitor ability to target folate enzymes. Biological evaluation of both series was performed against T. brucei and L. infantum and the toxicity against THP-1 human macrophages. Notably, the 5,6-dimethyl-2-guanidinobenzimidazole 2g resulted to be the most potent (Ki = 9 nM) and highly selective TbDHFR inhibitor, 6000-fold over TbPTR1 and 394-fold over hDHFR. The 5,6-dimethyl tricyclic analogue 1g, despite showing a lower potency and selectivity profile than 2g, shared a comparable antiparasitic activity against T. brucei in the low micromolar domain. The dichloro-substituted 2-guanidino benzimidazoles 2c and 2d revealed their potent and broad-spectrum antitrypanosomatid activity affecting the growth of T. brucei and L. infantum parasites. Therefore, both chemotypes could represent promising templates that could be valorized for further drug development.

Synthesis, molecular modeling simulations and anticancer activity of some new Imidazo[2,1-b]thiazole analogues as EGFR/HER2 and DHFR inhibitors.

New imidazo[2,1-b]thiazole analogs were designed, synthesized, and biologically evaluated as anticancer agents. In vitro biological evaluation of the anticancer properties of the compounds was performed against different cancer cell lines. Compounds 23 and 39 showed remarkable broad -spectrum cytotoxic potency on most of the tested cell lines. Compounds 23 and 39 exhibited potent activity against the MCF-7 breast cancer cell line, with IC50 values of 1.81 and 4.95 µM, respectively, compared to DOX and SOR (IC50 values of 4.17 and 7.26 µM, respectively). An enzyme inhibition assay was carried out to clarify the possible mode of action of the tested compounds. Compounds 23 and 39 were identified as possible EGFR, HER-2, and DHFR inhibitors. Cell cycle arrest results indicated that compound 23 caused cell cycle arrest at the G0/G1 phase in the MCF-7 cells and at the G2/M phase in the Hep G2 cells. Compound 39 induced cell cycle arrest at the G2/M phase in Hela cells. In vivo testing of the anticancer activity of the two most promising molecules in this study was conducted, and the results indicated that they possess considerable in vivo anticancer activity in mice. Data obtained from the molecular modeling simulation study were consistent with the biological evaluation results.

Synthesis, molecular docking study and biological evaluation of new pyrrole scaffolds as potential antitubercular agents for dual targeting of enoyl ACP reductase and dihydrofolate reductase.

In this study, new series of N'-(2-(substitutedphenoxy)acetyl)-4-(1H-pyrrol-1-yl)benzohydrazides (3a-j) 4-(2,5-dimethyl-1H-pyrrol-1-yl)-N'-(2-(substitutedphenoxy)acetyl)benzohydrazides (5a-j) were synthesized, characterized and assessed as inhibitors of enoyl ACP reductase and DHFR. Most of the compounds exhibited dual inhibition against the enzymes enoyl ACP reductase and DHFR. Several synthesized substances also demonstrated significant antibacterial and antitubercular properties. A molecular docking analysis was conducted in order to determine the potential mechanism of action of the synthesized compounds. The results indicated that there were binding interactions seen with the active sites of dihydrofolate reductase and enoyl ACP reductase. Additionally, important structural details were identified that play a critical role in sustaining the dual inhibitory activity. These findings were useful for the development of future dual inhibitors. Therefore, this study provided strong evidence that several synthesized molecules could exert their antitubercular properties at the cellular level through multi-target inhibition. By shedding light on the mechanisms through which these compounds exert their inhibitory effects, this research opens up promising avenues for the future development of dual inhibitors with enhanced antibacterial and antitubercular properties. The study's findings underscore the importance of multi-target approaches in drug design, providing a strong foundation for the design and optimization of novel compounds that can effectively target bacterial infections at the cellular level.

Synthesis and in vitro antitumor evaluation of new thieno[2,3-d]pyrimidine derivatives as EGFR and DHFR inhibitors.

New thienopyrimidine derivatives 2-16 have been synthesized and their in vitro cytotoxicity was evaluated against five different human cancer cell lines HCT-116, Hela, MDA-MB-231, MCF7 and PC3. Compounds 6e, 7a, 7b, 7d, 10c and 10e displayed the highest antitumor activity against all tested cell lines compared to Doxorubicin. Enzyme inhibition assay revealed that compounds 6e and 10e showed high inhibitory activity against EGFR-TK, with IC50 values of 0.133 and 0.151 µM, compared to Olmutinib (IC50 = 0.028 µM); while the highest DHFR inhibitory activity was shown by compounds 7d and 10e with IC50 values of 0.462 and 0.541 µM, compared to Methotrexate (IC50 = 0.117 µM). Cell cycle analysis following a flow cytometric study using colorectal HCT-116 cancer cell line proved that compound 6e induced cell cycle arrest in G0-G1 phase, while compound 10e arrested the cell cycle at both G0-G1 and S phases. Additionally, both compounds (6e and 10e) were potently able to induce apoptosis in HCT-116 cell line. Docking results of compounds 6e and 10e into the pocket of EGFR active site showed their similar main binding features with Olmutinib, while compounds 7d and 10e showed only moderate fitting into DHFR compared to methotrexate. In silico studies revealed that most of the tested compounds obeyed Lipinski's RO5 and showed positive drug likeness scores.

Direct radiolabeling of methotrexate and methotrexate micelles with Tc-99m using QbD approach.

The aim of the work presented in this manuscript was to radiolabel methotrexate and prepare radiolabeled methotrexate micelles, an antifolate drug with Tc-99m using QbD approach. The radiolabeling was executed using the experimental design and the radiolabeled drug was further encapsulated in micelles. The authors are of the view that the radiolabeled MTX could be used to target the folate receptor overexpressing cancers such as the kidney, colorectal, breast, brain etc thereby opening newer possibilities to the theranostic applications of the formed conjugate.

Quassinoids from Eurycoma longifolia as Potential Dihydrofolate Reductase Inhibitors: A Computational Study.

BACKGROUND: Quassinoids are degraded triterpene compounds that can be obtained from various species of the Simaroubaceae plant family, including Eurycoma longifolia. Quassinoids are the major compounds in E. longifolia, and they are known to have various medicinal potentials, such as anticancer and antimalarial properties. Dihydrofolate reductase (DHFR) was reported to be one of the important targets for certain anticancer and antimalarial drugs. Twelve quassinoids from E. longifolia were identified to have anticancer effects based on their IC50 values. This study aimed to evaluate the interactions of these twelve quassinoids with DHFR via Autodock 4.2 software and Biovia Discovery Studio Visualiser. METHODS: Twelve quassinoids from E. longifolia and their interactions with DHFR were evaluated via Autodock 4.2 software and Biovia Discovery Studio Visualiser. Their drug-likeness and pharmacokinetic properties were also assessed using the ADMETlab 2.0 program. RESULTS: The molecular docking results showed that eleven quassinoids showed better docking scores than methotrexate, in which the binding energy (BE) of these quassinoids ranged from - 7.87 to -9.58 kcal/mol. Their inhibition constant (Ki) ranged from 0.095 to 1.71 µM. At the same time, the BE and Ki values for methotrexate were -7.80 kcal/mol and 1.64 µM, respectively. CONCLUSION: From the analysis, 6-dehydrolongilactone and eurycomalide B are among the twelve compounds that showed great potential as hit-to-lead compounds based on the docking score on DHFR, drug-likeness, and ADMET properties. These results suggest a great potential to pursue validation studies via in vitro and in vivo models.

Repurposing FDA-approved drugs to target malaria through inhibition of dihydrofolate reductase in the folate biosynthesis pathway: A prospective approach.

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme that regulates the biosynthesis of tetrahydrofolate among various species of Plasmodium parasite. It is a validated target of the antifolate drug pyrimethamine (Pyr) in Plasmodium falciparum (Pf), but its clinical efficacy has been hampered due to the emergence of drug resistance. This has made the attempt to screen Food & Drug Administration-approved drugs against wild- and mutant PfDHFR by employing an in-silico pipeline to identify potent candidates. The current study has followed a virtual screening approach for identifying potential DHFR inhibitors from DrugBank database, based on a structure similarity search of candidates, followed by absorption, distribution, metabolism, and excretion estimation. The screened drugs were subjected to various parameters like docking, molecular mechanics with generalized born and surface area solvation calculations, and molecular simulations. We have thus identified two potential drug candidates, duloxetine and guanethidine, which can be repurposed to be tested for their efficacy against wild type and drug resistant falciparum malaria.

Methotrexate-based PROTACs as DHFR-specific chemical probes.

Methotrexate (MTX) is a tight-binding dihydrofolate reductase (DHFR) inhibitor, used as both an antineoplastic and immunosuppressant therapeutic. MTX, like folate undergoes folylpolyglutamate synthetase-mediated γ-glutamylation, which affects cellular retention and target specificity. Mechanisms of MTX resistance in cancers include a decrease in MTX poly-γ-glutamylation and an upregulation of DHFR. Here, we report a series of potent MTX-based proteolysis targeting chimeras (PROTACs) to investigate DHFR degradation pharmacology and one-carbon biochemistry. These on-target, cell-active PROTACs show proteasome- and E3 ligase-dependent activity, and selective degradation of DHFR in multiple cancer cell lines. By comparison, treatment with MTX increases cellular DHFR protein expression. Importantly, these PROTACs produced distinct, less-lethal phenotypes compared to MTX. The chemical probe set described here should complement conventional DHFR inhibitors and serve as useful tools for studying one-carbon biochemistry and dissecting complex polypharmacology of MTX and related drugs. Such compounds may also serve as leads for potential autoimmune and antineoplastic therapeutics.

Investigating the anti-cancer potential of pyrimethamine analogues through a modern chemical biology lens.

Pharmacological inhibition of dihydrofolate reductase (DHFR) is an established approach for treating a variety of human diseases, including foreign infections and cancer. However, treatment with classic DHFR inhibitors, such as methotrexate (MTX), are associated with negative side-effects and resistance mechanisms that have prompted the search for alternatives. The DHFR inhibitor pyrimethamine (Pyr) has compelling anti-cancer activity in in vivo models, but lacks potency compared to MTX, thereby requiring higher concentrations to induce therapeutic responses. The purpose of this work was to investigate structural analogues of Pyr to improve its in vitro and cellular activity. A series of 36 Pyr analogues were synthesized and tested in a sequence of in vitro and cell-based assays to monitor their DHFR inhibitory activity, cellular target engagement, and impact on breast cancer cell viability. Ten top compounds were identified, two of which stood out as potential lead candidates, 32 and 34. These functionalized Pyr analogues potently engaged DHFR in cells, at concentrations as low as 1 nM and represent promising DHFR inhibitors that could be further explored as potential anti-cancer agents.

Discovery of novel tumor-targeted near-infrared probes with 6-substituted pyrrolo[2,3-d]pyrimidines as targeting ligands.

Since the overexpression of folate receptors (FRs) in certain types of cancers, a variety of FR-targeted fluorescent probes for tumor detection have been developed. However, the reported probes almost all have the same targeting ligand of folic acid with various fluorophores and/or linkers. In the present study, a series of novel tumor-targeted near-infrared (NIR) molecular fluorescent probes were designed and synthesized based on previously reported 6-substituted pyrrolo[2,3-d]pyrimidine antifolates. All newly synthesized probes showed specific FR binding in vitro, whereas GT-NIR-4 and GT-NIR-5 with a benzene and a thiophene ring, respectively, on the side chain of pyrrolo[2,3-d]pyrimidine exhibited better FR binding affinity than that of GT-NIR-6 with folic acid as targeting ligand. GT-NIR-4 also showed high tumor uptake in KB tumor-bearing mice with good pharmacokinetic properties and biological safety. This work demonstrates the first attempt to replace folic acid with antifolates as targeting ligands for tumor-targeted NIR probes.

Structure-Based Design of Transport-Specific Multitargeted One-Carbon Metabolism Inhibitors in Cytosol and Mitochondria.

Multitargeted agents provide tumor selectivity with reduced drug resistance and dose-limiting toxicities. We previously described the multitargeted 6-substituted pyrrolo[3,2-d]pyrimidine antifolate 1 with activity against early- and late-stage pancreatic tumors with limited tumor selectivity. Structure-based design with our human serine hydroxymethyl transferase (SHMT) 2 and glycinamide ribonucleotide formyltransferase (GARFTase) structures, and published X-ray crystal structures of 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/inosine monophosphate cyclohydrolase (ATIC), SHMT1, and folate receptor (FR) α and ß afforded 11 analogues. Multitargeted inhibition and selective tumor transport were designed by providing promiscuous conformational flexibility in the molecules. Metabolite rescue identified mitochondrial C1 metabolism along with de novo purine biosynthesis as the targeted pathways. We identified analogues with tumor-selective transport via FRs and increased SHMT2, SHMT1, and GARFTase inhibition (28-, 21-, and 11-fold, respectively) compared to 1. These multitargeted agents represent an exciting new structural motif for targeted cancer therapy with substantial advantages of selectivity and potency over clinically used antifolates.

Kinetic Barrier to Enzyme Inhibition Is Manipulated by Dynamical Local Interactions in E. coli DHFR.

Dihydrofolate reductase (DHFR) is an important drug target and a highly studied model protein for understanding enzyme dynamics. DHFR's crucial role in folate synthesis renders it an ideal candidate to understand protein function and protein evolution mechanisms. In this study, to understand how a newly proposed DHFR inhibitor, 4'-deoxy methyl trimethoprim (4'-DTMP), alters evolutionary trajectories, we studied interactions that lead to its superior performance over that of trimethoprim (TMP). To elucidate the inhibition mechanism of 4'-DTMP, we first confirmed, both computationally and experimentally, that the relative binding free energy cost for the mutation of TMP and 4'-DTMP is the same, pointing the origin of the characteristic differences to be kinetic rather than thermodynamic. We then employed an interaction-based analysis by focusing first on the active site and then on the whole enzyme. We confirmed that the polar modification in 4'-DTMP induces additional local interactions with the enzyme, particularly, the M20 loop. These changes are propagated to the whole enzyme as shifts in the hydrogen bond networks. To shed light on the allosteric interactions, we support our analysis with network-based community analysis and show that segmentation of the loop domain of inhibitor-bound DHFR must be avoided by a successful inhibitor.

First-in-class multifunctional TYMS nonclassical antifolate inhibitor with potent in vivo activity that prolongs survival.

Although thymidylate synthase (TYMS) inhibitors have served as components of chemotherapy regimens, the currently available inhibitors induce TYMS overexpression or alter folate transport/metabolism feedback pathways that tumor cells exploit for drug resistance, limiting overall benefit. Here we report a small molecule TYMS inhibitor that i) exhibited enhanced antitumor activity as compared with current fluoropyrimidines and antifolates without inducing TYMS overexpression, ii) is structurally distinct from classical antifolates, iii) extended survival in both pancreatic xenograft tumor models and an hTS/Ink4a/Arf null genetically engineered mouse tumor model, and iv) is well tolerated with equal efficacy using either intraperitoneal or oral administration. Mechanistically, we verify the compound is a multifunctional nonclassical antifolate, and using a series of analogs, we identify structural features allowing direct TYMS inhibition while maintaining the ability to inhibit dihydrofolate reductase. Collectively, this work identifies nonclassical antifolate inhibitors that optimize inhibition of thymidylate biosynthesis with a favorable safety profile, highlighting the potential for enhanced cancer therapy.

Design, synthesis, anticancer, and antibacterial evaluation of some quinazolinone-based derivatives as DHFR inhibitors.

Two series of quinazolinone derivatives were designed and synthesized as dihydrofolate reductase (DHFR) inhibitors. All compounds were evaluated for their antibacterial and antitumor activities. Antibacterial activity was evaluated against three strains of Gram-positive and Gram-negative bacteria. Compound 3d exhibited the highest inhibitory activity against Staphylococcus aureus DHFR (SaDHFR) with IC50 of 0.769 ± 0.04 µM compared to 0.255 ± 0.014 µM for trimethoprim. Compound 3e was also more potent than trimethoprim against Escherichia coli DHFR (EcDHFR) with IC50 of 0.158 ± 0.01 µM and 0.226 ± 0.014 µM, respectively. Compound 3e exhibited a promising antiproliferative effect against most of the tested cancer cells. It also showed potent activity against leukemia (CCRF-CEM, and RPMI-8226); lung NCI-H522, and CNS U251 with GI% of 65.2, 63.22, 73.28, and 97.22, respectively. The cytotoxic activity of compound 3e was almost half the activity of doxorubicin against CCRF-CEM cell line with IC50 of 1.569 ± 0.06 µM and 0.822 ± 0.03 µM, respectively. In addition, compound 3e inhibited human DHFR with IC50 value of 0.527 ± 0.028 µM in comparison to methotrexate (IC50 = 0.118 ± 0.006 µM). Compound 3e caused an arrest of the cell cycle mainly at the S phase and caused a rise in the overall apoptotic percentage from 2.03% to 48.51%. (23.89-fold). Treatment of CCRF-CEM cells with compound 3e produced a significant increase in the active caspase-3 level by 6.25-fold compared to untreated cells. Molecular modeling studies were performed to evaluate the binding pattern of the most active compounds in the bacterial and human DHFR.

In silico screening of inhibitors against human dihydrofolate reductase to identify potential anticancer compounds.

In all species, dihydrofolate reductase (DHFR) is an essential enzyme that regulates the cellular amount of tetrahydrofolate. Human DHFR (hDHFR) activity inhibition results in tetrahydrofolate depletion and cell death. This property has made hDHFR a therapeutic target for cancer. Methotrexate is a well-known hDHFR inhibitor, but its administration has shown some light to severe adverse effects. Therefore, we aimed to find new potential hDHFR inhibitors using structure-based virtual screening, ADMET prediction, molecular docking, and molecular dynamics simulations. Here, we used the PubChem database to find all compounds with at least 90% structural similarity to known natural DHFR inhibitors. To explore their interaction pattern and estimate their binding affinities, the screened compounds (2023) were subjected to structure-based molecular docking against hDHFR. The fifteen compounds that showed higher binding affinity to the hDHFR than the reference compound (methotrexate) displayed important molecular orientation and interactions with key residues in the enzyme's active site. These compounds were subjected to Lipinski and ADMET prediction. PubChem CIDs: 46886812 and 638190 were identified as putative inhibitors. In addition, molecular dynamics simulations revealed that the binding of compounds (CIDs: 46886812 and 63819) stabilized the hDHFR structure and caused minor conformational changes. Our findings suggest that two compounds (CIDs: 46886812 and 63819) could be promising potential inhibitors of hDHFR in cancer therapy.Communicated by Ramaswamy H. Sarma.

Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies.

Dihydrofolate reductase (DHFR) is a crucial enzyme that maintains the levels of 5,6,7,8-tetrahydrofolate (THF) required for the biological synthesis of the building blocks of DNA, RNA, and proteins. Over-activation of DHFR results in the progression of multiple pathological conditions such as cancer, bacterial infection, and inflammation. Therefore, DHFR inhibition plays a major role in treating these illnesses. Sesquiterpenes of various types are prime metabolites derived from the marine sponge Dactylospongia elegans and have demonstrated antitumor, anti-inflammation, and antibacterial capacities. Here, we investigated the in silico potential inhibitory effects of 87 D. elegans metabolites on DHFR and predicted their ADMET properties. Compounds were prepared computationally for molecular docking into the selected crystal structure of DHFR (PDB: 1KMV). The docking scores of metabolites 34, 28, and 44 were the highest among this series (gscore values of -12.431, -11.502, and -10.62 kcal/mol, respectively), even above the co-crystallized inhibitor SRI-9662 score (-10.432 kcal/mol). The binding affinity and protein stability of these top three scored compounds were further estimated using molecular dynamic simulation. Compounds 34, 28, and 44 revealed high binding affinity to the enzyme and could be possible leads for DHFR inhibitors; however, further in vitro and in vivo investigations are required to validate their potential.

Design and development of novel N-(4-aminobenzoyl)- l-glutamic acid conjugated 1,3,5-triazine derivatives as Pf-DHFR inhibitor: An in-silico and in-vitro study.

In the present work, a library of 120 compounds was prepared using various aliphatic and aromatic amines. Finally, 10 compounds were selected through in silico screening carrying 4-aminobenzoyl-l-glutamic acid and 1,3,5-triazine moiety. The docking results of compounds 4d16 and 4d38 revealed higher binding interaction with amino acids Asp54 (-537.96 kcal/mol) and Asp54, Phe116 (-618.22 kcal/mol) against wild (1J3I) and quadruple mutant (1J3K) type of Pf-DHFR inhibitors and were comparable to standard WR99210. These compounds were developed by facile and microwave-assisted synthesis via nucleophilic substitution reaction and characterized by different spectroscopic methods. In vitro antimalarial assay results also suggested that these two compounds were having higher antimalarial activity against chloroquine-sensitive (3D7) and chloroquine-resistant (Dd2) strain out of the ten synthesized compounds with IC50 13.25 µM and 14.72 µM, respectively. These hybrid scaffolds might be useful in the lead discovery of a new class of Pf-DHFR inhibitors.

Key interactions of pyrimethamine derivatives specific to wild-type and mutant P. falciparum dihydrofolate reductase based on 3D-QSAR, MD simulations and quantum chemical calculations.

Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (PfDHFR-TS) is an important target enzyme in malarial chemotherapy. An understanding of how novel inhibitors interact with wild-type (wtPfDHFR), quadruple-mutant (qmPfDHFR), and human (hDHFR) enzymes is required for the development of these compounds as antimalarials. This study is focused on a series of des-Cl and m-Cl phenyl analogs of pyrimethamine with various flexible 6-substituents. The interactions of these compounds with DHFR enzymes were investigated by 3 D-QSAR, MD simulations, MM-PBSA, and DFT calculations. CoMFA and CoMSIA models were developed with good predictive abilities for wtPfDHFR and qmPfDHFR. For hDHFR, CoMSIA models combined with clogP descriptor were successfully derived. Binding free energy using MM-PBSA and comparison of per residue decomposition energy analyses with the DFT method at M06-2X/6-31G ++(d,p) level of theory indicated that Asp54 and Phe58 play important roles in the binding of the most potent compound in the series (compound 27) with both wtPfDHFR and qmPfDHFR, whereas Arg59 and Arg122 were additionally found to interact with this inhibitor in qmPfDHFR. For hDHFR, the residues Glu30 and Phe34 but not Arg70, equivalent to Asp54, Phe58, and Arg122 in PfDHFR, also play role in compound 27 binding through strong hydrophobic interactions (Phe34) and hydrogen bond network with Glu30, Ile7, and Val115. From the key interactions identified in the DHFR-inhibitor complexes, a general scheme is proposed for designing new inhibitors selective for PfDHFR that is important for the development of novel antifolate antimalarials.Communicated by Ramaswamy H. Sarma.