Interactions between cycloguanil derivatives and wild type and resistance-associated mutant Plasmodium falciparum dihydrofolate reductases

Abstract  Comparative molecular field analysis (CoMFA) and quantum chemical calculations were performed on cycloguanil (Cyc) derivatives of the wild type and the quadruple mutant (Asn51Ile, Cys59Arg, Ser108Asn, Ile164Leu) of Plasmodium falciparum dihydrofolate reductase (PfDHFR). The represented CoMFA models of wild type ( and r ² = 0.985) and mutant type ( and r ² = 0.979) can describe the differences of the Cyc structural requirements for the two types of PfDHFR enzymes and can be useful to guide the design of new inhibitors. Moreover, the obtained particular interaction energies between the Cyc and the surrounding residues in the binding pocket indicated that Asn108 of mutant enzyme was the cause of Cyc resistance by producing steric clash with p-Cl of Cyc. Consequently, comparing the energy contributions with the potent flexible WR99210 inhibitor, it was found that the key mutant residue, Asn108, demonstrates attractive interaction with this inhibitor and some residues, Leu46, Ile112, Pro113, Phe116, and Leu119, seem to perform as second binding site with WR99210. Therefore, quantum chemical calculations can be useful for investigating residue interactions to clarify the cause of drug resistance. Graphical Abstract  CoMFA steric contour maps of Cyc derivatives against the quadruple mutant PfDHFR. Electrostatic van der Waals surface interaction of Cyc and the key resistance residue Asn108.      

Selectivity analysis of 5-(arylthio)-2,4-diaminoquinazolines as inhibitors of Candida albicans dihydrofolate reductase by molecular dynamics simulations

A series of 5-(arylthio)-2,4-diaminoquinazolines are known as selective inhibitors of dihydrofolate reductase (DHFR) from Candida albicans. We have performed docking and molecular dynamics simulations of these inhibitors with C. albicans and human DHFR to understand the basis for selectivity of these agents. Study was performed on a selected set of 10 compounds with variation in structure and activity. Molecular dynamics simulations were performed at 300 K for 45 ps with equilibration for 10 ps. Trajectory data was analyzed on the basis of hydrogen bond interactions, energy of binding and conformational energy difference. The results indicate that hydrogen bonds formed between the compound and the active site residues are responsible for inhibition and higher potency. The selectivity index, i.e the ratio of I₅₀ against human DHFR to I₅₀ against fungal DHFR, is mainly determined by the conformation adapted by the compounds within the active site of two enzymes. Since the human DHFR active site is rigid, the compound is trapped in a higher energy conformation. This energy difference between the two conformations E mainly governs the selectivity against fungal DHFR. The information generated from this analysis of potency and selectivity should be useful for further work in the area of antifungal research.