Prediction of residues involved in inhibitor specificity in the dihydrofolate reductase family

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Mutational analysis confirms the presence of distal inhibitor-selectivity determining residues in B. stearothermophilus dihydrofolate reductase

Many antibacterial and antiparasitic drugs work by competitively inhibiting dihydrofolate reductase (DHFR), a vital enzyme in folate metabolism. The interactions between inhibitors and DHFR active site residues are known in many homologs but the contributions from distal residues are less understood. Identifying distal residues that aid in inhibitor binding can improve targeted drug development programs by accounting for distant influences that may be less conserved and subject to frequent resistance causing mutations. Previously, a novel, homology-based, computational approach that mines ligand inhibition data was used to predict residues involved in inhibitor selectivity in the DHFR family. Expectedly, some inhibitor selectivity determining residue positions were predicted to lie in the active site and coincide with experimentally known inhibitor selectivity determining positions. However, other residues that group spatially in clusters distal to the active site have not been previously investigated. In this study, the effect of introducing amino acid substitutions at one of these predicted clusters (His38-Ala39-Ile40) on the inhibitor selectivity profile in Bacillus stearothermophilus dihydrofolate reductase (Bs DHFR) was investigated. Mutations were introduced into these cluster positions to change sidechain chemistry and size. We determined kcat and KM values and measured KD values at equilibrium for two competitive DHFR inhibitors, trimethoprim (TMP) and pyrimethamine (PYR). Mutations in the His38-Ala39-Ile40 cluster significantly impacted inhibitor binding and TMP/PYR selectivity - seven out of nine mutations resulted in tighter binding to PYR when compared to TMP. These data suggest that the His38-Ala39-Ile40 cluster is a distal inhibitor selectivity determining region that favors PYR binding in Bs DHFR and, possibly, throughout the DHFR family.

Expression, purification, and inhibition profile of dihydrofolate reductase from the filarial nematode Wuchereria bancrofti

Filariasis is a tropical disease caused by the parasitic nematodes Wuchereria bancrofti and Brugia malayi. Known inhibitors of dihydrofolate reductase (DHFR) have been previously shown to kill Brugia malayi nematodes and to inhibit Brugia malayi DHFR (BmDHFR) at nanomolar concentrations. These data suggest that BmDHFR is a potential target for the treatment of filariasis. Here, protocols for cloning, expression and purification of Wuchereria bancrofti DHFR (WbDHFR) were developed. The Uniprot entry J9F199-1 predicts a 172 amino acid protein for WbDHFR but alignment of this sequence to the previously described BmDHFR shows that this WbDHFR sequence lacks a crucial, conserved 13 amino acid loop. The presence of the loop in WbDHFR is supported by a noncanonical splicing event and the loop sequence was therefore included in the gene design. Subsequently, the K_M for dihydrofolate (3.7 ± 2 μM), k_cat (7.4 ± 0.6 s^-1), and pH dependence of activity were determined. IC₅₀ values of methotrexate, trimethoprim, pyrimethamine, raltitrexed, aminopterin, (-)-epicatechin gallate, (-)-epicatechin, and vitexin were measured for WbDHFR and BmDHFR. Methotrexate and structurally related aminopterin were found to be effective inhibitors of WbDHFR, with an K_I of 1.2 ± 0.2 nM and 2.1 ± 0.5 nM, respectively, suggesting that repurposing of known antifolate compound may be an effective strategy to treating filariasis. Most compounds showed similar inhibition profiles toward both enzymes, suggesting that the two enzymes have important similarities in their active site environments and can be targeted with the same compound, once a successful inhibitor is identified.

Exploiting the 2-Amino-1,3,4-thiadiazole Scaffold To Inhibit Trypanosoma brucei Pteridine Reductase in Support of Early-Stage Drug Discovery

Pteridine reductase-1 (PTR1) is a promising drug target for the treatment of trypanosomiasis. We investigated the potential of a previously identified class of thiadiazole inhibitors of Leishmania major PTR1 for activity against Trypanosoma brucei (Tb). We solved crystal structures of several TbPTR1-inhibitor complexes to guide the structure-based design of new thiadiazole derivatives. Subsequent synthesis and enzyme- and cell-based assays confirm new, mid-micromolar inhibitors of TbPTR1 with low toxicity. In particular, compound 4m, a biphenyl-thiadiazole-2,5-diamine with IC₅₀ = 16 μM, was able to potentiate the antitrypanosomal activity of the dihydrofolate reductase inhibitor methotrexate (MTX) with a 4.1-fold decrease of the EC₅₀ value. In addition, the antiparasitic activity of the combination of 4m and MTX was reversed by addition of folic acid. By adopting an efficient hit discovery platform, we demonstrate, using the 2-amino-1,3,4-thiadiazole scaffold, how a promising tool for the development of anti-T. brucei agents can be obtained.

Expression, purification and enzymatic characterization of Brugia malayi dihydrofolate reductase

Brugia malayi (B. malayi) is one of the three causative agents of lymphatic filariasis, a neglected parasitic disease. Current literature suggests that dihydrofolate reductase is a potential drug target for the elimination of B. malayi. Here we report the recombinant expression and purification of a ∼20 kDa B. malayi dihydrofolate reductase (BmDHFR). A His6-tagged construct was expressed in E. coli and purified by affinity chromatography to yield active and homogeneous enzyme for steady-state kinetic characterization and inhibition studies. The catalytic activity kcat was found to be 1.4 ± 0.1 s(-1), the Michaelis Menten constant KM for dihydrofolate 14.7 ± 3.6 μM, and the equilibrium dissociation constant KD for NADPH 25 ± 24 nM. For BmDHFR, IC50 values for a six DHFR inhibitors were determined to be 3.1 ± 0.2 nM for methotrexate, 32 ± 22 μM for trimethoprim, 109 ± 34 μM for pyrimethamine, 154 ± 46 μM for 2,4-diaminoquinazoline, 771 ± 44 μM for cycloguanil, and >20,000 μM for 2,4-diaminopyrimidine. Our findings suggest that antifolate compounds can serve as inhibitors of BmDHFR.

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.

Homology, similarity, and identity in peptide epitope immunodefinition

The tendency to use the terms homology, similarity, and identity interchangeably persists in comparative biology. When translated to immunology, overlapping the concepts of homology, similarity, and identity complicates the exact definition of the self-nonself dichotomy and, in particular, affects immunopeptidomics, an emerging field aimed at cataloging and distinguishing immunoreactive peptide epitopes from silent nonreactive amino acid sequences. The definition of similar/dissimilar peptides in immunology is discussed with special attention to the analysis of immunological (dis)similarity between two or more protein sequences that equates to measuring sequence similarity with the use of a proper measurement unit such as a length determinant.