Interaction of pyrimethamine, cycloguanil, WR99210 and their analogues with Plasmodium falciparum dihydrofolate reductase: structural basis of antifolate resistance

An In Silico Study of the Interactions of Alkaloids from Cryptolepis sanguinolenta with Plasmodium falciparum Dihydrofolate Reductase and Dihydroorotate Dehydrogenase

The Plasmodium falciparum dihydrofolate reductase (PfDHFR) and dihydroorotate dehydrogenase (PfDHODH) are essential for Plasmodium falciparum growth and development, and have been validated as targets for the development of new antimalarial agents. Several alkaloids isolated from Cryptolepis sanguinolenta have been reported to have antiplasmodial activity, but their protein targets are unknown. Therefore, molecular docking and molecular dynamics simulations were used to investigate the interactions and stability of the alkaloids with PfDHFR and PfDHODH. Based on physicochemical characteristics, alkaloids were grouped as sterically bulky (sb) or planar (pg). Docking results revealed strong binding affinities (−6.0 to −13.4 kcal/mol) of the alkaloids against PfDHODH and various strains of PfDHFR while interacting with key residues such as Asp54 and Phe58 in PfDHFR. The pg alkaloids had high binding affinity and preference for the inhibitor binding domain over the flavin mononucleotide (FMN) binding domain in PfDHODH due to size considerations. From the molecular dynamics trajectories, protein-alkaloid complexes were stable throughout the simulation, with supporting evidence from root mean square deviations, root mean square fluctuations, radius of gyration, free binding energies, and other parameters. We report herein that biscryptolepine and cryptomisrine (sb class), as well as cryptolepinone, cryptoheptine, cryptolepine, and neocryptolepine (pg class), are capable of inhibiting PfDHFR effectively in pyrimethamine sensitive and resistant cells. Also, our results show that alkaloids of the pg class can inhibit PfDHODH as FMN decoys, as well as direct enzyme inhibitors, thereby halting crucial protein function.

Characterization of Human-malarial Parasite Species based on DHFR and GST Targets Resulting in Changes in Anti-malarial Drug Binding Conformations

BACKGROUND

In this study, we focused primarily on three anti-malarial drugs, namely chloroquine, mefloquine, and proguanil, and these were tested against two malarial targets DHFR and GST. The species Plasmodium falciparum, Plasmodium malariae, Plasmodium ovale, Plasmodium vivax were used for the study.

OBJECTIVE

The purpose of this study was to determine the sequence and structural similarity of the proteins DHFR and GST among four Plasmodium species as well as to discover the in silico interactions with the aforementioned drug candidates.

METHODS

Bioinformatics databases, such as PDB, UniProt, DrugBank, PubChem, and tools, and software like Phyre 2.0, Clustal O (1.2.4), AutoDock 4, AutoDock Vina, and Discovery Studio Visualizer were used to determine the evolutionary significance of the Plasmodium species.

RESULT

The variations showed a difference in the binding patterns of drugs with our target proteins. Our finding reveals the Plasmodium spp divergence or convergence as well as the structural and sequential similarity or dissimilarity features.

CONCLUSION

Our result suggests that due to the deviation in the sequences and structures, variations in protein-drug binding patterns have emerged.

Discovery of Traditional Chinese Medicine Derived Compounds as Wild Type and Mutant Plasmodium falciparum Dihydrofolate Reductase Inhibitors: Induced Fit Docking and ADME Studies

BACKGROUND

In a bid to come up with effective compounds as inhibitors for antimalarial treatment, we built a library of 2,000 traditional Chinese medicine(TCM)-derived compounds retrieved from TCM Database@Taiwan.

METHODS

The active sites of both the wild type and mutant Plasmodium falciparum dihydrofolatereductase (pfDHFR) were explored using computational tools. pfDHFR, one of the prime drug targets in the prevention of malaria infection induced by the female anopheles mosquito has continued to offer resistance to drugs (antifolates) due to mutation in some of the key amino acid residues crucial for its inhibition.

RESULTS

We utilized virtual throughput screening and glide XP docking to screen the compounds, and 8 compounds were found to have promising docking scores with both the wild type and mutant pfDHFR. They were further subjected to Induce Fit Docking (IFD) to affirm their inhibitory potency. The ADME properties and biological activity spectrum of the compounds were also considered. The inhibition profile of the compounds revealed that a number of compounds formed intermolecular interactions with ASP54, ILE14, LEU164, SER108/ASN108, ARG122 and ASP58. Most of the compounds can be considered as drug candidates due to their antiprotozoal activities and accordance with the Lipinski's Rule of Five (ROF).

CONCLUSION

The outcome of the present study should further be investigated to attest the efficacy of these compounds as better drug candidates than the antifolates.

Targeting Kinetoplastid and Apicomplexan Thymidylate Biosynthesis as an Antiprotozoal Strategy

Kinetoplastid and apicomplexan parasites comprise a group of protozoans responsible for human diseases, with a serious impact on human health and the socioeconomic growth of developing countries. Chemotherapy is the main option to control these pathogenic organisms and nucleotide metabolism is considered a promising area for the provision of antimicrobial therapeutic targets. Impairment of thymidylate (dTMP) biosynthesis severely diminishes the viability of parasitic protozoa and the absence of enzymatic activities specifically involved in the formation of dTMP (e.g. dUTPase, thymidylate synthase, dihydrofolate reductase or thymidine kinase) results in decreased deoxythymidine triphosphate (dTTP) levels and the so-called thymineless death. In this process, the ratio of deoxyuridine triphosphate (dUTP) versus dTTP in the cellular nucleotide pool has a crucial role. A high dUTP/dTTP ratio leads to uracil misincorporation into DNA, the activation of DNA repair pathways, DNA fragmentation and eventually cell death. The essential character of dTMP synthesis has stimulated interest in the identification and development of drugs that specifically block the biochemical steps involved in thymine nucleotide formation. Here, we review the available literature in relation to drug discovery studies targeting thymidylate biosynthesis in kinetoplastid (genera Trypanosoma and Leishmania) and apicomplexan (Plasmodium spp and Toxoplasma gondii) protozoans. The most relevant findings concerning novel inhibitory molecules with antiparasitic activity against these human pathogens are presented herein.

Design, synthesis and biological evaluation of 4-aminoquinoline-guanylthiourea derivatives as antimalarial agents

Guanylthiourea (GTU) has been identified as an important antifolate antimalarial pharmacophore unit, whereas, 4-amino quinolones are already known for antimalarial activity. In the present work molecules carrying 4-aminoquinoline and GTU moiety have been designed using molecular docking analysis with PfDHFR enzyme and heme unit. The docking results indicated that the necessary interactions (Asp54 and Ile14) and docking score (-9.63 to -7.36 kcal/mmol) were comparable to WR99210 (-9.89 kcal/mol). From these results nine molecules were selected for synthesis. In vitro analysis of these synthesized compounds reveal that out of the nine molecules, eight show antimalarial activity in the range of 0.61-7.55 μM for PfD6 strain and 0.43-8.04 μM for PfW2 strain. Further, molecular dynamics simulations were performed on the most active molecule to establish comparative binding interactions of these compounds and reference ligand with Plasmodium falciparum dihydrofolate reductase (PfDHFR).

Targeting plant DIHYDROFOLATE REDUCTASE with antifolates and mechanisms for genetic resistance

The folate biosynthetic pathway and its key enzyme dihydrofolate reductase (DHFR) is a popular target for drug development due to its essential role in the synthesis of DNA precursors and some amino acids. Despite its importance, little is known about plant DHFRs, which, like the enzymes from the malarial parasite Plasmodium, are bifunctional, possessing DHFR and thymidylate synthase (TS) domains. Here using genetic knockout lines we confirmed that either DHFR-TS1 or DHFR-TS2 (but not DHFR-TS3) was essential for seed development. Screening mutated Arabidopsis thaliana seeds for resistance to antimalarial DHFR-inhibitor drugs pyrimethamine and cycloguanil identified causal lesions in DHFR-TS1 and DHFR-TS2, respectively, near the predicted substrate-binding site. The different drug resistance profiles for the plants, enabled by the G137D mutation in DHFR-TS1 and the A71V mutation in DHFR-TS2, were consistent with biochemical studies using recombinant proteins and could be explained by structural models. These findings provide a great improvement in our understanding of plant DHFR-TS and suggest how plant-specific inhibitors might be developed, as DHFR is not currently targeted by commercial herbicides.

Synthesis, SAR and biological studies of sugar amino acid-based almiramide analogues: N-methylation leads the way

Leishmaniasis, caused by the protozoan parasites of the genus Leishmania, is one of the most neglected diseases endemic in many continents posing enormous global health threats and therefore the discovery of new antileishmanial compounds is of utmost urgency. The antileishmanial activities of a library of sugar amino acid-based linear lipopeptide analogues were examined with the aim to identify potential drug candidates to treat visceral leishmaniasis. It was found that among the synthesized analogues, most of the permethylated compounds exhibited more activity in in vitro studies against intra-macrophagic amastigotes than the non-methylated analogues. SAR and NMR studies revealed that introduction of the N-methyl groups inhibited the formation of any turn structure in these molecules, which led to their improved activities.

Design, synthesis and biological evaluation of antimalarial activity of new derivatives of 2,4,6-s-triazine

Dihydrofolate reductase (DHFR) is an important enzyme for de novo synthesis of nucleotides in Plasmodium falciparum and it is essential for cell proliferation. DHFR is a well known antimalarial target for drugs like cycloguanil and pyrimethamine which target its inhibition for their pharmacological actions. However, the clinical efficacies of these antimalarial drugs have been compromising due to multiple mutations occurring in DHFR that lead to drug resistance. In this background, we have designed 22 s -triazine compounds using the best five parameters based 3D-QSAR model built by using genetic function approximation. In-silico designed compounds were further filtered to 6 compounds based upon their ADME properties, docking studies and predicted minimum inhibitory concentrations (MIC). Out of 6 compounds, 3 compounds were synthesized in good yield over 95% and characterized using IR, ¹HNMR, ¹³CNMR and mass spectroscopic techniques. Parasitemia inhibition assay was used to evaluate the antimalarial activity of s -triazine compounds against 3D7 strain of P. falciparum. All the three compounds (7, 13 and 18) showed 30 times higher potency than cycloguanil (standard drug). It was observed that compound 18 was the most active while the compound 13 was the least active. On the closer inspection of physicochemical properties and SAR, it was observed that the presence of electron donating groups, number of hydrogen bond formation, lipophilicity of ligands and coulson charge of nitrogen atom present in the triazine ring enhances the DHFR inhibition significantly. This study will contribute to further endeavours of more potent DHFR inhibitors.

Theoretical Investigation of the Enantioselective Complexations between pfDHFR and Cycloguanil Derivatives

Point mutations in Plasmodium falciparum dihydrofolate reductase (pfDHFR), especially the double mutant variant (A16V + S108T), led to ineffective inhibiting by cycloguanil (Cyc). Cycloguanil derivatives showed good inhibiting properties against wild-type and mutant pfDHFR with an inhibition constant as low as the nanomolar level. However, there have been no reports on the stereochemistry of the compounds, and this is important because the pure enantiomeric form of a chiral drug can exert desirable, as well as non-desirable responses on the body or both. In this work, three-dimensional structures of Cyc derivatives in R and S configuration were constructed and optimized using Hartree-Fock/6-31G (d,p). Their structures were docked into the binding pocket of wild-type and double mutant (A16V + S108T) pfDHFR, complexed with nicotinamide adenine dinucleotide phosphate (NADPH). Results indicate that both wild-type and mutant pfDHFR are enantioselective towards enantiomeric Cyc derivatives (R and S configuration).

Diarylpyrazole Ligated Dihydropyrimidine Hybrids as Potent Non-Classical Antifolates and Their Efficacy Against Plasmodium falciparum

A series of diarylpyrazole clubbed dihydropyrimidine derivatives (J1-J30) was synthesized under microwave-assisted heating conditions by employing Biginelli reaction methodology and utilizing triethylammonium acetate both as a catalyst and as reaction medium, leading towards a greener reaction pathway. The synthesized entities were screened for their antimalarial efficacy against a Plasmodium falciparum strain in vitro. The active entities (J9, J15, J21, J25, and J27) obtained out of the in vitro screening were further evaluated for their enzyme inhibitory potency against the Pf-DHFR enzyme in vitro as well as in silico using Glide. Furthermore, the active scaffolds were tested for their cytotoxicity against Vero cells, proving their nontoxic behavior and selectivity. The ADME parameters were also evaluated and predicted in silico, indicating good oral bioavailability of the compounds.

Synthesis and evaluation of naphthyl bearing 1,2,3-triazole analogs as antiplasmodial agents, cytotoxicity and docking studies

Novel series of naphthyl bearing 1,2,3-triazoles (4a-t) were synthesized and evaluated for their in vitro antiplasmodial activity against pyrimethamine (Pyr)-sensitive and resistant strains of Plasmodium falciparum. The synthesized compounds were assessed for their cytotoxicity employing human embryonic kidney cell line (HEK-293), and none of them was found to be toxic. Among them 4j, 4k, 4l, 4m, 4n, 4t exhibited significant antiplasmodial activity in both strains, of which compounds 4m, 4n and 4t (∼3.0-fold) displayed superior activity to Pyr against resistant strain. Pyr and selected compounds (4n, 4p and 4t) that repressed parasite development also inhibited PfDHFR activity of the soluble parasite extract, suggesting that anti-parasitic activity of these compounds is a result of inhibition of the parasite DHFR. In silico studies suggest that activity of these compounds might be enhanced due to π-π stacking.

Microwave Assisted Synthesis of Pyrimidines in Ionic Liquid and Their Potency as Non-Classical Malarial Antifolates

Synthesis of pyrazole-linked triazolo-pyrimidine hybrids was achieved by employing Biginelli-type reaction methodology in an ionic liquid (triethylammonium acetate) under microwave irradiation. This method proved to be highly efficient and the ionic liquid employed was found recyclable for up to five consecutive cycles. The synthesized molecules were further screened for their antimalarial efficacy screening out the active scaffolds J15, J18, J21, J24, J27, and J30. The active molecules were evaluated in an enzyme inhibition study against the active Plasmodium falciparum dihydrofolate reductase (Pf-DHFR), computationally as well as in vitro, demonstrating their potency as DHFR inhibitors. The active entities were also investigated for their oral bioavailability by predicting ADME properties in silico, indicating good bioavailability.

A multilayer screening approach toward the discovery of novel Pf-DHFR inhibitors

A small yet diverse xanthone library was build and computationally docked against wild type Pf-DHFR by Molegro Virtual Docker (MolDock). For analysis of results an integrated approach based on re-ranking, scaling (based on heavy atom counts), pose clustering and visual inspection was implemented. Standard methods such as self-docking (for docking), EF analysis, average rank determinations (for size normalization), and cluster quality indices (for pose clustering) were used for validation of results. Three compounds X5, X113A and X164B displayed contact footprints similar to the known inhibitors with good scores. Finally, 16 compounds were extracted from ZINC data base by similarity based screening, docking score and drug/lead likeness. Out of these 16 compounds, 11 displayed very close contact footprints to experimentally known inhibitors, indicating there potential utility in further drug discovery efforts.

Design, synthesis and biological evaluation of 6-aryl-1,6-dihydro-1,3,5-triazine-2,4-diamines as antiplasmodial antifolates

One-pot multistep synthesis of 1,6-dihydro-1,3,5-triazines displaying potent antiplasmodial activity in vitro.

Molecular dynamics of interactions between rigid and flexible antifolates and dihydrofolate reductase from pyrimethamine-sensitive and pyrimethamine-resistant Plasmodium falciparum

Currently, the usefulness of antimalarials such as pyrimethamine (PYR) is drastically reduced due to the emergence of resistant Plasmodium falciparum (Pf) caused by its dihydrofolate reductase (PfDHFR) mutations, especially the quadruple N51I/C59R/S108N/I164L mutations. The resistance was due to the steric conflict of PYR with S108N. WR99210 (WR), a dihydrotriazine antifolate with a flexible side chain that can avoid such conflict, can overcome this resistance through tight binding with the mutant. To understand factors contributing to different binding affinities of PYR/WR to the wild type (WT) and quadruple mutant (QM), we performed simulations on WR-WT, WR-QM, PYR-WT, and PYR-QM complexes and found that Ile14 and Asp54 were crucial for PYR/WR binding to PfDHFR due to strong hydrogen bonds. The quadruple mutations cause PYR to form, on average, fewer hydrogen bonds with Ile14 and Leu164, and to be displaced from its optimal orientation for Asp54 interaction. The predicted binding affinity ranking (WR-QM ≈ WR-WT ≈ PYR-WT >> PYR-QM) reasonably agrees with the inhibition constant (K(i)) ranking. Our results reveal important residues for tight binding of PYR/WR to WT/QM, which may be used to evaluate the inhibition effectiveness of antimalarials and to provide fundamental information for designing new drugs effective against drug-resistant P. falciparum.

Accessible mutational trajectories for the evolution of pyrimethamine resistance in the malaria parasite Plasmodium vivax

Antifolate antimalarials, such as pyrimethamine, have experienced a dramatic reduction in therapeutic efficacy as resistance has evolved in multiple malaria species. We present evidence from one such species, Plasmodium vivax, which has experienced sustained selection for pyrimethamine resistance at the dihydrofolate reductase (DHFR) locus since the 1970s. Using a transgenic Saccharomyces cerevisiae model expressing the P. vivax DHFR enzyme, we assayed growth rate and resistance of all 16 combinations of four DHFR amino acid substitutions. These substitutions were selected based on their known association with drug resistance, both in natural isolates and in laboratory settings, in the related malaria species P. falciparum. We observed a strong correlation between the resistance phenotypes for these 16 P. vivax alleles and previously observed resistance data for P. falciparum, which was surprising since nucleotide diversity levels and common polymorphic variants of DHFR differ between the two species. Similar results were observed when we expressed the P. vivax alleles in a transgenic bacterial system. This suggests common constraints on enzyme evolution in the orthologous DHFR proteins. The interplay of negative trade-offs between the evolution of novel resistance and compromised endogenous function varies at different drug dosages, and so too do the major trajectories for DHFR evolution. In simulations, it is only at very high drug dosages that the most resistant quadruple mutant DHFR allele is favored by selection. This is in agreement with common polymorphic DHFR data in P. vivax, from which this quadruple mutant is missing. We propose that clinical dosages of pyrimethamine may have historically been too low to select for the most resistant allele, or that the fitness cost of the most resistant allele was untenable without a compensatory mutation elsewhere in the genome.

Discovery of a new class of natural product-inspired quinazolinone hybrid as potent antileishmanial agents

The high potential of quinazolinone containing natural products and their derivatives in medicinal chemistry led us to discover four novel series of 53 compounds of quinazolinone based on the concept of molecular hybridization. Most of the synthesized analogues exhibited potent leishmanicidal activity against intracellular amastigotes (IC50 from 0.65 ± 0.2 to 7.76 ± 2.1 μM) as compared to miltefosine (IC50 = 8.4 ± 2.1 μM) and nontoxic toward the J-774A.1 cell line and Vero cells. Moreover, activation of Th1 type and suppression of Th2 type immune responses and induction in nitric oxide generation proved that 8a and 8g induce murine macrophages to prevent survival of parasites. Compounds 8a and 8g exhibited significant in vivo inhibition of parasite 73.15 ± 12.69% and 80.93 ± 10.50% against Leishmania donovani /hamster model. Our results indicate that compounds 8a, 8g, and 9f represent a new structural lead for this serious and neglected disease.

Aminoalkyl derivatives of guanidine diaromatic minor groove binders with antiprotozoal activity

Considering the strong DNA minor groove binding observed for our previous series of diaromatic symmetric and asymmetric guanidinium and 2-aminoimidazolinium derivatives, we report now the synthesis of new aminoalkyl derivatives of diaromatic guanidines with potential as DNA minor groove binders and antiprotozoal activity. The preparation of these aminoalkyl derivatives (12a-e, 13a-e, 14a-c,e, 15a-e, 16a-e) is presented as well as their affinity for DNA which was evaluated by means of DNA thermal denaturation experiments. Finally, the antiprotozoal activity of most of these aminoalkyl minor groove binders was evaluated in vitro against Trypanosoma brucei rhodesiense (8 compounds) and Plasmodium falciparum (18 compounds). The O-linked derivatives 13c and 14c showed 100 nM activities against P. falciparum, whereas for T. b. rhodesiense all compounds tested showed micromolar activity. Some of the derivatives prepared seem to exert the antimalarial activity by binding to the DNA minor groove whereas other sets of compounds could exert this antimalarial activity by inhibiting the parasite dihydrofolate reductase, for example.

Dihydrofolate reductase as a therapeutic target for infectious diseases: opportunities and challenges

Infectious diseases caused by parasites continue to take a massive toll on human health in the poor regions of the world. Filling the anti-infective drug-discovery pipeline has never been as challenging as it is now. The organisms responsible for these diseases have interesting biology with many potential biochemical targets. Inhibition of metabolic enzymes has been established as an attractive strategy for anti-infectious drug development. In this field, dihydrofolate reductase (DHFR) is an important enzyme in nucleic and amino acid synthesis and an extensively studied drug target over the past 50 years. The challenges for novel DHFR inhibition-based chemotherapeutics for the treatment of infectious diseases are now focused on overcoming the resistance problem as well as cost–effectiveness. Each year, the large number of literature citations attest the continued popularity of DHFR. It becomes truly the ‘enzyme of choice for all seasons and almost all reasons’. Herein, we summarize the opportunities and challenges in developing novel lead based on this target.

Combined spatial limitation around residues 16 and 108 of Plasmodium falciparum dihydrofolate reductase explains resistance to cycloguanil

Natural mutations of Plasmodium falciparum dihydrofolate reductase (PfDHFR) at A16V and S108T specifically confer resistance to cycloguanil (CYC) but not to pyrimethamine (PYR). In order to understand the nature of CYC resistance, the effects of various mutations at A16 on substrate and inhibitor binding were examined. Three series of mutations at A16 with or without the S108T/N mutation were generated. Only three mutants with small side chains at residue 16 (G, C, and S) were viable from bacterial complementation assay in the S108 series, whereas these three and an additional four mutants (T, V, M, and I) with slightly larger side chains were viable with simultaneous S108T mutation. Among these combinations, the A16V+S108T mutant was the most CYC resistant, and all of the S108T series ranged from being highly to moderately sensitive to PYR. In the S108N series, a strict requirement for alanine was observed at position 16. Crystal structure analyses reveal that in PfDHFR-TS variant T9/94 (A16V+S108T) complexed with CYC, the ligand has substantial steric conflicts with the side chains of both A16V and S108T, whereas in the complex with PYR, the ligand only showed mild conflict with S108T. CYC analogs designed to avoid such conflicts improved the binding affinity of the mutant enzymes. These results show that there is greater spatial limitation around the S108T/N residue when combined with the limitation imposed by A16V. The limitation of mutation of this series provides opportunities for drug design and development against antifolate-resistant malaria.

2-Aminopyrimidine based 4-aminoquinoline anti-plasmodial agents. Synthesis, biological activity, structure-activity relationship and mode of action studies

2-Aminopyrimidine based 4-aminoquinolines were synthesized using an efficacious protocol. Some of the compounds showed in vitro anti-plasmodial activity against drug-sensitive CQ(S) (3D7) and drug-resistant CQ(R) (K1) strains of Plasmodium falciparum in the nM range. In particular, 5-isopropyloxycarbonyl-6-methyl-4-(2-nitrophenyl)-2-[(7-chloroquinolin-4-ylamino)butylamino] pyrimidine depicted the lowest IC(50) (3.6 nM) value (56-fold less than CQ) against CQ(R) strain. Structure-activity profile and binding with heme, μ-oxo-heme have been studied. Binding assays with DNA revealed better binding with target parasite type AT rich pUC18 DNA. Most compounds were somewhat cytotoxic, but especially cytostatic. Molecular docking analysis with Pf DHFR allowed identification of stabilizing interactions.

The evolutionary landscape of antifolate resistance in Plasmodium falciparum

Resistance to antifolates in Plasmodium falciparum is well described and has been observed in clinical settings for decades. At the molecular level, point mutations in the dhfr gene that lead to resistance have been identified, and the crystal structure of the wildtype and mutant dihydrofolate reductase enzymes have been solved in complex with native substrate and drugs. However, we are only beginning to understand the complexities of the evolutionary pressures that lead to the evolution of drug resistance in this system. Microbial systems that allow heterologous expression of malarial proteins provide a tractable way to investigate patterns of evolution that can inform our eventual understanding of the more complex factors that influence the evolution of drug resistance in clinical settings. In this paper we will review work in Escherichia coli and Saccharomyces cerevisiae expression systems that explore the fitness landscape of mutations implicated in drug resistance and show that (i) a limited number of evolutionary pathways to resistance are followed with high probability; (ii) fitness costs associated with the maintenance of high levels of resistance are modest; and (iii) different antifolates may exert opposing selective forces.

Binding modes of 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs to P. falciparum dihydrofolate reductase enzyme: Molecular docking studies

A molecular docking study was carried out on 28 compounds belonging to 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs using Glide, FlexX and GOLD programs and the X-ray crystallographic structures of the quadruple mutant (1J3K:pdb) and wild type (1J3I:pdb) Plasmodium falciparum dihydrofolate reductase enzyme. The experimental conformation the bound ligand WR99210 was precisely reproduced by the docking procedures as demonstrated by low (<2.00 Å) root-mean-square deviations. The results indicated that most of the compounds dock into the active sites of both the wild type and quadruple mutant P. falciparum dihydrofolate reductase enzymes. Visual inspection of the binding modes also demonstrated that most of the compounds could form H-bond interactions with the key amino acid residues (Asp54, Ile14 and Leu/Ile164) and with better docking scores than the bound compound (5). Their long side chains orient in the hydrophobic portion of the active site which is occupied by trichloro aryloxy side chain of WR99210 (5). Thus, avoid potential steric clashes with Asn108 (mutated from Ser108). Such a clash is known to be responsible for the resistance of the P. falciparum to pyrimethamine and cycloguanil.

Divalent N(I) compounds with two lone pairs on nitrogen

Carbon with the C(0) state has been reported recently, examples of which were known for the past decades. Silicon in the Si(0) state and phosphorus in the P(I) state are also known experimentally. This prompted us to search for divalent N(I) compounds, which resulted in the identification of ::N(←L)(2)(⊕) systems with bicoordinated nitrogen in the N(I) formal oxidation state. It was found that several biguanide derivatives (especially in their protonated state) belong to this class. Quantum chemical analysis provided the structural details, molecular orbitals, charge localization (vs delocalization) trends, etc. This class of compounds has been found to be characterized by two lone pairs on the central nitrogen, very similar to the central carbon in divalent C(0) compounds (::C(←L)(2)). The new bonding environment for nitrogen reported in this article, divalent nitrogen N(I), is clearly different from the nitrenium ions NR(2)(⊕). The electronic structure and reactivity of representative examples of this novel class of divalent nitrogen N(I) systems (::N(←L)(2)(⊕)) have been analyzed in detail, in terms of molecular orbitals, atomic charges, protonation energies, complexation energies with Lewis acids like BH(3), AlCl(3), and AuCl and compared with those of divalent C(0) systems.

Fitness trade-offs in the evolution of dihydrofolate reductase and drug resistance in Plasmodium falciparum

BACKGROUND

Patterns of emerging drug resistance reflect the underlying adaptive landscapes for specific drugs. In Plasmodium falciparum, the parasite that causes the most serious form of malaria, antifolate drugs inhibit the function of essential enzymes in the folate pathway. However, a handful of mutations in the gene coding for one such enzyme, dihydrofolate reductase, confer drug resistance. Understanding how evolution proceeds from drug susceptibility to drug resistance is critical if new antifolate treatments are to have sustained usefulness.

METHODOLOGY/PRINCIPAL FINDINGS

We use a transgenic yeast expression system to build on previous studies that described the adaptive landscape for the antifolate drug pyrimethamine, and we describe the most likely evolutionary trajectories for the evolution of drug resistance to the antifolate chlorcycloguanil. We find that the adaptive landscape for chlorcycloguanil is multi-peaked, not all highly resistant alleles are equally accessible by evolution, and there are both commonalities and differences in adaptive landscapes for chlorcycloguanil and pyrimethamine.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that cross-resistance between drugs targeting the same enzyme reflect the fitness landscapes associated with each particular drug and the position of the genotype on both landscapes. The possible public health implications of these findings are discussed.

A common feature-based 3D-pharmacophore model generation and virtual screening: identification of potential PfDHFR inhibitors

A four-feature 3D-pharmacophore model was built from a set of 24 compounds whose activities were reported against the V1/S strain of the Plasmodium falciparum dihydrofolate reductase (PfDHFR) enzyme. This is an enzyme harboring Asn51Ile + Cys59Arg + Ser108Asn + Ile164Leu mutations. The HipHop module of the Catalyst program was used to generate the model. Selection of the best model among the 10 hypotheses generated by HipHop was carried out based on rank and best-fit values or alignments of the training set compounds onto a particular hypothesis. The best model (hypo1) consisted of two H-bond donors, one hydrophobic aromatic, and one hydrophobic aliphatic features. Hypo1 was used as a query to virtually screen Maybridge2004 and NCI2000 databases. The hits obtained from the search were subsequently subjected to FlexX and Glide docking studies. Based on the binding scores and interactions in the active site of quadruple-mutant PfDHFR, a set of nine hits were identified as potential inhibitors.

Stepwise acquisition of pyrimethamine resistance in the malaria parasite

The spread of high-level pyrimethamine resistance in Africa threatens to curtail the therapeutic lifetime of antifolate antimalarials. We studied the possible evolutionary pathways in the evolution of pyrimethamine resistance using an approach in which all possible mutational intermediates were created by site-directed mutagenesis and assayed for their level of drug resistance. The coding sequence for dihydrofolate reductase (DHFR) from the malaria parasite Plasmodium falciparum was mutagenized, and tests were carried out in Escherichia coli under conditions in which the endogenous bacterial enzyme was selectively inhibited. We studied 4 key amino acid replacements implicated in pyrimethamine resistance: N51I, C59R, S108N, and I164L. Using empirical estimates of the mutational spectrum in P. falciparum and probabilities of fixation based on the relative levels of resistance, we found that the predicted favored pathways of drug resistance are consistent with those reported in previous kinetic studies, as well as DHFR polymorphisms observed in natural populations. We found that 3 pathways account for nearly 90% of the simulated realizations of the evolution of pyrimethamine resistance. The most frequent pathway (S108N and then C59R, N51I, and I164L) accounts for more than half of the simulated realizations. Our results also suggest an explanation for why I164L is detected in Southeast Asia and South America, but not at significant frequencies in Africa.

WISDOM-II: screening against multiple targets implicated in malaria using computational grid infrastructures

BACKGROUND

Despite continuous efforts of the international community to reduce the impact of malaria on developing countries, no significant progress has been made in the recent years and the discovery of new drugs is more than ever needed. Out of the many proteins involved in the metabolic activities of the Plasmodium parasite, some are promising targets to carry out rational drug discovery.

MOTIVATION

Recent years have witnessed the emergence of grids, which are highly distributed computing infrastructures particularly well fitted for embarrassingly parallel computations like docking. In 2005, a first attempt at using grids for large-scale virtual screening focused on plasmepsins and ended up in the identification of previously unknown scaffolds, which were confirmed in vitro to be active plasmepsin inhibitors. Following this success, a second deployment took place in the fall of 2006 focussing on one well known target, dihydrofolate reductase (DHFR), and on a new promising one, glutathione-S-transferase.

METHODS

In silico drug design, especially vHTS is a widely and well-accepted technology in lead identification and lead optimization. This approach, therefore builds, upon the progress made in computational chemistry to achieve more accurate in silico docking and in information technology to design and operate large scale grid infrastructures.

RESULTS

On the computational side, a sustained infrastructure has been developed: docking at large scale, using different strategies in result analysis, storing of the results on the fly into MySQL databases and application of molecular dynamics refinement are MM-PBSA and MM-GBSA rescoring. The modeling results obtained are very promising. Based on the modeling results, In vitro results are underway for all the targets against which screening is performed.

CONCLUSION

The current paper describes the rational drug discovery activity at large scale, especially molecular docking using FlexX software on computational grids in finding hits against three different targets (PfGST, PfDHFR, PvDHFR (wild type and mutant forms) implicated in malaria. Grid-enabled virtual screening approach is proposed to produce focus compound libraries for other biological targets relevant to fight the infectious diseases of the developing world.

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

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 (r(2) = 0.727 and r(2) = 0.985) and mutant type (r(2) = 0.786 and r(2) = 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.