Drug-resistant enzymes must balance catalytic function with inhibitor destabilization to provide a fitness advantage. This sensitive balance, often involving very subtle structural changes, must be achieved through a selection process involving a minimal number of eligible point mutations. As part of a program to design propargyl-linked antifolates (PLAs) against trimethoprim-resistant dihydrofolate reductase (DHFR) from Staphylococcus aureus, we have conducted a thorough study of several clinically observed chromosomal mutations in the enzyme at the cellular, biochemical, and structural levels. Through this work, we have identified a promising lead series that displays significantly greater activity against these mutant enzymes and strains than TMP. The best inhibitors have enzyme inhibition and MIC values near or below that of trimethoprim against wild-type S. aureus. Moreover, these studies employ a series of crystal structures of several mutant enzymes bound to the same inhibitor; analysis of the structures reveals a more detailed molecular understanding of drug resistance in this important enzyme.

Although classical, negatively charged antifolates such as methotrexate possess high affinity for the dihydrofolate reductase (DHFR) enzyme, they are unable to penetrate the bacterial cell wall, rendering them poor antibacterial agents. Herein, we report a new class of charged propargyl-linked antifolates that capture some of the key contacts common to the classical antifolates while maintaining the ability to passively diffuse across the bacterial cell wall. Eight synthesized compounds exhibit extraordinary potency against Gram-positive S. aureus with limited toxicity against mammalian cells and good metabolic profile. High resolution crystal structures of two of the compounds reveal extensive interactions between the carboxylate and active site residues through a highly organized water network.Keywords: antifolate; Dihydrofolate reductase; Escherichia coli; methotrexate; MRSA; Staphylococcus aureus; trimethoprim

Multidrug-resistant Enterobacteriaceae, notably Escherichia coli and Klebsiella pneumoniae, have become major health concerns worldwide. Resistance to effective therapeutics is often carried by class I and II integrons that can confer insensitivity to carbapenems, extended spectrum β-lactamases, the antifolate trimethoprim, fluoroquinolones, and aminoglycosides. Specifically of interest to the study here, a prevalent gene (dfrA1) coding for an insensitive dihydrofolate reductase (DHFR) confers 190- or 1000-fold resistance to trimethoprim for K. pneumoniae and E. coli, respectively. Attaining inhibition of both the wild-type and resistant forms of the enzyme is critical for new antifolates. For several years, we have been developing the propargyl-linked antifolates (PLAs) as effective inhibitors against trimethoprim-resistant DHFR enzymes. Here, we show that the PLAs are active against both the wild-type and DfrA1 DHFR proteins. We report two high-resolution crystal structures of DfrA1 bound to potent PLAs. The structure–activity relationships and crystal structures will be critical in driving the design of broadly active inhibitors against wild-type and resistant DHFR.Keywords: antifolate; DfrA1; Enterobacteriaceae; resistance; trimethoprim

While antifolates such as Bactrim (trimethoprim-sulfamethoxazole; TMP-SMX) continue to play an important role in treating community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA), resistance-conferring mutations, specifically F98Y of dihydrofolate reductase (DHFR), have arisen and compromise continued use. In an attempt to extend the lifetime of this important class, we have developed a class of propargyl-linked antifolates (PLAs) that exhibit potent inhibition of the enzyme and bacterial strains. Probing the role of the configuration at the single propargylic stereocenter in these inhibitors required us to develop a new approach to nonracemic 3-aryl-1-butyne building blocks by the pairwise use of asymmetric conjugate addition and aldehyde dehydration protocols. Using this new route, a series of nonracemic PLA inhibitors was prepared and shown to possess potent enzyme inhibition (IC50 values <50 nM), antibacterial effects (several with MIC values <1 μg/mL) and to form stable ternary complexes with both wild-type and resistant mutants. Unexpectedly, crystal structures of a pair of individual enantiomers in the wild-type DHFR revealed that the single change in configuration of the stereocenter drove the selection of an alternative NADPH cofactor, with the minor α-anomer appearing with R-27. Remarkably, this cofactor switching becomes much more prevalent when the F98Y mutation is present. The observation of cofactor site plasticity leads to a postulate for the structural basis of TMP resistance in DHFR and also suggests design strategies that can be used to target these resistant enzymes.

Methods to accurately predict potential drug target mutations in response to early-stage leads could drive the design of more resilient first generation drug candidates. In this study, a structure-based protein design algorithm (K* in the OSPREY suite) was used to prospectively identify single-nucleotide polymorphisms that confer resistance to an experimental inhibitor effective against dihydrofolate reductase (DHFR) from Staphylococcus aureus. Four of the top-ranked mutations in DHFR were found to be catalytically competent and resistant to the inhibitor. Selection of resistant bacteria in vitro reveals that two of the predicted mutations arise in the background of a compensatory mutation. Using enzyme kinetics, microbiology, and crystal structures of the complexes, we determined the fitness of the mutant enzymes and strains, the structural basis of resistance, and the compensatory relationship of the mutations. To our knowledge, this work illustrates the first application of protein design algorithms to prospectively predict viable resistance mutations that arise in bacteria under antibiotic pressure.

Species of Candida, primarily C. albicans and with increasing prevalence, C. glabrata, are responsible for the majority of fungal bloodstream infections that cause morbidity, especially among immune compromised patients. While the development of new antifungal agents that target the essential enzyme, dihydrofolate reductase (DHFR), in both Candida species would be ideal, previous attempts have resulted in antifolates that exhibit inconsistencies between enzyme inhibition and antifungal properties. In this article, we describe the evaluation of pairs of propargyl-linked antifolates that possess similar physicochemical properties but different shapes. All of these compounds are effective at inhibiting the fungal enzymes and the growth of C. glabrata; however, the inhibition of the growth of C. albicans is shape-dependent with extended para-linked compounds proving more effective than compact, meta-linked compounds. Using crystal structures of DHFR from C. albicans and C. glabrata bound to lead compounds, 13 new para-linked compounds designed to inhibit both species were synthesized. Eight of these compounds potently inhibit the growth of both fungal species with three compounds displaying dual MIC values less than 1 μg/mL. Analysis of the active compounds shows that shape and distribution of polar functionality is critical in achieving dual antifungal activity.

Thujaplicins are tropolone-derived natural products with antiproliferative properties. We recently reported that certain tropolones potently and selectively target histone deacetylases (HDAC) and inhibit the growth of hematological cell lines. Here, we investigated the mechanisms by which these compounds exert their antiproliferative activity in comparison with the pan-selective HDAC inhibitor, vorinostat, using Jurkat T-cell leukemia cells. The tropolones appear to work through a mechanism distinct from vorinostat. These studies suggest that tropolone derivatives may serve as selective epigenetic modulators of hematological cells with potential applications as anti-leukemic or anti-inflammatory agents.

Natural products have long been recognized as a rich source of potent therapeutics but further development is often limited by high structural complexity and high molecular weight. In contrast, at the core of the thujaplicins is a lead-like tropolone scaffold characterized by relatively low molecular weight, ample sites for diversification, and metal-binding functionality poised for targeting a range of metalloenzyme drug targets. Here, we describe the development of this underutilized scaffold for the discovery of tropolone derivatives that function as isozyme-selective inhibitors of the validated anticancer drug target, histone deacetylase (HDAC). Several monosubstituted tropolones display remarkable levels of selectivity for HDAC2 and potently inhibit the growth of T-cell lymphocyte cell lines. The tropolones represent a new chemotype of isozyme-selective HDAC inhibitors.Keywords: HDAC; isozyme-selectivity; metalloenzyme; T-lymphocyte cancer cell lines; thujaplicin; Tropolone;

The pursuit of antimicrobial drugs that target dihydrofolate reductase (DHFR) exploits differences in sequence and dynamics between the pathogenic and human enzymes. Here, we present five crystal structures of human DHFR bound to a new class of antimicrobial agents, the propargyl-linked antifolates (PLAs), with a range of potency (IC50 values of 0.045–1.07 μM) for human DHFR. These structures reveal that interactions between the ligands and Asn 64, Phe 31, and Phe 34 are important for increased affinity for human DHFR and that loop residues 58–64 undergo ligand-induced conformational changes. The utility of these structural studies was demonstrated through the design of three new ligands that reduce the number of contacts with Asn 64, Phe 31, and Phe 34. Synthesis and evaluation show that one of the designed inhibitors exhibits the lowest affinity for human DHFR of any of the PLAs (2.64 μM). Comparisons of structures of human and Staphylococcus aureus DHFR bound to the same PLA reveal a conformational change in the ligand that enhances interactions with residues Phe 92 (Val 115 in huDHFR) and Ile 50 (Ile 60 in huDHFR) in S. aureus DHFR, yielding selectivity. Likewise, comparisons of human and Candida glabrata DHFR bound to the same ligand show that hydrophobic interactions with residues Ile 121 and Phe 66 (Val 115 and Asn 64 in human DHFR) yield selective inhibitors. The identification of residue substitutions that are important for selectivity and the observation of active site flexibility will help guide antimicrobial antifolate development for the inhibition of pathogenic species.

A novel strategy for targeting the pathogenic organisms Candida albicans and Candida glabrata focuses on the development of potent and selective antifolates effective against dihydrofolate reductase. Crystal structure analysis suggested that an essential loop at the active site (Thr 58-Phe 66) differs from the analogous residues in the human enzyme, potentially providing a mechanism for achieving selectivity. In order to probe the role of this loop, we employed chemical synthesis, crystal structure determination and molecular dynamics simulations. The results of these analyses show that the loop residues undergo ligand-induced conformational changes that are similar among the fungal and human species.

Over the past six decades, the folate biosynthetic pathway has provided a rich source of drug targets for the treatment of proliferative diseases. Drugs targeting dihydrofolate reductase have been especially successful as anticancer (methotrexate), antibacterial (trimethoprim, TMP) and antiprotozoal (cycloguanil, pyrimethamine) therapeutics. While trimethoprim remains a clinically important antimicrobial DHFR inhibitor, resistance by point mutations in otherwise sensitive strains as well as the natural insensitivity of several species limits its use. In this review, an historical overview of the attempts to develop drugs that target the folate pathway is presented along with a discussion of the basis of insensitivity to trimethoprim. From this vantage, we have developed the propargyl-linked antifolates as potent inhibitors of TMP-insensitive enzymes and strains. The structural basis of the increased affinity is detailed to promote the development of further generations of antifolates.

Enzymes are often excellent drug targets. Yet drug pressure on an enzyme target often fosters the rise of cells with resistance-conferring mutations, some of which may compromise fitness and others that compensate to restore fitness. This review presents, first, a structural analysis of a diverse group of wild-type and mutant enzyme targets and, second, an in-depth analysis of five diverse targets to elucidate a broader perspective of the effects of resistance-conferring mutations on protein or organismal fitness. The structural analysis reveals that resistance-conferring mutations may introduce steric hindrance or eliminate critical interactions, as expected, but that they may also have indirect effects such as altering protein dynamics and enzyme kinetics. The structure-based development of the latest generation of inhibitors targeting HIV reverse transcriptase, P. falciparum and S. aureus dihydrofolate reductase, neuraminidase, and epithelial growth factor receptor (EGFR) tyrosine kinase, is highlighted to emphasize lessons that may be applied to future drug discovery to overcome mutation-induced resistance. Successful next-generation drugs tend to be more flexible and exploit a greater number of interactions mimicking those of the substrate with conserved residues.

Dihydrofolate reductase (DHFR) has been a well-recognized target for the development of therapeutics for human cancers for several decades. Classical inhibitors of DHFR use an active transport mechanism to gain access to the cell; disabling this mechanism creates a pathway for resistance. In response, recent research focuses on nonclassical lipid-soluble DHFR inhibitors that are designed to passively diffuse through the membrane. Here, a new series of propargyl-linked antifolates are investigated as potential nonclassical human DHFR inhibitors. Several of these compounds exhibit potent enzyme inhibition with 50% inhibition concentration values under 500 nM. Molecular docking investigations show that the compounds maintain conserved hydrogen bonds between the pyrimidine ring and the enzyme as well as form van der Waals interactions with critical residues in the active site. Interestingly, the most potent compound, 2,4-diamino-5-(3-(3,4,5-trimethoxyphenyl)prop-1-ynyl)-6-ethylpyrimidine (compound 35), is 3500-fold more potent than trimethoprim, a potent inhibitor of bacterial DHFR but weak inhibitor of human DHFR. The two structural differences between compound 35 and trimethoprim show that the propargyl linkage and the substitution at C6 of the pyrimidine ring are critical to the formation of contacts with Thr 56, Ser 59, Ile 60, Leu 22, Phe 31 and Phe 34 and hence, to enhancing potency. The propargyl-linked antifolates are efficient ligands with a high ratio of potency to the number of non-hydrogen atoms and represent a potentially fruitful avenue for future development of antineoplastic agents.Graphical abstractResearch highlights▶ Dihydrofolate reductase (DHFR) is a target for the development of new antineoplastic agents. ▶ A new series of nonclassical antifolates, the propargyl-linked inhibitors, are explored for their ability to inhibit human DHFR. ▶ Several of these compounds exhibit IC50 values under 500 nM. ▶ Molecular docking experiments show that interactions between the C6 substituent on the pyrimidine ring and residues Phe 31 and Phe 34 in the active site are important for activity, as are interactions between the aryl ring and Thr 56, Ser 59 and Ile 60. ▶ The most potent of these compounds (compound 35) has two critical substitutions relative to trimethoprim, a known weak inhibitor of human DHFR but potent inhibitor of bacterial species of DHFR. ▶ Docked complexes of compound 35 exhibit more extensive interactions with Thr 56, Ser 59, Ile 60, Leu 22, Phe 31 and Phe 34 compared to trimethoprim, most likely explaining its greater potency. ▶ The new series of propargyl-linked antifolates are highly efficient ligands and represent a new avenue for the development of antineoplastic agents.

Several antifolates, including trimethoprim (TMP) and a series of propargyl-linked analogues, bind dihydrofolate reductase from Bacillus anthracis (BaDHFR) with lower affinity than is typical in other bacterial species. To guide lead optimization for BaDHFR, we explored a new approach to determine structure−activity relationships whereby the enzyme is altered and the analogues remain constant, essentially reversing the standard experimental design. Active site mutants of the enzyme, Ba(F96I)DHFR and Ba(Y102F)DHFR, were created and evaluated with enzyme inhibition assays and crystal structures. The affinities of the antifolates increase up to 60-fold with the Y102F mutant, suggesting that interactions with Tyr 102 are critical for affinity. Crystal structures of the enzymes bound to TMP and propargyl-linked inhibitors reveal the basis of TMP resistance and illuminate the influence of Tyr 102 on the lipophilic linker between the pyrimidine and aryl rings. Two new inhibitors test and validate these conclusions and show the value of the technique for providing new directions during lead optimization.

Drug resistance resulting from mutations to the target is an unfortunate common phenomenon that limits the lifetime of many of the most successful drugs. In contrast to the investigation of mutations after clinical exposure, it would be powerful to be able to incorporate strategies early in the development process to predict and overcome the effects of possible resistance mutations. Here we present a unique prospective application of an ensemble-based protein design algorithm, K^∗, to predict potential resistance mutations in dihydrofolate reductase from Staphylococcus aureus using positive design to maintain catalytic function and negative design to interfere with binding of a lead inhibitor. Enzyme inhibition assays show that three of the four highly-ranked predicted mutants are active yet display lower affinity (18-, 9-, and 13-fold) for the inhibitor. A crystal structure of the top-ranked mutant enzyme validates the predicted conformations of the mutated residues and the structural basis of the loss of potency. The use of protein design algorithms to predict resistance mutations could be incorporated in a lead design strategy against any target that is susceptible to mutational resistance.

Ensembles of protein structures to simulate protein flexibility are widely used throughout several applications including virtual lead optimization where they have been shown to improve ligand ranking. Yet, there is no established convention for weighting individual scores generated from ensemble members. To investigate the best method for weighting ensemble scores for proper ligand ranking, a series of dihydrofolate reductase inhibitors was docked to ensembles of Candida albicans dihydrofolate reductase (CaDHFR) structures created from a molecular dynamics (MD) simulation. From a single MD simulation, two ensemble collections were generated, one of which was subjected to a minimization procedure to create a group of structures of equal probability. As expected, ligand ranking accuracy was significantly improved when Boltzmann weighting was applied to the energies of the ensemble without structural minimization (60%), relative to that achieved with averaging (36%). However, accuracy was further improved (72%) by averaging docking scores across a minimized ensemble. To examine whether this accuracy results from structural variation in the single trajectory versus the possibility that error is minimized by averaging, a third collection of receptor structures was created in which each member was taken from an independent molecular dynamics simulation after minimization. Comparison of the docking accuracy results from the single trajectory (72%) to this third collection (61%) showed decreased accuracy, suggesting that ligands are more accurately oriented and assessed when docked to the minimized ensemble from a single MD trajectory, an effect that is more than simply error minimization. Averaging docking scores over a minimized ensemble of another target, influenza A neuraminidase, yielded a ligand ranking accuracy of 83%, representing a 24% improvement over other methods tested.

In order to develop new antifungal agents effective against two species of Candida, we have designed a series of dihydrofolate reductase (DHFR) inhibitors. Here, we explore the structure–activity relationships of these inhibitors toward Candida albicans DHFR by evaluating enzyme inhibition, antifungal activity and toxicity to mammalian cells. Analysis of docked complexes of the enzyme and inhibitors yields the structural basis of relative potency. The meta-biphenyl series of this class exhibits the greatest enzyme inhibition, selectivity and antifungal activity.A comparison of the active sites of dihydrofolate reductase from Candida albicans (blue), Candida glabrata (green) and human (red) with a modeled bound inhibitor of C. albicans and C. glabrata.

There is a significant need for new therapeutics to treat infections caused by the biodefense agent Bacillus anthracis. In pursuit of drug discovery against this organism, we have developed novel propargyl-linked inhibitors that target the essential enzyme dihydrofolate reductase (DHFR) from B. anthracis. Previously, we reported an initial series of these inhibitors and a high-resolution crystal structure of the ternary complex of the enzyme bound to its cofactor and one of the most potent inhibitors, UCP120B [Beierlein, J., Frey, K., Bolstad, D., Pelphrey, P., Joska, T., Smith, A., Priestley, N., Wright, D., and Anderson, A. (2008) J. Med. Chem. 51, 7532−7540]. Herein, we describe a three-dimensional solution structure of the ternary complex as determined by NMR. A comparison of this solution structure to the crystal structure reveals a general conservation of the DHFR fold and cofactor interactions as well as differences in the location of an active site helix and specific ligand interactions. In addition to data for the fully assigned ternary complex, data for the binary (enzyme−cofactor) complex were collected, providing chemical shift comparisons and revealing perturbations in residues that accommodate ligand binding. Dynamics of the protein, measured using 15N T1 and T2 relaxation times and {1H}−15N heteronuclear NOEs, reveal residue flexibility at the active site that explains enzyme inhibition and structure−activity relationships for two different series of these propargyl-linked inhibitors. The information obtained from the solution structure regarding active site flexibility will be especially valuable in the design of inhibitors with increased potency.

Bacillus anthracis, the causative agent of anthrax, poses a significant biodefense danger. Serious limitations in approved therapeutics and the generation of resistance have produced a compelling need for new therapeutic agents against this organism. Bacillus anthracis is known to be insensitive to the clinically used antifolate, trimethoprim, because of a lack of potency against the dihydrofolate reductase enzyme. Herein, we describe a novel lead series of B. anthracis dihydrofolate reductase inhibitors characterized by an extended trimethoprim-like scaffold. The best lead compound adds only 22 Da to the molecular weight and is 82-fold more potent than trimethoprim. An X-ray crystal structure of this lead compound bound to B. anthracis dihydrofolate reductase in the presence of NADPH was determined to 2.25 Å resolution. The structure reveals several features that can be exploited for further development of this lead series.

Cryptosporidiosis is an emerging infectious disease that can be life-threatening in an immune-compromised individual and causes gastrointestinal distress lasting up to 2 weeks in an immune-competent individual. There are few therapeutics available for effectively treating this disease. We have been exploring dihydrofolate reductase (DHFR) as a potential target in Cryptosporidium. On the basis of the structure of the DHFR enzyme from C. hominis, we have developed a novel scaffold that led to the discovery of potent (38 nM) and efficient inhibitors of this enzyme. Recently, we have advanced these inhibitors to the next stage of development. Using the structures of both the protozoal and human enzymes, we have developed inhibitors with nanomolar potency (1.1 nM) against the pathogenic enzyme and high levels (1273-fold) of selectivity over the human enzyme.

Candida glabrata is a lethal fungal pathogen resistant to many antifungal agents and has emerged as a critical target for drug discovery. Over the past several years, we have been developing a class of propargyl-linked antifolates as antimicrobials and hypothesized that these compounds could be effective inhibitors of dihydrofolate reductase (DHFR) from C. glabrata. We initially screened a small collection of these inhibitors and found modest levels of potency. Subsequently, we determined the crystal structure of C. glabrata DHFR bound to a representative inhibitor with data to 1.6 Å resolution. Using this structure, we designed and synthesized second-generation inhibitors. These inhibitors bind the C. glabrata DHFR enzyme with subnanomolar potency, display greater than 2000-fold levels of selectivity over the human enzyme, and inhibit the growth of C. glabrata at levels observed with clinically employed therapeutics.

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 IC50 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.A trimethoprim derivative with an ethyl substitution at the C7 position forms van der Waals contacts with Cys 113 and shows increased potency.

Resistance to therapeutics such as trimethoprim-sulfamethoxazole has become an increasing problem in strains of methicillin-resistant Staphylococcus aureus (MRSA). Clinically isolated trimethoprim-resistant strains reveal a double mutation, H30N/F98Y, in dihydrofolate reductase (DHFR). In order to develop novel and effective therapeutics against these resistant strains, we evaluated a series of propargyl-linked antifolate lead compounds for inhibition of the mutant enzyme. For the propargyl-linked antifolates, the F98Y mutation generates minimal (between 1.2- and 6-fold) losses of affinity and the H30N mutation generates greater losses (between 2.4- and 48-fold). Conversely, trimethoprim affinity is largely diminished by the F98Y mutation (36-fold) and is not affected by the H30N mutation. In order to elucidate a mechanism of resistance, we determined a crystal structure of a complex of this double mutant with a lead propargyl-linked antifolate. This structure suggests a resistance mechanism consistent both for the propargyl-linked class of antifolates and for trimethoprim that is based on the loss of a conserved water-mediated hydrogen bond.

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.