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)     

A kinetic study of wild-type and mutant dihydrofolate reductases from Lactobacillus casei

A kinetic scheme is presented for Lactobacillus casei dihydrofolate reductase that predicts steady-state kinetic parameters. This scheme was derived from measuring association and dissociation rate constants and pre-steady-state transients by using stopped-flow fluorescence and absorbance spectroscopy. Two major features of this kinetic scheme are the following: (i) product dissociation is the rate-limiting step for steady-state turnover at low pH and follows a specific, preferred pathway in which tetrahydrofolate (H4F) dissociation occurs after NADPH replaces NADP+ in the ternary complex; (ii) the rate constant for hydride transfer from NADPH to dihydrofolate (H2F) is rapid (khyd = 430 s-1), favorable (Keq = 290), and pH dependent (pKa = 6.0), reflecting ionization of a single group. Not only is this scheme identical in form with the Escherichia coli kinetic scheme [Fierke et al. (1987) Biochemistry 26, 4085] but moreover none of the rate constants vary by more than 40-fold despite there being less than 30% amino acid homology between the two enzymes. This similarity is consistent with their overall structural congruence. The role of Trp-21 of L. casei dihydrofolate reductase in binding and catalysis was probed by amino acid substitution. Trp-21, a strictly conserved residue near both the folate and coenzyme binding sites, was replaced by leucine. Two major effects of this substitution are on (i) the rate constant for hydride transfer which decreases 100-fold, becoming the rate-limiting step in steady-state turnover, and (ii) the affinities for NADPH and NADP+ which decrease by approximately 3.5 and approximately 0.5 kcal mol-1, respectively.(ABSTRACT TRUNCATED AT 250 WORDS)     

Construction and evaluation of the kinetic scheme associated with dihydrofolate reductase from Escherichia coli

A kinetic scheme is presented for Escherichia coli dihydrofolate reductase that predicts steady-state kinetic parameters and full time course kinetics under a variety of substrate concentrations and pHs. This scheme was derived from measuring association and dissociation rate constants and pre-steady-state transients by using stopped-flow fluorescence and absorbance spectroscopy. The binding kinetics suggest that during steady-state turnover product dissociation follows a specific, preferred pathway in which tetrahydrofolate (H4F) dissociation occurs after NADPH replaces NADP+ in the ternary complex. This step, H4F dissociation from the E X NADPH X H4F ternary complex, is proposed to be the rate-limiting step for steady-state turnover at low pH because koff = VM. The rate constant for hydride transfer from NADPH to dihydrofolate (H2F), measured by pre-steady-state transients, has a deuterium isotope effect of 3 and is rapid, khyd = 950 s-1, essentially irreversible, Keq = 1700, and pH dependent, pKa = 6.5, reflecting ionization of a single group in the active site. This scheme accounts for the apparent pKa = 8.4 observed in the steady state as due to a change in the rate-determining step from product release at low pH to hydride transfer above pH 8.4. This kinetic scheme is a necessary background to analyze the effects of single amino acid substitutions on individual rate constants.     

Role of water in the catalytic cycle of E. coli dihydrofolate reductase

Dihydrofolate reductase (DHFR) catalyzes the nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reduction of 7,8-dihydrofolate (H2F) to 5,6,7,8-tetrahydrofolate (H4F). Because of the absence of any ionizable group in the vicinity of N5 of dihydrofolate it has been proposed that N5 could be protonated directly by a water molecule at the active site in the ternary complex of the Escherichia coli enzyme with cofactor and substrate. However, in the X-ray structures representing the Michaelis complex of the E. coli enzyme, a water molecule has never been observed in a position that could allow protonation of N5. In fact, the side chain of Met 20 blocks access to N5. Energy minimization reported here revealed that water could be placed in hydrogen bonding distance of N5 with only minor conformational changes. The r.m.s. deviation between the conformation of the M20 loop observed in the crystal structures of the ternary complexes and the conformation adopted after energy minimization was only 0.79 A. We performed molecular dynamics simulations to determine the accessibility by water of the active site of the Michaelis complex of DHFR. Water could access N5 relatively freely after an equilibration time of approximately 300 psec during which the side chain of Met 20 blocked water access. Protonation of N5 did not increase the accessibility by water. Surprisingly the number of near-attack conformations, in which the distance between the pro-R hydrogen of NADPH and C6 of dihydrofolate was less than 3.5 A and the angle between C4 and the pro-R hydrogen of NADPH and C6 of dihydrofolate was greater than 120 degrees, did not increase after protonation. However, when the hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the side chain of Met 20 moved away from N5 after approximately 100 psec thereby providing water access. The average time during which water was found in hydrogen bonding distance to N5 was significantly increased. These results suggest that hydride transfer might occur early to midway through the reaction followed by protonation. Such a mechanism is supported by the very close contact between C4 of NADP+ and C6 of folate observed in the crystal structures of the ternary enzyme complexes, when the M20 loop is in its closed conformation.     

Human dihydrofolate reductase: reduction of alternative substrates, pH effects, and inhibition by deazafolates.

The kinetics of the NADPH-dependent reduction of 7,8-dihydrofolate, folate, and 7,8-dihydrobiopterin by human dihydrofolate reductase have been examined over the pH range from 4.0 to 9.5. The V and V/K profiles obtained with the three substrates indicate that a single ionizing residue at the active site of the enzyme must be protonated for catalysis. Both the maximum velocity of the reactions and the rate of interaction of the substrates with the enzyme-NADPH complex decrease in the order dihydrofolate greater than dihydrobiopterin much greater than folate. From the pK values of the V/K profiles, it can be concluded that, while dihydrofolate behaves as a sticky substrate and dihydrobiopterin exhibits slight stickiness, folate is not a sticky substrate. Further support for this conclusion comes from the results of deuterium isotope effects. The pK values obtained from both the V and V/Kfolate profiles are similar to the intrinsic pK value of 5.6 for both the free enzyme and the enzyme-NADPH complex. The folate analogue, 5-deazafolate, is not a substrate, but it undergoes strong interaction with the enzyme. This interaction, which is enhanced by the presence of NADPH, is due to protonation of the bound ligand that does not involve the single ionizing group at the active center of the enzyme. Difference spectra yield evidence for the protonation of bound 5-deazafolate and show that, on binding to the enzyme-NADPH complex, the pK of the N-8 atom is raised to about 10 from a value of about 4 in solution. The results are in accord with those of a recent paper on the three-dimensional structure of the enzyme-5-deazafolate complex [Davies, J.F., Delcamp, T.J., Prendergast, N.J., Ashfors, V.A., Freisheim, J.H., & Kraut, J. (1990) Biochemistry 29, 9467-9479] which indicate that there is hydrogen bond formation between N-8 of the ligand and the carbonyl group of Ile-7. However, the present findings do not support the idea that bound 5-deazafolate resembles the transition-state complex for folate reduction. Quinazolines also interact strongly with the enzyme but in a pH-independent manner. The dissociation constants for the binary complexes are an order of magnitude lower than that for the binding to the enzyme of unprotonated 5-deazafolate. This difference reflects the hydrophobic nature of the amino acid residues at the active site that are near the N-5 and N-8 nitrogens of bound pterins.     

Impact on catalysis of secondary structural manipulation of the alpha C-helix of Escherichia coli dihydrofolate reductase

The alpha C-helix of Escherichia coli dihydrofolate reductase has been converted to its counterpart in Lactobacillus casei by a triple mutation in the helix (H45R, W47Y, and I50F). These changes result in a 2-fold increase in the steady-state reaction rate (kcat = 26 s-1) that is limited by an increased off rate for the release of tetrahydrofolate (koff = 40 s-1 versus 12 s-1). On the other hand the mutant protein exhibits a 10-fold increase in the KM value (6.8 microM) for dihydrofolate and a 10-fold decrease in the rate of hydride transfer (85 s-1) from NADPH to dihydrofolate. The elevated rate of tetrahydrofolate release upon the rebinding of NADPH, a characteristic of the wild-type enzyme-catalyzed reaction, is diminished. The intrinsic pKa (6.4) of the mutant enzyme binary complex with NADPH is similar to that of the wild type, but the pKa of the ternary complex is increased to 7.3, about on pH unit higher than the wild-type value. Further mutagenesis (G51P and an insertion of K52) was conducted to incorporate a hairpin turn unique to the C-terminus of the alpha C-helix of the L. casei enzyme in order to adjust a possible dislocation of the new helix. The resultant pentamutant enzyme shows restoration of many of the kinetic parameters, such as kcat (12 s-1), KM (1.1 microM for dihydrofolate), and khyd (526 s-1), to the wild-type values. The synergism in the product release is also largely restored. A substrate-induced conformational change responsible for the fine tuning of the catalytic process was found to be associated with the newly installed hairpin structure. The Asp27 residue of the mutant enzyme was found to be reprotonated before tetrahydrofolate release.     

Role of the conserved active site residue tryptophan-24 of human dihydrofolate reductase as revealed by mutagenesis

The active sites of all bacterial and vertebrate dihydrofolate reductases that have been examined have a tryptophan residue near the binding sites for NADPH and dihydrofolate. In cases where the three-dimensional structure has been determined by X-ray crystallography, this conserved tryptophan residue makes hydrophobic and van der Waals interactions with the nicotinamide moiety of bound NADPH, and its indole nitrogen interacts with the C4 oxygen of bound folate through a bridge provided by a bound water molecule. We have addressed the question of why even the very conservative replacement of this tryptophan by phenylalanine does not seem to occur naturally. Human dihydrofolate reductase with this replacement of tryptophan by phenylalanine has increased rate constants for dissociation of substrates and products and a considerably decreased rate of hydride transfer. These cause some changes in steady-state kinetic behavior, including substantial increases in Michaelis constants for NADPH and dihydrofolate, but there is also a 3-fold increase in kcat. The branched mechanistic pathway for this enzyme has been completely defined and is sufficiently different from that of wild-type enzyme to cause changes in some transient-state kinetics. The most critical changes resulting from the amino acid substitution appear to be a 50% decrease in stability and a decrease in efficiency from 69% to 21% under intracellular conditions.     

Analysis of hydride transfer and cofactor fluorescence decay in mutants of dihydrofolate reductase: possible evidence for participation of enzyme molecular motions in catalysis.

A remarkable correlation has been discovered between fluorescence lifetimes of bound NADPH and rates of hydride transfer among mutants of dihydrofolate reductase (DHFR) from Escherichia coli. Rates of hydride transfer from NADPH to dihydrofolate change by a factor of 1,000 for the series of mutant enzymes. Since binding constants for the initial complex between coenzyme and DHFR change by only a factor of 10, the major portion of the change in hydride transfer must be attributed to losses in transition-state stabilization. The time course of fluorescence decay for NADPH bound to DHFR is biphasic. Lifetimes ranging from 0.3 to 0.5 ns are attributed to a solvent-exposed dihydronicotinamide conformation of bound coenzyme which is presumably not active in catalysis, while decay times (tau 2) in the range of 1.3 to 2.3 ns are assigned to a more tightly bound species of NADPH in which dihydronicotinamide is sequestered from solvent. It is this slower component that is of interest. Ternary complexes with three different inhibitors, methotrexate, 5-deazafolate, and trimethoprim, were investigated, along with the holoenzyme complex; 3-acetylNADPH was also investigated. Fluorescence polarization decay, excitation polarization spectra, the temperature variation of fluorescence lifetimes, fluorescence amplitudes, and wavelength of absorbance maxima were measured. We suggest that dynamic quenching or internal conversion promotes decay of the excited state in NADPH-DHFR. When rates of hydride transfer are plotted against the fluorescence lifetime (tau 2) of tightly bound NADPH, an unusual correlation is observed. The fluorescence lifetime becomes longer as the rate of catalysis decreases for most mutants studied. However, the fluorescence lifetime is unchanged for those mutations that principally alter the binding of dihydrofolate while leaving most dihydronicotinamide interactions relatively undisturbed. The data are interpreted in terms of possible dynamic motions of a flexible loop region in DHFR which closes over both substrate and coenzyme binding sites. These motions could lead to faster rates of fluorescence decay in holoenzyme complexes and, when correlated over time, may be involved in other motions which give rise to enhanced rates of catalysis in DHFR.     

Catalytic mechanism of the dihydrofolate reductase reaction as determined by pH studies

The variation with pH of the kinetic parameters of the reaction catalyzed by dihydrofolate reductase from Escherichia coli has been determined with the aim of elucidating the chemical mechanism of the reaction. The (V/K)DHF and V profiles indicated that protonation enhances the observed rate of interaction of dihydrofolate (DHF) with the enzyme-NADPH complex as well as the maximum velocity of the reaction. The pKa value of 8.09 observed in the (V/K)DHF profile is similar to that of 7.9 observed in the Ki profile for 2,4-diamino-6,7-dimethylpteridine while the pKa value of the V profile is displaced to 8.4. From the magnitude of the pH-independent value for (V/K)DHF, it is concluded that unprotonated dihydrofolate must react, at neutral pH, with the protonated form of the enzyme. The D(V/K)DHF value is independent of pH and equal to unity whereas the DV value varies as a wave function of pH with limiting values of 1.5 and 1.0 at low and high pH, respectively. It is proposed that dihydrofolate reacts with the unprotonated enzyme-NADPH complex to form a dead-end complex and with the protonated form of the same complex to form a productive complex. Further, it is considered that the protonated carboxyl of Asp-27 at the active site of the enzyme is responsible for the protonation of the N-5 nitrogen of dihydrofolate and that this protonation precedes and facilitates hydride transfer.     

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.     

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase.

R67 is a Type II dihydrofolate reductase (DHFR) that catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate by facilitating the addition of a proton to N5 of DHF and the transfer of a hydride ion from NADPH to C6. Because this enzyme is a plasmid-encoded DHFR from trimethoprim-resistant bacteria, extensive studies on R67 with various methods have been performed to elucidate its reaction mechanism. Here, Raman difference measurements, conducted on the ternary complex of R67.NADP(+).DHF believed to be an accurate mimic of the productive DHFR.NADPH.DHF complex, show that the pK(a) of N5 in the complex is less than 4. This is in clear contrast to the behavior observed in Escherichia coli DHFR, a substantially more efficient enzyme, where the pK(a) of bound DHF at N5 is increased to 6.5 compared with its solution value of 2.6. A comparison of the ternary complexes in R67 and E. coli DHFRs suggests that enzymic raising of the pK(a) at N5 can significantly increase the catalytic efficiency of the hydride transfer step. However, R67 shows that even without such a strategy an effective DHFR can still be designed.     

A review of the binding-change mechanism for proton-translocating transhydrogenase

Proton-translocating transhydrogenase is found in the inner membranes of animal mitochondria, and in the cytoplasmic membranes of many bacteria. It catalyses hydride transfer from NADH to NADP(+) coupled to inward proton translocation. Evidence is reviewed suggesting the enzyme operates by a "binding-change" mechanism. Experiments with Escherichia coli transhydrogenase indicate the enzyme is driven between "open" and "occluded" states by protonation and deprotonation reactions associated with proton translocation. In the open states NADP(+)/NADPH can rapidly associate with, or dissociate from, the enzyme, and hydride transfer is prevented. In the occluded states bound NADP(+)/NADPH cannot dissociate, and hydride transfer is allowed. Crystal structures of a complex of the nucleotide-binding components of Rhodospirillum rubrum transhydrogenase show how hydride transfer is enabled and disabled at appropriate steps in catalysis, and how release of NADP(+)/NADPH is restricted in the occluded state. Thermodynamic and kinetic studies indicate that the equilibrium constant for hydride transfer on the enzyme is elevated as a consequence of the tight binding of NADPH relative to NADP(+). The protonation site in the translocation pathway must face the outside if NADP(+) is bound, the inside if NADPH is bound. Chemical shift changes detected by NMR may show where alterations in protein conformation resulting from NADP(+) reduction are initiated. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).     

A catalytically active complex formed from the recombinant dI protein of Rhodospirillum rubrum transhydrogenase, and the recombinant dIII protein of the human enzyme

Transhydrogenase is a proton pump. It has three components: dI and dIII protrude from the membrane and contain the binding sites for NAD(H) and NADP(H), respectively, and dII spans the membrane. We have expressed dIII from Homo sapiens transhydrogenase (hsdIII) in Escherichia coli. The purified protein was associated with stoichiometric amounts of NADP(H) bound to the catalytic site. The NADP+ and NADPH were released only slowly from the protein, supporting the suggestion that nucleotide-binding by dIII is regulated by the membrane-spanning dII. HsdIII formed a catalytically active complex with recombinant dI from Rhodospirillum rubrum (rrdI), even in the absence of dII. The rates of forward and reverse transhydrogenation catalysed by this complex are probably limited by slow release from dIII of NADPH and NADP+, respectively. The hybrid complex also catalysed high rates of 'cyclic' transhydrogenation, indicating that hydride transfer, and exchange of nucleotides with dI, are rapid. Stopped-flow experiments revealed a rapid, monoexponential, single-turnover burst of reverse transhydrogenation in pre-steady-state. The apparent first-order rate constant of the burst increased with the concentration of rrdI. A deuterium isotope effect (kH/kD approximately 2 at 27 degrees C) was observed when [4B-1H]NADPH was replaced with [4B-2H]NADPH. The characteristics of the burst of transhydrogenation with rrdI:hsdIII differed from those previously reported for rrdI:rrdIII (J.D. Venning et al., Eur. J. Biochem. 257 (1998) 202-209), but the differences are readily explained by a greater dissociation constant of the hybrid complex. The steady-state rate of reverse transhydrogenation by the rrdI:hsdIII complex was almost independent of pH, but there was a single apparent pKa ( approximately 9.1) associated with the cyclic reaction. The reactions of the dI:dIII complex probably proceed independently of those protonation/deprotonation reactions which, in the complete enzyme, are associated with H+ translocation.     

13C]Methylated ribonuclease A. 13C NMR studies of the interaction of lysine 41 with active site ligands

A unique resonance in the 13C NMR spectrum of [13C]methylated ribonuclease A has been assigned to a N epsilon, N-dimethylated active site residue, lysine 41. The chemical shift of this resonance was studied over the pH range 3 to 11, and the titration curve showed two inflection points, at pH 5.7 and 9.0. The higher pKa, designated pKa1, was assigned to the ionization of the lysyl residue itself while the pKa of 5.7, designated pKa2, was assigned on the basis of its pKa to the ionization of a histidyl residue which is somehow coupled to lysine 41. Both pKa values are measurably perturbed by the binding of active site ligands including nucleotides, nucleosides, phosphate, and sulfate. In most cases, the alterations in pKa values induced by the ligands were larger for pKa2. The ligand-induced perturbations in pKa2 generally paralleled those reported for histidine 12, another active site residue (Griffin, J. H., Schechter, A. N., and Cohen, J. S. (1973) Ann. N. Y. Acad. Sci. 222, 693-708). The sensitivity of the N epsilon, N-dimethylated lysine 41 resonance to the histidyl ionization may result from a conformational change in the active site region of ribonuclease which is coupled to the histidyl ionization. This coupling between lysine 41 and another ribonuclease residue, which has not been documented previously, offers new insight into the interrelationship between residues in the active site of this well characterized enzyme.     

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.     

Theoretical studies on the activation of the pterin cofactor in the catalytic mechanism of dihydrofolate reductase

Two mechanisms for facilitating hydride ion transfer from NADPH involving preprotonation of the pteridine rings of the dihydrofolate reductase substrates folate and dihydrofolate have been investigated by ab initio quantum mechanical methods. Protonation energies and effective solution pKas have been calculated for four protonated forms, three of which are nonpreferred in aqueous solution and therefore not directly accessible to experimental study. The pattern and degree of redistribution of the positive charge over the component rings of the N-heterobicyclic pi-system in these protonated forms have been analyzed in terms of changes in the electron populations of the ring atoms and total ring charges. The effects of such changes in promoting hydride ion transfer to C7 in folate and C6 in dihydrofolate have been evaluated by considering the extent of development of partial carbonium ion character at these carbon atoms and also the degree of electron deficiency in the pyrazine ring as a whole. The results illustrate that perturbations due, for instance, to protonation may be propagated by pi-electron coupling effects over medium-range distances of 4-6 A across the pteridine ring. The two mechanisms have been assessed in terms of the calculated absolute and relative pKas of the protonated species taking into account experimental information regarding possible stabilization of these forms in the enzyme active site and also the effectiveness of the various protonations in assisting the hydride ion transfer step. Judged against these criteria, the theoretical results favor the generally proposed mechanism involving preprotonation of N8 in folate and N5 in dihydrofolate. However, some support was also found for the alternative novel mechanism involving O4-protonation of both folate and dihydrofolate.     

Effect of enzyme and ligand protonation on the binding of folates to recombinant human dihydrofolate reductase: implications for the evolution of eukaryotic enzyme efficiency.

There is marked pH dependence of the rate constant (koff) for tetrahydrofolate (H4folate) dissociation from its ternary complex with human dihydrofolate reductase (hDHFR) and NADPH. Similar pH dependence of H4folate dissociation from the ternary complex of a variant of hDHFR with the substitution Phe31----Leu (F31L hDHFR) causes this dissociation to become rate limiting in the enzyme mechanism at pH approximately 5, and this accounts for the marked decrease in kcat for this variant as the pH is decreased from 7 to 5. This decreased kcat at low pH is not seen for most DHFRs. koff for dissociation of folate, dihydrofolate (H2folate), and H4folate from their binary complexes with hDHFR is similarly pH dependent. For all the complexes examined, the pH dependence of koff in the range pH 5-7 is well described by a pKa of about 6.2 and must be due to ionization of a group on the enzyme. In the higher pH range (7-10), koff increases further as the pH is raised, and this relation is governed by a second pKa which is close to the pKa for ionization of the amide group (HN3-C4O) of the respective ligands. Thus, ionization of the ligand amide group also increases koff. Evidence is presented that the dependence of pH on koff for hDHFR accounts for the shape of the kcat versus pH curve for both hDHFR as well as its F31L variant and contributes to the higher efficiency of hDHFR compared with bacterial DHFR.     

Probing the functional role of phenylalanine-31 of Escherichia coli dihydrofolate reductase by site-directed mutagenesis

The role of Phe-31 of Escherichia coli dihydrofolate reductase in binding and catalysis was probed by amino acid substitution. Phe-31, a strictly conserved residue located in a hydrophobic pocket and interacting with the pteroyl moiety of dihydrofolate (H2F), was replaced by Tyr and Val. The kinetic behavior of the mutant enzymes in general is similar to that of the wild type. The rate-limiting step for both mutant enzymes is the release of tetrahydrofolate (H4F) from the E X NADPH X H4F ternary complex as determined for the wild type. The 2-fold increase in V for the two mutant enzymes arises from faster dissociation of H4F from the enzyme-product complex. The quantitative effect of these mutations is to decrease the rate of hydride transfer, although not to the extent that this step becomes partially rate limiting, but to accelerate the dissociation rates of tetrahydrofolate from product complexes so that the opposing effects are nearly compensating.     

Hydrophobic interactions via mutants of Escherichia coli dihydrofolate reductase: separation of binding and catalysis

The strictly conserved residue leucine-54 of Escherichia coli dihydrofolate reductase forms part of the hydrophobic wall which binds the p-aminobenzoyl side chain of dihydrofolate. In addition to the previously reported glycine-54 mutant, isoleucine-54 and asparagine-54 substitutions have been constructed and characterized with regard to their effects on binding and catalysis. NADP+ and NADPH binding is virtually unaffected with the exception of a 15-fold decrease in NADPH dissociation from the Gly-54 mutant. The synergistic effect of NADPH on tetrahydrofolate dissociation seen in the wild-type enzyme is lost in the isoleucine-54 mutant: little acceleration is seen in tetrahydrofolate dissociation when cofactor is bound, and there is no discrimination between reduced and oxidized cofactor. The dissociation constants for dihydrofolate and methotrexate increase in the order Leu less than Ile less than Asn less than Gly, varying by a maximum factor of 1700 for dihydrofolate and 6300 for methotrexate. Despite these large changes in binding affinity, the hydride transfer rate of 950 s-1 in the wild-type enzyme is decreased by a constant factor of ca. 30 (2 kcal/mol) regardless of the mutant. Thus, the contributions of residue 54 to binding and catalysis appear to have been separated.     

Functional role for Tyr 31 in the catalytic cycle of chicken dihydrofolate reductase

Despite much work, many key aspects of the mechanism of the dihydrofolate reductase (DHFR) catalyzed reduction of dihydrofolate remain unresolved. In bacterial forms of DHFR both substrate and water access to the active site are controlled by the conformation of the mobile M20 loop. In vertebrate DHFRs only one conformation of the residues corresponding to the M20 loop has been observed. Access to the active site was proposed to be controlled by residue 31. MD simulations of chicken DHFR complexed with substrates and cofactor revealed a closing of the side chain of Tyr 31 over the active site on binding of dihydrofolate. This conformational change was dependent on the presence of glutamate on the para-aminobenzoylamide moiety of dihydrofolate. In its absence, the conformation remained open. Although water could enter the active site and hydrogen bond to N5 of dihydrofolate, indicating the feasibility of water as the proton donor, this was not controlled by the conformation of Tyr 31. The water accessibility of the active site was low for both conformations of Tyr 31. However, when hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the average time during which water was found in hydrogen bonding distance to N5 of dihydrofolate in the active site increased almost fivefold. These results indicated that water can serve as the Broensted acid for the protonation of N5 of dihydrofolate during the DHFR catalyzed reduction.