Tuning of the H-transfer coordinate in primitive versus well-evolved enzymes

Solvent Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects: Correlation to the Donor-Acceptor Distances

Protein structural effects on the temperature (T) dependence of kinetic isotope effects (KIEs) in H-tunneling reactions have recently been used to discuss about the role of enzyme thermal motions in catalysis. Frequently observed nearly T-independent KIEs in the wild-type enzymes and T-dependent KIEs in variants suggest that H-tunneling in the former is assisted by the naturally evolved protein constructive vibrations that help sample short donor-acceptor distances (DADs) needed. This explanation that correlates the T-dependence of KIEs with DAD sampling has been highly debated as simulations following other H-tunneling models sometimes gave alternative explanations. In this paper, solvent effects on the T-dependence of KIEs of two hydride tunneling reactions of NADH/NAD⁺ analogues (represented by ΔE_a = E_aD - E_aH) were determined in attempts to replicate the observations in enzymes and test the protein vibration-assisted DAD sampling concept. Effects of selected aprotic solvents on the DAD_PRC's of the productive reactant complexes (PRCs) and the DAD_TRS's of the activated tunneling ready states (TRSs) were obtained through computations and analyses of the kinetic data, including 2° KIEs, respectively. A weaker T-dependence of KIEs (i.e., smaller ΔE_a) was found in a more polar aprotic solvent in which the system has a shorter average DAD_PRC and DAD_TRS. Further results show that a charge-transfer (CT) complexation made of a stronger donor/acceptor gives rise to a smaller ΔE_a. Overall, the shorter and less broadly distributed DADs resulting from the stronger CT complexation vibrations give rise to a smaller ΔE_a. Our results appear to support the explanation that links the T-dependence of KIEs to the donor-acceptor rigidity in enzymes.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Tales of Dihydrofolate Binding to R67 Dihydrofolate Reductase

Homotetrameric R67 dihydrofolate reductase possesses 222 symmetry and a single active site pore. This situation results in a promiscuous binding site that accommodates either the substrate, dihydrofolate (DHF), or the cofactor, NADPH. NADPH interacts more directly with the protein as it is larger than the substrate. In contrast, the p-aminobenzoyl-glutamate tail of DHF, as monitored by nuclear magnetic resonance and crystallography, is disordered when bound. To explore whether smaller active site volumes (which should decrease the level of tail disorder by confinement effects) alter steady state rates, asymmetric mutations that decreased the half-pore volume by ∼35% were constructed. Only minor effects on k(cat) were observed. To continue exploring the role of tail disorder in catalysis, 1-ethyl-3-[3-(dimethylamino)propyl]carbodiimide-mediated cross-linking between R67 DHFR and folate was performed. A two-folate, one-tetramer complex results in the loss of enzyme activity where two symmetry-related K32 residues in the protein are cross-linked to the carboxylates of two bound folates. The tethered folate could be reduced, although with a ≤30-fold decreased rate, suggesting decreased dynamics and/or suboptimal positioning of the cross-linked folate for catalysis. Computer simulations that restrain the dihydrofolate tail near K32 indicate that cross-linking still allows movement of the p-aminobenzoyl ring, which allows the reaction to occur. Finally, a bis-ethylene-diamine-α,γ-amide folate adduct was synthesized; both negatively charged carboxylates in the glutamate tail were replaced with positively charged amines. The K(i) for this adduct was ∼9-fold higher than for folate. These various results indicate a balance between folate tail disorder, which helps the enzyme bind substrate while dynamics facilitates catalysis.

Role of long-range protein dynamics in different thymidylate synthase catalyzed reactions

Recent studies of Escherichia coli thymidylate synthase (ecTSase) showed that a highly conserved residue, Y209, that is located 8 Å away from the reaction site, plays a key role in the protein’s dynamics. Those crystallographic studies indicated that Y209W mutant is a structurally identical but dynamically altered relative to the wild type (WT) enzyme, and that its turnover catalytic rate governed by a slow hydride-transfer has been affected. The most challenging test of an examination of a fast chemical conversion that precedes the rate-limiting step has been achieved here. The physical nature of both fast and slow C-H bond activations have been compared between the WT and mutant by means of observed and intrinsic kinetic isotope effects (KIEs) and their temperature dependence. The findings indicate that the proton abstraction step has not been altered as much as the hydride transfer step. Additionally, the comparison indicated that other kinetic steps in the TSase catalyzed reaction were substantially affected, including the order of the substrate binding. Enigmatically, although Y209 is H-bonded to 3'-OH of 2'-deoxyuridine-5'-mono­phosphate (dUMP), its altered dynamics is more pronounced on the binding of the remote cofactor, (6R)-N⁵,N¹⁰-methylene-5,6,7,8-tetrahydrofolate (CH₂H₄folate), revealing the importance of long-range dynamics of the enzymatic complex and its catalytic function.

Effects of cavities at the nicotinamide binding site of liver alcohol dehydrogenase on structure, dynamics and catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.

Protein dynamics and enzymatic chemical barrier passage

After many decades of investigation, the manner in which enzymes increase the rate of chemical reactions, at times by a factor of 10(17) compared to the rate of the corresponding solution phase reaction, is still opaque. A topic of significant discussion in the literature of the past 5-10 years has been the importance of protein dynamics in this process. This Feature Article will discuss the authors' work on this still controversial topic with focus on both methodology and application to real systems. The end conclusion of this work has been that for specific enzymes under study protein dynamics on both rapid time scales of barrier crossing (termed promoting vibrations by the authors) and of conformational fluctuations are central to the function of biological catalysts. In another enzyme we will discuss, the results are far less clear. The manner of the coupling of chemistry to protein dynamics has deep implications for protein architecture, both natural and created, and recent results reinforce the complexity of the protein form that has evolved to support these functions.

Triple isotopic labeling and kinetic isotope effects: exposing H-transfer steps in enzymatic systems

Kinetic isotope effect (KIE) studies can provide insight into the mechanism and kinetics of specific chemical steps in complex catalytic cascades. Recent results from hydrogen KIE measurements have examined correlations between enzyme dynamics and catalytic function, leading to a surge of studies in this area. Unfortunately, most enzymatic H-transfer reactions are not rate limiting, and the observed KIEs do not reliably reflect the intrinsic KIEs on the chemical step of interest. Given their importance to understanding the chemical step under study, accurate determination of the intrinsic KIE from the observed data is essential. In 1975, Northrop developed an elegant method to assess intrinsic KIEs from their observed values [Northrop, D. B. (1975) Steady-state analysis of kinetic isotope effects in enzymic reactions. Biochemistry 14, 2644-2651]. The Northrop method involves KIE measurements using all three hydrogen isotopes, where one of them serves as the reference isotope. This method has been successfully used with different combinations of observed KIEs over the years, but criteria for a rational choice of reference isotope have never before been experimentally determined. Here we compare different reference isotopes (and hence distinct experimental designs) using the reduction of dihydrofolate and dihydrobiopterin by two dissimilar enzymes as model reactions. A number of isotopic labeling patterns have been applied to facilitate the comparative study of reference isotopes. The results demonstrate the versatility of the Northrop method and that such experiments are limited only by synthetic techniques, availability of starting materials, and the experimental error associated with the use of distinct combinations of isotopologues.

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

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

A twin-track approach has optimized proton and hydride transfer by dynamically coupled tunneling during the evolution of protochlorophyllide oxidoreductase

Protein dynamics are crucial for realizing the catalytic power of enzymes, but how enzymes have evolved to achieve catalysis is unknown. The light-activated enzyme protochlorophyllide oxidoreductase (POR) catalyzes sequential hydride and proton transfers in the photoexcited and ground states, respectively, and is an excellent system for relating the effects of motions to catalysis. Here, we have used the temperature dependence of isotope effects and solvent viscosity measurements to analyze the dynamics coupled to the hydride and proton transfer steps in three cyanobacterial PORs and a related plant enzyme. We have related the dynamic profiles of each enzyme to their evolutionary origin. Motions coupled to light-driven hydride transfer are conserved across all POR enzymes, but those linked to thermally activated proton transfer are variable. Cyanobacterial PORs require complex and solvent-coupled dynamic networks to optimize the proton donor-acceptor distance, but evolutionary pressures appear to have minimized such networks in plant PORs. POR from Gloeobacter violaceus has features of both the cyanobacterial and plant enzymes, suggesting that the dynamic properties have been optimized during the evolution of POR. We infer that the differing trajectories in optimizing a catalytic structure are related to the stringency of the chemistry catalyzed and define a functional adaptation in which active site chemistry is protected from the dynamic effects of distal mutations that might otherwise impact negatively on enzyme catalysis.

Characterizing the dynamics of functionally relevant complexes of formate dehydrogenase

The potential for femtosecond to picosecond time-scale motions to influence the rate of the intrinsic chemical step in enzyme-catalyzed reactions is a source of significant controversy. Among the central challenges in resolving this controversy is the difficulty of experimentally characterizing thermally activated motions at this time scale in functionally relevant enzyme complexes. We report a series of measurements to address this problem using two-dimensional infrared spectroscopy to characterize the time scales of active-site motions in complexes of formate dehydrogenase with the transition-state-analog inhibitor azide (N(3)(-)). We observe that the frequency-frequency time correlation functions (FFCF) for the ternary complexes with NAD(+) and NADH decay completely with slow time constants of 3.2 ps and 4.6 ps, respectively. This result suggests that in the vicinity of the transition state, the active-site enzyme structure samples a narrow and relatively rigid conformational distribution indicating that the transition-state structure is well organized for the reaction. In contrast, for the binary complex, we observe a significant static contribution to the FFCF similar to what is seen in other enzymes, indicating the presence of the slow motions that occur on time scales longer than our measurement window.

Hydrogen tunneling in adenosylcobalamin-dependent glutamate mutase: evidence from intrinsic kinetic isotope effects measured by intramolecular competition

Hydrogen atom transfer reactions between the substrate and coenzyme are key mechanistic features of all adenosylcobalamin-dependent enzymes. For one of these enzymes, glutamate mutase, we have investigated whether hydrogen tunneling makes a significant contribution to the mechanism by examining the temperature dependence of the deuterium kinetic isotope effect associated with the transfer of a hydrogen atom from methylaspartate to the coenzyme. To do this, we designed a novel intramolecular competition experiment that allowed us to measure the intrinsic kinetic isotope effect, even though hydrogen transfer may not be rate-determining. From the Arrhenius plot of the kinetic isotope effect, the ratio of the pre-exponential factors (A(H)/A(D)) was 0.17 +/- 0.04 and the isotope effect on the activation energy [DeltaE(a(D-H))] was 1.94 +/- 0.13 kcal/mol. The results imply that a significant degree of hydrogen tunneling occurs in glutamate mutase, even though the intrinsic kinetic isotope effects are well within the semiclassical limit and are much smaller than those measured for other AdoCbl enzymes and model reactions for which hydrogen tunneling has been implicated.

Examination of enzymatic H-tunneling through kinetics and dynamics

In recent years, kinetic measurements of isotope effects of enzyme-catalyzed reactions and their temperature dependence led to the development of theoretical models that were used to rationalize the findings. These models suggested that motions at the femto- to picosecond (fs to ps) time scale modulate the environment of the catalyzed reaction. Due to the fast nature of motions that directly affect the cleavage of a covalent bond, it is challenging to correlate the enzyme kinetics and dynamics related to that step. We report a study of formate dehydrogenase (FDH) that compares the temperature dependence of intrinsic kinetic isotope effects (KIEs) to measurements of the environmental dynamics at the fs-ps time scale (Bandaria et al. J. Am. Chem. Soc. 2008, 130, 22-23). The findings from this comparison of experimental kinetics and dynamics are consistent with models of environmentally coupled H-tunneling models, also known as Marcus-like models. Apparently, at tunneling ready conformations, the donor-acceptor distance, orientation, and fluctuations seems to be well tuned for H-transfer and are not affected by thermal fluctuations slower than 10 ps. This phenomenon has been suggested in the past to be quite general in enzymatic reactions. Here, the kinetics and the dynamics measurements on a single chemical step and on fs-ps time scale, respectively, provide new insight and support for the relevant theoretical models. Furthermore, this methodology could be applied to other systems and be used to examine mutants for which the organization of the donor and acceptor is not ideal, or enzymes with different rigidity and different temperature optimum.

Microscale synthesis and kinetic isotope effect analysis of (4R)-[Ad-14C, 4-2H] NADPH and (4R)-[Ad-3H,4-2H] NADPH

We present a one-pot chemo-enzymatic microscale synthesis of NADPH with two different patterns of isotopic labels: (4R)-[Ad-¹⁴C,4-²H] NADPH and (4R)-[Ad-³H,4-²H] NADPH. These co-factors are required by an enormous range of enzymes, and isotopically labeled nicotinamides are consequently of significant interest to researchers. In the current procedure, [Ad-¹⁴C] NAD⁺ and [Ad-³H] NAD⁺ were phosphorylated by NAD⁺ kinase to produce [Ad-¹⁴C] NADP⁺ and [Ad-³H] NADP⁺, respectively. Thermoanaerobium brockii alcohol dehydrogenase was then used to stereospecifically transfer deuterium from C2 of isopropanol to the re face of C4 of NADP⁺. After purification by HPLC, NMR analysis indicated the deuterium content at the 4R position is more than 99.7 %. The labeled cofactors were then used to successfully and sensitively measure kinetic isotope effects for R67 dihydrofolate reductase, providing strong evidence for the utility of this synthetic methodology.