Dihydrofolate Reductase Catalysis Research Articles

Single-Molecule Sampling of Dihydrofolate Reductase Shows Kinetic Pauses and an Endosteric Effect Linked to Catalysis.

The ability to sample multiple reactions on the same single enzyme is important to link rare intermediates with catalysis and to unravel the role of conformational changes. Despite decades of efforts, however, the single-molecule characterization of nonfluorogenic enzymes during multiple catalytic turnovers has been elusive. Here, we show that nanopore currents allow sampling the dynamic exchange between five structural intermediates during E. coli dihydrofolate reductase (DHFR) catalysis. We found that an endosteric effect promotes the binding of the substrate to the enzyme with a specific hierarchy. The chemical step then switched the enzyme from the closed to the occluded conformation, which in turn promotes the release of the reduced cofactor NADP+. Unexpectedly, only a few reactive complexes lead to catalysis. Furthermore, second-long catalytic pauses were observed, possibly reflecting an off-path conformation generated during the reaction. Finally, the free energy from multiple cofactor binding events were required to release the product and switch DHFR back to the reactive conformer. This catalytic fueled concerted mechanism is likely to have evolved to improve the catalytic efficiency of DHFR under the high concentrations of NADP+ in E. coli and might be a general feature for complex enzymatic reactions where the binding and release of the products must be tightly controlled.

Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

Electric Field Measurements Reveal the Pivotal Role of Cofactor-Substrate Interaction in Dihydrofolate Reductase Catalysis.

The contribution of ligand–ligand electrostatic interaction to transition state formation during enzyme catalysis has remained unexplored, even though electrostatic forces are known to play a major role in protein functions and have been investigated by the vibrational Stark effect (VSE). To monitor electrostatic changes along important steps during catalysis, we used a nitrile probe (T46C-CN) inserted proximal to the reaction center of three dihydrofolate reductases (DHFRs) with different biophysical properties, Escherichia coli DHFR (EcDHFR), its conformationally impaired variant (EcDHFR-S148P), and Geobacillus stearothermophilus DHFR (BsDHFR). Our combined experimental and computational approach revealed that the electric field projected by the substrate toward the probe negates those exerted by the cofactor when both are bound within the enzymes. This indicates that compared to previous models that focus exclusively on subdomain reorganization and protein–ligand contacts, ligand–ligand interactions are the key driving force to generate electrostatic environments conducive for catalysis.

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H2F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H2F in driving hydride transfer remains unclear. H2F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM simulations reveal that hydride transfer is triggered by polarization at the C6 of H2F. The σ Cβ–H bonds contribute to the buildup of the cationic character during the chemical transformation, and hyperconjugation influences the formation of the transition state. Our findings provide key insights into the hydride transfer mechanism of the DHFR-catalyzed reaction, which is a target for antiproliferative drugs and a paradigmatic model in mechanistic enzymology.

A Rapid Analysis of Variations in Conformational Behavior during Dihydrofolate Reductase Catalysis.

Protein flexibility is central to enzyme catalysis, yet it remains challenging both to predict conformational behavior on the basis of analysis of amino acid sequence and protein structure and to provide the necessary breadth of experimental support to any such predictions. Here a generic and rapid procedure for identifying conformational changes during dihydrofolate reductase (DHFR) catalysis is described. Using DHFR from Escherichia coli (EcDHFR), selective side-chain 13C labeling of methionine and tryptophan residues is shown to be sufficient to detect the closed-to-occluded conformational transition that follows the chemical step in the catalytic cycle, with clear chemical shift perturbations found for both methionine methyl and tryptophan indole groups. In contrast, no such perturbations are seen for the DHFR from the psychrophile Moritella profunda, where the equivalent conformational change is absent. Like EcDHFR, Salmonella enterica DHFR shows experimental evidence of a large-scale conformational change following hydride transfer that relies on conservation of a key hydrogen bonding interaction between the M20 and GH loops, directly comparable to the closed-to-occluded conformational change observed in EcDHFR. For the hyperthermophile Thermotoga maritima, no chemical shift perturbations were observed, suggesting that no major conformational change occurs during the catalytic cycle. In spite of their conserved tertiary structures, DHFRs display variations in conformational sampling that occurs concurrently with catalysis.

Dihydrofolate Reductase as a Model for Enzyme Catalysis

Dihydrofolate reductase (DHFR) has emerged as a model enzyme to investigate the role of dynamics in catalysis. This review begins with a description of the timing of the protonation of the dihydrofolate substrate and the subsequent hydride transfer reaction. This is followed by a summary of the results obtained from bioinformatics, theoretical and experimental studies of coevolving DHFR residues. Some of these residues appear to be dynamically coupled to each other and to the reaction coordinate. Importantly, the coupling appears to expand across the entire enzyme as many of the coevolved amino acids are remote from the active site. The importance of introducing mutations during evolution in the correct evolutionary order is also discussed. These studies involved a combination of steady state and pre-steady state rate measurements, kinetic isotope effects, NMR spectroscopy and computational investigations of the wild type DHFR from E. coli and several mutants. The cumulative results implicate a role of motions across the entire protein in DHFR catalysis, which may serve as a guide in studies of other enzymatic systems. Keywords: Dihydrofolate reductase, enzyme dynamics, enzyme evolution, global network, hydride transfer, kinetic isotope effect.

Network of remote and local protein dynamics in dihydrofolate reductase catalysis.

Molecular dynamics calculations and bionformatic studies of dihydrofolate reductase (DHFR) have suggested a network of coupled motions across the whole protein that is correlated to the reaction coordinate. Experimental studies demonstrated that distal residues G121, M42 and F125 in E. coli DHFR participate in that network. The missing link in our understanding of DHFR catalysis is the lack of a mechanism by which such remote residues can affect the catalyzed chemistry at the active site. Here, we present a study of the temperature dependence of intrinsic kinetic isotope effects (KIEs) that indicates synergism between a remote residue in that dynamic network, G121, and the active site's residue I14. The intrinsic KIEs for the I14A-G121V double mutant showed steeper temperature dependence (ΔEa(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔEa(T-H)G121V +ΔEa(T-H) I14A < ΔEa(T-H) double mutant, which indicates a synergism between the two residues. This finding links the remote residues in the network under investigation to the enzyme's active site, providing a mechanism by which these residues can be coupled to the catalyzed chemistry. This experimental evidence validates calculations proposing that both remote and active site residues constitute a network of coupled promoting motions correlated to the bond activation step (C-H→C hydride transfer in this case). Additionally, the effect of I14A and G121V mutations on single turnover rates was additive rather than synergistic. Although single turnover rate measurements are more readily available and thus more popular than assessing intrinsic kinetic isotope effects, the current finding demonstrates that for these rates, which in DHFR reflect several microscopic rate constants, can fall short of revealing the nature of the C-H bond activation per se.

Role of the occluded conformation in bacterial dihydrofolate reductases.

Dihydrofolate reductase (DHFR) from Escherichia coli (EcDHFR) adopts two major conformations, closed and occluded, and movement between these two conformations is important for progression through the catalytic cycle. DHFR from the cold-adapted organism Moritella profunda (MpDHFR) on the other hand is unable to form the two hydrogen bonds that stabilize the occluded conformation in EcDHFR and so remains in a closed conformation during catalysis. EcDHFR-S148P and MpDHFR-P150S were examined to explore the influence of the occluded conformation on catalysis by DHFR. Destabilization of the occluded conformation did not affect hydride transfer but altered the affinity for the oxidized form of nicotinamide adenine dinucleotide phosphate (NADP(+)) and changed the rate-determining step of the catalytic cycle for EcDHFR-S148P. Even in the absence of an occluded conformation, MpDHFR follows a kinetic pathway similar to that of EcDHFR with product release being the rate-limiting step in the steady state at pH 7, suggesting that MpDHFR uses a different strategy to modify its affinity for NADP(+). DHFRs from many organisms lack a hydrogen bond donor in the appropriate position and hence most likely do not form an occluded conformation. The link between conformational cycling between closed and occluded forms and progression through the catalytic cycle is specific to EcDHFR and not a general characteristic of prokaryotic DHFR catalysis.

First Three-Dimensional Structure of Toxoplasma gondii Thymidylate Synthase–Dihydrofolate Reductase: Insights for Catalysis, Interdomain Interactions, and Substrate Channeling

Most species, such as humans, have monofunctional forms of thymidylate synthase (TS) and dihydrofolate reductase (DHFR) that are key folate metabolism enzymes making critical folate components required for DNA synthesis. In contrast, several parasitic protozoa, including Toxoplasma gondii , contain a unique bifunctional thymidylate synthase-dihydrofolate reductase (TS-DHFR) having the catalytic activities contained on a single polypeptide chain. The prevalence of T. gondii infections across the world, especially for those immunocompromised, underscores the need to understand TS-DHFR enzyme function and to find new avenues to exploit for the design of novel antiparasitic drugs. As a first step, we have solved the first three-dimensional structures of T. gondii TS-DHFR at 3.7 Å and of a loop truncated TS-DHFR, removing several flexible surface loops in the DHFR domain, improving resolution to 2.2 Å. Distinct structural features of the TS-DHFR homodimer include a junctional region containing a kinked crossover helix between the DHFR domains of the two adjacent monomers, a long linker connecting the TS and DHFR domains, and a DHFR domain that is positively charged. The roles of these unique structural features were probed by site-directed mutagenesis coupled with presteady state and steady state kinetics. Mutational analysis of the crossover helix region combined with kinetic characterization established the importance of this region not only in DHFR catalysis but also in modulating the distal TS activity, suggesting a role for TS-DHFR interdomain interactions. Additional kinetic studies revealed that substrate channeling occurs in which dihydrofolate is directly transferred from the TS to DHFR active site without entering bulk solution. The crystal structure suggests that the positively charged DHFR domain governs this electrostatically mediated movement of dihydrofolate, preventing release from the enzyme. Taken together, these structural and kinetic studies reveal unique, functional regions on the T. gondii TS-DHFR enzyme that may be targeted for inhibition, thus paving the way for designing species specific inhibitors.

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis

Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy ((15)N, (13)C, (2)H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy enzyme reflects differences in environmental coupling to the hydride transfer step. Importantly, the barrier and contribution of quantum tunneling are not affected, indicating no significant role for "promoting motions" in driving tunneling or modulating the barrier. The chemical step is slower in the heavy enzyme because protein motions coupled to the reaction coordinate are slower. The fact that the heavy enzyme is only slightly less active than its light counterpart shows that protein dynamics have a small, but measurable, effect on the chemical reaction rate.

Temporally Overlapped but Uncoupled Motions in Dihydrofolate Reductase Catalysis

Temporal correlations between protein motions and enzymatic reactions are often interpreted as evidence for catalytically important motions. Using dihydrofolate reductase as a model system, we show that there are many protein motions that temporally overlapped with the chemical reaction, and yet they do not exhibit the same kinetic behaviors (KIE and pH dependence) as the catalyzed chemical reaction. Thus, despite the temporal correlation, these motions are not directly coupled to the chemical transformation, and they might represent a different part of the catalytic cycle or simply be the product of the intrinsic flexibility of the protein.

Effect of Dimerization on Dihydrofolate Reductase Catalysis

Dihydrofolate reductase (DHFR) from the hyperthermophile Thermotoga maritima (TmDHFR) forms a very stable homodimer, while DHFRs from other organisms are monomers. We investigated the effect of dimerization on DHFR catalysis by preparing a dimeric variant, Xet-3, of DHFR from Escherichia coli (EcDHFR). Introducing residues located at the TmDHFR dimer interface into EcDHFR increases the melting temperature to ∼60 °C, approximately 9 °C higher than that measured for EcDHFR. The steady-state and pre-steady-state rate constants measured for Xet-3 were similar to those of dimeric TmDHFR but significantly lower than those of the parent EcDHFR. This reduction in the degree of catalytic competence is likely a consequence of the loss of flexibility of catalytically important loop regions of EcDHFR on dimerization and a compromise of the electrostatic environment of the active site. In contrast, the reduced catalytic ability of TmDHFR relative to that of EcDHFR is not simply a consequence of reduced loop flexibility in the dimeric enzyme. Our studies demonstrate that EcDHFR is not a good model for understanding the properties of other DHFRs, including TmDHFR.

Evidence that a ‘dynamic knockout’ in Escherichia coli dihydrofolate reductase does not affect the chemical step of catalysis

The question of whether protein motions play a role in the chemical step of enzymatic catalysis has generated much controversy in recent years. Debate has recently reignited over possible dynamic contributions to catalysis in dihydrofolate reductase, following conflicting conclusions from studies of the N23PP/S148A variant of the Escherichia coli enzyme. By investigating the temperature dependence of kinetic isotope effects, we present evidence that the reduction in the hydride transfer rate constants in this variant is not a direct result of impairment of conformational fluctuations. Instead, the conformational state of the enzyme immediately before hydride transfer, which determines the electrostatic environment of the active site, affects the rate constant for the reaction. Although protein motions are clearly important for binding and release of substrates and products, there appears to be no detectable dynamic coupling of protein motions to the hydride transfer step itself.

Catalysis by dihydrofolate reductase and other enzymes arises from electrostatic preorganization, not conformational motions

The proposal that enzymatic catalysis is due to conformational fluctuations has been previously promoted by means of indirect considerations. However, recent works have focused on cases where the relevant motions have components toward distinct conformational regions, whose population could be manipulated by mutations. In particular, a recent work has claimed to provide direct experimental evidence for a dynamical contribution to catalysis in dihydrofolate reductase, where blocking a relevant conformational coordinate was related to the suppression of the motion toward the occluded conformation. The present work utilizes computer simulations to elucidate the true molecular basis for the experimentally observed effect. We start by reproducing the trend in the measured change in catalysis upon mutations (which was assumed to arise as a result of a "dynamical knockout" caused by the mutations). This analysis is performed by calculating the change in the corresponding activation barriers without the need to invoke dynamical effects. We then generate the catalytic landscape of the enzyme and demonstrate that motions in the conformational space do not help drive catalysis. We also discuss the role of flexibility and conformational dynamics in catalysis, once again demonstrating that their role is negligible and that the largest contribution to catalysis arises from electrostatic preorganization. Finally, we point out that the changes in the reaction potential surface modify the reorganization free energy (which includes entropic effects), and such changes in the surface also alter the corresponding motion. However, this motion is never the reason for catalysis, but rather simply a reflection of the shape of the reaction potential surface.

The Tail Wagging the Dog: Insights into Catalysis in R67 Dihydrofolate Reductase

Plasmid-encoded R67 dihydrofolate reductase (DHFR) catalyzes a hydride transfer reaction between substrate dihydrofolate (DHF) and its cofactor, nicotinamide adenine dinucleotide phosphate (NADPH). R67 DHFR is a homotetramer that exhibits numerous characteristics of a primitive enzyme, including promiscuity in binding of substrate and cofactor, formation of nonproductive complexes, and the absence of a conserved acid in its active site. Furthermore, R67's active site is a pore, which is mostly accessible by bulk solvent. This study uses a computational approach to characterize the mechanism of hydride transfer. Not surprisingly, NADPH remains fixed in one-half of the active site pore using numerous interactions with R67. Also, stacking between the nicotinamide ring of the cofactor and the pteridine ring of the substrate, DHF, at the hourglass center of the pore, holds the reactants in place. However, large movements of the p-aminobenzoylglutamate tail of DHF occur in the other half of the pore because of ion pair switching between symmetry-related K32 residues from two subunits. This computational result is supported by experimental results that the loss of these ion pair interactions (located >13 Å from the center of the pore) by addition of salt or in asymmetric K32M mutants leads to altered enzyme kinetics [Hicks, S. N., et al. (2003) Biochemistry 42, 10569-10578; Hicks, S. N., et al. (2004) J. Biol. Chem. 279, 46995-47002]. The tail movement at the edge of the active site, coupled with the fixed position of the pteridine ring in the center of the pore, leads to puckering of the pteridine ring and promotes formation of the transition state. Flexibility coupled to R67 function is unusual as it contrasts with the paradigm that enzymes use increased rigidity to facilitate attainment of their transition states. A comparison with chromosomal DHFR indicates a number of similarities, including puckering of the nicotinamide ring and changes in the DHF tail angle, accomplished by different elements of the dissimilar protein folds.

Disruption of the crossover helix impairs dihydrofolate reductase activity in the bifunctional enzyme TS–DHFR from Cryptosporidium hominis

In contrast with most species, including humans, which have monofunctional forms of the folate biosynthetic enzymes TS (thymidylate synthase) and DHFR (dihydrofolate reductase), several pathogenic protozoal parasites, including Cryptosporidium hominis, contain a bifunctional form of the enzymes on a single polypeptide chain having both catalytic activities. The crystal structure of the bifunctional enzyme TS-DHFR C. hominis reveals a dimer with a 'crossover helix', a swap domain between DHFR domains, unique in that this helical region from one monomer makes extensive contacts with the DHFR active site of the other monomer. In the present study, we used site-directed mutagenesis to probe the role of this crossover helix in DHFR catalysis. Mutations were made to the crossover helix: an 'alanine-face' enzyme in which the residues on the face of the helix close to the DHFR active site of the other subunit were mutated to alanine, a 'glycine-face' enzyme in which the same residues were mutated to glycine, and an 'all-alanine' helix in which all residues of the helix were mutated to alanine. These mutant enzymes were studied using a rapid transient kinetic approach. The mutations caused a dramatic decrease in the DHFR activity. The DHFR catalytic activity of the alanine-face mutant enzyme was 30 s(-1), the glycine-face mutant enzyme was 17 s(-1), and the all-alanine helix enzyme was 16 s(-1), all substantially impaired from the wild-type DHFR activity of 152 s(-1). It is clear that loss of helix interactions results in a marked decrease in DHFR activity, supporting a role for this swap domain in DHFR catalysis. The crossover helix provides a unique structural feature of C. hominis bifunctional TS-DHFR that could be exploited as a target for species-specific non-active site inhibitors.

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.

Origin of the Temperature Dependence of Isotope Effects in Enzymatic Reactions: The Case of Dihydrofolate Reductase

The origin of the temperature dependence of kinetic isotope effects (KIEs) in enzyme reactions is a problem of general interest and a major challenge for computational chemistry. The present work simulates the nuclear quantum mechanical (NQM) effects and the corresponding KIE in dihydrofolate reductase (DHFR) and two of its mutants by using the empirical valence bond (EVB) and the quantum classical path (QCP) centroid path integral approach. Our simulations reproduce the overall observed trend while using a fully microscopic rather than a phenomenological picture and provide an interesting insight. It appears that the KIE increases when the distance between the donor and acceptor increases, in a somewhat counter intuitive way. The temperature dependence of the KIE appears to reflect mainly the temperature dependence of the distance between the donor and acceptor. This trend is also obtained from a simplified vibronic treatment, but as demonstrated here, the vibronic treatment is not valid at short and medium distances, where it is essential to use the path integral or other approaches capable of moving seamlessly from the adiabatic to the diabatic limits. It is pointed out that although the NQM effects do not contribute to catalysis in DHFR, the observed temperature dependence can be used to refine the potential of mean force for the donor and acceptor distance and its change due to distanced mutations.

Ligand Binding and Circular Permutation Modify Residue Interaction Network in DHFR

Residue interaction networks and loop motions are important for catalysis in dihydrofolate reductase (DHFR). Here, we investigate the effects of ligand binding and chain connectivity on network communication in DHFR. We carry out systematic network analysis and molecular dynamics simulations of the native DHFR and 19 of its circularly permuted variants by breaking the chain connections in ten folding element regions and in nine nonfolding element regions as observed by experiment. Our studies suggest that chain cleavage in folding element areas may deactivate DHFR due to large perturbations in the network properties near the active site. The protein active site is near or coincides with residues through which the shortest paths in the residue interaction network tend to go. Further, our network analysis reveals that ligand binding has “network-bridging effects” on the DHFR structure. Our results suggest that ligand binding leads to a modification, with most of the interaction networks now passing through the cofactor, shortening the average shortest path. Ligand binding at the active site has profound effects on the network centrality, especially the closeness.

The Dynamic Energy Landscape of Dihydrofolate Reductase Catalysis

We used nuclear magnetic resonance relaxation dispersion to characterize higher energy conformational substates of Escherichia coli dihydrofolate reductase. Each intermediate in the catalytic cycle samples low-lying excited states whose conformations resemble the ground-state structures of preceding and following intermediates. Substrate and cofactor exchange occurs through these excited substates. The maximum hydride transfer and steady-state turnover rates are governed by the dynamics of transitions between ground and excited states of the intermediates. Thus, the modulation of the energy landscape by the bound ligands funnels the enzyme through its reaction cycle along a preferred kinetic path.