The coupling of structural fluctuations to hydride transfer in dihydrofolate reductase

The energy barrier for hydride transfer in wild-type G121V and G121S variants of Escherichia coli dihydrofolate reductase (DHFR) fluctuates in a time-dependent manner. This fluctuation may be attributed to structural changes in the protein that modulate the site of chemistry. Despite being far from the active site, mutations at position 121 of DHFR reduce the hydride transfer rate of the enzyme. This occurrence has been suggested to arise from modifications to the conformational ensemble of the protein. We elucidate the effects of the G121S and G121V mutations on the hydride transfer barrier by identifying structural changes in the protein that correlate with lowered barriers. The effect of these structural parameters on the hydride transfer barrier may be rationalized by simple considerations of the geometric constraints of the hydride transfer reaction. Fluctuations of these properties are associated with specific backbone dihedral angles of residues within the Methione-20 (M20) loop. The dihedral angle preferences are mediated by interactions with the region of the enzyme in the vicinity of residue 121 and are translated into distinct ligand conformations. We predict mutations within the M20 loop that may alter the conformational space explored by DHFR. Such mutational changes are anticipated to adjust the hydride transfer efficacy of DHFR by modifying equilibrium distributions of hydride transfer barriers found in the enzyme.     

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.     

The catalytic effect of dihydrofolate reductase and its mutants is determined by reorganization energies

The effect of distant mutations on the catalytic reaction of dihydrofolate reductase (DHFR) is reexamined by empirical valence bond simulations. The simulations reproduce for the first time the changes in the observed rate constants (without the use of adjustable parameters for this purpose) and show that the changes in activation barriers are strongly correlated with the corresponding changes in the reorganization energy. The preorganization of the polar groups of enzymes is the key catalytic factor, and anticatalytic mutations destroy this preorganization. Some anticatalytic mutations in DHFR also increase the distance between the donor and acceptor, but this effect is not directly related to catalysis since the native enzyme and the uncatalyzed reaction in water have similar average donor-acceptor distances. Insight into the effect of a mutation is provided by constructing the relevant free energy surfaces in terms of the generalized solute-solvent coordinates. It is shown how the mutations change the reaction coordinate and the activation barrier, and it is clarified that the corresponding changes do not reflect dynamical effects. It is also pointed out that all reactions in a condensed phase involve correlated motions (both in enzymes and in solution) and that the change of such motions upon mutations is a result of the change in the shape of the multidimensional reaction path on the solute-solvent surface, rather than the reason for the change in rate constant. Thus, as far as catalysis is concerned, the change in the activation barrier is due to the change in the electrostatic preorganization energy.     

Quantum-classical simulation methods for hydrogen transfer in enzymes: a case study of dihydrofolate reductase

A variety of theoretical approaches have been used to investigate hydrogen transfer in enzymatic reactions. The free energy barriers for hydrogen transfer in enzymes have been calculated using classical molecular dynamics simulations in conjunction with quantum mechanical/molecular mechanical and empirical valence bond potentials. Nuclear quantum effects have been included with several different approaches. Applications of these approaches to hydride transfer in dihydrofolate reductase are consistent with experimental measurements and provide significant insight into the protein conformational changes that facilitate the hydride transfer reaction.     

Dynamics of immobilized and native Escherichia coli dihydrofolate reductase by quasielastic neutron scattering

The internal dynamics of native and immobilized Escherichia coli dihydrofolate reductase (DHFR) have been examined using incoherent quasielastic neutron scattering. These results reveal no difference between the high frequency vibration mean-square displacement of the native and the immobilized E. coli DHFR. However, length-scale-dependent, picosecond dynamical changes are found. On longer length scales, the dynamics are comparable for both DHFR samples. On shorter length scales, the dynamics is dominated by local jump motions over potential barriers. The residence time for the protons to stay in a potential well is tau = 7.95 +/- 1.02 ps for the native DHFR and tau = 20.36 +/- 1.80 ps for the immobilized DHFR. The average height of the potential barrier to the local motions is increased in the immobilized DHFR, and may increase the activation energy for the activity reaction, decreasing the rate as observed experimentally. These results suggest that the local motions on the picosecond timescale may act as a lubricant for those associated with DHFR activity occurring on a slower millisecond timescale. Experiments indicate a significantly slower catalytic reaction rate for the immobilized E. coli DHFR. However, the immobilization of the DHFR is on the exterior of the enzyme and essentially distal to the active site, thus this phenomenon has broad implications for the action of drugs distal to the active site.     

Flexibility, diversity, and cooperativity: pillars of enzyme catalysis

This brief review discusses our current understanding of the molecular basis of enzyme catalysis. A historical development is presented, beginning with steady state kinetics and progressing through modern fast reaction methods, nuclear magnetic resonance, and single-molecule fluorescence techniques. Experimental results are summarized for ribonuclease, aspartate aminotransferase, and especially dihydrofolate reductase (DHFR). Multiple intermediates, multiple conformations, and cooperative conformational changes are shown to be an essential part of virtually all enzyme mechanisms. In the case of DHFR, theoretical investigations have provided detailed information about the movement of atoms within the enzyme-substrate complex as the reaction proceeds along the collective reaction coordinate for hydride transfer. A general mechanism is presented for enzyme catalysis that includes multiple intermediates and a complex, multidimensional standard free energy surface. Protein flexibility, diverse protein conformations, and cooperative conformational changes are important features of this model.     

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.     

The role of enzyme dynamics and tunnelling in catalysing hydride transfer: studies of distal mutants of dihydrofolate reductase

Residues M42 and G121 of Escherichia coli dihydrofolate reductase (ecDHFR) are on opposite sides of the catalytic centre (15 and 19 A away from it, respectively). Theoretical studies have suggested that these distal residues might be part of a dynamics network coupled to the reaction catalysed at the active site. The ecDHFR mutant G121V has been extensively studied and appeared to have a significant effect on rate, but only a mild effect on the nature of H-transfer. The present work examines the effect of M42W on the physical nature of the catalysed hydride transfer step. Intrinsic kinetic isotope effects (KIEs), their temperature dependence and activation parameters were studied. The findings presented here are in accordance with the environmentally coupled hydrogen tunnelling. In contrast to the wild-type (WT), fluctuations of the donor-acceptor distance were required, leading to a significant temperature dependence of KIEs and deflated intercepts. A comparison of M42W and G121V to the WT enzyme revealed that the reduced rates, the inflated primary KIEs and their temperature dependences resulted from an imperfect potential surface pre-arrangement relative to the WT enzyme. Apparently, the coupling of the enzyme's dynamics to the reaction coordinate was altered by the mutation, supporting the models in which dynamics of the whole protein is coupled to its catalysed chemistry.     

Impact of enzyme motion on activity

Experimental and theoretical data imply that enzyme motion plays an important role in enzymatic reactions. Enzyme motion can influence both the activation free energy barrier and the degree of barrier recrossing. A hybrid theoretical approach has been developed for the investigation of the relation between enzyme motion and activity. This approach includes both electronic and nuclear quantum effects. It distinguishes between thermally averaged promoting motions that influence the activation free energy barrier and dynamical motions that influence the barrier recrossings. Applications to hydride transfer in liver alcohol dehydrogenase and dihydrofolate reductase resulted in the identification and characterization of important enzyme motions. These applications have also led to the proposal of a network of coupled promoting motions in enzymatic reactions. These concepts have important implications for protein engineering and drug design.     

Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates

Recently, the participation in catalysis of residues spatially removed from the enzyme's active site has received considerable attention. The influence of the distal Gly-121 residue on the chemical step of hydride transfer in dihydrofolate reductase (DHFR) catalysis had been demonstrated previously [Cameron, C. E., and Benkovic, S. J. (1997) Biochemistry 36, 15792-15800]. In our continuing effort to identify catalytically important residues that are distal from the active site, we used sequence conservation information, kinetic data on site-directed mutants, dynamic motion information from NMR methods, and correlated motions from MD simulations to identify a subset of residues. Among them, the region spanning positions 41-45 is distal to the active site and was chosen as the focus for the mutagenesis and kinetic studies reported here. Specifically, the highly conserved Met-42 was selected for site-directed mutagenesis. While the reaction kinetics for the M42F mutant enzyme did not deviate from wild-type behavior, a 41-fold reduction in the forward hydride transfer rate was found for the M42W mutant. Given the established role of Gly-121 in the hydride transfer process, double mutant enzymes involving positions 42 and 121 were constructed and characterized. These double mutant enzymes generally showed little changes in substrate and cofactor binding but synergistic decreases in forward hydride transfer rates, while the decreases in reverse rates were additive. Along with supporting information from mixed quantum/classical MD simulations [Agarwal, P. K., Billeter, S. R., Rajagopalan, P. T., Benkovic, S. J., and Hammes-Schiffer, S. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 2794-2799], the data suggest that a coupled dynamic motion of these distal residues enhances crossing of the chemical reaction barrier and imply an expanded nonstatic role for the protein fold in catalysis.     

Role of aspartate 27 of dihydrofolate reductase from Escherichia coli in interconversion of active and inactive enzyme conformers and binding of NADPH

The apoenzyme of wild-type (WT) dihydrofolate reductase (DHRF) from Escherichia coli exists in two conformational states, Et and Ew, which differ in affinity for NADPH and in kinetic competence. Dissociation constants for the binary complex of NADPH with the two conformers differ by over 100-fold (KDt = 0.17 microM, KDw = 22 microM). Rate constants governing the interconversion of conformers are small (t1/2 for Ew----Et = 71 s), and since Ew is not catalytically competent, this conversion is accompanied by an increase in catalytic velocity. The equilibrium proportion of Et in the absence of ligands is 63%, but binding of NADPH greatly increases this proportion, and t1/2 for conversion of Ew.NADPH to Et.NADPH is 30 s. This conformational equilibrium has also been examined in mutant enzyme in which aspartate 27 is replaced by asparagine (D27N E. coli DHFR). Although ASp27 is an active site residue, it does not interact directly with bound NADPH, and in the mutant the rate constant for NADPH binding to Et is unchanged as are the dissociation constants for NADPH complexes with Et or Ew. However, for mutant apoenzyme, the proportion of Et is decreased to 18% in the absence of ligands so that the overall KD for NADPH is increased (0.15 microM for WT E. coli DHFR, 0.68 microM for D27N E. coli DHFR). The lower proportion of Et is due to a decreased rate for Ew----Et (t1/2 = 221 s) and an increased rate for Et----Ew (t1/2 = 50 s versus 120 s for WT E. coli DHFR).     

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.     

Mutagenesis Alters the Catalytic Mechanism of the Light-driven Enzyme Protochlorophyllide Oxidoreductase

The light-activated enzyme protochlorophyllide oxidoreductase (POR) catalyzes an essential step in the synthesis of the most abundant pigment on Earth, chlorophyll. This unique reaction involves the sequential addition of a hydride and proton across the C17=C18 double bond of protochlorophyllide (Pchlide) by dynamically coupled quantum tunneling and is an important model system for studying the mechanism of hydrogen transfer reactions. In the present work, we have combined site-directed mutagenesis studies with a variety of sensitive spectroscopic and kinetic measurements to provide new insights into the mechanistic role of three universally conserved Cys residues in POR. We show that mutation of Cys-226 dramatically alters the catalytic mechanism of the enzyme. In contrast to wild-type POR, the characteristic charge-transfer intermediate, formed upon hydride transfer from NADPH to the C17 position of Pchlide, is absent in C226S variant enzymes. This suggests a concerted hydrogen transfer mechanism where proton transfer only is rate-limiting. Moreover, Pchlide reduction does not require the network of solvent-coupled conformational changes that play a key role in the proton transfer step of wild-type POR. We conclude that this globally important enzyme is finely tuned to facilitate efficient photochemistry, and the removal of a key interaction with Pchlide in the C226S variants significantly affects the local active site structure in POR, resulting in a shorter donor-acceptor distance for proton transfer.     

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.     

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase.

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

Preservation of Protein Dynamics in Dihydrofolate Reductase Evolution

The hydride transfer reaction catalyzed by dihydrofolate reductase (DHFR) is a model for examining how protein dynamics contribute to enzymatic function. The relationship between functional motions and enzyme evolution has attracted significant attention. Recent studies on N23PP Escherichia coli DHFR (ecDHFR) mutant, designed to resemble parts of the human enzyme, indicated a reduced single turnover rate. NMR relaxation dispersion experiments with that enzyme showed rigidification of millisecond Met-20 loop motions (Bhabha, G., Lee, J., Ekiert, D. C., Gam, J., Wilson, I. A., Dyson, H. J., Benkovic, S. J., and Wright, P. E. (2011) Science 332, 234–238). A more recent study of this mutant, however, indicated that fast motions along the reaction coordinate are actually more dispersed than for wild-type ecDHFR (WT). Furthermore, a double mutant (N23PP/G51PEKN) that better mimics the human enzyme seems to restore both the single turnover rates and narrow distribution of fast dynamics (Liu, C. T., Hanoian, P., French, T. H., Hammes-Schiffer, S., and Benkovic, S. J. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 10159–11064). Here, we measured intrinsic kinetic isotope effects for both N23PP and N23PP/G51PEKN double mutant DHFRs over a temperature range. The findings indicate that although the C-H→C transfer and dynamics along the reaction coordinate are impaired in the altered N23PP mutant, both seem to be restored in the N23PP/G51PEKN double mutant. This indicates that the evolution of G51PEKN, although remote from the Met-20 loop, alleviated the loop rigidification that would have been caused by N23PP, enabling WT-like H-tunneling. The correlation between the calculated dynamics, the nature of C-H→C transfer, and a phylogenetic analysis of DHFR sequences are consistent with evolutionary preservation of the protein dynamics to enable H-tunneling from well reorganized active sites.     

An internal equilibrium preorganizes the enzyme-substrate complex for hydride tunneling in choline oxidase

The hydride transfer reaction catalyzed by choline oxidase under irreversible regime, i.e., at saturating oxygen, was shown in a recent study to occur quantum mechanically within a highly preorganized active site, with the reactive configuration for hydride tunneling being minimally affected by environmental vibrations of the reaction coordinate other than those affecting the distance between the alpha-carbon of the choline alkoxide substrate and the N(5) atom of the enzyme-bound flavin cofactor [Fan, F., and Gadda, G. (2005) J. Am. Chem. Soc. 127, 17954-17961]. In this study, we have determined the effects of pH and temperature on the substrate kinetic isotope effects with 1,2-[2H4]choline as substrate for choline oxidase at 0.2 mM oxygen to gain insights on the mechanism of hydride transfer under reversible catalytic regime. The data presented indicated that the kinetic complexity arising from the net flux through the reverse of the hydride transfer step changed with temperature, with the hydride transfer reaction becoming more reversible with increasing temperatures. After this kinetic complexity was accounted for, analyses of the kcat/Km and D(kcat/Km) values determined at 0.2 mM according to the Eyring and Arrhenius formalisms suggested that the quantum mechanical nature of the hydride transfer reaction is, not surprisingly, maintained during enzymatic catalysis under reversible regime. A comparison of the thermodynamic and kinetic parameters of the hydride transfer reaction under reversible and irreversible catalytic regimes showed that the enthalpies of activation (DeltaH++) were significantly larger in the reversible catalytic regime. This reflects the presence of an enthalpically unfavorable internal equilibrium of the enzyme-substrate Michaelis complex occurring prior to, and independently from, CH bond cleavage. Such an internal equilibrium is required to preorganize the enzyme-substrate complex for efficient quantum mechanical tunneling of the hydride ion from the substrate alpha-carbon to the flavin N(5) atom.     

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.     

A quantum mechanics/molecular mechanics study of the catalytic mechanism of the thymidylate synthase

A theoretical study of the molecular mechanism of the thymidylate synthase-catalyzed reaction has been carried out using hybrid quantum mechanics/molecular mechanics methods. We have examined all of the stationary points (reactants, intermediates, transition structures, and products) on the multidimensional potential energy surfaces for the multistep enzymatic process. The characterization of these relevant structures facilitates the gaining of insight into the role of the different residues in the active site. Furthermore, analysis of the full energy profile has revealed that the step corresponding to the reduction of the exocyclic methylene intermediate by hydride transfer from the 6S position of 5,6,7,8-tetrahydrofolate (H4folate), forming dTMP and 7,8-dihydrofolate (H2folate), is the rate-limiting step, in accordance with the experimental data. In this step, the hydride transfer and the scission of an overall conserved active site cysteine residue (Cys146 in Escherichia coli) take place in a concerted but very asynchronous way. These findings have also been tested with primary and secondary deuterium, tritium, and sulfur kinetic isotope effects, and the calculations have been compared to experimental data. Finally, the incorporation of high-level quantum mechanical corrections to the semiempirical AM1 Hamiltonian into our hybrid scheme has allowed us to obtain reasonable values of the energy barrier for the rate-limiting step. The resulting picture of the complete multistep enzyme mechanism that is obtained reveals several new features of substantial mechanistic interest.     

Biosensors based on binding-modulated donor-acceptor distances

The promising recognition characteristics exhibited by biomolecules have caused significant interest in biomolecule-based sensor strategies. Here we review several emerging biosensor designs that use modulated electron or energy transfer to a bio-specific ligand as the signaling mechanism. The efficiencies of both electron transfer and energy transfer are strongly dependent on donor-acceptor distance. When coupled with the large conformational changes sometimes associated with biomolecular recognition, these distance-dependent processes provide a robust means for generating optical and electronic signals.