Plasmodium falciparum dihydrofolate reductase (PfDHFR) is an important target for antimalarial chemotherapy. Unfortunately, the emergence of resistant parasites has significantly reduced the efficiency of classical antifolate drugs such as cycloguanil and pyrimethamine. In this study, an approach toward molecular docking of the structures contained in the Available Chemicals Directory (ACD) database to search for novel inhibitors of PfDHFR is described. Instead of docking the whole ACD database, specific 3D pharmacophores were used to reduce the number of molecules in the database by excluding a priori molecules lacking essential requisites for the interaction with the enzyme and potentially unable to bind to resistant mutant PfDHFRs. The molecules in the resulting "focused" database were then evaluated with regard to their fit into the PfDHFR active site. Twelve new compounds whose structures are completely unrelated to known antifolates were identified and found to inhibit, at the micromolar level, the wild-type and resistant mutant PfDHFRs harboring A16V, S108T, A16V + S108T, C59R + S108N + I164L, and N51I + C59R + S108N + I164L mutations. Depending on the functional groups interacting with key active site residues of the enzyme, these inhibitors were classified as N-hydroxyamidine, hydrazine, urea, and thiourea derivatives. The structures of the complexes of the most active inhibitors, as refined by molecular mechanics and molecular dynamics, provided insight into how these inhibitors bind to the enzyme and suggested prospects for these novel derivatives as potential leads for antimalarial development.

A 3D pharmacophore model able to quantitatively predict inhibition constants was derived for a series of inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR), a validated target for antimalarial therapy. The data set included 52 inhibitors, with 23 of these comprising the training set and 29 an external test set. The activity range, expressed as Ki, of the training set molecules was from 0.3 to 11 300 nM. The 3D pharmacophore, generated with the HypoGen module of Catalyst 4.7, consisted of two hydrogen bond donors, one positive ionizable feature, one hydrophobic aliphatic feature, and one hydrophobic aromatic feature and provided a 3D-QSAR model with a correlation coefficient of 0.954. Importantly, the type and spatial location of the chemical features encoded in the pharmacophore were in full agreement with the key binding interactions of PfDHFR inhibitors as previously established by molecular modeling and crystallography of enzyme-inhibitor complexes. The model was validated using several techniques, namely, Fisher's randomization test using CatScramble, leave-one-out test to ensure that the QSAR model is not strictly dependent on one particular compound of the training set, and activity prediction in an external test set of compounds. In addition, the pharmacophore was able to correctly classify as active and inactive the dihydrofolate reductase and aldose reductase inhibitors extracted from the MDDR database, respectively. This test was performed in order to challenge the predictive ability of the pharmacophore with two classes of inhibitors that target very different binding sites. Molecular diversity of the data sets was finally estimated by means of the Tanimoto approach. The results obtained provide confidence for the utility of the pharmacophore in the virtual screening of libraries and databases of compounds to discover novel PfDHFR inhibitors.

No abstract

No abstract