A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.     

Crystal structure of a type II dihydrofolate reductase catalytic ternary complex

Type II dihydrofolate reductase (DHFR) is a plasmid-encoded enzyme that confers resistance to bacterial DHFR-targeted antifolate drugs. It forms a symmetric homotetramer with a central pore which functions as the active site. Its unusual structure, which results in a promiscuous binding surface that accommodates either the dihydrofolate (DHF) substrate or the NADPH cofactor, has constituted a significant limitation to efforts to understand its substrate specificity and reaction mechanism. We describe here the first structure of a ternary R67 DHFR.DHF.NADP+ catalytic complex, resolved to 1.26 A. This structure provides the first clear picture of how this enzyme, which lacks the active site carboxyl residue that is ubiquitous in Type I DHFRs, is able to function. In the catalytic complex, the polar backbone atoms of two symmetry-related I68 residues provide recognition motifs that interact with the carboxamide on the nicotinamide ring, and the N3-O4 amide function on the pteridine ring. This set of interactions orients the aromatic rings of substrate and cofactor in a relative endo geometry in which the reactive centers are held in close proximity. Additionally, a central, hydrogen-bonded network consisting of two pairs of Y69-Q67-Q67'-Y69' residues provides an unusually tight interface, which appears to serve as a "molecular clamp" holding the substrates in place in an orientation conducive to hydride transfer. In addition to providing the first clear insight regarding how this extremely unusual enzyme is able to function, the structure of the ternary complex provides general insights into how a mutationally challenged enzyme, i.e., an enzyme whose evolution is restricted to four-residues-at-a-time active site mutations, overcomes this fundamental limitation.     

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.     

Effects of temperature and viscosity on R67 dihydrofolate reductase catalysis

R67 dihydrofolate reductase (DHFR) is a novel homotetrameric protein that possesses 222 symmetry and a single, voluminous active site pore. This symmetry poses numerous limitations on catalysis; for example, two dihydrofolate (DHF) molecules or two NADPH molecules, or one substrate plus one cofactor can bind. Only the latter combination leads to catalysis. To garner additional information on how this enzyme facilitates transition-state formation, the temperature dependence of binding and catalysis was monitored. The binding of NADPH and DHF is enthalpy-driven. Previous primary isotope effect studies indicate hydride transfer is at least partially rate-determining. Accordingly, the activation energy associated with transition-state formation was measured and is found to be 6.9 kcal/mol (DeltaH(++)(25) = 6.3 kcal/mol). A large entropic component is also found associated with catalysis, TDeltaS(++)(25) = -11.3 kcal/mol. The poor substrate, dihydropteroate, binds more weakly than dihydrofolate (DeltaDeltaG = 1.4 kcal/mol) and displays a large loss in the binding enthalpy value (DeltaDeltaH = 3.8 kcal/mol). The k(cat) value for dihydropteroate reduction is decreased 1600-fold compared to DHF usage. This effect appears to derive mostly from the DeltaDeltaH difference in binding, demonstrating that the glutamate tail is important for catalysis. This result is surprising, as the para-aminobenzoyl-glutamate tail of DHF has been previously shown to be disordered by both NMR and crystallography studies. Viscosity studies were also performed and confirmed that the hydride transfer rate is not sensitive to sucrose addition. Surprisingly, binding of DHF, by both K(m) and K(d) determination, was found to be sensitive to added viscogens, suggesting a role for water in DHF binding.     

NMR studies of the interaction of a type II dihydrofolate reductase with pyridine nucleotides reveal unexpected phosphatase and reductase activity

The interaction of type II R67 dihydrofolate reductase (DHFR) with its cofactor nicotinamide adenine dinucleotide phosphate (NADP(+)) has been studied using nuclear magnetic resonance (NMR). Doubly labeled [U-(13)C,(15)N]DHFR was obtained from Escherichia coli grown on a medium containing [U-(13)C]-D-glucose and (15)NH(4)Cl, and the 16 disordered N-terminal amino acids were removed by treatment with chymotrypsin. Backbone and side chain NMR assignments were made using triple-resonance experiments. The degeneracy of the amide (1)H and (15)N shifts of the tetrameric DHFR was preserved upon addition of NADP(+), consistent with kinetic averaging among equivalent binding sites. Analysis of the more titration-sensitive DHFR amide resonances as a function of added NADP(+) gave a K(D) of 131 +/- 50 microM, consistent with previous determinations using other methodology. We have found that the (1)H spectrum of NADP(+) in the presence of the R67 DHFR changes as a function of time. Comparison with standard samples and mass spectrometric analysis indicates a slow conversion of NADP(+) to NAD(+), i.e., an apparent NADP(+) phosphatase activity. Studies of this activity in the presence of folate and a folate analogue support the conclusion that this activity results from an interaction with the DHFR rather than a contaminating phosphatase. (1)H NMR studies of a mixture of NADP(+) and NADPH in the presence of the enzyme reveal that a ternary complex forms in which the N-4A and N-4B nuclei of the NADPH are in the proximity of the N-4 and N-5 nuclei of NADP(+). Studies using the NADP(+) analogue acetylpyridine adenosine dinucleotide phosphate (APADP(+)) demonstrated a low level of enzyme-catalyzed hydride transfer from NADPH. Analysis of DHFR backbone dynamics revealed little change upon binding of NADP(+). These additional catalytic activities and dynamic behavior are in marked contrast to those of type I DHFR.     

Conformational changes in the active site loops of dihydrofolate reductase during the catalytic cycle

Escherichia coli dihydrofolate reductase (DHFR) has several flexible loops surrounding the active site that play a functional role in substrate and cofactor binding and in catalysis. We have used heteronuclear NMR methods to probe the loop conformations in solution in complexes of DHFR formed during the catalytic cycle. To facilitate the NMR analysis, the enzyme was labeled selectively with [(15)N]alanine. The 13 alanine resonances provide a fingerprint of the protein structure and report on the active site loop conformations and binding of substrate, product, and cofactor. Spectra were recorded for binary and ternary complexes of wild-type DHFR bound to the substrate dihydrofolate (DHF), the product tetrahydrofolate (THF), the pseudosubstrate folate, reduced and oxidized NADPH cofactor, and the inactive cofactor analogue 5,6-dihydroNADPH. The data show that DHFR exists in solution in two dominant conformational states, with the active site loops adopting conformations that closely approximate the occluded or closed conformations identified in earlier X-ray crystallographic analyses. A minor population of a third conformer of unknown structure was observed for the apoenzyme and for the disordered binary complex with 5,6-dihydroNADPH. The reactive Michaelis complex, with both DHF and NADPH bound to the enzyme, could not be studied directly but was modeled by the ternary folate:NADP(+) and dihydrofolate:NADP(+) complexes. From the NMR data, we are able to characterize the active site loop conformation and the occupancy of the substrate and cofactor binding sites in all intermediates formed in the extended catalytic cycle. In the dominant kinetic pathway under steady-state conditions, only the holoenzyme (the binary NADPH complex) and the Michaelis complex adopt the closed loop conformation, and all product complexes are occluded. The catalytic cycle thus involves obligatory conformational transitions between the closed and occluded states. Parallel studies on the catalytically impaired G121V mutant DHFR show that formation of the closed state, in which the nicotinamide ring of the cofactor is inserted into the active site, is energetically disfavored. The G121V mutation, at a position distant from the active site, interferes with coupled loop movements and appears to impair catalysis by destabilizing the closed Michaelis complex and introducing an extra step into the kinetic pathway.     

Breaking symmetry: mutations engineered into R67 dihydrofolate reductase,a D2 symmetric homotetramer possessing a single active site pore

R67 dihydrofolate reductase (DHFR) is an enzyme, encoded by an R-plasmid, that confers resistance to the antibacterial agent trimethoprim. This homotetramer possesses a single active site pore and exact 222 symmetry. The symmetry imposes constraints on the ability of the enzyme to optimize binding of the substrate, dihydrofolate (DHF), and the cofactor, NADPH, resulting in a "one site fits both ligands" approach. This approach allows formation of either a NADPH.NADPH, dihydrofolate.dihydrofolate, or NADPH.dihydrofolate complex. The first two complexes are nonproductive, while the third is the productive catalytic species. To break the symmetry of the active site, a tandem array of four R67 DHFR genes has been linked in frame, allowing individual manipulation of each gene copy. Various numbers and combinations of asymmetric Q67H mutations have been engineered into the tandem gene array. The Q67H mutation was chosen for investigation as it was previously found to tighten binding to both dihydrofolate and NADPH by approximately 100-fold in homotetrameric R67 DHFR [Park, H., Bradrick, T. D., and Howell, E. E. (1997) Protein Eng. 10, 1415-1424]. Nonadditive effects on ligand binding are observed when one to four mutations are inserted, indicating either conformational changes in the protein or different cooperativity patterns in the ligand-ligand interactions. From steady state kinetics, addition of Q67H mutations does not drastically affect formation of the NADPH.dihydrofolate complex; however, a large energy difference between the productive and nonproductive complexes is no longer maintained. A role for Q67 in discriminating between these various states is proposed. Since theories of protein evolution suggest gene duplication followed by accumulation of mutations can lead to divergence of activity, this study is a first step toward asking if introduction of asymmetric mutations in the quadrupled R67 DHFR gene can lead to optimization of ligand binding sites.     

Gammaherpesviruses encode functional dihydrofolate reductase activity

We overexpressed and purified from Escherichia coli the dihydrofolate reductase (DHFR) of the gammaherpesviruses human herpesvirus 8 (HHV-8), herpesvirus saimiri (HVS), and rhesus rhadinovirus (RRV). All three enzymes proved catalytically active. The K(m) value of HHV-8 DHFR for dihydrofolate (DHF) was 2.02+/-0.44 microM, that of HVS DHFR was 4.31+/-0.56 microM, and that of RRV DHFR is 7.09+/-0.11 microM. These values are approximately 5-15-fold higher than the K(m) value reported for the human DHFR. The K(m) value of HHV-8 DHFR for NADPH was 1.31+/-0.23 microM, that of HVS DHFR was 3.78+/-0.61 microM, and that of RRV DHFR was 7.47+/-0.59 microM. These values are similar or slightly higher than the corresponding K(m) value of the human enzyme. Methotrexate, aminopterin, trimethoprim, pyrimethamine, and N(alpha)-(4-amino-4-deoxypteroyl)-N(delta)-hemiphthaloyl-L-ornithine (PT523), all well-known folate antagonists, inhibited the DHFR activity of the three gammaherpesviruses competitively with respect to DHF but proved markedly less inhibitory to the viral than towards the human enzyme.     

Structure and hydride transfer mechanism of a moderate thermophilic dihydrofolate reductase from Bacillus stearothermophilus and comparison to its mesophilic and hyperthermophilic homologues

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.     

Diagnostic chemical shift markers for loop conformation and substrate and cofactor binding in dihydrofolate reductase complexes

Heteronuclear NMR methods have been used to probe the conformation of four complexes of Escherichia coli dihydrofolate reductase (DHFR) in solution. (1)H(N), (15)N, and (13)C(alpha) resonance assignments have been made for the ternary complex with folate and oxidized NADP(+) cofactor and the ternary complex with folate and a reduced cofactor analog, 5,6-dihydroNADPH. The backbone chemical shifts have been compared with those of the binary complex of DHFR with the substrate analog folate and the binary complex with NADPH (the holoenzyme). Analysis of (1)H(N) and (15)N chemical shifts has led to the identification of marker resonances that report on the active site conformation of the enzyme. Other backbone amide resonances report on the presence of ligands in the pterin binding pocket and in the adenosine and nicotinamide-ribose binding sites of the NADPH cofactor. The chemical shift data indicate that the enzyme populates two dominant structural states in solution, with the active site loops in either the closed or occluded conformations defined by X-ray crystallography; there is no evidence that the open conformation observed in some X-ray structures of E. coli DHFR are populated in solution.     

Role of S65, Q67, I68, and Y69 residues in homotetrameric R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) shares no sequence or structural homology with chromosomal DHFRs. This enzyme arose recently in response to the clinical use of the antibacterial drug trimethoprim. R67 DHFR is a homotetramer possessing a single active site pore. A high-resolution crystal structure shows the homotetramer possesses exact 222 symmetry [Narayana, N., et al. (1995) Nat. Struct. Biol. 2, 1018-1025]. This symmetry dictates four symmetry-related binding sites must exist for each substrate as well as each cofactor. Isothermal titration calorimetry studies, however, indicate only two molecules bind: either two dihydrofolate molecules, two NADPH molecules, or one substrate and one cofactor [Bradrick, T. D., et al. (1996) Biochemistry 35, 11414-11424]. The latter is the productive ternary complex. To evaluate the role of S65, Q67, I68, and Y69 residues, located near the center of the active site pore, site-directed mutagenesis was performed. One mutation in the gene creates four mutations per active site pore which typically result in large cumulative effects. Steady state kinetic data indicate the mutants have altered K(m) values for both cofactor and substrate. For example, the Y69F R67 DHFR displays an 8-fold increase in the K(m) for dihydrofolate and a 20-fold increase in the K(m) for NADPH. Residues involved in ligand binding in R67 DHFR display very little, if any, specificity, consistent with their possessing dual roles in binding. These results support a model where R67 DHFR utilizes an unusual "hot spot" binding surface capable of binding both ligands and indicate this enzyme has adopted a novel yet simple approach to catalysis.     

Insight into the molecular mechanism about lowered dihydrofolate binding affinity to dihydrofolate reductase-like 1 (DHFRL1)

Human dihydrofolate reductase-like 1 (DHFRL1) has been identified as a second human dihydrofolate reductase (DHFR) enzyme. Although DHFRL1 have high sequence homology with human DHFR, dihydrofolate (DHF) exhibits a lowered binding affinity to DHFRL1 and the corresponding molecular mechanism is still unknown. To address this question, we studied the binding of DHF to DHFRL1 and DHFR by using molecular dynamics simulation. Moreover, to investigate the role the 24th residue of DHFR/DHFRL1 plays in DHF binding, R24W DHFRL1 mutant was also studied. The van der Waals interaction are more crucial for the total DHF binding energies, while the difference between the DHF binding energies of human DHFR and DHFRL1 can be attributed to the electrostatic interaction and the polar desolvation free energy. More specifically, lower DHF affinity to DHFRL1 can be mainly attributed to the reduction of net electrostatic interactions of residues Arg32 and Gln35 of DHFRL1 with DHF as being affected by Arg24. The side chain of Arg24 in DHFRL1 can extend deeply into the binding sites of DHF and NADPH, and disturb the DHF binding by steric effect, which rarely happens in human DHFR and R24W DHFRL1 mutant. Additionally, the conformation of loop I in DHFRL1 was also studied in this work. Interestingly, the loop conformation resemble to normal closed state of Escherichia coli DHFR other than the closed state of human DHFR. We hope this work will be useful to understand the general characteristics of DHFRL1.     

R67, the other dihydrofolate reductase: rational design of an alternate active site configuration

R67 dihydrofolate reductase (DHFR) bears no sequence or structural homologies with chromosomal DHFRs. The gene for this enzyme produces subunits that are 78 amino acids long, which assemble into a homotetramer possessing 222 symmetry. More recently, a tandem array of four gene copies linked in-frame was constructed, which produces a monomer containing 312 amino acids named Quad3. Asymmetric mutations in Quad3 have also been constructed to probe the role of Q67 and K32 residues in catalysis. This present study mixes and matches mutations to determine if the Q67H mutation, which tightens binding approximately 100-fold to both dihydrofolate (DHF) and NADPH, can help rescue the K32M mutation. While the latter mutation weakens DHF binding over 60-fold, it concurrently increases kcat by a factor of 5. Two Q67H mutations were added to gene copies 1 and 4 in conjunction with the K32M mutation in gene copies 1 and 3. Addition of these Q67H mutations tightens binding 40-fold, and the catalytic efficiency (kcat/Km(DHF)) of the resulting protein is similar to that of Quad3. Since these Q67H mutations can mostly compensate for the K32M lesion, K32 must not be necessary for DHF binding. Another multimutant combines the K32M mutation in gene copies 1 and 3 with the Q67H mutation in all gene copies. This mutant is inhibited by DHF but not NADPH, indicating that NADPH binds only to the wild type half of the pore, while DHF can bind to either the wild type or mutant half of the pore. This inhibition pattern contrasts with the mutant containing only the Q67H substitution in all four gene copies, which is severely inhibited by both NADPH and substrate. Since gene duplication and divergence are evolutionary tools for gaining function, these constructs are a first step toward building preferences for NADPH and DHF in each half of the active site pore of this primitive enzyme.     

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ～5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.     

The kinetic mechanism of wild-type and mutant mouse dihydrofolate reductases

A kinetic mechanism is presented for mouse dihydrofolate reductase that predicts all the steady-state parameters and full time-course kinetics. This mechanism was derived from association and dissociation rate constants and pre-steady-state transients by using stopped-flow fluorescence and absorbance measurements. The major features of this kinetic mechanism are as follows: (1) the two native enzyme conformers, E1 and E2, bind ligands with varying affinities although only one conformer, E1, can support catalysis in the forward direction, (2) tetrahydrofolate dissociation is the rate-limiting step under steady-state turnover at low pH, and (3) the pH-independent rate of hydride transfer from NADPH to dihydrofolate is fast (khyd = 9000 s-1) and favorable (Keq = 100). The overall mechanism is similar in form to the Escherichia coli kinetic scheme (Fierke et al., 1987), although several differences are observed: (1) substrates and products predominantly bind the same form of the E. coli enzyme, and (2) the hydride transfer rate from NADPH to either folate or dihydrofolate is considerably faster for the mouse enzyme. The role of Glu-30 (Asp-27 in E. coli) in mouse DHFR has also been examined by using site-directed mutagenesis as a potential source of these differences. While aspartic acid is strictly conserved in all bacterial DHFRs, glutamic acid is conserved in all known eucaryotes. The two major effects of substituting Asp for Glu-30 in the mouse enzyme are (1) a decreased rate of folate reduction and (2) an increased rate of hydride transfer from NADPH to dihydrofolate.(ABSTRACT TRUNCATED AT 250 WORDS)     

Role of Lys-32 residues in R67 dihydrofolate reductase probed by asymmetric mutations

R67 dihydrofolate reductase (R67 DHFR) is a novel protein encoded by an R-plasmid that confers resistance to the antibiotic, trimethoprim. This homotetrameric enzyme possesses 222 symmetry, which imposes numerous constraints on the single active site pore, including a "one-site-fits-both" strategy for binding its ligands, dihydrofolate (DHF) and NADPH. Previous studies uncovered salt effects on binding and catalysis (Hicks, S. N., Smiley, R. D., Hamilton, J. B., and Howell, E. E. (2003) Biochemistry 42, 10569-10578), however the one or more residues that participate in ionic contacts with the negatively charged tail of DHF as well as the phosphate groups in NADPH were not identified. Several studies predict that Lys-32 residues were involved, however mutations at this residue destabilize the R67 DHFR homotetramer. To study the role of Lys-32 in binding and catalysis, asymmetric K32M mutations have been utilized. To create asymmetry, individual mutations were added to a tandem array of four in-frame gene copies. These studies show one K32M mutation is tolerated quite well, whereas addition of two mutations has variable effects. Two double mutants, K32M:1+2 and K32M: 1+4, which place the mutations on opposite sides of the pore, reduce kcat. However a third double mutant, K32M: 1+3, that places two mutations on the same half pore, enhances kcat 4- to 5-fold compared with the parent enzyme, albeit at the expense of weaker binding of ligands. Because the kcat/Km values for this double mutant series are similar, these mutations appear to have uncovered some degree of non-productive binding. This non-productive binding mode likely arises from formation of an ionic interaction that must be broken to allow access to the transition state. The K32M:1+3 mutant data suggest this interaction is an ionic interaction between Lys-32 and the charged tail of dihydrofolate. This unusual catalytic scenario arises from the 222 symmetry imposed on the single active site pore.     

Calorimetric studies of ligand binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) is a novel bacterial protein that possesses 222 symmetry and a single active site pore. Although the 222 symmetry implies that four symmetry-related binding sites must exist for each substrate as well as for each cofactor, various studies indicate only two molecules bind. Three possible combinations include two dihydrofolate molecules, two NADPH molecules, or one substrate plus one cofactor. The latter is the productive ternary complex. To explore the role of various ligand substituents during binding, numerous analogues, inhibitors, and fragments of NADPH and/or folate were used in both isothermal titration calorimetry (ITC) and K(i) studies. Not surprisingly, as the length of the molecule is shortened, affinity is lost, indicating that ligand connectivity is important in binding. The observed enthalpy change in ITC measurements arises from all components involved in the binding process, including proton uptake. As a buffer dependence for binding of folate was observed, this likely correlates with perturbation of the bound N3 pK(a), such that a neutral pteridine ring is preferred for pairwise interaction with the protein. Of interest, there is no enthalpic signal for binding of folate fragments such as dihydrobiopterin where the p-aminobenzoylglutamate tail has been removed, pointing to the tail as providing most of the enthalpic signal. For binding of NADPH and its analogues, the nicotinamide carboxamide is quite important. Differences between binary (binding of two identical ligands) and ternary complex formation are observed, indicating interligand pairing preferences. For example, while aminopterin and methotrexate both form binary complexes, albeit weakly, neither readily forms ternary complexes with the cofactor. These observations suggest a role for the O4 atom of folate in a pairing preference with NADPH, which ultimately facilitates catalysis.     

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.     

Hydride transfer by dihydrofolate reductase. Causes and consequences of the wide range of rates exhibited by bacterial and vertebrate enzymes

Transient and steady-state kinetics have been examined for dihydrofolate reductase (DHFR) from a number of sources. Rates of hydride transfer at pH 7.65 cover a wide range, from 7 s-1 for DHFR from a strain of Lactobacillus casei (LCDHFR1) to 3000 s-1 for recombinant human DHFR (rHDHFR). In all cases as the pH is increased from 7 to 10, Vmax for the steady-state reaction decreases, and DVmax, the primary isotope effect, increases. This indicates a decrease in the rate of hydride transfer with increasing pH. The cross-over points, at which rates of product release and hydride transfer become equal, were calculated to occur at DVmax = 2.34. The higher the rate of hydride transfer at pH 7.65, the higher the pH of the cross-over point. For LCDHFR1 the low rate of hydride transfer results in this process being partially rate-limiting for the steady-state reaction even at pH 5, with a cross-over point at about pH 7. At pH 7.65 the burst phase associated with the initial conversion of enzyme-bound substrates to enzyme-bound products has an isotope effect of 3 or higher for LCDHFR and for DHFR from Escherichia coli (ECDHFR). In contrast, the vertebrate DHFRs (bovine, BDHFR; chicken, CDHFR; and rHDHFR) exhibit a burst of product formation which is only partially limited by hydride transfer at this pH (Dkb: 2.3, 2.2, and 2.1, respectively). An obligatory isomerization of the ternary substrate complex or of the ternary product complex is postulated to be partially rate-limiting for the vertebrate enzymes. At pH 5 LCDHFR1 and ECDHFR also exhibit evidence of such a rate-limiting obligatory conformational transition of the substrate or product ternary complex.