Fast protein motions are coupled to enzyme H-transfer reactions

Study of the Effects of Remote Heavy Group Vibrations on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions

It has recently been observed that the temperature(T)-dependence of KIEs in H-tunneling reactions, characterized by isotopic activation energy difference (ΔEa = EaD - EaH), is correlated to the rigidity of the tunneling ready states (TRSs) in enzymes. A more rigid system with narrowly distributed H-donor-acceptor distances (DADs) at the TRSs gives rise to a weaker T-dependence of KIEs (i.e., a smaller ΔEa). Theoreticians have attempted to develop new H-tunneling models to explain this, but none has been universally accepted. In order to further understand the observations in enzymes and provide useful data to build new theoretical models, we have studied the electronic and solvent effects on ΔEa's for the hydride-tunneling reactions of NADH/NAD+ analogues. We found that a tighter charge-transfer (CT) complex system gives rises to a smaller ΔEa, consistent with the enzyme observations. In this paper, we use the remote heavy group (R) vibrational effects to mediate the system rigidity to study the rigidity-ΔEa relationship. The specific hypothesis is that slower vibrations of a heavier remote group would broaden the DAD distributions and increase the ΔEa value. Four NADH/NAD+ systems were studied in acetonitrile but most of such heavy group vibrations do not appear to significantly increase the ΔEa. The remote heavy group vibrations in these systems may have not affected the CT complexation rigidity to a degree that can significantly increase the DADs, and further, the ΔEa values.

Allosteric activation unveils protein-mass modulation of ATP phosphoribosyltransferase product release

Heavy-isotope substitution into enzymes slows down bond vibrations and may alter transition-state barrier crossing probability if this is coupled to fast protein motions. ATP phosphoribosyltransferase from Acinetobacter baumannii is a multi-protein complex where the regulatory protein HisZ allosterically enhances catalysis by the catalytic protein HisGS. This is accompanied by a shift in rate-limiting step from chemistry to product release. Here we report that isotope-labelling of HisGS has no effect on the nonactivated reaction, which involves negative activation heat capacity, while HisZ-activated HisGS catalytic rate decreases in a strictly mass-dependent fashion across five different HisGS masses, at low temperatures. Surprisingly, the effect is not linked to the chemical step, but to fast motions governing product release in the activated enzyme. Disruption of a specific enzyme-product interaction abolishes the isotope effects. Results highlight how altered protein mass perturbs allosterically modulated thermal motions relevant to the catalytic cycle beyond the chemical step.

Phonon-assisted electron-proton transfer in [FeFe] hydrogenases: Topological role of clusters

[FeFe] hydrogenases are enzymes that have acquired a unique capacity to synthesize or consume molecular hydrogen (H2). This function relies on a complex catalytic mechanism involving the active site and two distinct electron and proton transfer networks working in concert. By an analysis based on terahertz vibrations of [FeFe] hydrogenase structure, we are able to predict and identify the existence of rate-promoting vibrations at the catalytic site and the coupling with functional residues involved in reported electron and proton transfer networks. Our findings suggest that the positioning of the cluster is influenced by the response of the scaffold to thermal fluctuations, which in turn drives the formation of networks for electron transfer through phonon-assisted mechanisms. Thus, we address the problem of linking the molecular structure to the catalytic function through picosecond dynamics, while raising the functional gain brought by the cofactors or clusters, using the concept of fold-encoded localized vibrations.

Employing deuterium kinetic isotope effects to uncover the mechanism of (R)-3-hydroxybutyrate dehydrogenase

Short-chain dehydrogenases/reductases (SDR) form a large enzyme superfamily playing important roles in health and disease. Furthermore, they are useful tools in biocatalysis. Unveiling the nature of the transition state for hydride transfer is a crucial undertaking toward defining the physicochemical underpinnings of catalysis by SDR enzymes, including possible contributions from quantum mechanical tunneling. Primary deuterium kinetic isotope effects can uncover the contribution from chemistry to the rate-limiting step and potentially provide detailed information on the hydride-transfer transition state in SDR-catalyzed reactions. For the latter, however, one needs to determine the intrinsic isotope effect: that which would be measured if hydride transfer were rate determining. Alas, as is the case for many other enzymatic reactions, those catalyzed by SDRs are often limited by the rate of isotope-insensitive steps, such as product release and conformational changes, which masks the expression of the intrinsic isotope effect. This can be overcome by the powerful yet underexplored method of Palfey and Fagan via which intrinsic kinetic isotope effects can be extracted from pre-steady-state kinetics data. SDRs are ideal systems to which this method can be applied. We have employed this approach to elucidate the transition states for hydride transfer catalyzed by NADH-dependent cold- and warm-adapted (R)-3-hydroxybutyrate dehydrogenase. Experimental conditions which simplify the analysis are discussed.

Allosteric rescue of catalytically impaired ATP phosphoribosyltransferase variants links protein dynamics to active-site electrostatic preorganisation

ATP phosphoribosyltransferase catalyses the first step of histidine biosynthesis and is controlled via a complex allosteric mechanism where the regulatory protein HisZ enhances catalysis by the catalytic protein HisG_S while mediating allosteric inhibition by histidine. Activation by HisZ was proposed to position HisG_S Arg56 to stabilise departure of the pyrophosphate leaving group. Here we report active-site mutants of HisG_S with impaired reaction chemistry which can be allosterically restored by HisZ despite the HisZ:HisG_S interface lying ~20 Å away from the active site. MD simulations indicate HisZ binding constrains the dynamics of HisG_S to favour a preorganised active site where both Arg56 and Arg32 are poised to stabilise leaving-group departure in WT-HisG_S. In the Arg56Ala-HisG_S mutant, HisZ modulates Arg32 dynamics so that it can partially compensate for the absence of Arg56. These results illustrate how remote protein-protein interactions translate into catalytic resilience by restoring damaged electrostatic preorganisation at the active site.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Inverse heavy enzyme isotope effects in methylthioadenosine nucleosidases

Heavy enzyme isotope effects occur in proteins substituted with 2H-, 13C-, and 15N-enriched amino acids. Mass alterations perturb femtosecond protein motions and have been used to study the linkage between fast motions and transition-state barrier crossing. Heavy enzymes typically show slower rates for their chemical steps. Heavy bacterial methylthioadenosine nucleosidases (MTANs from Helicobactor pylori and Escherichia coli) gave normal isotope effects in steady-state kinetics, with slower rates for the heavy enzymes. However, both enzymes revealed rare inverse isotope effects on their chemical steps, with faster chemical steps in the heavy enzymes. Computational transition-path sampling studies of H. pylori and E. coli MTANs indicated closer enzyme-reactant interactions in the heavy MTANs at times near the transition state, resulting in an improved reaction coordinate geometry. Specific catalytic interactions more favorable for heavy MTANs include improved contacts to the catalytic water nucleophile and to the adenine leaving group. Heavy bacterial MTANs depart from other heavy enzymes as slowed vibrational modes from the heavy isotope substitution caused improved barrier-crossing efficiency. Improved sampling frequency and reactant coordinate distances are highlighted as key factors in MTAN transition-state stabilization.

Monoclonal antibody stability can be usefully monitored using the excitation-energy-dependent fluorescence edge-shift

Among the major challenges in the development of biopharmaceuticals are structural heterogeneity and aggregation. The development of a successful therapeutic monoclonal antibody (mAb) requires both a highly active and also stable molecule. Whilst a range of experimental (biophysical) approaches exist to track changes in stability of proteins, routine prediction of stability remains challenging. The fluorescence red edge excitation shift (REES) phenomenon is sensitive to a range of changes in protein structure. Based on recent work, we have found that quantifying the REES effect is extremely sensitive to changes in protein conformational state and dynamics. Given the extreme sensitivity, potentially this tool could provide a ‘fingerprint’ of the structure and stability of a protein. Such a tool would be useful in the discovery and development of biopharamceuticals and so we have explored our hypothesis with a panel of therapeutic mAbs. We demonstrate that the quantified REES data show remarkable sensitivity, being able to discern between structurally identical antibodies and showing sensitivity to unfolding and aggregation. The approach works across a broad concentration range (µg–mg/ml) and is highly consistent. We show that the approach can be applied alongside traditional characterisation testing within the context of a forced degradation study (FDS). Most importantly, we demonstrate the approach is able to predict the stability of mAbs both in the short (hours), medium (days) and long-term (months). The quantified REES data will find immediate use in the biopharmaceutical industry in quality assurance, formulation and development. The approach benefits from low technical complexity, is rapid and uses instrumentation which exists in most biochemistry laboratories without modification.

Protein Mass-Modulated Effects in Alkaline Phosphatase

Recent experimental studies engaging isotopically substituted protein (heavy protein) have revealed that many, but not all, enzymatic systems exhibit altered chemical steps in response to an altered mass. The results have been interpreted as femtosecond protein dynamics at the active site being linked (or not) to transition-state barrier crossing. An altered enzyme mass can influence several kinetic parameters (kcat, Km, and kchem) in amounts of ≤30% relative to light enzymes. An early report on deuterium-labeled Escherichia coli alkaline phosphatase (AP) showed an unusually large enzyme kinetic isotope effect on kcat. We examined steady-state and chemical step properties of native AP, [2H]AP, and [2H,13C,15N]AP to characterize the role of heavy enzyme protein dynamics in reactions catalyzed by AP. Both [2H]- and [2H,13C,15N]APs showed unaltered steady-state and single-turnover rate constants. These findings characterize AP as one of the enzymes in which mass-dependent catalytic site dynamics is dominated by reactant-linked atomic motions. Two catalytic site zinc ions activate the oxygen nucleophiles in the catalytic site of AP. The mass of the zinc ions is unchanged in light and heavy APs. They are essentially linked to catalysis and provide a possible explanation for the loss of linkage between catalysis and protein mass in these enzymes.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Role of Protein Motions in Catalysis by Formate Dehydrogenase

We have analyzed the reaction catalyzed by formate dehydrogenase using transition path sampling. This system has recently received experimental attention using infrared spectroscopy and heavy-enzyme studies. Some of the experimental results point to the possible importance of protein motions that are coupled to the chemical step. We found that the residue Val123 that lies behind the nicotinamide ring occasionally comes into van der Waals contact with the acceptor and that in all reactive trajectories, the barrier-crossing event is preceded by this contact, meaning that the motion of Val123 is part of the reaction coordinate. Experimental results have been interpreted with a two-dimensional formula for the chemical rate, which cannot capture effects such as the one we describe.

Directed Evolution's Influence on Rapid Density Fluctuations Illustrates How Protein Dynamics Can Become Coupled to Chemistry

Protein engineering is a growing field with a variety of experimental techniques available for altering protein function. However, creating an enzyme de novo is still in its infancy, so far yielding enzymes of modest catalytic efficiency. In this study, a system of artificial retro-aldolase enzymes found to have chemistry coupled to protein dynamics was examined. The original design was created computationally, and this protein was then subjected to directed evolution to improve the initial low catalytic efficiency. We found that this re-engineering of the enzyme resulted in rapid density fluctuations throughout the enzyme being reshaped via alterations in the hydrogen bonding network. This work also led to the discovery of a second important motion which aids in the release of an intermediate product. These results provide compelling evidence that to engineer efficient protein catalysts, fast protein dynamics need to be considered in the design.

Ultrafast Vibrational Energy Transfer between Protein and Cofactor in a Flavoenzyme

Protein motions and enzyme catalysis are often linked. It is hypothesized that ultrafast vibrations (femtosecond-picosecond) enhance the rate of hydride transfer catalyzed by members of the old yellow enzyme (OYE) family of ene-reductases. Here, we use time-resolved infrared (TRIR) spectroscopy in combination with stable "heavy" isotopic labeling (2H, 13C, 15N) of protein and/or cofactor to probe the vibrational energy transfer (VET) between pentaerythritol tetranitrate reductase (a member of the OYE family) and its noncovalently bound flavin mononucleotide (FMN) cofactor. We show that when the FMN cofactor is photoexcited with visible light, vibrational energy is transferred from the flavin to the surrounding protein environment on the picosecond timescale. This finding expands the scope of VET investigation in proteins, which are limited by suitable intrinsic probes, and may have implications in the understanding of the mechanism of recently discovered photoactive flavoenzymes.

What are the signatures of tunnelling in enzyme-catalysed reactions?

Computed tunnelling contributions and correlations between apparent activation enthalpy and entropy are explored for the interpretation of enzyme-catalysed H-transfer reactions.

Hydrogen Tunnelling as a Probe of the Involvement of Water Vibrational Dynamics in Aqueous Chemistry?

Our study of tunnelling in proton-coupled electron transfer (PCET) oxidation of ascorbate with hexacyanoferrate(III) follows the insights obtained from ultrafast 2D IR spectroscopy and theoretical studies of the vibrational water dynamics that led to the proposal of the involvement of collective intermolecular excitonic vibrational water dynamics in aqueous chemistry. To test the proposal, the hydrogen tunnelling modulation observed in the PCET reaction studied in the presence of low concentrations of various partial hydrophobic solutes in the water reaction system has been analyzed in terms of the proposed involvement of the collective intermolecular vibrational water dynamics in activation process in the case. The strongly linear correlation between common tunnelling signatures, isotopic values of Arrhenius prefactor ratios ln AH/AD and isotopic differences in activation enthalpies ΔΔH‡ (H,D) observed in the process in fairly diluted water solutions containing various partial hydrophobic solutes (such as dioxane, acetonitrile, ethanol, and quaternary ammonium ions) points to the common physical origin of the phenomenon in all the cases. It is suggested that the phenomenon can be rooted in an interplay of delocalized collective intermolecular vibrational dynamics of water correlated with vibrations of the coupled transition configuration, where the donor-acceptor oscillations, the motions being to some degree along the reaction coordinate, lead to modulation of hydrogen tunnelling in the reaction.

Oscillatory Active-site Motions Correlate with Kinetic Isotope Effects in Formate Dehydrogenase

Thermal motions of enzymes have been invoked to explain the temperature dependence of kinetic isotope effects (KIE) in enzyme-catalyzed hydride transfers. Formate dehydrogenase (FDH) from Candida boidinii exhibits a temperature independent KIE that becomes temperature dependent upon mutation of hydrophobic residues in the active site. Ternary complexes of FDH that mimic the transition state structure allow investigation of how these mutations influence active-site dynamics. A combination of X-ray crystallography, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamic simulations characterize the structure and dynamics of the active site. FDH exhibits oscillatory frequency fluctuations on the picosecond timescale, and the amplitude of these fluctuations correlates with the temperature dependence of the KIE. Both the kinetic and dynamic phenomena can be reproduced computationally. These results provide experimental evidence for a connection between the temperature dependence of KIEs and motions of the active site in an enzyme-catalyzed reaction consistent with activated tunneling models of the hydride transfer reaction.

Comparative kinetic isotope effects on first- and second-order rate constants of soybean lipoxygenase variants uncover a substrate-binding network

Lipoxygenases are widespread enzymes found in virtually all eukaryotes, including fungi, and, more recently, in prokaryotes. These enzymes act on long-chain polyunsaturated fatty acid substrates (C18 to C20), raising questions regarding how the substrate threads its way from solvent to the active site. Herein, we report a comparison of the temperature dependence of isotope effects on first- and second-order rate constants among single-site variants of the prototypic plant enzyme soybean lipoxygenase-1 substituted at amino acid residues inferred to impact substrate binding. We created 10 protein variants including four amino acid positions, Val-750, Ile-552, Ile-839, and Trp-500, located within a previously proposed substrate portal. The conversion of these bulky hydrophobic side chains to smaller side chains is concluded to increase the mobility of flanking helices, giving rise to increased off rates for substrate dissociation from the enzyme. In this manner, we identified a specific "binding network" that can regulate movement of the substrate from the solvent to the active site. Taken together with our previous findings on C-H and O2 activation of soybean lipoxygenase-1, these results support the emergence of multiple complementary networks within a single protein scaffold that modulate different steps along the enzymatic reaction coordinate.

Heavy Enzymes and the Rational Redesign of Protein Catalysts

An unsolved mystery in biology concerns the link between enzyme catalysis and protein motions. Comparison between isotopically labelled "heavy" dihydrofolate reductases and their natural-abundance counterparts has suggested that the coupling of protein motions to the chemistry of the catalysed reaction is minimised in the case of hydride transfer. In alcohol dehydrogenases, unnatural, bulky substrates that induce additional electrostatic rearrangements of the active site enhance coupled motions. This finding could provide a new route to engineering enzymes with altered substrate specificity, because amino acid residues responsible for dynamic coupling with a given substrate present as hotspots for mutagenesis. Detailed understanding of the biophysics of enzyme catalysis based on insights gained from analysis of "heavy" enzymes might eventually allow routine engineering of enzymes to catalyse reactions of choice.

Towards the competent conformation for catalysis in the ferredoxin-NADP+ reductase from the Brucella ovis pathogen

Brucella ovis encodes a bacterial subclass 1 ferredoxin-NADP(H) reductase (BoFPR) that, by similarity with other FPRs, is expected either to deliver electrons from NADPH to the redox-based metabolism and/or to oxidize NADPH to regulate the soxRS regulon that protects bacteria against oxidative damage. Such potential roles for the pathogen survival under infection conditions make of interest to understand and to act on the BoFPR mechanism. Here, we investigate the NADP⁺/H interaction and NADPH oxidation by hydride transfer (HT) to BoFPR. Crystal structures of BoFPR in free and in complex with NADP⁺ hardly differ. The latter shows binding of the NADP⁺ adenosine moiety, while its redox-reactive nicotinamide protrudes towards the solvent. Nonetheless, pre-steady-state kinetics show formation of a charge-transfer complex (CTC-1) prior to the hydride transfer, as well as conversion of CTC-1 into a second charge-transfer complex (CTC-2) concomitantly with the HT event. Thus, during catalysis nicotinamide and flavin reacting rings stack. Kinetic data also identify the HT itself as the rate limiting step in the reduction of BoFPR by NADPH, as well as product release limiting the overall reaction. Using all-atom molecular dynamics simulations with a thermal effect approach we are able to visualise a potential transient catalytically competent interaction of the reacting rings. Simulations indicate that the architecture of the FAD folded conformation in BoFPR might be key in catalysis, pointing to its adenine as an element to orient the reactive atoms in conformations competent for HT.

Crystal Structure and Biophysical Analysis of Furfural-Detoxifying Aldehyde Reductase from Clostridium beijerinckii

Many aldehydes, such as furfural, are present in high quantities in lignocellulose lysates and are fermentation inhibitors, which makes biofuel production from this abundant carbon source extremely challenging. Cbei_3974 has recently been identified as an aldo-keto reductase responsible for partial furfural resistance in Clostridium beijerinckii Rational engineering of this enzyme could enhance the furfural tolerance of this organism, thereby improving biofuel yields. We report an extensive characterization of Cbei_3974 and a single-crystal X-ray structure of Cbei_3974 in complex with NADPH at a resolution of 1.75 Å. Docking studies identified residues involved in substrate binding, and an activity screen revealed the substrate tolerance of the enzyme. Hydride transfer, which is partially rate limiting under physiological conditions, occurs from the pro-R hydrogen of NADPH. Enzyme isotope labeling revealed a temperature-independent enzyme isotope effect of unity, indicating that the enzyme does not use dynamic coupling for catalysis and suggesting that the active site of the enzyme is optimally configured for catalysis with the substrate tested.IMPORTANCE Here we report the crystal structure and biophysical properties of an aldehyde reductase that can detoxify furfural, a common inhibitor of biofuel fermentation found in lignocellulose lysates. The data contained here will serve as a guide for protein engineers to develop improved enzyme variants that would impart furfural resistance to the microorganisms used in biofuel production and thus lead to enhanced biofuel yields from this sustainable resource.

Why are some Enzymes Dimers? Flexibility and Catalysis in Thermotoga Maritima Dihydrofolate Reductase

Dihydrofolate reductase from Thermotoga maritima (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue Escherichia coli DHFR (EcDHFR). Using QM/MM potentials, we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme: it protects its structure against denaturation by reducing thermal fluctuations and it provides a less negative activation entropy, toning down the increase of the activation free energy with temperature. Our molecular picture is confirmed by the analysis of the temperature dependence of enzyme kinetic isotope effects in different DHFR enzymes.

Optimization of the Turnover in Artificial Enzymes via Directed Evolution Results in the Coupling of Protein Dynamics to Chemistry

The design of artificial enzymes is an emerging field of research. Although progress has been made, the catalytic proficiency of many designed enzymes is low compared to natural enzymes. Nevertheless, recently Hilvert et al. ( Nat. Chem. 2017, 9, 50-56) created a series of five artificial retro-aldolase enzymes via directed evolution, with the final variant exhibiting a rate comparable to the naturally occurring enzyme fructose 1,6 bisphosphate aldolase. We present a study of this system in atomistic detail that elucidates the effects of mutational changes on the chemical step. Transition path sampling is used to create ensembles of reactive trajectories, and committor analysis is used to identify the stochastic separatrix of each ensemble. The application of committor distribution analysis to constrained trajectories allows the identification of changes in important protein motions coupled to reaction across the generated series of the artificial retro-aldolases. We observed two different reaction mechanisms and analyzed the role of the residues participating in the reaction coordinate of each enzyme. However, only in the most evolved variant we identified a fast motion that promotes catalysis, suggesting that this rate promoting vibration was introduced during directed evolution. This study provides further evidence that protein dynamics must be taken into account in designing efficient artificial enzymes.

Discovery, Characterisation, Engineering and Applications of Ene Reductases for Industrial Biocatalysis

Recent studies of multiple enzyme families collectively referred to as ene-reductases (ERs) have highlighted potential industrial application of these biocatalysts in the production of fine and speciality chemicals. Processes have been developed whereby ERs contribute to synthetic routes as isolated enzymes, components of multi-enzyme cascades, and more recently in metabolic engineering and synthetic biology programmes using microbial cell factories to support chemicals production. The discovery of ERs from previously untapped sources and the expansion of directed evolution screening programmes, coupled to deeper mechanistic understanding of ER reactions, have driven their use in natural product and chemicals synthesis. Here we review developments, challenges and opportunities for the use of ERs in fine and speciality chemicals manufacture. The ER research field is rapidly expanding and the focus of this review is on developments that have emerged predominantly over the last 4 years.

Selectivity through discriminatory induced fit enables switching of NAD(P)H coenzyme specificity in Old Yellow Enzyme ene-reductases

Most ene‐reductases belong to the Old Yellow Enzyme (OYE) family of flavin‐dependent oxidoreductases. OYEs use nicotinamide coenzymes as hydride donors to catalyze the reduction of alkenes that contain an electron‐withdrawing group. There have been many investigations of the structures and catalytic mechanisms of OYEs. However, the origin of coenzyme specificity in the OYE family is unknown. Structural NMR and X‐ray crystallographic data were used to rationally design variants of two OYEs, pentaerythritol tetranitrate reductase (PETNR) and morphinone reductase (MR), to discover the basis of coenzyme selectivity. PETNR has dual‐specificity and reacts with NADH and NADPH; MR accepts only NADH as hydride donor. Variants of a β‐hairpin motif in an active site loop of both these enzymes were studied using stopped‐flow spectroscopy. Specific attention was placed on the potential role of arginine residues within the β‐hairpin motif. Mutagenesis demonstrated that Arg130 governs the preference of PETNR for NADPH, and that Arg142 interacts with the coenzyme pyrophosphate group. These observations were used to switch coenzyme specificity in MR by replacing either Glu134 or Leu146 with arginine residues. These variants had increased (~15‐fold) affinity for NADH. Mutagenesis enabled MR to accept NADPH as a hydride donor, with E134R MR showing a significant (55‐fold) increase in efficiency in the reductive half‐reaction, when compared to the essentially unreactive wild‐type enzyme. Insight into the question of coenzyme selectivity in OYEs has therefore been addressed through rational redesign. This should enable coenzyme selectivity to be improved and switched in other OYEs.

Active-Site Glu165 Activation in Triosephosphate Isomerase and Its Deprotonation Kinetics

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) via an enediol(ate) intermediate. The active-site residue Glu165 serves as the catalytic base during catalysis. It abstracts a proton from C1 carbon of DHAP to form the reaction intermediate and donates a proton to C2 carbon of the intermediate to form product GAP. Our difference Fourier transform infrared spectroscopy studies on the yeast TIM (YeTIM)/phosphate complex revealed a C═O stretch band at 1706 cm^-1 from the protonated Glu165 carboxyl group at pH 7.5, indicating that the p K_a of the catalytic base is increased by >3.0 pH units upon phosphate binding, and that the Glu165 carboxyl environment in the complex is still hydrophilic in spite of the increased p K_a. Hence, the results show that the binding of the phosphodianion group is part of the activation mechanism which involves the p K_a elevation of the catalytic base Glu165. The deprotonation kinetics of Glu165 in the μs to ms time range were determined via infrared (IR) T-jump studies on the YeTIM/phosphate and ("heavy enzyme") [U-¹³C,-¹⁵N]YeTIM/phosphate complexes. The slower deprotonation kinetics in the ms time scale is due to phosphate dissociation modulated by the loop motion, which slows down by enzyme mass increase to show a normal heavy enzyme kinetic isotope effect (KIE) ∼1.2 (i.e., slower rate in the heavy enzyme). The faster deprotonation kinetics in the tens of μs time scale is assigned to temperature-induced p K_a decrease, while phosphate is still bound, and it shows an inverse heavy enzyme KIE ∼0.89 (faster rate in the heavy enzyme). The IR static and T-jump spectroscopy provides atomic-level resolution of the catalytic mechanism because of its ability to directly observe the bond breaking/forming process.

Isotopically labeled flavoenzymes and their uses in probing reaction mechanisms

The incorporation of stable isotopes into proteins is beneficial or essential for a range of experiments, including NMR, neutron scattering and reflectometry, proteomic mass spectrometry, vibrational spectroscopy and "heavy" enzyme kinetic isotope effect (KIE) measurements. Here, we present detailed protocols for the stable isotopic labeling of pentaerythritol tetranitrate reductase (PETNR) via recombinant expression in E. coli. PETNR is an ene-reductase belonging to the Old Yellow Enzyme (OYE) family of flavoenzymes, and is regarded as a model system for studying hydride transfer reactions. Included is a discussion of how efficient back-exchange of amide protons in the protein core can be achieved and how the intrinsic flavin mononucleotide (FMN) cofactor can be exchanged, allowing the production of isotopologues with differentially labeled protein and cofactor. In addition to a thorough description of labeling strategies, we briefly exemplify how data analysis and interpretation of "heavy" enzyme KIEs can be performed.

Nonequivalence of Second Sphere "Noncatalytic" Residues in Pentaerythritol Tetranitrate Reductase in Relation to Local Dynamics Linked to H-Transfer in Reactions with NADH and NADPH Coenzymes

Many enzymes that catalyze hydride transfer reactions work via a mechanism dominated by quantum mechanical tunneling. The involvement of fast vibrational modes of the reactive complex is often inferred in these reactions, as in the case of the NAD(P)H-dependent pentaerythritol tetranitrate reductase (PETNR). Herein, we interrogated the H-transfer mechanism in PETNR by designing conservative (L25I and I107L) and side chain shortening (L25A and I107A) PETNR variants and using a combination of experimental approaches (stopped-flow rapid kinetics, X-ray crystallography, isotope/temperature dependence studies of H-transfer and NMR spectroscopy). X-ray data show subtle changes in the local environment of the targeted side chains but no major structural perturbation caused by mutagenesis of these two second sphere active site residues. However, temperature dependence studies of H-transfer revealed a coenzyme-specific and complex thermodynamic equilibrium between different reactive configurations in PETNR–coenzyme complexes. We find that mutagenesis of these second sphere “noncatalytic” residues affects differently the reactivity of PETNR with NADPH and NADH coenzymes. We attribute this to subtle, dynamic structural changes in the PETNR active site, the effects of which impact differently in the nonequivalent reactive geometries of PETNR−NADH and PETNR−NADPH complexes. This inference is confirmed through changes observed in the NMR chemical shift data for PETNR complexes with unreactive 1,4,5,6-tetrahydro-NAD(P) analogues. We show that H-transfer rates can (to some extent) be buffered through entropy–enthalpy compensation, but that use of integrated experimental tools reveals hidden complexities that implicate a role for dynamics in this relatively simple H-transfer reaction. Similar approaches are likely to be informative in other enzymes to understand the relative importance of (distal) hydrophobic side chains and dynamics in controlling the rates of enzymatic H-transfer.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.

Promoting Vibrations and the Function of Enzymes. Emerging Theoretical and Experimental Convergence

A complete understanding of enzyme catalysis requires knowledge of both transition state features and the detailed motions of atoms that cause reactant molecules to form and traverse the transition state. The seeming intractability of the problem arises from the femtosecond lifetime of chemical transition states, preventing most experimental access. Computational chemistry is admirably suited to short time scale analysis but can be misled by inappropriate starting points or by biased assumptions. Kinetic isotope effects provide an experimental approach to transition state structure and a method for obtaining transition state analogues but, alone, do not inform how that transition state is reached. Enzyme structures with transition state analogues provide computational starting points near the transition state geometry. These well-conditioned starting points, combined with the unbiased computational method of transition path sampling, provide realistic atomistic motions involved in transition state formation and passage. In many, but not all, enzymatic systems, femtosecond local protein motions near the catalytic site are linked to transition state formation. These motions are not inherently revealed by most approaches of transition state theory, because transition state theory replaces dynamics with the statistics of the transition state. Experimental and theoretical convergence of the link between local catalytic site vibrational modes and catalysis comes from heavy atom ("Born-Oppenheimer") enzymes. Fully labeled and catalytic site local heavy atom labels perturb the probability of finding enzymatic transition states in ways that can be analyzed and predicted by transition path sampling. Recent applications of these experimental and computational approaches reveal how subpicosecond local catalytic site protein modes play important roles in creating the transition state.

1H, 15N and 13C backbone resonance assignments of pentaerythritol tetranitrate reductase from Enterobacter cloacae PB2

Pentaerythritol tetranitrate reductase (PETNR) is a flavoenzyme possessing a broad substrate specificity and is a member of the Old Yellow Enzyme family of oxidoreductases. As well as having high potential as an industrial biocatalyst, PETNR is an excellent model system for studying hydrogen transfer reactions. Mechanistic studies performed with PETNR using stopped-flow methods have shown that tunneling contributes towards hydride transfer from the NAD(P)H coenzyme to the flavin mononucleotide (FMN) cofactor and fast protein dynamics have been inferred to facilitate this catalytic step. Herein, we report the near-complete ¹H, ¹⁵N and ¹³C backbone resonance assignments of PETNR in a stoichiometric complex with the FMN cofactor in its native oxidized form, which were obtained using heteronuclear multidimensional NMR spectroscopy. A total of 97% of all backbone resonances were assigned, with 333 out of a possible 344 residues assigned in the ¹H-¹⁵N TROSY spectrum. This is the first report of an NMR structural study of a flavoenzyme from the Old Yellow Enzyme family and it lays the foundation for future investigations of functional dynamics in hydride transfer catalytic mechanism.

Directed Evolution as a Probe of Rate Promoting Vibrations Introduced via Mutational Change

In this article, we study with transition path sampling and reaction coordinate analysis how directed evolution in the Kemp eliminase family of artificial enzymes makes differential use of rapid rate promoting vibrations as a component of their chemical mechanism. Even though this family was initially created by placing the expected active site in a fixed protein matrix, we find a shift from largely static to more dynamic active sites that make use of donor-acceptor compression as the evolutionary process proceeds. We see that this introduction of dynamics significantly shifts the order of processes in the reaction. We also suggest that the lack of "design for dynamics" may help explain the relatively low proficiency of such designed enzymes.

Isotope Substitution of Promiscuous Alcohol Dehydrogenase Reveals the Origin of Substrate Preference in the Transition State

The origin of substrate preference in promiscuous enzymes was investigated by enzyme isotope labelling of the alcohol dehydrogenase from Geobacillus stearothermophilus (BsADH). At physiological temperature, protein dynamic coupling to the reaction coordinate was insignificant. However, the extent of dynamic coupling was highly substrate-dependent at lower temperatures. For benzyl alcohol, an enzyme isotope effect larger than unity was observed, whereas the enzyme isotope effect was close to unity for isopropanol. Frequency motion analysis on the transition states revealed that residues surrounding the active site undergo substantial displacement during catalysis for sterically bulky alcohols. BsADH prefers smaller substrates, which cause less protein friction along the reaction coordinate and reduced frequencies of dynamic recrossing. This hypothesis allows a prediction of the trend of enzyme isotope effects for a wide variety of substrates.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Protein Mass Effects on Formate Dehydrogenase

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128β. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

Evidence for proton tunneling and a transient covalent flavin-substrate adduct in choline oxidase S101A

The effect of temperature on the reaction of alcohol oxidation catalyzed by choline oxidase was investigated with the S101A variant of choline oxidase. Anaerobic enzyme reduction in a stopped-flow spectrophotometer was biphasic using either choline or 1,2-[²H₄]-choline as a substrate. The limiting rate constants k_lim1 and k_lim2 at saturating substrate were well separated (k_lim1/k_lim2>9), and were >15-fold slower than for wild-type choline oxidase. Solvent deuterium kinetic isotope effects (KIEs) ~4 established that k_lim1 probes the proton transfer from the substrate hydroxyl to a catalytic base. Primary substrate deuterium KIEs ≥7 demonstrated that k_lim2 reports on hydride transfer from the choline alkoxide to the flavin. Between 15°C and 39°C the k_lim1 and k_lim2 values increased with increasing temperature, allowing for the analyses of H⁺ and H^- transfers using Eyring and Arrhenius formalisms. Temperature-independent KIE on the k_lim1 value (^H2Ok_lim1/^D2Ok_lim1) suggests that proton transfer occurs within a highly reorganized tunneling-ready-state with a narrow distribution of donor-acceptor distances. Eyring analysis of the k_lim2 value gave lines with the slope_(choline)>slope_(D-choline), suggesting kinetic complexity. Spectral evidence for the transient occurrence of a covalent flavin-substrate adduct during the first phase of the anaerobic reaction of S101A CHO with choline is presented, supporting the notion that an important role of amino acid residues in the active site of flavin-dependent enzymes is to eliminate alternative reactions of the versatile enzyme-bound flavin for the reaction that needs to be catalyzed.

Chemical Ligation and Isotope Labeling to Locate Dynamic Effects

Heavy isotope labeling of enzymes slows protein motions without disturbing the electrostatics and can therefore be used to probe the role of dynamics in enzyme catalysis. To identify the structural elements responsible for dynamic effects, individual segments of an enzyme can be labeled and the resulting effect on the kinetics of the reaction can be measured. Such hybrid isotopomers can be constructed by expressed protein ligation, in which complementary labeled and unlabeled peptide segments are prepared by recombinant gene expression and linked by means of chemical ligation. The construction of such hybrid isotopomers is exemplified here with the paradigmatic enzyme dihydrofolate reductase (DHFR) from Escherichia coli.

Measurement of Enzyme Isotope Effects

Enzyme isotope effects, or the kinetic effects of "heavy" enzymes, refer to the effect of isotopically labeled protein residues on the enzyme's activity or physical properties. These effects are increasingly employed in the examination of the possible contributions of protein dynamics to enzyme catalysis. One hypothesis assumed that isotopic substitution of all ¹²C, ¹⁴N, and nonexchangeable ¹H by ¹³C, ¹⁵N, and ²H, would slow down protein picosecond to femtosecond dynamics without any effect on the system's electrostatics following the Born-Oppenheimer approximation. It was suggested that reduced reaction rates reported for several "heavy" enzymes accords with that hypothesis. However, numerous deviations from the predictions of that hypothesis were also reported. Current studies also attempt to test the role of individual residues by site-specific labeling or by labeling a pattern of residues on activity. It appears that in several systems the protein's fast dynamics are indeed reduced in "heavy" enzymes in a way that reduces the probability of barrier crossing of its chemical step. Other observations, however, indicated that slower protein dynamics are electrostatically altered in isotopically labeled enzymes. Interestingly, these effects appear to be system dependent, thus it might be premature to suggest a general role of "heavy" enzymes' effect on catalysis.

A complete thermodynamic analysis of enzyme turnover links the free energy landscape to enzyme catalysis

Our understanding of how enzymes work is coloured by static structure depictions where the enzyme scaffold is presented as either immobile, or in equilibrium between well-defined static conformations. Proteins, however, exhibit a large degree of motion over a broad range of timescales and magnitudes and this is defined thermodynamically by the enzyme free energy landscape (FEL). The role and importance of enzyme motion is extremely contentious. Much of the challenge is in the experimental detection of so called 'conformational sampling' involved in enzyme turnover. Herein we apply combined pressure and temperature kinetics studies to elucidate the full suite of thermodynamic parameters defining an enzyme FEL as it relates to enzyme turnover. We find that the key thermodynamic parameters governing vibrational modes related to enzyme turnover are the isobaric expansivity term and the change in heat capacity for enzyme catalysis. Variation in the enzyme FEL affects these terms. Our analysis is supported by a range of biophysical and computational approaches that specifically capture information on protein vibrational modes and the FEL (all atom flexibility calculations, red edge excitation shift spectroscopy and viscosity studies) that provide independent evidence for our findings. Our data suggest that restricting the enzyme FEL may be a powerful strategy when attempting to rationally engineer enzymes, particularly to alter thermal activity. Moreover, we demonstrate how rational predictions can be made with a rapid computational approach.

Incorporating Fast Protein Dynamics into Enzyme Design: A Proposed Mutant Aromatic Amine Dehydrogenase

In recent years, there has been encouraging progress in the engineering of enzymes that are designed to catalyze reactions not accelerated by natural enzymes. We tested the possibility of reengineering an existing enzyme by introducing a fast protein motion that couples to the reaction. Aromatic amine dehydrogenase is a system that has been shown to use a fast substrate motion as part of the reaction mechanism. We identified a mutation that preserves this fast motion but also introduces a favorable fast motion near the active site that did not exist in the native enzyme. Transition path sampling was used for the analysis of the atomic details of the mechanism.

Triple Isotope Effects Support Concerted Hydride and Proton Transfer and Promoting Vibrations in Human Heart Lactate Dehydrogenase

Transition path sampling simulations have proposed that human heart lactate dehydrogenase (LDH) employs protein promoting vibrations (PPVs) on the femtosecond (fs) to picosecond (ps) time scale to promote crossing of the chemical barrier. This chemical barrier involves both hydride and proton transfers to pyruvate to form l-lactate, using reduced nicotinamide adenine dinucleotide (NADH) as the cofactor. Here we report experimental evidence from three types of isotope effect experiments that support coupling of the promoting vibrations to barrier crossing and the coincidence of hydride and proton transfer. We prepared the native (light) LDH and a heavy LDH labeled with ¹³C, ¹⁵N, and nonexchangeable ²H (D) to perturb the predicted PPVs. Heavy LDH has slowed chemistry in single turnover experiments, supporting a contribution of PPVs to transition state formation. Both the [4-²H]NADH (NADD) kinetic isotope effect and the D₂O solvent isotope effect were increased in dual-label experiments combining both NADD and D₂O, a pattern maintained with both light and heavy LDHs. These isotope effects support concerted hydride and proton transfer for both light and heavy LDHs. Although the transition state barrier-crossing probability is reduced in heavy LDH, the concerted mechanism of the hydride-proton transfer reaction is not altered. This study takes advantage of triple isotope effects to resolve the chemical mechanism of LDH and establish the coupling of fs-ps protein dynamics to barrier crossing.

Untangling Heavy Protein and Cofactor Isotope Effects on Enzyme-Catalyzed Hydride Transfer

"Heavy" (isotopically labeled) enzyme isotope effects offer a direct experimental probe of the role of protein vibrations on enzyme-catalyzed reactions. Here we have developed a strategy to generate isotopologues of the flavoenzyme pentaerythritol tetranitrate reductase (PETNR) where the protein and/or intrinsic flavin mononucleotide (FMN) cofactor are isotopically labeled with ²H, ¹⁵N, and ¹³C. Both the protein and cofactor contribute to the enzyme isotope effect on the reductive hydride transfer reaction, but their contributions are not additive and may partially cancel each other out. However, the isotope effect specifically arising from the FMN suggests that vibrations local to the active site play a role in the hydride transfer chemistry, while the protein-only "heavy enzyme" effect demonstrates that protein vibrations contribute to catalysis in PETNR. In all cases, enthalpy-entropy compensation plays a major role in minimizing the magnitude of "heavy enzyme" isotope effects. Fluorescence lifetime measurements of the intrinsic flavin mononucleotide show marked differences between "light" and "heavy" enzymes on the nanosecond-picosecond time scale, suggesting relevant time scale(s) for those vibrations implicated in the "heavy enzyme" isotope effect on the PETNR reaction.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.