Dihydrofolate reductase as a therapeutic target

The folate antagonists are an important class of therapeutic compounds, as evidenced by their use as antiinfective, antineoplastic, and antiinflammatory drugs. Thus far, all of the clinically useful drugs of this class have been inhibitors of dihydrofolate reductase (DHFR), a key enzyme in the synthesis of thymidylate, and therefore, of DNA. The basis of the antiinfective selectivity of these compounds is clear; the antifolates trimethoprim and pyrimethamine are potent inhibitors of bacterial and protozoal DHFRs, respectively, but are only weak inhibitors of mammalian DHFRs. These species-selective agents apparently exploit the differences in the active site regions of the parasite and host enzymes. Methotrexate is the DHFR inhibitor used most often in a clinical setting as an anticancer drug and as an antiinflammatory and immunosuppressive agent. Considerable progress has been made recently in understanding the biochemical basis for the selectivity of this drug and the biochemical mechanism (or mechanisms) responsible for the development of resistance to treatment with the drug. This understanding has led to a new generation of DHFR inhibitors that are now in clinical trials.     

Dihydrofolate reductase gene variations in susceptibility to disease and treatment outcomes

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate to tetrahydrofolate (THF). THF is needed for the action of folate-dependent enzymes and is thus essential for DNA synthesis and methylation. The importance of this reaction is demonstrated by the effectiveness of antifolate medications used to treat cancer by inhibiting DHFR, thereby depleting THF and slowing DNA synthesis and cell proliferation. Due to the pivotal role that DHFR plays in folate metabolism and cancer treatment, changes in the level of DHFR expression can affect susceptibility to a variety of diseases dependent on folate status such as spina bifida and cancer. Likewise, variability in DHFR expression can affect sensitivity to anti-cancer drugs such as the folate antagonist methotrexate. Alterations in DHFR expression can be due to polymorphisms in the DHFR gene. Several variations have recently been described in DHFR, including promoter polymorphisms, the 19-bp deletion allele and variations in 3'UTR. These polymorphisms seem to be functional, affecting mRNA levels through various interesting mechanisms, including regulation through RNA interference. Several groups have assessed the association of these polymorphisms with folate levels, risk of cancer and spina bifida as well as the outcome of diseases treated with MTX. The latter may lead to different treatment schedules, improving treatment efficacy and/or allowing for a reduction in drug side effects. This review will summarize present knowledge regarding the predictive potential of DHFR polymorphisms in disease and treatment.     

Molecular modeling of wild-type and antifolate resistant mutant Plasmodium falciparum DHFR

The development of drug resistance is reducing the efficiency of antifolates as antimalarials. This phenomenon has been linked to the occurrence of mutations in the parasite's dihydrofolate reductase (DHFR). In this way, the resistance to pyrimethamine and cycloguanil, two potent inhibitors of P. falciparum DHFR, is mainly related to mutations (single and crossed) at residues 16, 51, 59, 108 and 164 of the enzyme. In this work, we have refined a recently proposed homology-model of P. falciparum DHFR, and the resulting structure was used to obtain models for 14 mutant enzymes, employing molecular modeling. Ternary complexes of the mutant enzymes with these inhibitors have been superimposed to equivalent ternary complexes of the wild-type enzyme, allowing the proposition of hypotheses for the role of each mutation in drug resistance. Based on these results, possible reasons for antifolate resistance have been proposed.     

Three-dimensional structure of M. tuberculosis dihydrofolate reductase reveals opportunities for the design of novel tuberculosis drugs

Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate to tetrahydrofolate and is essential for the synthesis of thymidylate, purines and several amino acids. Inhibition of the enzyme's activity leads to arrest of DNA synthesis and cell death. The enzyme has been studied extensively as a drug target for bacterial, protozoal and fungal infections, and also for neoplastic and autoimmune diseases. Here, we report the crystal structure of dihydrofolate reductase from Mycobacterium tuberculosis, a human pathogen responsible for the death of millions of human beings per year. Three crystal structures of ternary complexes of M. tuberculosis DHFR with NADP and different inhibitors have been determined, as well as the binary complex with NADP, with resolutions ranging from 1.7 to 2.0 A. The three DHFR inhibitors are the anticancer drug methotrexate, the antimicrobial trimethoprim and Br-WR99210, an analogue of the antimalarial agent WR99210. Structural comparison of these complexes with human dihydrofolate reductase indicates that the overall protein folds are similar, despite only 26 % sequence identity, but that the environments of both NADP and of the inhibitors contain interesting differences between the enzymes from host and pathogen. Specifically, residues Ala101 and Leu102 near the N6 of NADP are distinctly more hydrophobic in the M. tuberculosis than in the human enzyme. Another striking difference occurs in a region near atoms N1 and N8 of methotrexate, which is also near atom N1 of trimethoprim, and near the N1 and two methyl groups of Br-WR99210. A glycerol molecule binds here in a pocket of the M. tuberculosis DHFR:MTX complex, while this pocket is essentially filled with hydrophobic side-chains in the human enzyme. These differences between the enzymes from pathogen and host provide opportunities for designing new selective inhibitors of M. tuberculosis DHFR.     

Thymidylate synthase-dihydrofolate reductase in protozoa

In protozoa, thymidylate synthase (TS) and dihydrofolate reductase (DHFR) exist on the same polypeptide. The DHFR domain is on the amino terminus, TS is on the carboxy terminus, and the domains are separated by a junction peptide of varying size depending on the source. The native protein is a dimer of two such subunits and is 110-140 kDa. Most studies of bifunctional TS-DHFR have been performed with the protein from anti-folate resistant strains of Leishmania major, which show amplification of the TS-DHFR gene and overproduction of the bifunctional protein. The Leishmania TS-DHFR has also been highly expressed in heterologous systems. There is extensive communication between domains, and channeling of the H2folate product of TS to DHFR. Anti-folates commonly used to treat microbial infections are poor inhibitors of L. major DHFR. However, selective inhibitors of L. major vs human DHFR have been found. The TS-DHFR from Plasmodium falciparum has also been cloned and sequenced. Interestingly, pyrimethamine-resistant strains of P. falciparum have a common point mutation in the DHFR coding sequence which causes decreased binding of the folate analog. A detailed knowledge of the structure and function of protozoan TS-DHFRs will soon be available.     

Analysis of Coxiela burnetti dihydrofolate reductase via in silico docking with inhibitors and molecular dynamics simulation

Coxiella burnetii is a gram-negative bacterium able to infect several eukaryotic cells, mainly monocytes and macrophages. It is found widely in nature with ticks, birds, and mammals as major hosts. C. burnetii is also the biological warfare agent that causes Q fever, a disease that has no vaccine or proven chemotherapy available. Considering the current geopolitical context, this fact reinforces the need for discovering new treatments and molecular targets for drug design against C. burnetii. Among the main molecular targets against bacterial diseases reported, the enzyme dihydrofolate reductase (DHFR) has been investigated for several infectious diseases. In the present work, we applied molecular modeling techniques to evaluate the interactions of known DHFR inhibitors in the active sites of human and C. burnetii DHFR (HssDHFR and CbDHFR) in order to investigate their potential as selective inhibitors of CbDHFR. Results showed that most of the ligands studied compete for the binding site of the substrate more effectively than the reference drug trimethoprim. Also the most promising compounds were proposed as leads for the drug design of potential CbDHFR inhibitors.     

Mobility of the spin-labeled side chains of some novel antifolate inhibitors in their complexes with dihydrofolate reductase

Four spin-labeled inhibitors of dihydrofolate reductase (DHFR) have been synthesized, each of which has the 2,2,6,6-tetramethyl-1-piperidinyloxy (TEMPO) reporting group at a different distance from the 2,4-diaminopyrimidine moiety by which the inhibitors are anchored and oriented in the active site. Inhibitors in which the TEMPO group is attached by a short side chain are weakly bound to DHFR from bacteria (Streptococcus faecium and Lactobacillus casei), to the bovine enzyme and to recombinant human DHFR. However, binding is sufficiently tight, especially in the ternary complexes with NADPH, for recording of the EPR spectra of the bound ligands. The spectra indicate that when these inhibitors are bound to the enzyme the TEMPO group is highly immobilized with correlation time, tau c, 4-20ns. Inhibitors that have the reporter group attached to the glutamate moiety of methotrexate bind to all four DHFRs more tightly than the inhibitors with shorter side chains by factors of up to 10(6). However, in most complexes formed by the inhibitors with longer side chains immobilization of the TEMPO group is slight (tau c 0.2-4 ns). These results are in general agreement with predictions from X-ray crystallographic results including thermal factors but there are some unanticipated differences between some results for bacterial and eukaryotic enzymes. Three of the splin-labeled inhibitors would provide good probes for distance measurements in and around the active site of mammalian DHFR.     

Characterization of dihydrofolate reductase of Pneumocystis carinii and Toxoplasma gondii

Pneumocystis carinii and Toxoplasma gondii are opportunistic pathogens of immunosuppressed patients that are susceptible to therapy with inhibitors of dihydrofolate reductase (DHFR). The DHFR of these two organisms was characterized to facilitate the identification of more selective inhibitors. Similar to all reported protozoa, T. gondii has a bifunctional enzyme, of 120,000 Da, that possesses both DHFR and thymidylate synthase (TS) activity. Unexpectedly, P. carinii DHFR activity was present on a small molecule, of 26,000 Da. T. gondii DHFR and TS activity coeluted during affinity chromatography using a methotrexate-Sepharose column, whereas P. carinii DHFR and TS activity could be separated by affinity chromatography using the same column. P. carinii DHFR could be easily distinguished from rat DHFR, which is similar in size, by the differences in Km for dihydrofolate (P. carinii, 17.6 +/- 3.9 microM; rat, 4.0 +/- 2.2 microM). Since all protozoa reported have a large molecular weight, bifunctional DHFR, these studies support the classification of P. carinii as a fungus. These studies also provide a basis for the development of more effective therapeutic agents for these pathogens.     

Crystal Structures of Wild-type and Mutant Methicillin-resistant Staphylococcus aureus Dihydrofolate Reductase Reveal an Alternate Conformation of NADPH That May Be Linked to Trimethoprim Resistance

Both hospital- and community-acquired Staphylococcus aureus infections have become major health concerns in terms of morbidity, suffering and cost. Trimethoprim-sulfamethoxazole (TMP-SMZ) is an alternative treatment for methicillin-resistant S. aureus (MRSA) infections. However, TMP-resistant strains have arisen with point mutations in dihydrofolate reductase (DHFR), the target for TMP. A single point mutation, F98Y, has been shown biochemically to confer the majority of this resistance to TMP. Using a structure-based approach, we have designed a series of novel propargyl-linked DHFR inhibitors that are active against several trimethoprim-resistant enzymes. We screened this series against wild-type and mutant (F98Y) S. aureus DHFR and found that several are active against both enzymes and specifically that the meta-biphenyl class of these inhibitors is the most potent. In order to understand the structural basis of this potency, we determined eight high-resolution crystal structures: four each of the wild-type and mutant DHFR enzymes bound to various propargyl-linked DHFR inhibitors. In addition to explaining the structure-activity relationships, several of the structures reveal a novel conformation for the cofactor, NADPH. In this new conformation that is predominantly associated with the mutant enzyme, the nicotinamide ring is displaced from its conserved location and three water molecules complete a network of hydrogen bonds between the nicotinamide ring and the protein. In this new position, NADPH has reduced interactions with the inhibitor. An equilibrium between the two conformations of NADPH, implied by their occupancies in the eight crystal structures, is influenced both by the ligand and the F98Y mutation. The mutation induced equilibrium between two NADPH-binding conformations may contribute to decrease TMP binding and thus may be responsible for TMP resistance.     

Virtual screening reveals allosteric inhibitors of the Toxoplasma gondii thymidylate synthase–dihydrofolate reductase

The parasite Toxoplasma gondii can lead to toxoplasmosis in those who are immunocompromised. To combat the infection, the enzyme responsible for nucleotide synthesis thymidylate synthase–dihydrofolate reductase (TS–DHFR) is a suitable drug target. We have used virtual screening to determine novel allosteric inhibitors at the interface between the two TS domains. Selected compounds from virtual screening inhibited TS activity. Thus, these results show that allosteric inhibition by small drug-like molecules can occur in T. gondii TS–DHFR and pave the way for new and potent species-specific inhibitors.     

Bifunctional thymidylate synthase-dihydrofolate reductase in protozoa

Protozoa contain thymidylate synthase (TS) and dihydrofolate reductase (DHFR) on the same polypeptide. In the bifunctional protein, the DHFR domain is on the amino terminus, TS is on the carboxyl terminus, and the two domains are separated by a junction peptide of varying size depending on the source. The native protein is composed of a dimer of two such subunits and is 110-140 kDa. Most studies of the bifunctional TS-DHFR have been performed with the protein from anti-folate resistant strains of Leishmania major, which show amplification of the TS-DHFR gene and overproduction of the bifunctional protein. The Leishmania TS-DHFR has also been highly expressed in heterologous systems. There appears to be extensive communication among domains and channeling of the H2folate product of TS to DHFR. Anti-folates commonly used to treat microbial infections are poor inhibitors of L. major DHFR. However, selective inhibition of L. major vs. human DHFR does not appear difficult to achieve, and selective inhibitors are known. The TS-DHFR from Plasmodium falciparum has also been cloned and has recently been expressed in Escherichia coli, albeit in small amounts. Interestingly, pyrimethamine-resistant strains of P. falciparum all have a common point mutation in the DHFR coding sequence (Thr/Ser 108 to Asn), which causes decreased binding of the folate analog. It is suggested that if an appropriate inhibitor of the pyrimethamine-resistant P. falciparum DHFRs can be found, it may serve in combination with pyrimethamine as an antimalarial regimen with low propensity for the development of resistance. In the future, we project that we will have a detailed knowledge of the structure and function of TS-DHFRs, and have the essential tools necessary for a molecular-based approach to drug design.     

Much more than you expected: The non-DHFR-mediated effects of methotrexate

BackgroundFor decades, methotrexate (MTX; amethopterin) has been known as an antifolate inhibitor of dihydrofolate reductase (DHFR), and it is widely used for the treatment of various malignancies and autoimmune diseases. Although the inclusion of MTX in various therapeutic regimens is based on its ability to inhibit DHFR and consequently to suppress the synthesis of pyrimidine and purine precursors, recent studies have shown that MTX is also able to target other intracellular pathways that are independent of folate metabolism.Scope of reviewThe main aim of this review is to summarize the most important, up-to-date findings of studies regarding the non-DHFR-mediated mechanisms of MTX action.Major conclusionsThe effectiveness of MTX is undoubtedly caused by its capability to affect various intracellular pathways at many levels. Although the most important therapeutic mechanism of MTX is strongly based on the inhibition of DHFR, many other effects of this compound have been described and new studies bring new insights into the pharmacology of MTX every year.General significanceIdentification of these new targets for MTX is especially important for a better understanding of MTX action in new protocols of combination therapy.     

Structure-based approach to pharmacophore identification, in silico screening, and three-dimensional quantitative structure-activity relationship studies for inhibitors of Trypanosoma cruzi dihydrofolate reductase function

We have employed a structure-based three-dimensional quantitative structure-activity relationship (3D-QSAR) approach to predict the biochemical activity for inhibitors of T. cruzi dihydrofolate reductase-thymidylate synthase (DHFR-TS). Crystal structures of complexes of the enzyme with eight different inhibitors of the DHFR activity together with the structure in the substrate-free state (DHFR domain) were used to validate and refine docking poses of ligands that constitute likely active conformations. Structural information from these complexes formed the basis for the structure-based alignment used as input for the QSAR study. Contrary to indirect ligand-based approaches the strategy described here employs a direct receptor-based approach. The goal is to generate a library of selective lead inhibitors for further development as antiparasitic agents. 3D-QSAR models were obtained for T. cruzi DHFR-TS (30 inhibitors in learning set) and human DHFR (36 inhibitors in learning set) that show a very good agreement between experimental and predicted enzyme inhibition data. For crossvalidation of the QSAR model(s), we have used the 10% leave-one-out method. The derived 3D-QSAR models were tested against a few selected compounds (a small test set of six inhibitors for each enzyme) with known activity, which were not part of the learning set, and the quality of prediction of the initial 3D-QSAR models demonstrated that such studies are feasible. Further refinement of the models through integration of additional activity data and optimization of reliable docking poses is expected to lead to an improved predictive ability.     

Kinetics of the inhibition of bovine liver dihydrofolate reductase by tea catechins: origin of slow-binding inhibition and pH studies

Dihydrofolate reductase (DHFR) is the subject of intensive investigation since it appears to be the primary target enzyme for "antifolate" drugs, such as methotrexate and trimethoprim. Fluorescence quenching and stopped-flow fluorimetry show that the ester bond-containing tea polyphenols (-)-epigallocatechin gallate (EGCG) and (-)-epicatechin gallate (ECG) are potent and specific inhibitors of DHFR with inhibition constants (K(I)) of 120 and 82 nM, respectively. Both tea compounds showed the characteristics of slow-binding inhibitors of bovine liver DHFR. In this work, we have determined a complete kinetic scheme to explain the slow-binding inhibition and the pH effects observed during the inhibition of bovine liver DHFR by these tea polyphenols. Experimental data, based on fluorimetric titrations, and transient phase and steady-state kinetic studies confirm that EGCG and ECG are competitive inhibitors with respect to 7,8-dihydrofolate, which bind preferentially to the free form of the enzyme. The origin of their slow-binding inhibition is proposed to be the formation of a slow dissociation ternary complex by the reaction of NADPH with the enzyme-inhibitor complex. The pH controls both the ionization of critical catalytic residues of the enzyme and the protonation state of the inhibitors. At acidic pH, EGCG and ECG are mainly present as protonated species, whereas near neutrality, they evolve toward deprotonated species due to ionization of the ester-bonded gallate moiety (pK = 7.8). Although DHFR exhibits different affinities for the protonated and deprotonated forms of EGCG and ECG, it appears that the ionization state of Glu-30 in DHFR is critical for its inhibition. The physiological implications of these pH dependencies are also discussed.     

Novel non-classical C9-methyl-5-substituted-2,4-diaminopyrrolo[2,3-d]pyrimidines as potential inhibitors of dihydrofolate reductase and as anti-opportunistic agents

Six novel C9-methyl-5-substituted-2,4-diaminopyrrolo[2,3-d]pyrimidines 18-23 were synthesized as potential inhibitors of dihydrofolate reductase (DHFR) and as anti-opportunistic agents. These compounds represent the only examples of 9-methyl substitution in the carbon-carbon bridge of 2,4-diaminopyrrolo[2,3-d]pyrimidines. The analogs 18-23 were synthesized in a concise eight-step procedure starting from the appropriate commercially available aromatic methyl ketones. The key step involved a Michael addition reaction of 2,4,6-triaminopyrimidine to the appropriate 1-nitroalkene, followed by ring closure of the nitro adducts via a Nef reaction. The compounds were evaluated as inhibitors of DHFR from Pneumocystis carinii (pc), Toxoplasma gondii (tg), Mycobacterium avium (ma) and rat liver (rl). The biological result indicated that some of these analogs are potent inhibitors of DHFR and some have selectivity for pathogen DHFR. Compound 23 was a two digit nanomolar inhibitor of tgDHFR with 9.6-fold selectivity for tgDHFR.     

N-terminal amino acid sequence of the chromosomal dihydrofolate reductase purified from trimethoprim-resistant Staphylococcus aureus

The existence of two distinct dihydrofolate reductases (DHFR) in highly trimethoprim-resistant clinical isolates has been unequivocally demonstrated. The enzymes have been characterized with regard to the affinity for substrates and sensitivity to inhibitors. The chromosomal, trimethoprim-sensitive DHFR was purified to homogeneity by a new simple two-step procedure. Its N-terminal amino acid sequence, determined up to the first 35 amino acids, showed 69% homology with the Escherichia coli DHFR.     

The influence of glycerol on the binding of methotrexate to Mycobacterium tuberculosis dihydrofolate reductase: a molecular modelling study

The enzyme, dihydrofolate reductase (DHFR), from Mycobacterium tuberculosis (mt-DHFR) is believed to be a potential drug target for the treatment of tuberculosis. The co-crystal structure of mt-DHFR bound with glycerol (GOL), NAPDH and methotrexate (MTX) reveals a GOL binding site on the enzyme. This GOL binding site could be very important for the design of novel, selective mt-DHFR inhibitors, because this binding site is absent on human DHFR (h-DHFR). We have performed molecular dynamic simulations and free energy calculations to evaluate the binding affinity of GOL and its free energy contribution to the binding of MTX to mt-DHFR. The results showed that GOL does not bind tightly to mt-DHFR. Although GOL itself contributed free energy on MTX binding to mt-DHFR, GOL also increased the flexibilities of MTX, so that MTX cannot maintain strong electronic interactions with ARG32 and ARG60, which caused the total binding free energy to decrease. These data suggest that GOL binding is weak and it could be expelled from the binding site, to allow inhibitors containing appropriate side chains to bind. This observation can be used to inform future drug design studies, especially those aimed at improving drug selectivity against h-DHFR.     

Thermodynamic and NMR analysis of inhibitor binding to dihydrofolate reductase

Isothermal titration calorimetry (ITC) was used to determine the thermodynamic driving force for inhibitor binding to the enzyme dihydrofolate reductase (DHFR) from Escherichia coli. 1,4-Bis-{[N-(1-imino-1-guanidino-methyl)]sulfanylmethyl}-3,6-dimethyl-benzene (1) binds DHFR:NADPH with a K(d) of 13±5 nM while the related inhibitor 1-{[N-(1-imino-guanidino-methyl)]sulfanylmethyl}-3-trifluoromethyl-benzene (2) binds DHFR:NADPH with a K(d) of 3.2±2.2 μM. The binding of these inhibitors has both a favorable entropy and enthalpy of binding. Additionally, we observe positive binding cooperativity between both 1 and 2 and the cofactor NADPH. Binding of compound 1 to DHFR is 285-fold tighter in the presence of the NADPH cofactor than in its absence. We did not detect binding of 2 to DHFR in the absence of NADPH. The backbone amide (1)H and (15)N NMR resonances of DHFR:NADPH and both DHFR:NADPH inhibitor complexes were assigned in order to better understand the binding of these inhibitors in solution. The chemical shift perturbations observed with the binding of 1 were greatest at residues closest to the binding site, but significant perturbations also occur away from the inhibitor location at amino acids in the vicinity of residue 58 and in the GH loop. The pattern of chemical shift changes observed with the binding of 2 is similar to that seen with 1. The main differences in chemical shift perturbation between the two inhibitors are in the Met20 loop and in residues at the interface between the inhibitor and NADPH.     

Structure-based design of new DHFR-based antibacterial agents: 7-aryl-2,4-diaminoquinazolines

Dihydrofolate reductase (DHFR) inhibitors such as trimethoprim (TMP) have long played a significant role in the treatment of bacterial infections. Not surprisingly, after decades of use there is now bacterial resistance to TMP and therefore a need to develop novel antibacterial agents with expanded spectrum including these resistant strains. In this study, we investigated the optimization of 2,4-diamnoquinazolines for antibacterial potency and selectivity. Using structure-based drug design, several 7-aryl-2,4-diaminoquinazolines were discovered that have excellent sub-100 picomolar potency against bacterial DHFR. These compounds have good antibacterial activity especially on gram-positive pathogens including TMP-resistant strains.     

Computer modeling studies of the structural role of NADPH binding to active site mutants of human dihydrofolate reductase in complex with piritrexim

Dihydrofolate reductase (DHFR, EC 1.5.1.3) is one of the enzymes active in the folate cycle which plays an important role in DNA synthesis. Inhibition of DHFR is a key element in the treatment of many diseases, including cancer and AIDS related infections. A search for new selective inhibitors is motivated by the resistance to common drugs observed in the course of treatment. In this paper, results of a detailed computer analysis of human DHFR interactions with the lipophilic inhibitor piritrexim (PTX) are presented. It was found that the NADPH cofactor contributes 30% of the total PTX-enzyme interaction energy. Substitution of the highly conserved Glu30 with alanine does not lead to the release of the inhibitor from the hDHFR pocket. The important L22F point mutation does affect PTX orientation but does not changethe binding energy. Simulations of the dynamics of binary hDHFR-PTX complexes were performed with the use of Extensible Systematic Force Field (ESFF) and the results indicate structural changes in the enzyme induced by NADPH binding.