Microsecond subdomain folding in dihydrofolate reductase

The characterization of microsecond dynamics in the folding of multisubdomain proteins has been a major challenge in understanding their often complex folding mechanisms. Using a continuous-flow mixing device coupled with fluorescence lifetime detection, we report the microsecond folding dynamics of dihydrofolate reductase (DHFR), a two-subdomain α/β/α sandwich protein known to begin folding in this time range. The global dimensions of early intermediates were monitored by Förster resonance energy transfer, and the dynamic properties of the local Trp environments were monitored by fluorescence lifetime detection. We found that substantial collapse occurs in both the locally connected adenosine binding subdomain and the discontinuous loop subdomain within 35 μs of initiation of folding from the urea unfolded state. During the fastest observable ∼ 550 μs phase, the discontinuous loop subdomain further contracts, concomitant with the burial of Trp residue(s), as both subdomains achieve a similar degree of compactness. Taken together with previous studies in the millisecond time range, a hierarchical assembly of DHFR—in which each subdomain independently folds, subsequently docks, and then anneals into the native conformation after an initial heterogeneous global collapse—emerges. The progressive acquisition of structure, beginning with a continuously connected subdomain and spreading to distal regions, shows that chain entropy is a significant organizing principle in the folding of multisubdomain proteins and single-domain proteins. Subdomain folding also provides a rationale for the complex kinetics often observed.     

An examination of the involvement of proline peptide isomerization in protein folding.

The proline peptide isomerization model of protein folding predicts that the fraction of denatured polypeptide chains which rapidly fold should be quantitatively related to the numbers of cis and trans proline residues in the folded polypeptide conformation. However, we find that neither the comparative nor the absolute kinetic pattern for folding of the homologous proteins, tuna heart and horse heart ferricytochrome c which differ by one proline residue, is compatible with the quantitative predictions of the proline peptide isomerization model.     

The coordination of the isomerization of a conserved non-prolyl cis peptide bond with the rate-limiting steps in the folding of dihydrofolate reductase

The propensity for peptide bonds to adopt the trans configuration in native and unfolded proteins, and the relatively slow rates of cis-trans isomerization reactions, imply that the formation of cis peptide bonds in native conformations are likely to limit folding reactions. The role of the conserved cis Gly95-Gly96 peptide bond in dihydrofolate reductase (DHFR) from Escherichia coli was examined by replacing Gly95 with alanine. The introduction of a beta carbon at position 95 is expected to increase the propensity for the trans isomer and perturb the isomerization reaction required to reach the native conformation. Although G95A DHFR is 1.30 kcal mol(-1) less stable than the wild-type protein, it adopts a well-folded structure that can be chemically denatured in a cooperative fashion. The mutant protein also retains the complex refolding kinetic pattern attributed to a parallel-channel mechanism in wild-type DHFR. The spectroscopic response upon refolding monitored by Trp fluorescence and the absence of a Trp/Trp exciton coupling apparent in the far-UV CD spectrum of the wild-type protein, however, indicated that the tertiary structure of the folded state for G95A DHFR is altered. The addition of methotrexate (MTX), a tight-binding inhibitor, to folded G95A DHFR restored the exciton coupling and the fluorescence properties through five slow kinetic events whose relaxation times are independent of the ligand and the denaturant concentrations. The results were interpreted to mean that MTX-binding drives the formation of the cis isomer of the peptide bond between Ala95 and Gly96 in five compact and stable but not wild-type-like conformations that contain the trans isomer. Folding studies in the presence of MTX for both wild-type and G95A DHFR support the notion that the cis peptide bond between Gly95 and Gly96 in the wild-type protein forms during four parallel rate-limiting steps, which are primarily controlled by folding reactions, and lead directly to a set of native, or native-like, conformers. The isomerization of the cis peptide bond is not a source of the parallel channels that characterize the complex folding mechanism for DHFR.     

A hydrophobic cluster forms early in the folding of dihydrofolate reductase

The rapid kinetic phase that leads from unfolded species to transient folding intermediates in dihydrofolate reductase from Escherichia coli was examined by site-directed mutagenesis and by physicochemical means. The absence of this fluorescence-detected phase in the refolding of the Trp-74Phe mutant protein strongly implies that this early phase in refolding can be assigned to just one of the five Trp residues in the protein, Trp-74. In addition, water-soluble fluorescence quenching agents, iodide and cesium, have a much less significant effect on this early step in refolding than on the slower phases that lead to native and native-like conformers. These and other data imply that an important early event in the folding of dihydrofolate reductase is the formation of a hydrophobic cluster which protects Trp-74 from solvent.     

Stability and folding of dihydrofolate reductase from the hyperthermophilic bacterium Thermotoga maritima.

Dihydrofolate reductase (DHFR) has been a well-established model system for protein folding. The enzyme DHFR from the hyperthermophilic bacterium Thermotoga maritima (TmDHFR) displays distinct adaptations toward high temperatures at the level of both structure and stability. The enzyme represents an extremely stable dimer; no isolated structured monomers could be detected in equilibrium or during unfolding. The equilibrium unfolding strictly follows the two-state model for a dimer (N(2) right harpoon over left harpoon 2U), with a free energy of stabilization of DeltaG = -142 +/- 10 kJ/mol at 15 degrees C. The two-state model is applicable over the whole temperature range (5-70 degrees C), yielding a DeltaG vs T profile with maximum stability at around 35 degrees C. There is no flattening of the stability profile. Instead, the enhanced thermostability is characterized by shifts toward higher overall stability and higher temperature of maximum stability. TmDHFR unfolds in a highly cooperative manner via a nativelike transition state without intermediates. The unfolding reaction is much slower (ca. 10(8) times) compared to DHFR from Escherichia coli (EcDHFR). In contrast to EcDHFR, no evidence for heterogeneity of the native state is detectable. Refolding proceeds via at least two intermediates and a burst-phase of rather low amplitude. Reassociation of monomeric intermediates is not rate-limiting on the folding pathway due to the high association constant of the dimer.     

Kinetic mechanism and catalysis of a native-state prolyl isomerization reaction

Folding of tendamistat is a rapid two-state process for the majority of the unfolded molecules. In fluorescence-monitored refolding kinetics about 8% of the unfolded molecules fold slowly (lambda=0.083s(-1)), limited by peptidyl-prolyl cis-trans isomerization. This is significantly less than expected from the presence of three trans prolyl-peptide bonds in the native state. In interrupted refolding experiments we detected an additional very slow folding reaction (lambda=0.008s(-1) at pH 2) with an amplitude of about 12%. This reaction is caused by the interconversion of a highly structured intermediate to native tendamistat. The intermediate has essentially native spectroscopic properties and about 2% of it remain populated in equilibrium after folding is complete. Catalysis by human cyclophilin 18 identifies this very slow reaction as a prolyl isomerization reaction. This shows that prolyl-isomerases are able to efficiently catalyze native state isomerization reactions, which allows them to influence biologically important regulatory conformational transitions. Folding kinetics of the proline variants P7A, P9A, P50A and P7A/P9A show that the very slow reaction is due to isomerization of the Glu6-Pro7 and Ala8-Pro9 peptide bonds, which are located in a region that makes strong backbone and side-chain interactions to both beta-sheets. In the P50A variant the very slow isomerization reaction is still present but native state heterogeneity is not observed any more, indicating a long-range destabilizing effect on the alternative native state relative to N. These results enable us to include all prolyl and non-prolyl peptide bond isomerization reactions in the folding mechanism of tendamistat and to characterize the kinetic mechanism and the energetics of a native-state prolyl isomerization reaction.     

Conformational relaxation following hydride transfer plays a limiting role in dihydrofolate reductase catalysis

The catalytic cycle of an enzyme is frequently associated with conformational changes that may limit maximum catalytic throughput. In Escherichia coli dihydrofolate reductase, release of the tetrahydrofolate (THF) product is the rate-determining step under physiological conditions and is associated with an "occluded" to "closed" conformational change. In this study, we demonstrate that in dihydrofolate reductase the closed to occluded conformational change in the product ternary complex (E.THF.NADP (+)) also gates progression through the catalytic cycle. Using NMR relaxation dispersion, we have measured the temperature and pH dependence of microsecond to millisecond time scale backbone dynamics of the occluded E.THF.NADP (+) complex. Our studies indicate the presence of three independent dynamic regions, associated with the active-site loops, the cofactor binding cleft, and the C-terminus and an adjacent loop, which fluctuate into discrete conformational substates with different kinetic and thermodynamic parameters. The dynamics of the C-terminally associated region is pH-dependent (p K a < 6), but the dynamics of the active-site loops and cofactor binding cleft are pH-independent. The active-site loop dynamics access a closed conformation, and the accompanying closed to occluded rate constant is comparable to the maximum pH-independent hydride transfer rate constant. Together, these results strongly suggest that the closed to occluded conformational transition in the product ternary complex is a prerequisite for progression through the catalytic cycle and that the rate of this process places an effective limit on the maximum rate of the hydride transfer step.     

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

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

Coupled kinetic traps in cytochrome c folding: His-heme misligation and proline isomerization

The effect of His-heme misligation on folding has been investigated for a triple mutant of yeast iso-2 cytochrome c (N26H,H33N,H39K iso-2). The variant contains a single misligating His residue at position 26, a location at which His residues are found in several cytochrome c homologues, including horse, tuna, and yeast iso-1. The amplitude for fast phase folding exhibits a strong initial pH dependence. For GdnHCl unfolded protein at an initial pH<5, the observed refolding at final pH 6 is dominated by a fast phase (tau(2f)=20 ms, alpha(2f)=90 %) that represents folding in the absence of misligation. For unfolded protein at initial pH 6, folding at final pH 6 occurs in a fast phase of reduced amplitude (alpha(2f) approximately 20 %) but the same rate (tau(2f)=20 ms), and in two slower phases (tau(m)=6-8 seconds, alpha(m) approximately 45 %; and tau(1b)=16-20 seconds, alpha(1b) approximately 35 %). Double jump experiments show that the initial pH dependence of the folding amplitudes results from a slow pH-dependent equilibrium between fast and slow folding species present in the unfolded protein. The slow equilibrium arises from coupling of the His protonation equilibrium to His-heme misligation and proline isomerization. Specifically, Pro25 is predominantly in trans in the unligated low-pH unfolded protein, but is constrained in a non-native cis isomerization state by His26-heme misligation near neutral pH. Refolding from the misligated unfolded form proceeds slowly due to the large energetic barrier required for proline isomerization and displacement of the misligated His26-heme ligand.     

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.     

Kinetics of the formation and isomerization of methotrexate complexes of recombinant human dihydrofolate reductase

The kinetics of inhibitor binding to highly purified recombinant human dihydrofolate reductase (rHDHFR) have been examined. Methotrexate (MTX) binds rapidly (kon = 1.0 x 10(8) M-1 s-1) and tightly (koff/kon = 210 pM) to the preformed complex of rHDHFR with NADPH. The initial association reaction between rHDHFR.NADPH and MTX is followed by an isomerization of the resulting complex (kiso = 0.4 s-1) leading to a new conformer in which MTX is bound even more tightly (Ki = 3.4 pM). Similar results have been obtained with a major metabolite of MTX having four additional glutamate residues for which Ki = 1.4 pM. 7-HydroxyMTX, another major metabolite of MTX, is a weak inhibitor of rHDHFR (Ki = 8.9 nM), and a polyglutamate form of this metabolite is an equally weak inhibitor (Ki = 9.9 nM), so that the addition of glutamate residues to MTX or 7-hydroxyMTX has little effect on their binding. It follows that the significance of MTX polyglutamate formation relates to other roles such as increasing the cytotoxicity of MTX by prolonging intracellular retention of the drug. Another antifolate, trimethoprim, binds tightly to dihydrofolate reductases from bacterial sources, but weakly to rHDHFR in the ternary complex (KD = 0.5 microM). Although the association step is rapid (kon = 0.4 x 10(8) M-1 s-1), the dissociation rate is also rapid (koff = 15 s-1). Furthermore, there is no isomerization of the ternary complex of trimethoprim with rHDHFR, in contrast to the known isomerization of complexes of trimethoprim with bacterial dihydrofolate reductases.     

Folding of chymotrypsin inhibitor 2. 2. Influence of proline isomerization on the folding kinetics and thermodynamic characterization of the transition state of folding

The refolding of chymotrypsin inhibitor 2 (CI2) is, at least, a triphasic process. The rate constants are 53 s-1 for the major phase (77% of the total amplitude) and 0.43 and 0.024 s-1 for the slower phases (23% of the total amplitude) at 25 degrees C and pH 6.3. The multiphase nature of the refolding reaction results from heterogeneity in the denatured state because of proline isomerization. The fast phase corresponds to the refolding of the fraction of protein that has all its prolines in a native trans conformation in the denatured state. It is not catalyzed by peptidyl-prolyl isomerase. The rate-limiting step of folding for the slower phases, however, is proline isomerization, and they are both catalyzed by peptidyl-prolyl isomerase. The slowest phase has properties consistent with a process involving proline isomerization in a denatured state. In particular, the activation enthalpy is large, 16 kcal mol-1 K-1, and the rate is independent of guanidinium chloride concentration ([GdnHCl]). In comparison, the intermediate phase shows properties consistent with a process involving proline isomerization in a partially structured state. The activation enthalpy is small, 8 kcal mol-1 K-1, and the rate has a strong dependence on [GdnHCl]. Temperature dependences of the rate constants for unfolding and for the fast refolding phase, both in the absence and in the presence of GdnHCl, were used to characterize the thermodynamic nature of the transition state and its relative exposure to solvent. The Eyring plot for unfolding is linear, indicating that there is relatively little change in heat capacity between native state and transition state.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Transient intermediates in the folding of dihydrofolate reductase as detected by far-ultraviolet circular dichroism spectroscopy

The kinetics of the reversible folding and unfolding of Escherichia coli dihydrofolate reductase have been studied by stopped-flow circular dichroism in the peptide region at pH 7.8 and 15 degrees C. The reactions were induced by concentration jumps of a denaturant, urea. The method can detect various intermediates transiently populated in the reactions although the equilibrium unfolding of the protein is apparently approximated by a two-state reaction. The results can be summarized as follows. (1) From transient circular dichroism spectra measured as soon as the refolding is started, a substantial amount of secondary structure is formed in the burst phase, i.e., within the dead time of stopped-flow mixing (18 ms). (2) The kinetics from this burst-phase intermediate to the native state are multiphasic, consisting of five phases designated as tau 1, tau 2, tau 3, tau 4, and tau 5 in increasing order of the reaction rate. Measurements of the kinetics at various wavelengths have provided kinetic difference circular dichroism spectra for the individual phases. (3) The tau 5 phase shows a kinetic difference spectrum consistent with an exciton contribution of two aromatic residues in the peptide CD region. The absence of the tau 5 phase in a mutant protein, in which Trp 74 is replaced by leucine, suggests that Trp 74 is involved in the exciton pair and that the tau 5 phase reflects the formation of a hydrophobic cluster around Trp 74. From the similarity of the kinetic difference spectrum to the difference between the native spectra of the mutant and wild-type proteins, it appears that Trp 47 is the partner in the exciton pair and that the structure formed in the tau 5 phase persists during the later stages of folding. (4) The later stages of folding show kinetic difference spectra that can be interpreted by rearrangement of secondary structure, particularly the central beta sheet of the protein. The pairwise similarities in the spectrum between the tau 3 and tau 4 phases, and between the tau 1 and tau 2 phases, also suggest the presence of two parallel folding channels for refolding. (5) The unfolding kinetics show three to four phases and are interpreted in terms of the presence of multiple native species. The total ellipticity change in kinetic unfolding reaction, however, agrees with the ellipticity difference between the native and unfolding states, indicating the absence of the burst phase in unfolding.(ABSTRACT TRUNCATED AT 400 WORDS)     

Probing the interactions between the folding elements early in the folding of Escherichia coli dihydrofolate reductase by systematic sequence perturbation analysis

One of the necessary conditions for a protein to be foldable is the presence of a complete set of folding elements (FEs) that are short contiguous peptide segments distributed over an amino acid sequence. Previous studies indicated the FE assembly model of protein folding, in which the FEs interact with each other and coalesce to form an intermediate(s) early in the folding reaction. This suggests that a clue to the understanding of the determinants of protein foldability can be found by investigating how the FEs interact with each other early in the folding and thereby elucidating roles of the FEs in protein folding. To reveal the formation process of FE-FE interactions, we studied the early folding events of Escherichia coli dihydrofolate reductase (DHFR) utilizing systematic sequence perturbation analysis. Here, systematic single amino acid substitutions were introduced inside of the FEs (W30X in FE2, V40X in FE3, N59X in FE4, and I155X in FE10; X refers to various amino acid residues), and their kinetic refolding reactions were measured by stopped-flow circular dichroism and fluorescence. We show that the interactions around Trp30 and Ile155 are formed in the burst phase intermediate, while those around Val40 and Asn59 are formed in the transition state of the subsequent folding phase (tau5-phase) and in much later processes, respectively. These and previous results suggest that FE2 and FE10, and also FE1 and FE7, involved in the loop subdomain of DHFR, interact with each other within a millisecond time range, while the stable FE3-FE4 interactions are formed in the later processes. This may highlight the important roles of the FEs mainly inside of the loop subdomain in formation of the burst phase intermediate having a hydrophobic cluster and native-like overall topology and in acquisition of the foldability of DHFR.     

Role of ionic interactions in ligand binding and catalysis of R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), which catalyzes the NADPH dependent reduction of dihydrofolate to tetrahydrofolate, belongs to a type II family of R-plasmid encoded DHFRs that confer resistance to the antibacterial drug trimethoprim. Crystal structure data reveals this enzyme is a homotetramer that possesses a single active site pore. Only two charged residues in each monomer are located near the pore, K32 and K33. Site-directed mutants were constructed to probe the role of these residues in ligand binding and/or catalysis. As a result of the 222 symmetry of this enzyme, mutagenesis of one residue results in modification at four related sites. All mutants at K32 affected the quaternary structure, producing an inactive dimer. The K33M mutant shows only a 2-4-fold effect on K(m) values. Salt effects on ligand binding and catalysis for K33M and wildtype R67 DHFRs were investigated to determine if these lysines are involved in forming ionic interactions with the negatively charged substrates, dihydrofolate (overall charge of -2) and NADPH (overall charge of -3). Binding studies indicate that two ionic interactions occur between NADPH and R67 DHFR. In contrast, the binding of folate, a poor substrate, to R67 DHFR.NADPH appears weak as a titration in enthalpy is lost at low ionic strength. Steady-state kinetic studies for both wild type (wt) and K33M R67 DHFRs also support a strong electrostatic interaction between NADPH and the enzyme. Interestingly, quantitation of the observed salt effects by measuring the slopes of the log of ionic strength versus the log of k(cat)/K(m) plots indicates that only one ionic interaction is involved in forming the transition state. These data support a model where two ionic interactions are formed between NADPH and symmetry related K32 residues in the ground state. To reach the transition state, an ionic interaction between K32 and the pyrophosphate bridge is broken. This unusual scenario likely arises from the constraints imposed by the 222 symmetry of the enzyme.     

Interloop contacts modulate ligand cycling during catalysis by Escherichia coli dihydrofolate reductase

As a continuation to our studies on the importance of interloop interactions in the Escherichia coli DHFR catalytic cycle, we have investigated the role of the betaG-betaH loop in modulating the closed and occluded conformations of the Met20 loop during the DHFR catalytic cycle. Specifically, to assess the importance of the hydrogen bond formed between Ser148 in the betaG-betaH loop and the Met20 loop, Ser148 was independently substituted with aspartic acid, alanine, and lysine. Moreover, the betaG-betaH loop was deleted entirely to yield the Delta(146-148) DHFR mutant. Steady-state turnover rates for all mutants were at most 3-fold lower than the wild-type rate. Lack of an isotope effect on this rate indicated the chemistry step does not contribute to the steady-state turnover. Consistent with this finding, hydride transfer rates for the DHFR mutants were at least 10-fold greater than the observed steady-state rates. The values ranged from a 30% decrease (Ser148Ala and Ser148Lys) to a 50% increase (Ser148Asp) in rate relative to that of the wild type. Modifications of the betaG-betaH loop enhanced the affinity for the cofactor and decreased the affinity for pterin, as determined by the K(D) values of the mutant proteins. Further analysis of Ser148Ala and Delta(146-148) DHFRs indicated these effects were manifest mainly in ligand off rates, although in some cases the on rate was affected. The Ser148Asp and Delta(146-148) mutations perturbed the preferred catalytic cycle through the introduction of branching at key intermediates. Rather than following the single WT pathway which involves loss of NADP(+) and rebinding of NADPH to precede loss of the product H4F (negative cooperativity), the mutants can reenter the catalytic cycle through different pathways. These findings suggest that the role of the interloop interaction between the betaG-betaH loop and the Met20 loop is to modulate ligand off rates allowing for proper cycling through the preferred kinetic pathway.     

Conformation coupled enzyme catalysis: single-molecule and transient kinetics investigation of dihydrofolate reductase

Ensemble kinetics and single-molecule fluorescence microscopy were used to study conformational transitions associated with enzyme catalysis by dihydrofolate reductase (DHFR). The active site loop of DHFR was labeled with a fluorescence quencher, QSY35, at amino acid position 17, and the fluorescent probe, Alexa555, at amino acid 37, by introducing cysteines at these sites with site-specific mutagenesis. The distance between the probes was such that approximately 50% fluorescence resonance energy transfer (FRET) occurred. The double-labeled enzyme retained essentially full catalytic activity, and stopped-flow studies of both the forward and reverse reactions revealed that the distance between probes increased prior to hydride transfer. A fluctuation in fluorescence intensity of single molecules of DHFR was observed in an equilibrium mixture of substrates but not in their absence. Ensemble rate constants were derived from the distributions of lifetimes observed and attributed to a reversible conformational change. Studies were carried out with both NADPH and NADPD as substrates, with no measurable isotope effect. Similar studies with a G121V mutant DHFR resulted in smaller rate constants. This mutant DHFR has reduced catalytic activity, so that the collective data for the conformational change suggest that the conformational change being observed is associated with catalysis and probably represents a conformational change prior to hydride transfer. If the change in fluorescence is attributed to a change in FRET, the distance change associated with the conformational change is approximately 1-2 A. These results are correlated with other measurements related to conformation coupled catalysis.     

Methotrexate, a high-affinity pseudosubstrate of dihydrofolate reductase.

Investigations have been made of the slow, tight-binding inhibition by methotrexate of the reaction catalyzed by dihydrofolate reductase from Streptococcus faecium A. Quantitative analysis has shown that progress curve data are in accord with a mechanism that involves the rapid formation of an enzyme-NADPH-methotrexate complex that subsequently undergoes a relatively slow, reversible isomerization reaction. From the Ki value for the dissociation of methotrexate from the E-NADPH-methotrexate complex (23 nM) and values of 5.1 and 0.013 min-1 for the forward and reverse rate constants of the isomerization reaction, the overall inhibition constant for methotrexate was calculated to be 58 pM. The formation of an enzyme-methotrexate complex was demonstrated by means of fluorescence quenching, and a value of 0.36 muM was determined for its dissociation constant. The same technique was used to determine dissociation constants for the reaction of methotrexate with the E-NADP and E-NADPH complexes. The results indicate that in the presence of either NADPH or NADP there is enhancement of the binding of methotrexate to the enzyme. It is proposed that methotrexate behaves as a pseudosubstrate for dihydrofolate reductase.