Protonated state of methotrexate, trimethoprim, and pyrimethamine bound to dihydrofolate reductase

13C nuclear magnetic resonance (NMR) of methotrexate, trimethoprim, and pyrimethamine enriched 90% with 13C at C2 has provided a sensitive means of detecting the state of protonation of the heterocyclic rings of these inhibitors. In each case, protonation of N1 causes an upfield movement of the chemical shift of C2 by more than 6 ppm. By this method it has been shown that, at pH values up to 9.2, methotrexate is bound to bovine liver dihydrofolate reductase with N1 of the inhibitor protonated, just as in the case of the complex with reductase from Streptococcus faecium and Lactobacillus casei. Furthermore, trimethoprim bound to reductase from any of the three sources, and pyrimethamine bound to either of the bacterial reductases also have N1 protonated even at pH values up to 10. This implies that in all cases there is a strong interaction between protonated N1 of the inhibitor and the carboxylate group of the active site aspartate or glutamate. In every case pKa of the bound inhibitor is increased by several units, a finding in accord with crystallographic evidence that inhibitor bound to L. casei reductase is in a hydrophobic environment and that N1 is not hydrogen-bonded to water. It was confirmed by titration of protein fluorescence that trimethoprim has greater affinity for bacterial reductase than for vertebrate (bovine) reductase, and that this selectivity is more marked in ternary complexes in which NADPH is also bound to the active site. However, the data cited above indicate that this difference in affinities is not due to a weaker ionic interaction between protonated N1 of trimethoprim and the bovine enzyme. Instead, binding of the trimethoprim side chain to hydrophobic sites on the enzyme must provide less binding energy in the case of the mammalian enzyme.     

Folding of dihydrofolate reductase from Escherichia coli

The urea-induced equilibrium unfolding transition of dihydrofolate reductase from Escherichia coli was monitored by UV difference, circular dichroism (CD), and fluorescence spectroscopy. Each of these data sets were well described by a two-state unfolding model involving only native and unfolded forms. The free energy of folding in the absence of urea at pH 7.8, 15 degrees C is 6.13 +/- 0.36 kcal mol-1 by difference UV, 5.32 +/- 0.67 kcal mol-1 by CD, and 5.42 +/- 1.04 kcal mol-1 by fluorescence spectroscopy. The midpoints for the difference UV, CD, and fluorescence transitions are 3.12, 3.08, and 3.18 M urea, respectively. The near-coincidence of the unfolding transitions monitored by these three techniques also supports the assignment of a two-state model for the equilibrium results. Kinetic studies of the unfolding and refolding reactions show that the process is complex and therefore that additional species must be present. Unfolding jumps in the absence of potassium chloride revealed two slow phases which account for all of the amplitude predicted by equilibrium experiments. Unfolding in the presence of 400 mM KCl results in the selective loss of the slower phase, implying that there are two native forms present in equilibrium prior to unfolding. Five reactions were observed in refolding: two slow phases designated tau 1 and tau 2 that correspond to the slow phases in unfolding and three faster reactions designated tau 3, tau 4, and tau 5 that were followed by stopped-flow techniques. The kinetics of the recovery of the native form was monitored by following the binding of methotrexate, a tight-binding inhibitor of dihydrofolate reductase, at 380 nm.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate

We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.     

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.     

The crystal structure of Escherichia coli ketopantoate reductase with NADP+ bound

The NADPH-dependent reduction of ketopantoate to pantoate, catalyzed by ketopantoate reductase (KPR; EC 1.1.1.169), is essential for the biosynthesis of pantothenate (vitamin B(5)). Here we present the crystal structure of Escherichia coli KPR with NADP(+) bound, solved to 2.1 A resolution. The cofactor is bound in the active site cleft between the N-terminal Rossmann-fold domain and the C-terminal alpha-helical domain. The thermodynamics of cofactor and substrate binding were characterized by isothermal titration calorimetry. The dissociation constant for NADP(+) was found to be 6.5 muM, 20-fold larger than that for NADPH (0.34 muM). The difference is primarily due to the entropic term, suggesting favorable hydrophobic interactions of the more lipophilic nicotinamide ring in NADPH. Comparison of this binary complex structure with the previously studied apoenzyme reveals no evidence for large domain movements on cofactor binding. This observation is further supported both by molecular dynamics and by calorimetric analysis. A model of the ternary complex, based on the structure presented here, provides novel insights into the molecular mechanism of enzyme catalysis. We propose a conformational switch of the essential Lys176 from the "resting" state observed in our structure to an "active" state, to bind ketopantoate. Additionally, we identify the importance of Asn98 for substrate binding and enzyme catalysis.     

Site-specific mutagenesis of dihydrofolate reductase from Escherichia coli

Two site-specific mutations of dihydrofolate reductase from Escherichia coli based on the x-ray crystallographic structure were constructed. The first mutation (His-45----Gln) is aimed at assessing the interaction between the imidazole moiety and the pyrophosphate backbone of NADPH. The second (Thr-113----Val) is part of a hydrogen bonding network that contacts the dihydrofolate substrate and may be involved in proton delivery to the N5-C6 imine undergoing reduction. The first mutation was shown to alter both the association and dissociation rate constants for the cofactor so that the dissociation constant was increased 6-40-fold. A corresponding but smaller (fourfold) effect was noted in V/K but not in V compared to the wild-type enzyme. The second was demonstrated to increase the dissociation rate constant for methotrexate 20-30-fold, and presumably dihydrofolate also, with a corresponding 20-30-fold increase in the dissociation constant. In this case an identical effect was noted on V/K but not in V relative to the native enzyme. Thus, in both mutant enzymes the decrease in binding has not been translated into a loss of catalytic efficiency.     

Characterization of the cloned Escherichia coli dihydrofolate reductase

Dihydrofolate reductase (5,6,7,8-tetrahydrofolate: NADP+ oxidoreductase, EC 1.5.1.3) was purified from Escherichia coli strains that carried derivatives of the multicopy recombinant plasmid, pJFM8. The results of enzyme kinetic and two-dimensional gel electrophoresis experiments showed that the cloned enzyme is indistinguishable from the chromosomal enzyme. Therefore it can be concluded that these strains are ideal for use as a source of enzyme for further studies on the biochemistry and regulation of this important enzyme. The plasmid derivatives were constructed by recloning experiments that utilized several restriction endonucleases. From the analysis both of these plasmids and the purified dihydrofolate reductase enzymes it was possible to deduce the location and orientation of the dihydrofolate reductase structural gene on the parent plasmid, pJFM8.     

Purification and properties of Escherichia coli dihydrofolate reductase.

Dihydrofolate reductase has been purified 40-fold to apparent homogeneity from a trimethoprim-resistant strain of Escherichia coli (RT 500) using a procedure that includes methotrexate affinity column chromatography. Determinations of the molecular weight of the enzyme based on its amino acid composition, sedimentation velocity, and sodium dodecyl sulfate gel electrophoresis gave values of 17680, 17470 and 18300, respectively. An aggregated form of the enzyme with a low specific activity can be separated from the monomer by gel filtration; treatment of the aggregate with mercaptoethanol or dithiothreitol results in an increase in enzymic activity and a regeneration of the monomer. Also, multiple molecular forms of the monomer have been detected by polyacrylamide gel electrophoresis. The unresolved enzyme exhibits two pH optima (pH 4.5 and pH 7.0) with dihydrofolate as a substrate. Highest activities are observed in buffers containing large organic cations. In 100 mM imidazolium chloride (pH 7), the specific activity is 47 mumol of dihydrofolate reduced per min per mg at 30 degrees. Folic acid also serves as a substrate with a single pH optimum of pH 4.5. At this pH the Km for folate is 16 muM, and the Vmax is 1/1000 of the rate observed with dihydrofolate as the substrate. Monovalent cations (Na+, K+, Rb+, and Cs+) inhibit dihydrofolate reductase; at a given ionic strength the degree of inhibition is a function of the ionic radius of the cation. Divalent cations are more potent inhibitors; the I50 of BaCl2 is 250 muM, as compared to 125 mM for KCl. Anions neither inhibit nor activate the enzyme.     

Design of dihydrofolate reductase inhibitors from X-ray crystal structures

Dihydrofolate reductase (DHFR) is an important therapeutic target for treatment of cancer and microbial disease. Its species specificity has resulted in the sequencing of a number of vertebrate and bacterial DHFRs, and the three-dimensional structure of isozymes from Escherichia coli, Lactobacillus casei, and chicken liver has been elucidated, in the presence of the coenzyme NADPH and of a number of inhibitors. This information has enabled scientists to try to design improved and more selective inhibitors, based on the known coordinates of the enzyme features. Simple use of computer graphics or wire models has resulted in the design of inhibitors with 50 times the activity of trimethoprim, an antibacterial DHFR inhibitor, by making use of an unused ionic binding site. However, in a number of instances this approach was completely unsuccessful because hydrophobic sites of interaction were preferred. More sophisticated techniques involve energy minimization of the small molecule-macromolecule interactions to optimize the geometry. In this paper I describe the use of a molecular mechanics program, AMBER, for predicting the geometry and relative energetics of binding. Very encouraging results have been obtained for a closely related series of compounds. Where differing entropic and solvent effects are involved, predictions may be poor. The use of super computers and molecular dynamics methods should increase this capability in the near future.     

Purification of dihydrofolate reductase amplified in Escherichia coli K-12

The Escherichia coli strain carrying pTP 6-10 which was constructed in our previous work (Iwakura, M., et al. (1983) J. Biochem. 93, 927-930) produces more than 400-fold dihydrofolate reductase as compared with the strain without the plasmid. Dihydrofolate reductase was highly purified from the cell-free extract of the plasmid strain simply by two steps; ammonium sulfate fractionation and ion-exchange chromatography. By 10-fold purification, the enzyme was essentially homogeneous as judged by polyacrylamide gel electrophoresis in the presence of sodium dodecyl sulfate. The restriction map of pTP 6-10 was also determined and the plasmid was shown to have an Ava I, an EcoR I, a Pst I, a Pvu I, and a Pvu II site. Our results indicate that the plasmid strain is suitable as a source of the enzyme and that plasmid pTP 6-10 is promising as a versatile plasmid vector for efficiently yielding the product of the cloned gene.     

Kinetic analysis of the mechanism of Escherichia coli dihydrofolate reductase

A kinetic mechanism is presented for Escherichia coli dihydrofolate reductase which describes the full time course of the enzymatic reaction over a wide range of substrate and enzyme concentrations at pH 7.2 and 20 degrees C. Specific rate constants were estimated by computer simulation of the full time course of single turnover, burst, and steady-state experiments using both nondeuterated and deuterated NADPH. The mechanism involves the random addition of substrates, but the substrates and enzyme are not at equilibrium prior to the chemical transformation step. The rate-limiting step follows the chemical transformation, and the maximum velocity of the reaction is limited by the release of the product tetrahydrofolate. The full time course of the reaction is markedly affected by the formation of the enzyme-NADPH-tetrahydrofolate abortive complex, but not by the enzyme-NADP-dihydrofolate abortive complex.     

Substrate-induced hysteresis in the activity of Escherichia coli dihydrofolate reductase

Full time course studies of the kinetic activity of Escherichia coli dihydrofolate reductase show that there is an increase in activity with time. The half-time for this hysteretic behavior is about 9 s. Preincubation of the enzyme with either of the substrates abolishes the lag and results in initial velocities which are 2-2.3-fold faster than those observed for the non-preincubated enzyme. The kinetic properties of the activated and nonactivated forms of the enzyme appear to be similar as measured by the full time course of the reaction. The results are consistent with observations for NADPH binding studies that the enzyme exists in two interconvertible forms, one of which is incapable of binding NADPH (Cayley, P. J., Dunn, S. M. J., and King, R. W. (1981) Biochemistry 20, 874-879).     

Lanthanide ions inhibit the activity of dihydrofolate reductase from chicken liver

The influences of mono-, bi- and trivalent metal ions (as chloride salts) on the activity of dihydrofolate reductase (DHFR) from chicken liver have been studied to elucidate the mechanism of ion-activation of this enzyme. The results show that monovalent ions (Na⁺ and K⁺) activate DHFR at low concentration reaching a maximum activation of about 2.5 fold at 0.4–0.5 M and declining at higher concentrations. Ca²⁺ shows similar activation but at lower concentration, reaching a maximum at 0.1 M; activity declines with further increases in concentration. At very high concentration (>0.4 M), Ca²⁺ is inhibitory. The trivalent lanthanide ions, however, show a dramatic inhibition of activity of DHFR even at very low concentration. The activity of DHFR declines to 50% of that of the control at 0.02 mM EuCl₃. Intrinsic fluorescence measurements show that the ion-dependent activation in the presence of mono- and bivalent metal ions is due to the conformational changes in the protein. Energy transfer phenomenon suggests that the specific interaction of Eu³⁺ with Trp24 located in a loop at the active site of DHFR is responsible for the strong inhibition. The possible mechanism for the ion-inhibition is proposed and discussed.     

The pH-dependence of the binding of dihydrofolate and substrate analogues to dihydrofolate reductase from Escherichia coli

The interaction of dihydrofolate reductase (EC 1.5.1.3) from Escherichia coli with dihydrofolate and folate analogues has been studied by means of binding and spectroscopic experiments. The aim of the investigation was to determine the number and identity of the binary complexes that can form, as well as pKa values for groups on the ligand and enzyme that are involved with complex formation. The results obtained by ultraviolet difference spectroscopy indicate that, when bound to the enzyme, methotrexate and 2,4-diamino-6,7-dimethylpteridine exist in their protonated forms and exhibit pKa values for their N-1 nitrogens of above 10.0. These values are about five pH units higher than those for the compounds in free solution. The binding data suggest that both folate analogues interact with the enzyme to yield a protonated complex which may be formed by reaction of ionized enzyme with protonated ligand and/or protonated enzyme with unprotonated ligand. The protonated complex formed with 2,4-diamino-6,7-dimethylpteridine can undergo further protonation to form a protonated enzyme-protonated ligand complex, while that formed with methotrexate can ionize to give an unprotonated complex. A group on the enzyme with a pKa value of about 6.3 is involved with the interactions. However, the ionization state of this group has little effect on the binding of dihydrofolate to the enzyme. For the formation of an enzyme-dihydrofolate complex it is essential that the N-3/C-4 amide of the pteridine ring of the substrate be in its neutral form. It appears that dihydrofolate is not protonated in the binary complex.     

Effects of pressure on enzyme function of Escherichia coli dihydrofolate reductase

To elucidate the effects of pressure on the function of Escherichia coli dihydrofolate reductase (DHFR), the enzyme activity and the dissociation constants of substrates and cofactors were measured at pressures up to 250 MPa at 25 degrees C and pH 7.0. The enzyme activity decreased with increasing pressure, accompanying the activation volume of 7.8 ml mol(-1). The values of the Michaelis constant (K(m)) for dihydrofolate and NADPH were slightly higher at 200 MPa than at atmospheric pressure. The hydride-transfer step was insensitive to pressure, as monitored by the effects of the deuterium isotope of NADPH on the reaction velocity. The dissociation constants of substrates and cofactors increased with pressure, producing volume reductions from 6.5 ml mol(-1) (tetrahydrofolate) to 33.5 ml mol(-1) (NADPH). However, the changes in Gibbs free energy with dissociation of many ligands showed different pressure dependences below and above 50 MPa, suggesting conformational changes of the enzyme at high pressure. The enzyme function at high pressure is discussed based on the volume levels of the intermediates and the candidates for the rate-limiting process.     

Escherichia coli dihydrofolate reductase: isolation and characterization of two isozymes.

A combination of affinity column chromatography and preparative gel electrophoresis has been used to purify to homogeneity the two isozymes of dihydrofolate reductase from a trimethoprim-resistant strain of Escherichia coli B (RT 500). These enzyme forms are noninterconvertible and are present in crude cell lysates, but other electrophoretic species can be generated durng purification if sulfhydryl-protecting agents, such as dithiothreitol, are not present. The two isozymes, numbered form 1 and form 2 with respect to their decreasing electrophoretic mobilities, have similar molecular weights (18 500), molecular radii (21 A), and apparent Km values for reduced nico inamide adenin- dinucleotide (NADH) and NADH phosphate (NADPH). Both forms contain 2 mol of sulfhydryl/mol of enzyme which can be oxidized to intramolecular disulfide bonds. However, forms 1 and 2 differ physically in their electrophoretic mobility and isoelectric point and kinetically in their pH-activity profile, specific activity, Km for dihydrofolate, and their affinity toward a number of inhibitors.     

FolM, a new chromosomally encoded dihydrofolate reductase in Escherichia coli

Escherichia coli (thyA DeltafolA) mutants are viable and can grow in minimal medium when supplemented with thymidine alone. Here we present evidence from in vivo and in vitro studies that the ydgB gene determines an alternative dihydrofolate reductase that is related to the trypanosomatid pteridine reductases. We propose to rename this gene folM.