The structure of Pneumocystis carinii dihydrofolate reductase to 1.9 A resolution

Deep mutational scanning of Pneumocystis jirovecii dihydrofolate reductase reveals allosteric mechanism of resistance to an antifolate

Pneumocystis jirovecii is a fungal pathogen that causes pneumocystis pneumonia, a disease that mainly affects immunocompromised individuals. This fungus has historically been hard to study because of our inability to grow it in vitro. One of the main drug targets in P. jirovecii is its dihydrofolate reductase (PjDHFR). Here, by using functional complementation of the baker's yeast ortholog, we show that PjDHFR can be inhibited by the antifolate methotrexate in a dose-dependent manner. Using deep mutational scanning of PjDHFR, we identify mutations conferring resistance to methotrexate. Thirty-one sites spanning the protein have at least one mutation that leads to resistance, for a total of 355 high-confidence resistance mutations. Most resistance-inducing mutations are found inside the active site, and many are structurally equivalent to mutations known to lead to resistance to different antifolates in other organisms. Some sites show specific resistance mutations, where only a single substitution confers resistance, whereas others are more permissive, as several substitutions at these sites confer resistance. Surprisingly, one of the permissive sites (F199) is without direct contact to either ligand or cofactor, suggesting that it acts through an allosteric mechanism. Modeling changes in binding energy between F199 mutants and drug shows that most mutations destabilize interactions between the protein and the drug. This evidence points towards a more important role of this position in resistance than previously estimated and highlights potential unknown allosteric mechanisms of resistance to antifolate in DHFRs. Our results offer unprecedented resources for the interpretation of mutation effects in the main drug target of an uncultivable fungal pathogen.

Deciphering Antifungal Drug Resistance in Pneumocystis jirovecii DHFR with Molecular Dynamics and Machine Learning

Drug resistance impacts the effectiveness of many new therapeutics. Mutations in the therapeutic target confer resistance; however, deciphering which mutations, often remote from the enzyme active site, drive resistance is challenging. In a series of Pneumocystis jirovecii dihydrofolate reductase variants, we elucidate which interactions are key bellwethers to confer resistance to trimethoprim using homology modeling, molecular dynamics, and machine learning. Six molecular features involving mainly residues that did not vary were the best indicators of resistance.

Pharmacoinformatic Study on the Selective Inhibition of the Protozoan Dihydrofolate Reductase Enzymes

Dihydrofolate reductase (DHFR) is an essential enzyme of the folate metabolic pathway in protozoa and it is a validated, potential drug target in many infectious diseases. Information about unique conserved residues of the DHFR enzyme is required to understand residual selectivity of the protozoan DHFR enzyme. The three dimensional crystal structures are not available for all the protozoan DHFR enzymes. Enzyme-substrate/inhibitor interaction information is required for the binding mode characterization in protozoan DHFR for selective inhibitor design. In this work, multiple sequence analysis was carried out in all the studied species. Homology models were built for protozoan DHFR enzymes, for which 3D structures are not available in PDB. The molecular docking and Prime-MMGBSA calculations of the natural substrate (dihydrofolate, DHF) and classical DHFR inhibitor (methotrexate, MTX) were performed in protozoan DHFR enzymes. Comparative sequence analysis showed that an overall sequence identity between the studied species ranging from 22.94 % (CfDHFR-BgDHFR) to 94.61 % (LdDHFR-LmDHFR). Interestingly, it was observed that most of the active site residues were conserved in all the cases and all the enzymes exhibit similar key binding interactions with DHF and MTX in molecular docking analysis, but there are a few key binding residues which differ in protozoan species that makes it suitable for target selectivity. This information can be used to design selective and potent protozoan DHFR enzyme inhibitors.

Trimethoprim resistance of dihydrofolate reductase variants from clinical isolates of Pneumocystis jirovecii

Pneumocystis jirovecii is an opportunistic pathogen that causes serious pneumonia in immunosuppressed patients. Standard therapy and prophylaxis include trimethoprim (TMP)-sulfamethoxazole; trimethoprim in this combination targets dihydrofolate reductase (DHFR). Fourteen clinically observed variants of P. jirovecii DHFR were produced recombinantly to allow exploration of the causes of clinically observed failure of therapy and prophylaxis that includes trimethoprim. Six DHFR variants (S31F, F36C, L65P, A67V, V79I, and I158V) showed resistance to inhibition by trimethoprim, with Ki values for trimethoprim 4-fold to 100-fold higher than those for the wild-type P. jirovecii DHFR. An experimental antifolate with more conformational flexibility than trimethoprim showed strong activity against one trimethoprim-resistant variant. The two variants that were most resistant to trimethoprim (F36C and L65P) also had increased Km values for dihydrofolic acid (DHFA). The catalytic rate constant (kcat) was unchanged for most variant forms of P. jirovecii DHFR but was significantly lowered in F36C protein; one naturally occurring variant with two amino acid substitutions (S106P and E127G) showed a doubling of kcat, as well as a Km for NADPH half that of the wild type. The strongest resistance to trimethoprim occurred with amino acid changes in the binding pocket for DHFA or trimethoprim, and the strongest effect on binding of NADPH was linked to a mutation involved in binding the phosphate group of the cofactor. This study marks the first confirmation that naturally occurring mutations in the gene for DHFR from P. jirovecii produce variant forms of DHFR that are resistant to trimethoprim and may contribute to clinically observed failures of standard therapy or prophylaxis.

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

Design, synthesis, and molecular modeling of novel pyrido[2,3-d]pyrimidine analogues as antifolates; application of Buchwald-Hartwig aminations of heterocycles

Opportunistic infections caused by Pneumocystis jirovecii (P. jirovecii, pj), Toxoplasma gondii (T. gondii, tg), and Mycobacterium avium (M. avium, ma) are the principal causes of morbidity and mortality in patients with acquired immunodeficiency syndrome (AIDS). The absence of any animal models for human Pneumocystis jirovecii pneumonia and the lack of crystal structures of pjDHFR and tgDHFR make the design of inhibitors challenging. A novel series of pyrido[2,3-d]pyrimidines as selective and potent DHFR inhibitors against these opportunistic infections are presented. Buchwald-Hartwig coupling reaction of substituted anilines with pivaloyl protected 2,4-diamino-6-bromo-pyrido[2,3-d]pyrimidine was successfully explored to synthesize these analogues. Compound 26 was the most selective inhibitor with excellent potency against pjDHFR. Molecular modeling studies with a pjDHFR homology model explained the potency and selectivity of 26. Structural data are also reported for 26 with pcDHFR and 16 and 22 with variants of pcDHFR.

Kinetic and structural analysis for potent antifolate inhibition of Pneumocystis jirovecii, Pneumocystis carinii, and human dihydrofolate reductases and their active-site variants

A major concern of immunocompromised patients, in particular those with AIDS, is susceptibility to infection caused by opportunistic pathogens such as Pneumocystis jirovecii, which is a leading cause of pneumonia in immunocompromised patients. We report the first kinetic and structural data for 2,4-diamino-6-[(2',5'-dichloro anilino)methyl]pyrido[2,3-d]pyrimidine (OAAG324), a potent inhibitor of dihydrofolate reductase (DHFR) from P. jirovecii (pjDHFR), and also for trimethoprim (TMP) and methotrexate (MTX) with pjDHFR, Pneumocystis carinii DHFR (pcDHFR), and human DHFR (hDHFR). OAAG324 shows a 9.0-fold selectivity for pjDHFR (Ki, 2.7 nM) compared to its selectivity for hDHFR (Ki, 24.4 nM), whereas there is only a 2.3-fold selectivity for pcDHFR (Ki, 6.3 nM). In order to understand the determinants of inhibitory potency, active-site mutations of pj-, pc-, and hDHFR were explored to make these enzymes more like each other. The most unexpected observations were that the variant pcDHFR forms with K37Q and K37Q/F69N mutations, which made the enzyme more like the human form, also made these enzymes more sensitive to the inhibitory activity of OAAG324, with Ki values of 0.26 and 0.71 nM, respectively. A similar gain in sensitivity was also observed for the hDHFR N64F variant, which showed a lower Ki value (0.58 nM) than native hDHFR, pcDHFR, or pjDHFR. Structural data are reported for complexes of OAAG324 with hDHFR and its Q35K and Q35S/N64F variants and for the complex of the K37S/F69N variant of pcDHFR with TMP. These results provide useful insight into the role of these residues in the optimization of highly selective inhibitors of DHFR against the opportunistic pathogen P. jirovecii.

Structure-activity relationship for enantiomers of potent inhibitors of B. anthracis dihydrofolate reductase

BACKGROUND

Bacterial resistance to antibiotic therapies is increasing and new treatment options are badly needed. There is an overlap between these resistant bacteria and organisms classified as likely bioterror weapons. For example, Bacillus anthracis is innately resistant to the anti-folate trimethoprim due to sequence changes found in the dihydrofolate reductase enzyme. Development of new inhibitors provides an opportunity to enhance the current arsenal of anti-folate antibiotics while also expanding the coverage of the anti-folate class.

METHODS

We have characterized inhibitors of B. anthracis dihydrofolate reductase by measuring the K(i) and MIC values and calculating the energetics of binding. This series contains a core diaminopyrimidine ring, a central dimethoxybenzyl ring, and a dihydrophthalazine moiety. We have altered the chemical groups extended from a chiral center on the dihydropyridazine ring of the phthalazine moiety. The interactions for the most potent compounds were visualized by X-ray structure determination.

RESULTS

We find that the potency of individual enantiomers is divergent with clear preference for the S-enantiomer, while maintaining a high conservation of contacts within the binding site. The preference for enantiomers seems to be predicated largely by differential interactions with protein residues Leu29, Gln30 and Arg53.

CONCLUSIONS

These studies have clarified the activity of modifications and of individual enantiomers, and highlighted the role of the less-active R-enantiomer in effectively diluting the more active S-enantiomer in racemic solutions. This directly contributes to the development of new antimicrobials, combating trimethoprim resistance, and treatment options for potential bioterrorism agents.

Protein pharmacophore selection using hydration-site analysis

Virtual screening using pharmacophore models is an efficient method to identify potential lead compounds for target proteins. Pharmacophore models based on protein structures are advantageous because a priori knowledge of active ligands is not required and the models are not biased by the chemical space of previously identified actives. However, in order to capture most potential interactions between all potentially binding ligands and the protein, the size of the pharmacophore model, i.e. number of pharmacophore elements, is typically quite large and therefore reduces the efficiency of pharmacophore based screening. We have developed a new method to select important pharmacophore elements using hydration-site information. The basic premise is that ligand functional groups that replace water molecules in the apo protein contribute strongly to the overall binding affinity of the ligand, due to the additional free energy gained from releasing the water molecule into the bulk solvent. We computed the free energy of water released from the binding site for each hydration site using thermodynamic analysis of molecular dynamics (MD) simulations. Pharmacophores which are colocalized with hydration sites with estimated favorable contributions to the free energy of binding are selected to generate a reduced pharmacophore model. We constructed reduced pharmacophore models for three protein systems and demonstrated good enrichment quality combined with high efficiency. The reduction in pharmacophore model size reduces the required screening time by a factor of 200-500 compared to using all protein pharmacophore elements. We also describe a training process using a small set of known actives to reliably select the optimal set of criteria for pharmacophore selection for each protein system.

Crystallographic analysis reveals a novel second binding site for trimethoprim in active site double mutants of human dihydrofolate reductase

In order to produce a more potent replacement for trimethoprim (TMP) used as a therapy for Pneumocystis pneumonia and targets dihydrofolate reductase from Pneumocystis jirovecii (pjDHFR), it is necessary to understand the determinants of potency and selectivity against DHFR from the mammalian host and fungal pathogen cells. To this end, active site residues in human (h) DHFR were replaced with those from pjDHFR. Structural data are reported for two complexes of TMP with the double mutants Gln35Ser/Asn64Phe (Q35S/N64F) and Gln35Lys/Asn64Phe (Q35K/N64F) of hDHFR that unexpectedly show evidence for the binding of two molecules of TMP: one molecule that binds in the normal folate binding site and the second molecule that binds in a novel subpocket site such that the mutated residue Phe64 is involved in van der Waals contacts to the trimethoxyphenyl ring of the second TMP molecule. Kinetic data for the binding of TMP to hDHFR and pjDHFR reveal an 84-fold selectivity of TMP against pjDHFR (K(i) 49 nM) compared to hDHFR (K(i) 4093 nM). Two mutants that contain one substitution from pj--and one from the closely related Pneumocystis carinii DHFR (pcDHFR) (Q35K/N64F and Q35S/N64F) show K(i) values of 593 and 617 nM, respectively; these K(i) values are well above both the K(i) for pjDHFR and are similar to pcDHFR (Q35K/N64F and Q35S/N64F) (305nM). These results suggest that active site residues 35 and 64 play key roles in determining selectivity for pneumocystis DHFR, but that other residues contribute to the unique binding of inhibitors to these enzymes.

Structural analysis of Pneumocystis carinii and human DHFR complexes with NADPH and a series of five potent 6-[5'-(ω-carboxyalkoxy)benzyl]pyrido[2,3-d]pyrimidine derivatives

Structural data are reported for five antifolates, namely 2,4-diamino-6-[5'-(5-carboxypentyloxy)-2'-methoxybenzyl]-5-methylpyrido[2,3-d]pyrimidine, (1), and the 5'-[3-(ethoxycarbonyl)propoxy]-, (2), 5'-[3-(ethoxycarbonyl)butoxy]-, (3), 5'-[3-(ethoxycarbonyl)pentyloxy]-, (4), and 5'-benzyloxy-, (5), derivatives, which are potent and selective for Pneumocystis carinii dihydrofolate reductase (pcDHFR). Crystal structures are reported for their ternary complexes with NADPH and pcDHFR refined to between 1.4 and 2.0 Å resolution and for that of 3 with human DHFR (hDHFR) to 1.8 Å resolution. These data reveal that the carboxylate of the ω-carboxyalkoxy side chain of 1, the most potent inhibitor in this series, forms ionic interactions with the conserved Arg75 in the substrate-binding pocket of pcDHFR, whereas the less potent ethyl esters of 2-4 bind with variable side-chain conformations. The benzyloxy side chain of 5 makes no contact with Arg75 and is the least active inhibitor in this series. These structural results suggest that the weaker binding of this series compared with that of their pyrimidine homologs in part arises from the flexibility observed in their side-chain conformations, which do not optimize intermolecular contact to Arg75. Structural data for the binding of 3 to both hDHFR and pcDHFR reveals that the inhibitor binds in two different conformations, one similar to each of the two conformations observed for the parent pyrido[2,3-d]pyrimidine, piritrexim (PTX), bound to hDHFR. The structure of the pcDHFR complex of 4 reveals disorder in the side-chain orientation; one orientation has the ω-carboxyalkoxy side chain positioned in the folate-binding pocket similar to the others in this series, while the second orientation occupies a new site near the nicotinamide ring of NADPH. This alternate binding site has not been observed in other DHFR structures. Structural data for the pcDHFR complex of 5 show that its benzyl side chain forms intermolecular van der Waals interactions with Phe69 in the binding pocket that could account for its enhanced binding selectivity compared with the other analogs in this series.

Correlations of inhibitor kinetics for Pneumocystis jirovecii and human dihydrofolate reductase with structural data for human active site mutant enzyme complexes

To understand the role of specific active site residues in conferring selective dihydrofolate reductase (DHFR) inhibition from pathogenic organisms such as Pneumocystis carinii (pc) or Pneumocystis jirovecii (pj), the causative agent in AIDS pneumonia, it is necessary to evaluate the role of these residues in the human enzyme. We report the first kinetic parameters for DHFR from pjDHFR and pcDHFR with methotrexate (MTX), trimethoprim (TMP), and its potent analogue, PY957. We also report the mutagenesis and kinetic analysis of active site mutant proteins at positions 35 and 64 of human (h) DHFR and the crystal structure determinations of hDHFR ternary complexes of NADPH and PY957 with the wild-type DHFR enzyme, the single mutant protein, Gln35Lys, and two double mutant proteins, Gln35Ser/Asn64Ser and Gln35Ser/Asn64Phe. These substitutions place into human DHFR amino acids found at those sites in the opportunistic pathogens pcDHFR (Q35K/N64F) and pjDHFR (Q35S/N64S). The K(i) inhibition constant for PY957 showed greatest potency of the compound for the N64F single mutant protein (5.2 nM), followed by wild-type pcDHFR (K(i) 22 nM) and then wild-type hDHFR enzyme (K(i) 230 nM). Structural data reveal significant conformational changes in the binding interactions of PY957 in the hDHFR Q35S/N64F mutant protein complex compared to the other hDHFR mutant protein complexes and the pcDHFR ternary complex. The conformation of PY957 in the wild-type DHFR is similar to that observed for the single mutant protein. These data support the hypothesis that the enhanced selectivity of PY957 for pcDHFR is in part due to the contributions at positions 37 and 69 (pcDHFR numbering). This insight will help in the design of more selective inhibitors that target these opportunistic pathogens.

Crystal structure of Bacillus anthracis dihydrofolate reductase with the dihydrophthalazine-based trimethoprim derivative RAB1 provides a structural explanation of potency and selectivity

Bacillus anthracis possesses an innate resistance to the antibiotic trimethoprim due to poor binding to dihydrofolate reductase (DHFR); currently, there are no commercial antibacterials that target this enzyme in B. anthracis. We have previously reported a series of dihydrophthalazine-based trimethoprim derivatives that are inhibitors for this target. In the present work, we have synthesized one compound (RAB1) displaying favorable 50% inhibitory concentration (54 nM) and MIC (< or =12.8 microg/ml) values. RAB1 was cocrystallized with the B. anthracis DHFR in the space group P2(1)2(1)2(1), and X-ray diffraction data were collected to a 2.3-A resolution. Binding of RAB1 causes a conformational change of the side chain of Arg58 and Met37 to accommodate the dihydrophthalazine moiety. Unlike the natural substrate or trimethoprim, the dihydrophthalazine group provides a large hydrophobic anchor that embeds within the DHFR active site and accounts for its selective inhibitory activity against B. anthracis.

In pursuit of virtual lead optimization: pruning ensembles of receptor structures for increased efficiency and accuracy during docking

Representing receptors as ensembles of protein conformations during docking is a powerful method to approximate protein flexibility and increase the accuracy of the resulting ranked list of compounds. Unfortunately, docking compounds against a large number of ensemble members can increase computational cost and time investment. In this article, we present an efficient method to evaluate and select the most contributive ensemble members prior to docking for targets with a conserved core of residues that bind a ligand moiety. We observed that ensemble members that preserve the geometry of the active site core are most likely to place ligands in the active site with a conserved orientation, generally rank ligands correctly and increase interactions with the receptor. A relative distance approach is used to quantify the preservation of the three-dimensional interatomic distances of the conserved ligand-binding atoms and prune large ensembles quickly. In this study, we investigate dihydrofolate reductase as an example of a protein with a conserved core; however, this method for accurately selecting relevant ensemble members a priori can be applied to any system with a conserved ligand-binding core, including HIV-1 protease, kinases, and acetylcholinesterase. Representing a drug target as a pruned ensemble during in silico screening should increase the accuracy and efficiency of high-throughput analyses of lead analogs.

Mutational 'hot-spots' in mammalian, bacterial and protozoal dihydrofolate reductases associated with antifolate resistance: sequence and structural comparison

Human dihydrofolate reductase (DHFR) is a primary target for antifolate drugs in cancer treatment, while DHFRs from Plasmodium falciparum, Plasmodium vivax and various bacterial species are primary targets in the treatment of malaria and bacterial infections. Mutations in each of these DHFRs can result in resistance towards clinically relevant antifolates. We review the structural and functional impact of active-site mutations with respect to enzyme activity and antifolate resistance of DHFRs from mammals, protozoa and bacteria. The high structural homology between DHFRs results in a number of cross-species, active-site 'hot-spots' for broad-based antifolate resistance. In addition, we identify mutations that confer species-specific resistance, or antifolate-specific resistance. This comparative review of antifolate binding in diverse species provides new insights into the relationship between antifolate design and the development of mutational resistance. It also presents avenues for designing antifolate-resistant mammalian DHFRs as chemoprotective agents.

In pursuit of virtual lead optimization: the role of the receptor structure and ensembles in accurate docking

Accurate ranking during in silico lead optimization is critical to drive the generation of new ligands with higher affinity, yet it is especially difficult because of the subtle changes between analogs. In order to assess the role of the structure of the receptor in delivering accurate lead ranking results, we docked a set of forty related inhibitors to structures of one species of dihydrofolate reductase (DHFR) derived from crystallographic, NMR solution data, and homology models. In this study, the crystal structures yielded the superior results: the compounds were placed in the active site in the conserved orientation and the docking scores for 80% percent of the compounds clustered into the same bins as the measured affinity. Single receptor structures derived from NMR data or homology models did not serve as accurate docking receptors. To our knowledge, these are the first experiments that assess ranking of homologous lead compounds using a variety of receptor structures. We then extended the study to investigate whether ensembles, either computationally or experimentally derived, of all of the single starting structures aid, hinder or have no effect on the performance of the starting template. Impressively, when ensembles of receptor structures derived from NMR data or homology models were employed, docking accuracy improved to a level equal to that of the high resolution crystal structures. The same experiments using a second species of DHFR and set of ligands confirm the results. A comparison of the structures of the individual ensemble members to the starting structures shows that the effect of the ensembles can be ascribed to protein flexibility in addition to absorption of computational error.

N9-substituted 2,4-diaminoquinazolines: synthesis and biological evaluation of lipophilic inhibitors of pneumocystis carinii and toxoplasma gondii dihydrofolate reductase

N9-substituted 2,4-diaminoquinazolines were synthesized and evaluated as inhibitors of Pneumocystis carinii (pc) and Toxoplasma gondii (tg) dihydrofolate reductase (DHFR). Reduction of commercially available 2,4-diamino-6-nitroquinazoline 14 with Raney nickel afforded 2,4,6-triaminoquinazoline 15. Reductive amination of 15 with the appropriate benzaldehydes or naphthaldehydes, followed by N9-alkylation, afforded the target compounds 5- 13. In the 2,5-dimethoxybenzylamino substituted quinazoline analogues, replacement of the N9-CH 3 group of 4 with the N9-C2H5 group of 8 resulted in a 9- and 8-fold increase in potency against pcDHFR and tgDHFR, respectively. The N9-C2H5 substituted compound 8 was highly potent, with IC50 values of 9.9 and 3.7 nM against pcDHFR and tgDHFR, respectively. N9-propyl and N9-cyclopropyl methyl substitutions did not afford further increases in potency. This study indicates that the N9-ethyl substitution is optimum for inhibitory activity against pcDHFR and tgDHFR for the 2,4-diaminoquinazolines. Selectivity was unaffected by N9 substitution.

The effect of 5-alkyl modification on the biological activity of pyrrolo[2,3-d]pyrimidine containing classical and nonclassical antifolates as inhibitors of dihydrofolate reductase and as antitumor and/or antiopportunistic infection agents

Novel classical antifolates (3 and 4) and 17 nonclassical antifolates (11-27) were synthesized as antitumor and/or antiopportunistic infection agents. Intermediates for the synthesis of 3, 4, and 11-27 were 2,4-diamino-5-alkylsubstituted-7H-pyrrolo[2,3-d]pyrimidines, 31 and 38, prepared by a ring transformation/ring annulation sequence of 2-amino-3-cyano-4-alkyl furans to which various aryl thiols were attached at the 6-position via an oxidative addition reaction using I2. The condensation of alpha-hydroxy ketones with malonodinitrile afforded the furans. For the classical analogues 3 and 4, the ester precursors were deprotected, coupled with diethyl-L-glutamate, and saponified. Compounds 3 (IC50 = 60 nM) and 4 (IC50 = 90 nM) were potent inhibitors of human DHFR. Compound 3 inhibited tumor cells in culture with GI50 500-fold selectivity over human DHFR. Analogue 17 was 50-fold more potent than trimethoprim and about twice as selective against T. gondii DHFR.

Towards in silico lead optimization: scores from ensembles of protein/ligand conformations reliably correlate with biological activity

Accurately ranking protein/ligand interactions and distinguishing subtle differences between homologous compounds in a virtual focused library in silico is essential in a structure-based drug discovery program. In order to establish a predictive model to design novel inhibitors of dihydrofolate reductase (DHFR) from the parasitic protozoa, Cryptosporidium hominis, we docked a series of 30 DHFR inhibitors with measured inhibition constants against the crystal structure of the protein. By including protein flexibility and averaging the energies of the 25 lowest protein/ligand conformers we obtained more accurate total nonbonded energies from which we calculated a predicted biological activity. The calculated and measured biological activities showed reliable correlations of 72.9%. Additionally, visual analysis of the ensemble of protein/ligand conformations revealed alternative ligand binding pockets in the active site. Using the same principles we then created a homology model of DHFR from Toxoplasma gondii and docked 11 inhibitors. A correlation of 50.2% between docking score and activity validates both the method and the model. The correlations presented here are particularly compelling considering the high structural similarity of the ligands and the fact that we have used structures derived from crystallographic data and homology modeling. These docking principles may be useful in any lead optimization study where accurate ranking of similar compounds is desired.

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.

New insights into DHFR interactions: analysis of Pneumocystis carinii and mouse DHFR complexes with NADPH and two highly potent 5-(omega-carboxy(alkyloxy) trimethoprim derivatives reveals conformational correlations with activity and novel parallel ring stacking interactions

Structural data are reported for two highly potent antifolates, 2,4-diamino-5-[3',4'-dimethoxy-5'-(5-carboxy-1-pentynyl)]benzylpyrimidine (PY1011), with 5000-fold selectivity for Pneumocystis carinii dihydrofolate reductase (pcDHFR), relative to rat liver DHFR, and 2,4-diamino-5-[2-methoxy-5-(4-carboxybutyloxy)benzyl]pyrimidine (PY957), that has 80-fold selectivity for pcDHFR. Crystal structures are reported for NADPH ternary complexes with PY957 and pcDHFR, refined to 2.2 A resolution; with PY1011 and pcDHFR, refined to 2.0 A resolution; and with PY1011 and mouse DHFR (mDHFR), refined to 2.2 A resolution. These results reveal that the carboxylate of the omega-carboxyalkyloxy side chain of these inhibitors form ionic interactions with the conserved Arg in the substrate binding pocket of DHFR. These data suggest that the enhanced inhibitory activity of PY1011 compared with PY957 is, in part, due to the favorable contacts with Phe69 of pcDHFR by the methylene carbons of the inhibitor side chain that are oriented by the triple bond of the 1-pentynyl side chain. These contacts are not present in the PY957 pcDHFR complex, or in the PY1011 mDHFR complex. In the structure of mDHFR the site of Phe69 in pcDHFR is occupied by Asn64. These data also revealed a preference for an unusual parallel ring stacking interaction between Tyr35 of the active site helix and Phe199 of the C-terminal beta sheet in pcDHFR and by Tyr33 and Phe179 in mDHFR that is independent of bound ligand. A unique His174-His187 parallel ring stacking interaction was also observed only in the structure of pcDHFR. These ring stacking interactions are rarely found in any other protein families and may serve to enhance protein stability.

Highly Efficient Ligands for Dihydrofolate Reductase from Cryptosporidium hominis and Toxoplasma gondii Inspired by Structural Analysis

The search for effective therapeutics for cryptosporidiosis and toxoplasmosis has led to the discovery of novel inhibitors of dihydrofolate reductase (DHFR) that possess high ligand efficiency: compounds with high potency and low molecular weight. Detailed analysis of the crystal structure of dihydrofolate reductase-thymidylate synthase from Cryptosporidium hominis and a homology model of DHFR from Toxoplasma gondii inspired the synthesis of a new series of compounds with a propargyl-based linker between a substituted 2,4-diaminopyrimidine and a trimethoxyphenyl ring. An enantiomerically pure compound in this series exhibits IC50 values of 38 and 1 nM against C. hominis and T. gondii DHFR, respectively. Improvements of 368-fold or 5714-fold (C. hominis and T. gondii) relative to trimethoprim were generated by synthesizing just 14 new analogues and by adding only a total of 52 Da to the mass of the parent compound, creating an efficient ligand as an excellent candidate for further study.

The crystal structure of the bifunctional deaminase/reductase RibD of the riboflavin biosynthetic pathway in Escherichia coli: implications for the reductive mechanism

We have determined the crystal structure of the bi-functional deaminase/reductase enzyme from Escherichia coli (EcRibD) that catalyzes two consecutive reactions during riboflavin biosynthesis. The polypeptide chain of EcRibD is folded into two domains where the 3D structure of the N-terminal domain (1-145) is similar to cytosine deaminase and the C-terminal domain (146-367) is similar to dihydrofolate reductase. We showed that EcRibD is dimeric and compared our structure to tetrameric RibG, an ortholog from Bacillus subtilis (BsRibG). We have also determined the structure of EcRibD in two binary complexes with the oxidized cofactor (NADP(+)) and with the substrate analogue ribose-5-phosphate (RP5) and superposed these two in order to mimic the ternary complex. Based on this superposition we propose that the invariant Asp200 initiates the reductive reaction by abstracting a proton from the bound substrate and that the pro-R proton from C4 of the cofactor is transferred to C1 of the substrate. A highly flexible loop is found in the reductase active site (159-173) that appears to control cofactor and substrate binding to the reductase active site and was therefore compared to the corresponding Met20 loop of E. coli dihydrofolate reductase (EcDHFR). Lys152, identified by comparing substrate analogue (RP5) coordination in the reductase active site of EcRibD with the homologous reductase from Methanocaldococcus jannaschii (MjaRED), is invariant among bacterial RibD enzymes and could contribute to the various pathways taken during riboflavin biosynthesis in bacteria and yeast.

Analysis of complexes of inhibitors with Cryptosporidium hominis DHFR leads to a new trimethoprim derivative

Cryptosporidiosis, an opportunistic infection affecting immunocompromised patients, the elderly, and children, is still an untreatable disease since the causative agent, Cryptosporidium hominis, is essentially resistant to all clinically used antimicrobial agents. In order to accelerate the design of new potent and selective inhibitors targeting dihydrofolate reductase of C. hominis (ChDHFR), we determined the structural basis for the potency of existing DHFR inhibitors using superpositions of the structure of ChDHFR with other species and analysis of active site complexes of ChDHFR bound to ligands exhibiting a wide range of IC(50) values. This information was used to develop an accurate docking model capable of identifying potent inhibitors in silico. A series of C7-trimethoprim derivatives, designed to exploit a unique pocket in ChDHFR, was synthesized and evaluated; 7-ethyl TMP has four times higher activity than TMP against ChDHFR.

Structures of Leishmania major pteridine reductase complexes reveal the active site features important for ligand binding and to guide inhibitor design

Pteridine reductase (PTR1) is an NADPH-dependent short-chain reductase found in parasitic trypanosomatid protozoans. The enzyme participates in the salvage of pterins and represents a target for the development of improved therapies for infections caused by these parasites. A series of crystallographic analyses of Leishmania major PTR1 are reported. Structures of the enzyme in a binary complex with the cofactor NADPH, and ternary complexes with cofactor and biopterin, 5,6-dihydrobiopterin, and 5,6,7,8-tetrahydrobiopterin reveal that PTR1 does not undergo any major conformational changes to accomplish binding and processing of substrates, and confirm that these molecules bind in a single orientation at the catalytic center suitable for two distinct reductions. Ternary complexes with cofactor and CB3717 and trimethoprim (TOP), potent inhibitors of thymidylate synthase and dihydrofolate reductase, respectively, have been characterized. The structure with CB3717 reveals that the quinazoline moiety binds in similar fashion to the pterin substrates/products and dominates interactions with the enzyme. In the complex with TOP, steric restrictions enforced on the trimethoxyphenyl substituent prevent the 2,4-diaminopyrimidine moiety from adopting the pterin mode of binding observed in dihydrofolate reductase, and explain the inhibition properties of a range of pyrimidine derivates. The molecular detail provided by these complex structures identifies the important interactions necessary to assist the structure-based development of novel enzyme inhibitors of potential therapeutic value.

Design, synthesis, and antifolate activity of new analogues of piritrexim and other diaminopyrimidine dihydrofolate reductase inhibitors with omega-carboxyalkoxy or omega-carboxy-1-alkynyl substitution in the side chain

As part of a search for dihydrofolate reductase (DHFR) inhibitors combining the high potency of piritrexim (PTX) with the high antiparasitic vs mammalian selectivity of trimethoprim (TMP), the heretofore undescribed 2,4-diamino-6-(2',5'-disubstituted benzyl)pyrido[2,3-d]pyrimidines 6-14 with O-(omega-carboxyalkyl) or omega-carboxy-1-alkynyl groups on the benzyl moiety were synthesized and tested against Pneumocystis carinii, Toxoplasma gondii, and Mycobacterium avium DHFR vs rat DHFR. Three N-(2,4-diaminopteridin-6-yl)methyl)-2'-(omega-carboxy-1-alkynyl)dibenz[b,f]azepines (19-21) were also synthesized and tested. The pyridopyrimidine with the best combination of potency and selectivity was 2,4-diamino-5-methyl-6-[2'-(5-carboxy-1-butynyl)-5'-methoxy]benzyl]pyrimidine (13), with an IC(50) value of 0.65 nM against P. carinii DHFR, 0.57 nM against M. avium DHFR, and 55 nM against rat DHFR. The potency of 13 against P. carinii DHFR was 20-fold greater than that of PTX (IC(50) = 13 nM), and its selectivity index (SI) relative to rat DHFR was 85, whereas PTX was nonselective. The activity of 13 against P. carinii DHFR was 20 000 times greater than that of TMP, with an SI of 96, whereas that of TMP was only 14. However 13 was no more potent than PTX against M. avium DHFR, and its SI was no better than that of TMP. Molecular modeling dynamics studies using compounds 10 and 13 indicated a slight binding preference for the latter, in qualitative agreement with the IC(50) data. Among the pteridines, the most potent against P. carinii DHFR and M. avium DHFR was the 2'-(5-carboxy-1-butynyl)dibenz[b,f]azepinyl derivative 20 (IC(50) = 2.9 nM), whereas the most selective was the 2'-(5-carboxy-1-pentynyl) analogue 21, with SI values of >100 against both P. carinii and M. avium DHFR relative to rat DHFR. The final compound, 2,4-diamino-5-[3'-(4-carboxy-1-butynyl)-4'-bromo-5'-methoxybenzyl]pyrimidine (22), was both potent and selective against M. avium DHFR (IC(50) = 0.47 nM, SI = 1300) but was not potent or selective against either P. carinii or T. gondii DHFR.

Mutations of Pneumocystis jirovecii dihydrofolate reductase associated with failure of prophylaxis

Most drugs used for prevention and treatment of Pneumocystis jirovecii pneumonia target enzymes involved in the biosynthesis of folic acid, i.e., dihydropteroate synthase (DHPS) and dihydrofolate reductase (DHFR). Emergence of P. jirovecii drug resistance has been suggested by the association between failure of prophylaxis with sulfa drugs and mutations in DHPS. However, data on the occurrence of mutations in DHFR, the target of trimethoprim and pyrimethamine, are scarce. We examined polymorphisms in P. jirovecii DHFR from 33 patients diagnosed with P. jirovecii pneumonia who were receiving prophylaxis with a DHFR inhibitor (n = 15), prophylaxis without a DHFR inhibitor (n = 11), or no prophylaxis (n = 7). Compared to the wild-type sequence present in GenBank, 19 DHFR nucleotide substitution sites were found in 18 patients with 3 synonymous and 16 nonsynonymous mutations. Of 16 amino acid changes, 6 were located in positions conserved among distant organisms, and five of these six positions are probably involved in the putative active sites of the enzyme. Patients with failure of prophylaxis, including a DHFR inhibitor, were more likely to harbor nonsynonymous DHFR mutations than those who did not receive such prophylaxis (9 of 15 patients versus 2 of 18; P = 0.008). Analysis of the rate of nonsynonymous versus synonymous mutations was consistent with selection of amino acid substitutions in patients with failure of prophylaxis including a DHFR inhibitor. The results suggest that P. jirovecii populations may evolve under selective pressure from DHFR inhibitors, in particular pyrimethamine, and that DHFR mutations may contribute to P. jirovecii drug resistance.

Design, Synthesis, and Biological Evaluation of 2,4-Diamino-5-methyl-6-substituted-pyrrolo[2,3-d]pyrimidines as Dihydrofolate Reductase Inhibitors

2,4-diamino-5-methyl-6-(substituted-phenyl)thiopyrrolo[2,3-d]pyrimidines 4-11 were synthesized as dihydrofolate reductase (DHFR) inhibitors against opportunistic pathogens that afflict patients with AIDS. Synthesis was achieved from 2,4-diamino-5-methypyrrolo[2,3-d]pyrimidine and substituted phenylthiols under modified conditions reported by Gangjee et al. Some of these compounds were potent and selective against DHFR from both Toxoplasma gondii and Mycobacterium avium compared to mammalian DHFR. Compound 11 with a 1-naphthyl substituent is 16-fold more potent and equally selective against Toxoplasma gondii DHFR as the clinically used trimethoprim.

Design, synthesis, and computational affinity prediction of ester soft drugs as inhibitors of dihydrofolate reductase from Pneumocystis carinii

A series of dihydrofolate reductase (DHFR) inhibitors, where the methylenamino-bridge of non-classical inhibitors was replaced with an ester function, have been prepared as potential soft drugs intended for inhalation against Pneumocystis carinii pneumonia (PCP). Several of the new ester-based inhibitors that should serve as good substrates for the ubiquitous esterases and possibly constitute safer alternatives to metabolically stable DHFR inhibitors administered orally, were found to be potent inhibitors of P. carinii DHFR (pcDHFR). Although the objectives of the present program is to achieve a favorable toxicity profile by applying the soft drug concept, a high preference for inhibition of the fungal DHFR versus the mammalian DHFR is still desirable to suppress host toxicity at the site of administration. Compounds with a slight preference for the fungal enzyme were identified. The selection of the target compounds for synthesis was partly guided by an automated docking and scoring procedure as well as molecular dynamics simulations. The modest selectivity of the synthesized inhibitors was reasonably well predicted, although a correct ranking of the relative affinities was not successful in all cases.

Phylogenetic classification of protozoa based on the structure of the linker domain in the bifunctional enzyme, dihydrofolate reductase-thymidylate synthase

We have determined the crystal structure of dihydrofolate reductase-thymidylate synthase (DHFR-TS) from Cryptosporidium hominis, revealing a unique linker domain containing an 11-residue alpha-helix that has extensive interactions with the opposite DHFR-TS monomer of the homodimeric enzyme. Analysis of the structure of DHFR-TS from C. hominis and of previously solved structures of DHFR-TS from Plasmodium falciparum and Leishmania major reveals that the linker domain primarily controls the relative orientation of the DHFR and TS domains. Using the tertiary structure of the linker domains, we have been able to place a number of protozoa in two distinct and dissimilar structural families corresponding to two evolutionary families and provide the first structural evidence validating the use of DHFR-TS as a tool of phylogenetic classification. Furthermore, the structure of C. hominis DHFR-TS calls into question surface electrostatic channeling as the universal means of dihydrofolate transport between TS and DHFR in the bifunctional enzyme.

Insights into antifolate resistance from malarial DHFR-TS structures

Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (PfDHFR-TS) is an important target of antimalarial drugs. The efficacy of this class of DHFR-inhibitor drugs is now compromised because of mutations that prevent drug binding yet retain enzyme activity. The crystal structures of PfDHFR-TS from the wild type (TM4/8.2) and the quadruple drug-resistant mutant (V1/S) strains, in complex with a potent inhibitor WR99210, as well as the resistant double mutant (K1 CB1) with the antimalarial pyrimethamine, reveal features for overcoming resistance. In contrast to pyrimethamine, the flexible side chain of WR99210 can adopt a conformation that fits well in the active site, thereby contributing to binding. The single-chain bifunctional PfDHFR-TS has a helical insert between the DHFR and TS domains that is involved in dimerization and domain organization. Moreover, positively charged grooves on the surface of the dimer suggest a function in channeling of substrate from TS to DHFR active sites. These features provide possible approaches for the design of new drugs to overcome antifolate resistance.

Synthesis of classical and nonclassical, partially restricted, linear, tricyclic 5-deaza antifolates

Seven novel 2,4-diamino-5-deaza-6,7,8,9-tetrahydropyrido[3,4-g]pteridine derivatives 3-9 with different benzyl and a benzoyl substitution at the N7 position were designed and synthesized, as classical and nonclassical, partially restricted, linear tricyclic 5-deaza antifolates. The purpose was to investigate the effect of conformational restriction of the C6-C9 (tau(1)) and C9-N10 (tau(2)) bonds via an ethyl bridge from the N10 to the C7 position of 5-deaza methotrexate (MTX) on the inhibitory potency against dihydrofolate reductase (DHFR) from different sources and on antitumor activity. The synthetic methodology for most of the target compounds was a concise five-step total synthesis to construct the tricyclic nucleus, 2,4-diamino-5-deaza-7H-6,7,8,9-tetrahydropyrido[3,4-g]pteridine (23), followed by regioselective alkylation of the N7 nitrogen. Biological results indicated that this partial conformational modification for the classical analogue N-[4-[(2,4-diamino-5-deaza-6,7,8,9-tetrahydropyrido[3,4-g]pteridin-7-yl)methyl]benzoyl]-L-glutamic acid 3 was detrimental to DHFR inhibitory activity as well as to antitumor activity compared to MTX or 5-deaza MTX. However, the classical analogue 3 was a better substrate for folypolyglutamate synthetase (FPGS) than MTX. These results show that a classical 5-deaza folate partially restricted via a bridge between the N10 and C7 positions retains FPGS substrate activity and that the antitumor activity of classical tricyclic analogues such as 3 would be influenced by FPGS levels in tumor systems. Interestingly, the nonclassical analogues 4-9 showed moderate to good selectivity against DHFR from pathogenic microbes compared to recombinant human DHFR. These results support the idea that removal of the 5-methyl group of piritrexim along with restriction of tau(1) and tau(2) can translate into selectivity for DHFR from pathogens.

Structural studies on bioactive compounds. Part 36: design, synthesis and biological evaluation of pyrimethamine-based antifolates against Pneumocystis carinii

As part of a research effort to improve the quality of current chemotherapy of Pneumocystis carinii pneumonia, we report a structure-based design project to optimise activity, species selectivity and pharmaceutical properties of the triazenyl-pyrimethamine TAB (4) (IC(50)=0.17 microM; rat liver DHFR IC(50)/P. carinii DHFR IC(50)=114). This has led us to design, synthesise and evaluate four new series of pyrimethamine derivatives bearing triazole, triazolium, triazinium and amino moieties at the 3'-position of the p-chlorophenyl ring. Such stabilised 'triazene' derivatives address the potentially compromised pharmaceutical profile of TAB and the 3'-amine substituted agents afford conformationally flexible substitutes. The benzylamino-pyrimethamine derivative (24a) (IC(50)=0.12 microM, rat liver DHFR IC(50)/P. carinii DHFR IC(50): 5.26) was the most potent and the only P. carinii-selective antifolate of the new series.

A structurally conserved water molecule in Rossmann dinucleotide-binding domains

A computational comparison of 102 high-resolution (</=1.90 A) enzyme-dinucleotide (NAD, NADP, FAD) complexes was performed to investigate the role of solvent in dinucleotide recognition by Rossmann fold domains. The typical binding site contains about 9-12 water molecules, and about 30% of the hydrogen bonds between the protein and the dinucleotide are water mediated. Detailed inspection of the structures reveals a structurally conserved water molecule bridging dinucleotides with the well-known glycine-rich phosphate-binding loop. This water molecule displays a conserved hydrogen-bonding pattern. It forms hydrogen bonds to the dinucleotide pyrophosphate, two of the three conserved glycine residues of the phosphate-binding loop, and a residue at the C-terminus of strand four of the Rossmann fold. The conserved water molecule is also present in high-resolution structures of apo enzymes. However, the water molecule is not present in structures displaying significant deviations from the classic Rossmann fold motif, such as having nonstandard topology, containing a very short phosphate-binding loop, or having alpha-helix "A" oriented perpendicular to the beta-sheet. Thus, the conserved water molecule appears to be an inherent structural feature of the classic Rossmann dinucleotide-binding domain.

Atomic structures of human dihydrofolate reductase complexed with NADPH and two lipophilic antifolates at 1.09 a and 1.05 a resolution

The crystal structures of two human dihydrofolate reductase (hDHFR) ternary complexes, each with bound NADPH cofactor and a lipophilic antifolate inhibitor, have been determined at atomic resolution. The potent inhibitors 6-([5-quinolylamino]methyl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9439) and (Z)-6-(2-[2,5-dimethoxyphenyl]ethen-1-yl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9662) were developed at Southern Research Institute against Toxoplasma gondii DHFR-thymidylate synthase. The 5-deazapteridine ring of each inhibitor adopts an unusual puckered conformation that enables the formation of identical contacts in the active site. Conversely, the quinoline and dimethoxybenzene moieties exhibit distinct binding characteristics that account for the differences in inhibitory activity. In both structures, a salt-bridge is formed between Arg70 in the active site and Glu44 from a symmetry-related molecule in the crystal lattice that mimics the binding of methotrexate to DHFR.

Inhibition of Pneumocystis carinii, Toxoplasma gondii, and Mycobacterium avium dihydrofolate reductases by 2,4-diamino-5-[2-methoxy-5-(omega-carboxyalkyloxy)benzyl]pyrimidines: marked improvement in potency relative to trimethoprim and species selectivity relative to piritrexim

A series of previously undescribed 2,4-diamino-5-[2-methoxy-5-alkoxybenzyl]pyrimidines (3a-e) and 2,4-diamino-5-[2-methoxy-5-(omega-carboxyalkyloxy)benzyl]pyrimidines (3f-k) with up to eight CH2 groups in the alkoxy or omega-carboxyalkyloxy side chain were synthesized and tested for the ability to inhibit partially purified dihydrofolate reductase (DHFR) from Pneumocystis carinii (Pc), Toxoplasma gondii (Tg), Mycobacterium avium (Ma), and rat liver in comparison with two standard inhibitors, trimethoprim (1) and piritrexim (2). The latter drug is known to be extremely potent but shows a marked preference for binding to mammalian DHFR, whereas the former is very selective for the parasite enzymes but is a much weaker inhibitor. The underlying strategy for the synthesis of compounds 3a-k was that a hybrid structure embodying some features of both 1 and 2 might possess a more favorable combination of potency and selectivity than either parent drug. The choice of analogues 3f-k was based on the idea that the acidic omega-carboxyl group might interact preferentially with a basic center in the active site of DHFR from any of the parasite species relative to the active site of mammalian DHFR. In addition, the omega-carboxyl group was expected to improve water solubility relative to 1 or 2. In standardized spectrophotometric assays with dihydrofolate as the substrate and NADPH as the cofactor, 2,4-diamino-5-[(2-methoxy-4-carboxybutyloxy)benzyl]pyrimidine (3g) inhibited Pc DHFR with an IC(50) of 0.049 microM and rat DHFR with IC(50) of 3.9 microM. Its potency against Pc DHFR was 140-fold greater than that of 1 and close to that of 2, and its selectivity index, defined as the ratio IC(50)(rat liver)/IC(50)(P. carinii), was 8-fold higher than that of 1 and >10(4)-fold higher than that of 2. Although it was less potent and less selective against Tg than Pc DHFR, it was very potent as well as highly selective against Ma DHFR, with an IC(50) of 0.0058 microM and an IC(50)(rat liver)/IC(50)(M. avium) ratio of >600. Because of this favorable combination of potency and selectivity relative to 1 and 2, compound 3g may be viewed as a promising lead in the search for new antifolates with potential clinical activity against P. carinii and other opportunistic pathogens in patients with AIDS.

Isolation of rat dihydrofolate reductase gene and characterization of recombinant enzyme

While assays of many antifolate inhibitors for dihydrofolate reductase (DHFR) have been performed using rat DHFR as a target, neither the sequence nor the structure of rat DHFR is known. Here, we report the isolation of the rat DHFR gene through screening of a rat liver cDNA library. The rat liver DHFR gene has an open reading frame of 561 bp encoding a protein of 187 amino acids. Comparisons of the rat enzyme with those from other species indicate a high level of conservation at the primary sequence level and more so for the amino acid residues comprising the active site of the enzyme. Expression of the rat DHFR gene in bacteria produced a recombinant protein with high enzymatic activity. The recombinant protein also paralleled the human enzyme with respect to the inhibition by most of the antifolates tested with PT652 and PT653 showing a reversal in their patterns. Our results indicated that rat DHFR can be used as a model to study antifolate compounds as potential drug candidates. However, variations between rat and human DHFR enzymes, coupled with unique features in the inhibitors, could lead to the observed differences in enzyme sensitivity and selectivity.

Structural studies on bioactive compounds. 34. Design, synthesis, and biological evaluation of triazenyl-substituted pyrimethamine inhibitors of Pneumocystis carinii dihydrofolate reductase

The triazenyl-pyrimethamine derivative 3a (TAB), a potent and selective inhibitor of Pneumocystis carinii DHFR, was selected as the starting point for a lead optimization study. Molecular modeling studies, corroborated by a recent crystal structure determination of the ternary complex of P. carinii DHFR--NADPH bound to TAB, predicted that modifications to the acetoxy residue of the lead inhibitor could exploit binding opportunities in the vicinity of an active site pocket bounded by residues Ile33, Lys37, and Leu72. Substitutions in the benzyl moiety with electron-donating and electron-withdrawing groups were predicted to probe face-edge interactions with amino acid Phe69 unique to the P. carinii enzyme. New triazenes 10a--v and 12a--f were prepared by coupling the diazonium tetrafluoroborate salt 6b of aminopyrimethamine with substituted benzylamines or phenethylamines. The most potent of the new inhibitors against P. carinii DHFR was the naphthylmethyl-substituted triazene 10t (IC(50): 0.053 microM), but a more substantial increase in potency against the rat liver DHFR led to a reduction in selectivity (ratio rat liver DHFR IC(50)/P. carinii DHFR IC(50): 5.36) compared to the original lead structure 3a (ratio rat liver DHFR IC(50)/P. carinii DHFR IC(50): 114).