The relationship between protein dynamics and function is a subject of considerable contemporary interest. Although protein motions are frequently observed during ligand binding and release steps, the contribution of protein motions to the catalysis of bond making/breaking processes is more difficult to probe and verify. Here, we show how the quantum mechanical hydrogen tunneling associated with enzymatic C-H bond cleavage provides a unique window into the necessity of protein dynamics for achieving optimal catalysis. Experimental findings support a hierarchy of thermodynamically equilibrated motions that control the H-donor and -acceptor distance and active-site electrostatics, creating an ensemble of conformations suitable for H-tunneling. A possible extension of this view to methyl transfer and other catalyzed reactions is also presented. The impact of understanding these dynamics on the conceptual framework for enzyme activity, inhibitor/drug design, and biomimetic catalyst design is likely to be substantial.

The role of the enzyme’s dynamic motions in catalysis is at the center of heated contemporary debates among both theoreticians and experimentalists. Resolving these apparent disputes is of both intellectual and practical importance: incorporation of enzyme dynamics could be critical for any calculation of enzymatic function and may have profound implications for structure-based drug design and the design of biomimetic catalysts. Analysis of the literature suggests that while part of the dispute may reflect substantial differences between theoretical approaches, much of the debate is semantic. For example, the term “protein dynamics” is often used by some researchers when addressing motions that are in thermal equilibrium with their environment, while other researchers only use this term for nonequilibrium events. The last cases are those in which thermal energy is “stored” in a specific protein mode and “used” for catalysis before it can dissipate to its environment (i.e., “nonstatistical dynamics”). This terminology issue aside, a debate has arisen among theoreticians around the roles of nonstatistical vs statistical dynamics in catalysis. However, the author knows of no experimental findings available today that examined this question in enzyme catalyzed reactions. Another source of perhaps nonsubstantial argument might stem from the varying time scales of enzymatic motions, which range from seconds to femtoseconds. Motions at different time scales play different roles in the many events along the catalytic cascade (reactant binding, reprotonation of reactants, structural rearrangement toward the transition state, product release, etc.). In several cases, when various experimental tools have been used to probe catalytic events at differing time scales, illusory contradictions seem to have emerged. In this Account, recent attempts to sort the merits of those questions are discussed along with possible future directions. A possible summary of current studies could be that enzyme, substrate, and solvent dynamics contribute to enzyme catalyzed reactions in several ways: first via mutual “induced-fit” shifting of their conformational ensemble upon binding; then via thermal search of the conformational space toward the reaction’s transition-state (TS) and the rare event of the barrier crossing toward products, which is likely to be on faster time scales then the first and following events; and finally via the dynamics associated with products release, which are rate-limiting for many enzymatic reactions. From a chemical perspective, close to the TS, enzymatic systems seem to stiffen, restricting motions orthogonal to the chemical coordinate and enabling dynamics along the reaction coordinate to occur selectively. Studies of how enzymes evolved to support those efficient dynamics at various time scales are still in their infancy, and further experiments and calculations are needed to reveal these phenomena in both enzymes and uncatalyzed reactions.

This paper draws attention to selected experiments on enzyme-catalyzed reactions that show convex Arrhenius plots, which are very rare, and points out that Tolman's interpretation of the activation energy places a fundamental model-independent constraint on any detailed explanation of these reactions. The analysis presented here shows that in such systems, the rate coefficient as a function of energy is not just increasing more slowly than expected, it is actually decreasing. This interpretation of the data provides a constraint on proposed microscopic models, i.e., it requires that any successful model of a reaction with a convex Arrhenius plot should be consistent with the microcanonical rate coefficient being a decreasing function of energy. The implications and limitations of this analysis to interpreting enzyme mechanisms are discussed. This model-independent conclusion has broad applicability to all fields of kinetics, and we also draw attention to an analogy with diffusion in metastable fluids and glasses.

One of the most intriguing questions in modern enzymology is whether enzyme dynamics evolved to enhance the catalyzed chemical transformation. In this study, dihydrofolate reductase, a small monomeric protein that catalyzes a single C-H–C transfer, is used as a model system to address this question. Experimental and computational studies have proposed a dynamic network that includes two residues remote from the active site (G121 and M42). The current study compares the nature of the H-transfer step of the WT enzyme, two single mutants, and their double mutant. The contribution of quantum mechanical tunneling and enzyme dynamics to the H-transfer step was examined by determining intrinsic kinetic isotope effects, their temperature dependence, and activation parameters. Different patterns of environmentally coupled tunneling were found for these four enzymes. The findings indicate that the naturally evolved WT dihydrofolate reductase requires no donor–acceptor distance fluctuations (no gating). Both single mutations affect the rearrangement of the system before tunneling, so some gating is required, but the overall nature of the environmentally coupled tunneling appears similar to that of the WT enzyme. The double mutation, on the other hand, seems to cause a major change in the nature of H transfer, leading to poor reorganization and substantial gating. These findings support the suggestion that these distal residues synergistically affect the H transfer at the active site of the enzyme. This observation is in accordance with the notion that these remote residues are part of a dynamic network that is coupled to the catalyzed chemistry.