Inhibition of dihydrofolate reductase from bacterial and vertebrate sources by folate, aminopterin, methotrexate and their 5-deaza analogues

Chemical space of Escherichia coli dihydrofolate reductase inhibitors: New approaches for discovering novel drugs for old bugs

Escherichia coli Dihydrofolate reductase is an important enzyme that is essential for the survival of the Gram-negative microorganism. Inhibitors designed against this enzyme have demonstrated application as antibiotics. However, either because of poor bioavailability of the small-molecules resulting from their inability to cross the double membrane in Gram-negative bacteria or because the microorganism develops resistance to the antibiotics by mutating the DHFR target, discovery of new antibiotics against the enzyme is mandatory to overcome drug-resistance. This review summarizes the field of DHFR inhibition with special focus on recent efforts to effectively interface computational and experimental efforts to discover novel classes of inhibitors that target allosteric and active-sites in drug-resistant variants of EcDHFR.

Ligand binding studies, preliminary structure-activity relationship and detailed mechanistic characterization of 1-phenyl-6,6-dimethyl-1,3,5-triazine-2,4-diamine derivatives as inhibitors of Escherichia coli dihydrofolate reductase

Gram-negative bacteria are implicated in the causation of life-threatening hospital-acquired infections. They acquire rapid resistance to multiple drugs and available antibiotics. Hence, there is the need to discover new antibacterial agents with novel scaffolds. For the first time, this study explores the 1,3,5-triazine-2,4-diamine and 1,2,4-triazine-2,4-diamine group of compounds as potential inhibitors of Escherichia coli DHFR, a pivotal enzyme in the thymidine and purine synthesis pathway. Using differential scanning fluorimetry, DSF, fifteen compounds with various substitutions on either the 3rd or 4th positions on the benzene group of 6,6-dimethyl-1-(benzene)-1,3,5-triazine-2,4-diamine were shown to bind to the enzyme with varying affinities. Then, the dose dependence of inhibition by these compounds was determined. Preliminary quantitative structure-activity relationship analysis and docking studies implicate the alkyl linker group and the sulfonyl fluoride group in increasing the potency of inhibition. 4-[4-[3-(4,6-diamino-2,2-dimethyl-1,3,5-triazin-1-yl)phenyl]butyl]benzenesulfonyl fluoride (NSC120927), the best hit from the study and a molecule with no reported inhibition of E. coli DHFR, potently inhibits the enzyme with a Ki value of 42.50 ± 5.34 nM, followed by 4-[6-[4-(4,6-diamino-2,2-dimethyl-1,3,5-triazin-1-yl)phenyl]hexyl]benzenesulfonyl fluoride (NSC132279), with a Ki value of 100.9 ± 12.7 nM. Detailed kinetic characterization of the inhibition brought about by five small-molecule hits shows that these inhibitors bind to the dihydrofolate binding site with preferential binding to the NADPH-bound binary form of the enzyme. Furthermore, in search of novel diaminotriazine scaffolds, it is shown that lamotrigine, a 1,2,4-triazine-3,5-diamine and a sodium-ion channel blocker class of antiepileptic drug, also inhibits E. coli DHFR. This is the first comprehensive study on the binding and inhibition brought about by diaminotriazines of a gram-negative prokaryotic enzyme and provides valuable insights into the SAR as an aid to the discovery of novel antibiotics.

PoLi: A Virtual Screening Pipeline Based on Template Pocket and Ligand Similarity

Often in pharmaceutical research the goal is to identify small molecules that can interact with and appropriately modify the biological behavior of a new protein target. Unfortunately, most proteins lack both known structures and small molecule binders, prerequisites of many virtual screening, VS, approaches. For such proteins, ligand homology modeling, LHM, that copies ligands from homologous and perhaps evolutionarily distant template proteins, has been shown to be a powerful VS approach to identify possible binding ligands. However, if we want to target a specific pocket for which there is no homologous holo template protein structure, then LHM will not work. To address this issue, in a new pocket-based approach, PoLi, we generalize LHM by exploiting the fact that the number of distinct small molecule ligand-binding pockets in proteins is small. PoLi identifies similar ligand-binding pockets in a holo template protein library, selectively copies relevant parts of template ligands, and uses them for VS. In practice, PoLi is a hybrid structure and ligand-based VS algorithm that integrates 2D fingerprint-based and 3D shape-based similarity metrics for improved virtual screening performance. On standard DUD and DUD-E benchmark databases, using modeled receptor structures, PoLi achieves an average enrichment factor of 13.4 and 9.6, respectively, in the top 1% of the screened library. In contrast, traditional docking-based VS using AutoDock Vina and homology-based VS using FINDSITE(filt) have an average enrichment of 1.6 (3.0) and 9.0 (7.9) on the DUD (DUD-E) sets, respectively. Experimental validation of PoLi predictions on dihydrofolate reductase, DHFR, using differential scanning fluorimetry, DSF, identifies multiple ligands with diverse molecular scaffolds, thus demonstrating the advantage of PoLi over current state-of-the-art VS methods.

Insights into the slow-onset tight-binding inhibition of Escherichia coli dihydrofolate reductase: detailed mechanistic characterization of pyrrolo [3,2-f] quinazoline-1,3-diamine and its derivatives as novel tight-binding inhibitors

Dihydrofolate reductase (DHFR) is a pivotal enzyme involved in the de novo pathway of purine synthesis, and hence, represents an attractive target to disrupt systems that require rapid DNA turnover. The enzyme acquires resistance to available drugs by various molecular mechanisms, which necessitates the continuous discovery of novel antifolates. Previously, we identified a set of novel molecules that showed binding to E. coli DHFR by means of a thermal shift without establishing whether they inhibited the enzyme. Here, we show that a fraction of those molecules represent potent and novel inhibitors of DHFR activity. 7-[(4-aminophenyl)methyl]-7H-pyrrolo [3,2-f] quinazoline-1,3-diamine, a molecule with no reported inhibition of DHFR, potently inhibits the enzyme with a Ki value of 7.42 ± 0.92 nm by competitive displacement of the substrate dihydrofolic acid. It shows uncompetitive inhibition vis-à-vis NADPH, indicating that the inhibitor has markedly increased affinity for the NADPH-bound form of the enzyme. Further, we demonstrate that the mode of binding of the inhibitor to the enzyme-NADPH binary complex conforms to the slow-onset, tight-binding model. By contrast, mechanistic characterization of the parent molecule 7H-pyrrolo [3,2-f] quinazoline-1,3-diamine shows that lack of (4-aminophenyl)-methyl group at the seventh position abolishes the slow onset of inhibition. This finding provides novel insights into the role of substitutions on inhibitors of E. coli DHFR and represents the first detailed kinetic investigation of a novel diaminopyrroloquinazoline derivative on a prokaryotic DHFR. Furthermore, marked differences in the potency of inhibition for E. coli and human DHFR makes this molecule a promising candidate for development as an antibiotic.

Replacing sulfa drugs with novel DHPS inhibitors

More research effort needs to be invested in antimicrobial drug development to address the increasing threat of multidrug-resistant organisms. The enzyme DHPS has been a validated drug target for over 70 years as the target for the highly successful sulfa drugs. The use of sulfa drugs has been compromised by the widespread presence of resistant organisms and the adverse side effects associated with their use. Despite the large amount of structural information available for DHPS, few recent publications address the possibility of using this knowledge for novel drug design. This article reviews the relevant papers and patents that report promising new small-molecule inhibitors of DHPS, and discuss these data in light of new insights into the DHPS catalytic mechanism and recently determined crystal structures of DHPS bound to potent small-molecule inhibitors. This new functional understanding confirms that DHPS deserves further consideration as an antimicrobial drug target.

Three-dimensional quantitative structure activity relationship (QSAR) of cytotoxic active 3,5-diaryl-4,5-dihydropyrazole analogs: a comparative molecular field analysis (CoMFA) revisited study

In vitro antitumor evaluation of the synthesized 46 compounds of 3,5-diaryl-4,5-dihydropyrazoles against EAC cell lines and 3D QSAR study using pharmacophore and Comparative Molecular Field Analysis (CoMFA) methods were described. CoMFA derived QSAR model shows a good conventional squared correlation coefficient r² and cross validated correlation coefficient r²cv 0.896 and 0.568 respectively. In this analysis steric and electrostatic field contribute to the QSAR equation by 70% and 30% respectively, suggesting that variation in biological activity of the compounds is dominated by differences in steric (van der Waals) interactions. To visualize the CoMFA steric and electrostatic field from partial least squares (PLS) analysis, contour maps are plotted as percentage contribution to the QSAR equation and are associated with the differences in biological activity.

Synthesis and anti-tumor activities of some new pyridines and pyrazolo[1,5-a]pyrimidines

Cyclocondensation of cyanoacetamide, cyanothioacetamide and 3-aminopyrazols with sodium salt of 3-hydroxy-1-(2-naphthyl)prop-2-en-1-one gives pyridin-2-one, pyridin-2(1H)-thione, and pyrazolo[1,5-a]pyrimidine derivatives. These derivatives showed potent anti-tumor cytotoxic activity in vitro using different human cancer cell lines.

Critical evaluation of methods to incorporate entropy loss upon binding in high-throughput docking

Proper accounting of the positional/orientational/conformational entropy loss associated with protein-ligand binding is important to obtain reliable predictions of binding affinity. Herein, we critically examine two simplified statistical mechanics-based approaches, namely a constant penalty per rotor method, and a more rigorous method, referred to here as the partition function-based scoring (PFS) method, to account for such entropy losses in high-throughput docking calculations. Our results on the estrogen receptor beta and dihydrofolate reductase proteins demonstrate that, while the constant penalty method over-penalizes molecules for their conformational flexibility, the PFS method behaves in a more "DeltaG-like" manner by penalizing different rotors differently depending on their residual entropy in the bound state. Furthermore, in contrast to no entropic penalty or the constant penalty approximation, the PFS method does not exhibit any bias towards either rigid or flexible molecules in the hit list. Preliminary enrichment studies using a lead-like random molecular database suggest that an accurate representation of the "true" energy landscape of the protein-ligand complex is critical for reliable predictions of relative binding affinities by the PFS method.

Candidate target mechanisms of the growth inhibitor cyromazine: studies of phenylalanine hydroxylase, puparial amino acids, and dihydrofolate reductase in dipteran insects

Cyromazine, an insect growth regulator, affects larval and pupal cuticles in dipterans and some other insects. The mode of action of this aminotriazine is not known yet, though it has been shown not to inhibit the synthesis of chitin and cuticular proteins. Cyromazine may, however, act on some step(s) of sclerotization of the cuticle. In the present study, we have analyzed the key enzyme for the production of sclerotization agents, phenylalanine hydroxylase (PAH), using the enzyme from Drosophila, a cyromazine-sensitive insect. PAH was studied in vitro with cyromazine and three biologically less active derivatives at concentrations ranging from 1 microM to 1 mM. None of the compounds did significantly affect PAH activity. Nor did cyromazine, fed to last instar larvae of Musca domestica, change the relative content of phenylalanine and tyrosine, or the total amount and profile of amino acids of puparial cuticles, which showed a larviform shape typical for cyromazine intoxication. Taken together, this study does not support the hypothesis that phenylalanine hydroxylase represents a target site of cyromazine. In additional studies, the conflicting results, as reported by others, on in vitro inhibition of dihydrofolate reductase (DHFR) by cyromazine were re-examined using the enzymes from larvae of the blowfly, Protophormia terraenovae, and from hen liver. There was no significant inhibition of either DHFR at 100 microM by cyromazine as well as by dicylanil, a pyrimidine analog that is biologically more active than cyromazine. In conclusion, the mode of action of cyromazine remains completely open.

Textilinins from Pseudonaja textilis textilis. Characterization of two plasmin inhibitors that reduce bleeding in an animal model

The incidence of vein-graft occlusion associated with myocardial infarction and thrombosis following the use of the plasmin inhibitor, aprotinin, to reduce blood loss during vascular surgery has prompted the isolation of an alternative kinetically distinct inhibitor of plasmin from the venom of Pseudonaja textilis. This inhibitor has been called textilinin (Txln) and two distinct forms have been isolated from the Brown-snake venom (molecular weight, 6688 and 6692). A comparison of plasmin inhibitor constants for aprotinin and the Txlns 1 and 2 indicated that the former bound very tightly (inhibitor constant, Ki approximately 10(-11) mol/l), while both of the latter bound less tightly (Ki approximately 10(-9) mol/l). Homogeneity of Txlns 1 and 2 was confirmed by sodium dodecyl sulphate-polyacrylamide gel electrophoresis and mass spectrometry. A sequence difference of six amino acids was observed between the two forms of Txln. Txln 1 and 2 showed, respectively, 45 and 43% homology with aprotinin, while there was 58 and 55% homology, respectively, with a plasmin inhibitor from the venom of eastern Taipan, Oxyuranus scutellatus. Both Txlns have six cysteines, like other inhibitors of this group, and homology was determined by alignment of these cysteines. Both have been shown to reduce blood loss by about 60% in a murine tail vein bleeding model. It is proposed that the kinetic profiles of Txln 1 and 2 for plasmin allow the arrest of haemorrhage without the possible threat of thrombosis.

Expression and characterization of the Trypanosoma cruzi dihydrofolate reductase domain

We have cloned and expressed in Escherichia coli a 702-base pair gene coding for the dihydrofolate reductase (DHFR) domain of the bifunctional dihydrofolate reductase-thymidylate synthase (DHFR-TS) from Trypanosoma cruzi. The DHFR domain was purified to homogeneity by methotrexate-Sepharose chromatography followed by an anion-exchange chromatography step in a mono Q column, and displayed a single 27-kDa band on SDS-PAGE. Gel filtration showed that the catalytic domain was expressed as a monomer. Kinetic parameters were similar to those reported for the wild-type bifunctional enzyme with Km values of 0.75 microM for dihydrofolate and 16 microM for NADPH and a kcat value of 16.5 s-1. T. cruzi DHFR is poorly inhibited by trimethoprim and pyrimethamine and the inhibition constants were always lower for the bifunctional enzyme. The binding of methotrexate was characteristic of a class of inhibitors that form an initial complex which isomerizes slowly to a tighter complex and are referred to as 'slow, tight-binding' inhibitors. While the slow-binding step of inhibition was apparently unaffected in the individually expressed DHFR domain, the overall inhibition constant was two-fold higher as a consequence of the superior inhibition constant value obtained for the initial inhibitory complex.

Collapse of parallel folding channels in dihydrofolate reductase from Escherichia coli by site-directed mutagenesis

The rate-limiting steps in the folding of dihydrofolate reductase from Escherichia coli have been shown to involve the conversion of a set of four intermediates to a corresponding set of native conformers via four parallel channels [Jennings et al. (1993) Biochemistry 32, 3783-3789]. Fluorescence and absorbance studies of the unfolding and refolding of the C85S/C152E double mutant at various final urea concentrations reveal two slow folding reactions, two fewer than observed in the wild-type protein. Refolding in the presence of substoichiometric levels of the inhibitor methotrexate shows that the two remaining slow reactions correspond to two parallel channels which lead to a pair of native conformers capable of binding the inhibitor. A combination of stopped-flow circular dichroism and cofactor binding studies confirms that the four parallel channels observed in the wild-type protein have collapsed into two channels in the mutant. Kinetic and equilibrium studies of the single cysteine mutants suggest that replacements of Cysteine-85 which perturb the hydrophobic core containing this side chain are responsible for the simplification of the kinetic mechanism. These results demonstrate that at least two of the parallel folding channels in dihydrofolate reductase arise when tertiary structure develops and are not dependent upon cis/trans isomerization at prolyl peptide bonds.

Computer-aided drug design: a free energy perturbation study on the binding of methyl-substituted pterins and N5-deazapterins to dihydrofolate reductase

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of 8-methylpterins and 8-methyl-N5-deazapterins to dihydrofolate reductase (DHFR). Methyl-substitution at the 5, 6 and 7 positions in the N-heterocyclic ring gives rise to a variety of ring substituent patterns and biological activity: several of these methyl derivatives of the 8-methyl parent compounds (8-methylpterin and 8-methyl-N5-deazapterin) have been identified as substrates or inhibitors of vertebrate DHFR in previous work. The calculated free energy differences reveal that the methyl-substituted compounds are thermodynamically more stable than the primary compounds (8-methylpterin and 8-methyl-N5-deazapterin) when bound to the enzyme, due largely to hydrophobic hydration phenomena. Methyl substitution at the 5 and/or 7 positions in the 6-methyl-substituted compounds has only a small effect on the stability of ligand binding. Furthermore, repulsive interactions between the 6-methyl substituent and DHFR are minimal, suggesting that the 6-methyl position is optimal for binding. The results also show that similarly substituted 8-methylpterins and 8-methyl-N5-deazapterins have very similar affinities for binding to DHFR. The computer simulation predictions are in broad agreement with experimental data obtained from kinetic studies, i.e. 6,8-dimethylpterin is a more efficient substrate than 8-methylpterin and 6,8-dimethyl-N5-deazapterin is a better inhibitor than 8-methyl-N5-deazapterin.

Evidence for two interconverting protein isomers in the methotrexate complex of dihydrofolate reductase from Escherichia coli

Two-dimensional 1H NMR methods and a knowledge of the X-ray crystal structure have been used to make resonance assignments for the amino acid side chains of dihydrofolate reductase from Escherichia coli complexed with methotrexate. The H7 proton on the pteridine ring of methotrexate was found to have NOEs to the methyl protons of Leu-28 which were assigned by using the L28F mutant. These NOEs indicated that the orientation of the methotrexate pteridine ring is similar in both solution and crystal structures. During the initial assignment process, it became evident that many of the resonances in this complex, unlike those of the folate complex, are severely broadened or doubled. The observation of two distinct sets of resonances in a ratio of approximately 2:1 was attributed to the presence of two protein isomers. At 303 K, NOESY spectra with mixing times of 100 ms did not show interconversion between these isomers. However, exchange cross-peaks were observed in a 700-ms NOESY spectrum at 323 K which demonstrated that these isomers are interconverting slowly on the NMR time scale. Many of the side chains with clearly doubled resonances were located in the beta-sheet and the active site. Preliminary studies on the apoprotein also revealed doubled resonances in the absence of the inhibitor, indicating the existence of the protein isomers prior to methotrexate binding. In contrast to the methotrexate complex, the binary complex with folate and the ternary MTX-NADPH-DHFR complex presented a single enzyme form. These results are proposed to reflect the ability of folate and NADPH to bind predominantly to one protein isomer.

Crystal structures of recombinant human dihydrofolate reductase complexed with folate and 5-deazafolate

The 2.3-A crystal structure of recombinant human dihydrofolate reductase (EC 1.5.1.3, DHFR) has been solved as a binary complex with folate (a poor substrate at neutral pH) and also as a binary complex with an inhibitor, 5-deazafolate. The inhibitor appears to be protonated at N8 on binding, whereas folate is not. Rotation of the peptide plane joining I7 and V8 from its position in the folate complex permits hydrogen bonding of 5-deazafolate's protonated N8 to the backbone carbonyl of I7, thus contributing to the enzyme's greater affinity for 5-deazafolate than for folate. In this respect it is likely that bound 5-deazafolate furnishes a model for 7,8-dihydrofolate binding and, in addition, resembles the transition state for folate reduction. A hypothetical transition-state model for folate reduction, generated by superposition of the DHFR binary complexes human.5-deazafolate and chicken liver.NADPH, reveals a 1-A overlap of the binding sites for folate's pteridine ring and the dihydronicotinamide ring of NADPH. It is proposed that this binding-site overlap accelerates the reduction of both folate and 7,8-dihydrofolate by simultaneously binding substrate and cofactor with a sub van der Waals separation that is optimal for hydride transfer.

Monitoring the effect of triamterene and hydrochlorothiazide on dihydrofolate reductase activity using a new spectrophotometric method

A new spectrophotometric method is developed and applied for the study of the inhibitory effect of triamterene, hydrochlorothiazide and their combinations on the in vitro activity of dihydrofolate reductase enzyme. The method is based on incubating the drug (0.1-1.0 microM) or a buffer control with a solution containing reduced nicotinamide adenine dinucleotide phosphate (0.5 mM), magnesium chloride (1.29 mM), and folic acid as a substrate (0.01-0.1 mM) with the dihydrofolate reductase (0.25 unit). The resulting tetrahydrofolic acid is determined by first hydrolysing it by a methanol-hydrochloric acid mixture to produce p-aminobenzoyl glutamic acid, then adding p-dimethylaminocinnamic aldehyde reagent to form a stable pink coloured product. The colour is found to develop within 5 min and is stable over 12 h, with a maximum absorption at 545 nm. A linear calibration curve is formed by using standard solutions of tetrahydrofolic acid. The presence of the studied drugs did not interfere with the determination. Lineweaver-Burk plots of the reaction kinetics, in the presence of triamterene and/or hydrochlorothiazide showed a competitive inhibition of the dihydrofolate reductase in the presence of triamterene with or without hydrochlorothiazide. A 100% inhibition is obtained by 1 microM solution of triamterene at a folic acid concentration of 0.01 mM. No measurable effect of hydrochlorothiazide at the studied concentration range is demonstrated.

Nuclear magnetic resonance study of the state of protonation of inhibitors bound to mutant dihydrofolate reductase lacking the active-site carboxyl

13C nuclear magnetic resonance spectra have been obtained for complexes of [2-13C]methotrexate and [2-13C]trimethoprim with wild-type dihydrofolate reductase (DHFR) from Escherichia coli and with two mutant enzymes in which aspartic acid-27 is replaced by asparagine and by serine, respectively. In both the wild-type and mutated enzymes, exchange between the free inhibitor and the enzyme-complexed inhibitor is slow on the NMR time scale; hence, despite the considerably increased dissociation constants for binary complexes with the enzymes, the dissociation rate remains small relative to the frequency separation of the resonances. In all cases but one, the pKa of an inhibitor that is complexed to enzyme differs greatly from that of the free inhibitor. However, while the pKa of both inhibitors in complexes with the wild-type enzyme is elevated to above 10, the pKa of the inhibitors complexed with the Asn-27 and Ser-27 enzymes is lowered to a value below 4. Exact determinations of bound pKa values are limited by the solubility of the enzyme and the dissociation constants of the complexes. The single exception to these general conclusions is the ternary complex of the Ser-27 DHFR with trimethoprim and NADPH. In this complex, both free and enzyme-complexed trimethoprim exhibit similar pKa values (approximately equal to 7.6). However, both the exchange between free and enzyme-complexed inhibitor and the protonation of the enzyme-complexed inhibitor are slow in the NMR time scale, so that the spectra reveal three resonances corresponding to free inhibitor, to protonated enzyme-complexed inhibitor, and to unprotonated enzyme-complexed inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)

Role of isomerization of initial complexes in the binding of inhibitors to dihydrofolate reductase

Stopped-flow measurements of protein fluorescence quenching when methotrexate (MTX) binds to dihydrofolate reductase (isoenzyme II) of Streptococcus faecium (SFDHFR II) analyze as the sum of two differentials: a rapid binding phase and a second phase for which the observed rate constant is independent of methotrexate concentration. Analysis of variation of the ratio of the amplitude of the fast and slow phases with methotrexate concentration indicates that the second phase is an isomerization of the initial binary complex. At pH 7.3, the equilibrium constant for this isomerization is 21.9, and the forward and reverse rate constants are 0.57 and 0.026 s-1, respectively. Similar results were obtained for binding of 3-deazamethotrexate to SFDHFR II, but the forward rate constant is greater (2.9 s-1 at pH 7.3). The equilibrium constants for these isomerizations are pH independent, but the rate constants decrease as the pH is raised, probably due to deprotonation of one or more groups on the enzyme. Analysis of progress curves obtained by the development of inhibition when SFDHFR II is added last to reaction mixtures containing dihydrofolate, NADPH, and MTX gives an association constant for initial reactions of 4.3 X 10(7) M-1. Since a preliminary estimate of the association constant for the binding reaction is 7.6 X 10(5) M-1, this suggests an isomerization of the ternary complex(es) with an equilibrium constant of about 56. In addition, analysis of the progress of development of inhibition indicates a further very slow isomerization with equilibrium constant 419 and forward rate constant 2.6 min-1.(ABSTRACT TRUNCATED AT 250 WORDS)