Structural characteristics of antifolate dihydrofolate reductase enzyme interactions

Oxidative degradation of dihydrofolate reductase increases CD38-mediated ferroptosis susceptibility

High expression of CD38 in tissues is a characteristic of aging, resulting in a decline in nicotinamide adenine dinucleotide (NAD) and increasing cellular reactive oxygen species (ROS). However, whether CD38 increases susceptibility to ferroptosis remains largely unexplored. Our previous study showed that CD38 overexpression decreased dihydrofolate reductase (DHFR). In the present study, we confirmed that high expression of CD38 increased ROS levels and induced DHFR degradation, which was prevented by nicotinamide mononucleotide (NMN) replenishment. We further revealed that ROS-mediated sulfonation on Cys7 of DHFR induced its degradation via the autophagy and non-canonical proteasome pathways. Mutation of Cys7 to alanine abolished ROS-induced DHFR degradation. Moreover, oxidative degradation of DHFR was responsible for the increased ferroptosis susceptibility of cells in which CD38 was highly expressed. We also found that CD38 expression was higher in bone-marrow-derived macrophages (BMDMs) from aged mice than those from young mice, while the DHFR level was lower. Consequently, we demonstrated that BMDMs from aged mice were more susceptible to ferroptosis that can be reverted by NMN replenishment, suggesting that CD38 high expression rendered cells more susceptible to ferroptosis. Taken together, these results indicated that CD38-mediated NAD⁺ decline promoted DHFR oxidative degradation, thus resulting in increased cellular susceptibility to ferroptosis and suggesting that NMN replenishment may protect macrophages from ferroptosis in aged mice.

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.

Synthesis, pharmacological evaluation and Molecular modelling studies of pregnenolone derivatives as inhibitors of human dihydrofolate reductase

In current study, we synthesized chalcone derivatives (13a-c) via base-catalyzed Claisen-Schmidt condensation reaction. We further treated diamino compounds with synthesized chalcones to produce 3,4-dihydropyrimidin-2(1H)-one (18a-c), 3,4-dihydropyrimidin-2(1H)-thione (19a-c) and 2-aminopyrimidine (20a-c) derivatives of pregnenolone by cyclization reaction. Cell viability test of synthesized steroidal chalcones and their pyrimidine and thiopyrimidine derivatives against human breast (MCF-7), human lung (A549) and human prostate (PC-3) cancer cell lines was performed using (4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), assay. Compounds were further evaluated for their inhibition potential against recombinant human DHFR (rhDHFR). All compounds showed activity from low micromolar to submicromolar range. Compound 20b with IC₅₀ value of 180 nM emerged as most potent compound against rhDHFR. Interaction of the newly synthesized pregnenolone derivatives with hDHFR and estrogen receptor alpha (ERα) were also explored via docking simulations. The overall results of hDHFR inhibition have shown that these analogues can be further optimized and developed as potent anticancer agents.

Exploring the ability of dihydropyrimidine-5-carboxamide and 5-benzyl-2,4-diaminopyrimidine-based analogues for the selective inhibition of L. major dihydrofolate reductase

To tackle leishmaniasis, search for efficient therapeutic drug targets should be pursued. Dihydrofolate reductase (DHFR) is considered as a key target for the treatment of leishmaniasis. In current study, we are interested in the design and synthesis of selective antifolates targeting DHFR from L. major. We focused on the development of new antifolates based on 3,4-dihydropyrimidine-2-one and 5-(3,5-dimethoxybenzyl)pyrimidine-2,4-diamine motif. Structure activity relationship (SAR) studies were performed on 4-phenyl ring of dihydropyrimidine (26-30) template. While for 5-(3,5-dimethoxybenzyl)pyrimidine-2,4-diamine, the impact of different amino acids (valine, tryptophan, phenylalanine, and glutamic acid) and two carbon linkers were explored (52-59). The synthesized compounds were assayed against LmDHFR. Compound 59 with the IC₅₀ value of 0.10 μM appeared as potent inhibitors of L. major. Selectivity for parasite DHFR over human DHFR was also determined. Derivatives 55-59 demonstrated excellent selectivity for LmDHFR. Highest selectivity for LmDHFR was shown by compounds 56 (SI = 84.5) and 58 (SI = 87.5). Compounds Antileishmanial activity against L. major and L. donovani promastigotes was also performed. To explore the interaction pattern of the synthesized compounds with biological macromolecules, the docking studies were carried out against homology modelled LmDHFR and hDHFR targets.

Targeting EGFR tyrosine kinase: Synthesis, in vitro antitumor evaluation, and molecular modeling studies of benzothiazole-based derivatives

New benzothiazole-based derivatives were synthesized in the present work with the aim of evaluating their antitumor activity. They were in vitro tested against hepatocellular carcinoma (HepG2), colorectal carcinoma (HCT-116), mammary gland cancer (MCF-7), prostate cancer (PC-3), and epithelioid carcinoma (HeLa). The results of the in vitro antitumor evaluation revealed that the most active compounds were 39, 40, 51, 56, and 61 exhibiting IC₅₀ values comparable to the reference drug lapatinib. The most active compounds were further subjected to EGFR inhibitory activity assay to rationalize their potency mode. Notably, the most active antitumor compounds 39 and 40 represented the most potent inhibitors to EGFR with IC₅₀ values of 24.58 and 30.42 nM respectively in comparison with 17.38 nM for lapatinib as a standard drug. Molecular modeling studies were also conducted for the synthesized compounds, including docking into EGFR active site and surface mapping. Results proved the superior binding of the hydrazone derivatives 39 and 40 with EGFR suggesting them as good candidates for targeted antitumor therapy through EGFR kinase inhibition.

Allostery and Epistasis: Emergent Properties of Anisotropic Networks

Understanding the underlying mechanisms behind protein allostery and non-additivity of substitution outcomes (i.e., epistasis) is critical when attempting to predict the functional impact of mutations, particularly at non-conserved sites. In an effort to model these two biological properties, we extend the framework of our metric to calculate dynamic coupling between residues, the Dynamic Coupling Index (DCI) to two new metrics: (i) EpiScore, which quantifies the difference between the residue fluctuation response of a functional site when two other positions are perturbed with random Brownian kicks simultaneously versus individually to capture the degree of cooperativity of these two other positions in modulating the dynamics of the functional site and (ii) DCI_asym, which measures the degree of asymmetry between the residue fluctuation response of two sites when one or the other is perturbed with a random force. Applied to four independent systems, we successfully show that EpiScore and DCI_asym can capture important biophysical properties in dual mutant substitution outcomes. We propose that allosteric regulation and the mechanisms underlying non-additive amino acid substitution outcomes (i.e., epistasis) can be understood as emergent properties of an anisotropic network of interactions where the inclusion of the full network of interactions is critical for accurate modeling. Consequently, mutations which drive towards a new function may require a fine balance between functional site asymmetry and strength of dynamic coupling with the functional sites. These two tools will provide mechanistic insight into both understanding and predicting the outcome of dual mutations.

Sensitivity of Functional Loop Conformations on Long-Range Electrostatics: Implications for M20 Loop Dynamics in E. coli Dihydrofolate Reductase

In E. coli dihydrofolate reductase, unusual conformational motions of a functional M20 loop that interacts with substrate and coenzyme have been construed as evidence for dynamical effects in enzyme catalysis. By computing this loop's conformational free energies in the apoenzyme, we show that it is sensitive to the treatment of long-range electrostatic interactions and the solvation box size in modeling/simulations. These results provide important guidelines for computing reaction/binding free energy profiles of proteins with functional loops.

Towards breast cancer targeting: Synthesis of tetrahydroindolocarbazoles, antibreast cancer evaluation, uPA inhibition, molecular genetic and molecular modelling studies

A series of some new tetrahydroindolocarbazole derivatives has been synthesized. The structure of the synthesized compounds has been confirmed by different spectroscopic techniques such as IR, NMR, elemental analysis and mass spectrometry. The target compounds were evaluated for their antitumor activity against breast cancer cell line MCF-7, their GI% and their LC₅₀ have been determined. Six of the synthesized compounds exhibited GI% values against MCF-7 cell lines exceeding 70% ranging from 71.9 to 85.0% in addition that compound 11 expressed GI% values of 99.9% and considered the most active derivatives among the synthesized ones. Compound 11 showed a remarkable decrease of u PA level to 3.5 ng/ml compared to DOX. Compound 5, 11 and 15 showed significant decrease in expression of MTAP and CDKN2A, in addition to a remarkable decrease in DNA damage comet assay method. Molecular modeling studies were performed to interpretate the behavior of active ligands as uPA inhibitors.

Medicines associated with folate-homocysteine-methionine pathway disruption

Folate is vital for cell development and growth. It is involved in one-carbon transfer reactions essential for the synthesis of purines and pyrimidines. It also acts in conjunction with cobalamin (vitamin B₁₂) as a fundamental cofactor in the remethylation cycle that converts homocysteine to methionine. A deficiency in folate or vitamin B₁₂ can lead to elevated homocysteine level, which has been identified as an independent risk factor in several health-related conditions. Adequate folate levels are essential in women of childbearing age and in pregnant women, and folate deficiency is associated with several congenital malformations. Low folate levels can be caused by dietary deficiencies, a genetic predisposition or treatment with medicines that affect folate concentration. Women who are pregnant or of child-bearing age commonly use medicines, so it is important to identify the basic biochemical mechanisms by which medicines interfere with the folate-homocysteine-methionine pathway. This review focuses on prescription medicines associated with folate disruption. It also summarizes their undesirable/toxic effects. Recommendations regarding folate supplementation during medical therapy are also reviewed.

Tetrahydroindolocarbazoles (THICZs) as new class of urokinase (uPA) inhibitors: Synthesis, anticancer evaluation, DNA-damage determination, and molecular modelling study

Tetrahydroindolocarbazoles (THICZs) with versatile substituents, have been designed, synthesized, structure characterized, then investigated for their in-vitro anticancer screening, urokinase inhibition (uPA) evaluated, DNA-damage determination was further explored. Compounds 5, 8, 10 and 17 displayed the most promising antitumor activities against the breast cancer cell line as compared to the standard drug, doxorubicin with IC₅₀ = 5.24 ± 0.37, 4.00 ± 0.52, 7.20 ± 0.90 and 9.60 ± 1.10 µg/ml (versus 3.30 ± 0.48 µg/ml for doxorubicin). Compounds 5, 8, 10 and 17 represents the most significant uPA inhibitors of our study with IC₅₀ of 3.80, 2.70. 4.75, 10.80 (ng/ml) respectively. The expression levels of CDKN2A gene were decreased in 8, 10 and 17 cell lines as compared to those in positive control samples. Cell lines treated with 5, 8, 10 and 17 clearly observed a high score of damaged DNA cells. A deeper examination revealed that our hetroaromatics showed an extensive hydrogen bonding interactions that is required in the S pocket which is important for activity Arg 217, Gly 219, Gly 216, Lys 143 and Ser 190. So we present THICZs as promising uPA inhibitors expected as significant promise for further development as anti-invasiveness drugs.

Targeting species specific amino acid residues: Design, synthesis and biological evaluation of 6-substituted pyrrolo[2,3-d]pyrimidines as dihydrofolate reductase inhibitors and potential anti-opportunistic infection agents

To combine the potency of trimetrexate (TMQ) or piritrexim (PTX) with the species selectivity of trimethoprim (TMP), target based design was carried out with the X-ray crystal structure of human dihydrofolate reductase (hDHFR) and the homology model of Pneumocystis jirovecii DHFR (pjDHFR). Using variation of amino acids such as Met33/Phe31 (in pjDHFR/hDHFR) that affect the binding of inhibitors due to their distinct positive or negative steric effect at the active binding site of the inhibitor, we designed a series of substituted-pyrrolo[2,3-d]pyrimidines. The best analogs displayed better potency (IC₅₀) than PTX and high selectivity for pjDHFR versus hDHFR, with 4 exhibiting a selectivity for pjDHFR of 24-fold.

Tyrosine kinase inhibition effects of novel Pyrazolo[1,5-a]pyrimidines and Pyrido[2,3-d]pyrimidines ligand: Synthesis, biological screening and molecular modeling studies

Tyrosine kinases are one of the most critical mediators in the signaling path way. Late studies have proved the part of tyrosine kinases in the pathophysiology of cancer diseases. This current research paper has focused on investigating the novel Pyrazolo[1,5-a]pyrimidines and Pyrido[2,3-d]pyrimidines as a small molecules that can inhibit tyrosine kinase in cancer cells. NCI protocol was applied to test the antitumor activity of such compounds. Leukemia and renal cancer cell lines proved to be sensitive to some derivatives such as 6b-d, 9a and 11 with GI% values ranging from 30.4 to 41.3%. In addition, compound 11 proved to be the most active against MCF-7 with GI% 62.5. The synthesized compounds were also evaluated for their inhibitory effects against EGFR kinase enzyme. Compound 9b proved to be the most active one among the synthesized series with inhibition % value of 81.72 at 25 nM concentration and IC₅₀ 8.4 nM which is very close to the reference drug Sorafenib. In vitro cytotoxicity test was also performed using the MCF-7 breast cell line. Computer modeling using the active site of tyrosine kinase as a template and the most active tyrosine kinase inhibitors were calculated. Docking studies of the synthesized compounds into the active site of EGFR kinase domain showed good agreement with the obtained biological results.

Synthesis, biological evaluation and molecular modeling study of some new methoxylated 2-benzylthio-quinazoline-4(3H)-ones as nonclassical antifolates

A new series of 2,3,6-substituted-quinazolin-4-ones was designed, synthesized, and evaluated for their in vitro DHFR inhibition, antimicrobial, and antitumor activities. Compounds 28 and 61 proved to be active DHFR inhibitors with IC50 0.02 and 0.01μM, respectively. Molecular modeling studies concluded that recognition with the key amino acid Phe34 is essential for binding and hence DHFR inhibition. Compounds 34, 56 and 66 showed broad spectrum antimicrobial activity comparable to Gentamicin and Ciprofloxacin. Compounds 40 and 64 showed broad spectrum antitumor activity toward several tumor cell lines and proved to be 10 fold more active than 5-FU, with GI50 MG-MID values of 2.2 and 2.4μM, respectively.

Structural analysis of dihydrofolate reductases enables rationalization of antifolate binding affinities and suggests repurposing possibilities

Antifolates are competitive inhibitors of dihydrofolate reductase (DHFR), a conserved enzyme that is central to metabolism and widely targeted in pathogenic diseases, cancer and autoimmune disorders. Although most clinically used antifolates are known to be target specific, some display a fair degree of cross-reactivity with DHFRs from other species. A method that enables identification of determinants of affinity and specificity in target DHFRs from different species and provides guidelines for the design of antifolates is currently lacking. To address this, we first captured the potential druggable space of a DHFR in a substructure called the 'supersite' and classified supersites of DHFRs from 56 species into 16 'site-types' based on pairwise structural similarity. Analysis of supersites across these site-types revealed that DHFRs exhibit varying extents of dissimilarity at structurally equivalent positions in and around the binding site. We were able to explain the pattern of affinities towards chemically diverse antifolates exhibited by DHFRs of different site-types based on these structural differences. We then generated an antifolate-DHFR network by mapping known high-affinity antifolates to their respective supersites and used this to identify antifolates that can be repurposed based on similarity between supersites or antifolates. Thus, we identified 177 human-specific and 458 pathogen-specific antifolates, a large number of which are supported by available experimental data. Thus, in the light of the clinical importance of DHFR, we present a novel approach to identifying differences in the druggable space of DHFRs that can be utilized for rational design of antifolates.

Utility of the Biosynthetic Folate Pathway for Targets in Antimicrobial Discovery

The need for new antimicrobials is great in face of a growing pool of resistant pathogenic organisms. This review will address the potential for antimicrobial therapy based on polypharmacological activities within the currently utilized bacterial biosynthetic folate pathway. The folate metabolic pathway leads to synthesis of required precursors for cellular function and contains a critical node, dihydrofolate reductase (DHFR), which is shared between prokaryotes and eukaryotes. The DHFR enzyme is currently targeted by methotrexate in anti-cancer therapies, by trimethoprim for antibacterial uses, and by pyrimethamine for anti-protozoal applications. An additional anti-folate target is dihyropteroate synthase (DHPS), which is unique to prokaryotes as they cannot acquire folate through dietary means. It has been demonstrated as a primary target for the longest standing antibiotic class, the sulfonamides, which act synergistically with DHFR inhibitors. Investigations have revealed most DHPS enzymes possess the ability to utilize sulfa drugs metabolically, producing alternate products that presumably inhibit downstream enzymes requiring the produced dihydropteroate. Recent work has established an off-target effect of sulfonamide antibiotics on a eukaryotic enzyme, sepiapterin reductase, causing alterations in neurotransmitter synthesis. Given that inhibitors of both DHFR and DHPS are designed to mimic their cognate substrate, which contain shared substructures, it is reasonable to expect such “off-target” effects. These inhibitors are also likely to interact with the enzymatic neighbors in the folate pathway that bind products of the DHFR or DHPS enzymes and/or substrates of similar substructure. Computational studies designed to assess polypharmacology reiterate these conclusions. This leads to hypotheses exploring the vast utility of multiple members of the folate pathway for modulating cellular metabolism, and includes an appealing capacity for prokaryotic-specific polypharmacology for antimicrobial applications.

Nonclassical antifolates, part 4. 5-(2-aminothiazol-4-yl)-4-phenyl-4H-1,2,4-triazole-3-thiols as a new class of DHFR inhibitors: synthesis, biological evaluation and molecular modeling study

A new series of compounds possessing 5-(2-aminothiazol-4-yl)-4-phenyl-4H-1,2,4-triazole-3-thiol skeleton was designed, synthesized, and evaluated for their in vitro DHFR inhibition, antimicrobial, antitumor and schistosomicidal activities. Four active compounds were allocated, the antibacterial 22 (comparable to gentamicin and ciprofloxacin), the schistosomicidal 29 (comparable to praziquantel), the DHFR inhibitor 34 (IC₅₀ 0.03 μM, 2.7 fold more active than MTX), and the antitumor 36 (comparable to doxorubicin). Molecular modeling studies concluded that recognition with key amino acid Leu4 and Val1 is essential for DHFR binding. Flexible alignment and surface mapping revealed that the obtained model could be useful for the development of new class of DHFR inhibitors.

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.

Structural analysis of Pneumocystis carinii dihydrofolate reductase complexed with NADPH and 2,4-diamino-6-[2-(5-carboxypent-1-yn-1-yl)-5-methoxybenzyl]-5-methylpyrido[2,3-d]pyrimidine

Structural data are reported for 2,4-diamino-6-[2-(5-carboxypent-1-yn-1-yl)-5-methoxybenzyl]-5-methylpyrido[2,3-d]pyrimidine (PY1014) complexed with Pneumocystis carinii dihydrofolate reductase (pcDHFR) refined to 1.8 Å resolution. These data reveal that the carboxylate of the ω-carboxyalkynyl side chain of PY1014, the most pcDHFR-selective analog in this series, forms ionic interactions with the conserved Arg75 in the substrate-binding pocket of pcDHFR. The reversal of the 2',5'-substitution pattern of this analog compared with the highly selective diaminopyrimidine analog PY1011 (i.e. the 5'-pentynylcarboxy-5'-methoxy pattern of PY1014 versus the 3',4'-dimethoxy-5'-pentynylcarboxy pattern of PY1011) is necessary to achieve optimal interaction with Arg75 as observed in other structures. The larger diaminopyrido[2,3-d]pyrimidine ring of PY1014 places the 5'-methoxy group closer to Leu25 and Ser64 than does the diaminopyrimidine ring of PY1011. The 5'-methoxy O atom forms a hydrogen bond to the amide of Leu25 (O···N, 2.7 Å) and the 5'-methoxy methyl group makes a hydrophobic contact of 3.1 Å with C(β) of Ser64. Although the IC(50) values of PY1014 and PY1011 are similar, inhibition data show that the selectivity of PY1011 for pcDHFR is significantly greater. The greater selectivity for pcDHFR compared with mammalian DHFR of these inhibitors is also influenced by the enhanced hydrophobic interactions of the side-chain methylene atoms with Phe69 of pcDHFR compared with Asn64 of mammalian DHFR.

Structural analysis of human dihydrofolate reductase as a binary complex with the potent and selective inhibitor 2,4-diamino-6-{2'-O-(3-carboxypropyl)oxydibenz[b,f]-azepin-5-yl}methylpteridine reveals an unusual binding mode

In order to understand the structure-activity profile observed for a series of substituted dibenz[b,f]azepine antifolates, the crystal structure of the binary complex of human dihydrofolate reductase (hDHFR) with the potent and selective inhibitor 2,4-diamino-6-{2'-O-(3-carboxypropyl)oxydibenz[b,f]-azepin-5-yl}methylpteridine (PT684) was determined to 1.8 Å resolution. These data revealed that the carboxylate side chain of PT684 occupies two alternate positions, neither of which interacts with the conserved Arg70 in the active-site pocket, which in turn hydrogen bonds to water. These observations are in contrast to those reported for the ternary complex of mouse DHFR (mDHFR) with NADPH [Cody et al. (2008), Acta Cryst. D64, 977-984], in which the 3-carboxypropyl side chain of PT684 was hydrolyzed to its hydroxyl derivative, PT684a. The crystallization conditions differed for the human and mouse DHFR crystals (100 mM K2HPO4 pH 6.9, 30% ammonium sulfate for hDHFR; 15 mM Tris pH 8.3, 75 mM sodium cacodylate, PEG 4K for mDHFR). Additionally, the side chains of Phe31 and Gln35 in the hDHFR complex have a single conformation, whereas in the mDHFR complex they occupied two alternative conformations. These data show that the hDHFR complex has a decreased active-site volume compared with the mDHFR complex, as reflected in a relative shift of helix C (residues 59-64) of 1.2 Å, and a shift of 1.5 Å compared with the ternary complex of Pneumocystis carinii DHFR (pcDHFR) with the parent dibenz[b,f]azepine PT653. These data suggest that the greater inhibitory potency of PT684 against pcDHFR is consistent with the larger active-site volume of pcDHFR and the predicted interactions of the carboxylate side chain with Arg75.

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.

Targeted mutations of Bacillus anthracis dihydrofolate reductase condense complex structure−activity relationships

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.

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.

Structural comparison of chromosomal and exogenous dihydrofolate reductase from Staphylococcus aureus in complex with the potent inhibitor trimethoprim

Dihydrofolate reductase (DHFR) is the enzyme responsible for the NADPH-dependent reduction of 5,6-dihydrofolate to 5,6,7,8-tetrahydrofolate, an essential cofactor in the synthesis of purines, thymidylate, methionine, and other key metabolites. Because of its importance in multiple cellular functions, DHFR has been the subject of much research targeting the enzyme with anticancer, antibacterial, and antimicrobial agents. Clinically used compounds targeting DHFR include methotrexate for the treatment of cancer and diaminopyrimidines (DAPs) such as trimethoprim (TMP) for the treatment of bacterial infections. DAP inhibitors of DHFR have been used clinically for >30 years and resistance to these agents has become widespread. Methicillin-resistant Staphylococcus aureus (MRSA), the causative agent of many serious nosocomial and community acquired infections, and other gram-positive organisms can show resistance to DAPs through mutation of the chromosomal gene or acquisition of an alternative DHFR termed "S1 DHFR." To develop new therapies for health threats such as MRSA, it is important to understand the molecular basis of DAP resistance. Here, we report the crystal structure of the wild-type chromosomal DHFR from S. aureus in complex with NADPH and TMP. We have also solved the structure of the exogenous, TMP resistant S1 DHFR, apo and in complex with TMP. The structural and thermodynamic data point to important molecular differences between the two enzymes that lead to dramatically reduced affinity of DAPs to S1 DHFR. These differences in enzyme binding affinity translate into reduced antibacterial activity against strains of S. aureus that express S1 DHFR.

Synthetic and crystallographic studies of a new inhibitor series targeting Bacillus anthracis dihydrofolate reductase

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 A resolution. The structure reveals several features that can be exploited for further development of this lead series.

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

Comparative molecular field analysis (CoMFA) and quantum chemical calculations were performed on cycloguanil (Cyc) derivatives of the wild type and the quadruple mutant (Asn51Ile, Cys59Arg, Ser108Asn, Ile164Leu) of Plasmodium falciparum dihydrofolate reductase (PfDHFR). The represented CoMFA models of wild type (r(2) = 0.727 and r(2) = 0.985) and mutant type (r(2) = 0.786 and r(2) = 0.979) can describe the differences of the Cyc structural requirements for the two types of PfDHFR enzymes and can be useful to guide the design of new inhibitors. Moreover, the obtained particular interaction energies between the Cyc and the surrounding residues in the binding pocket indicated that Asn108 of mutant enzyme was the cause of Cyc resistance by producing steric clash with p-Cl of Cyc. Consequently, comparing the energy contributions with the potent flexible WR99210 inhibitor, it was found that the key mutant residue, Asn108, demonstrates attractive interaction with this inhibitor and some residues, Leu46, Ile112, Pro113, Phe116, and Leu119, seem to perform as second binding site with WR99210. Therefore, quantum chemical calculations can be useful for investigating residue interactions to clarify the cause of drug resistance.

Recombinant bovine dihydrofolate reductase produced by mutagenesis and nested PCR of murine dihydrofolate reductase cDNA

Recent reports of the slow-tight binding inhibition of bovine liver dihydrofolate reductase (bDHFR) in the presence of polyphenols isolated from green tea leaves has spurred renewed interest in the biochemical properties of bDHFR. Earlier studies were done with native bDHFR but in order to validate models of polyphenol binding to bDHFR, larger quantities of bDHFR are necessary to support structural studies. Bovine DHFR differs from its closest sequence homologue, murine DHFR, by 19 amino acids. To obtain the bDHFR cDNA, murineDHFR cDNA was transformed by a series of nested PCRs to reproduce the amino acid coding sequence for bovine DHFR. The bovine liver DHFR cDNA has an open reading frame of 561 base pairs encoding a protein of 187 amino acids that has a high level of conservation at the primary sequence level with other DHFR enzymes, and more so for the amino acid residues in the active site of the mammalian DHFR enzymes. Expression of the bovine DHFR cDNA in bacterial cells produced a stable recombinant protein with high enzymatic activity and kinetic properties similar to those previously reported for the native protein.

Regulation of human dihydrofolate reductase activity and expression

Dihydrofolate reductase (DHFR) enzyme catalyzes tetrahydrofolate regeneration by reduction of dihydrofolate using NADPH as a cofactor. Tetrahydrofolate and its one carbon adducts are required for de novo synthesis of purines and thymidylate, as well as glycine, methionine and serine. DHFR inhibition causes disruption of purine and thymidylate biosynthesis and DNA replication, leading to cell death. Therefore, DHFR has been an attractive target for chemotherapy of many diseases including cancer. Over the following years, in order to develop better antifolates, a detailed understanding of DHFR at every level has been undertaken such as structure-functional analysis, mechanisms of action, transcriptional and translation regulation of DHFR using a wide range of technologies. Because of this wealth of information created, DHFR has been used extensively as a model system for enzyme catalysis, investigating the relations between structure in-silico structure-based drug design, transcription from TATA-less promoters, regulation of transcription through the cell cycle, and translational autoregulation. In this review, the current understanding of human DHFR in terms of structure, function and regulation is summarized.