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

Transitioning enzyme catalysis towards photocatalysis

Enzyme biocatalysis is being industrialized at a phenomenal rate. Biocatalysis offers routes to chemical transformations that avoid the use of expensive metal catalysts, high temperatures and pressures, while providing impressive enantio-, regio- and chemo-selectivities. Working individually or as cascades, in live cells or cell-free preparations, to manufacture everyday chemicals, materials, healthcare products, fuels and pharmaceuticals and in diagnostic and industrial sensing applications, enzymes are key enablers in a circular bioeconomy. An ability to exploit and tailor biocatalysts rapidly and predictably requires knowledge of structure-mechanism relationships and the physical chemistry of enzyme action. This knowledge has advanced since our millennium article on this topic (Sutcliffe and Scrutton Phil Trans R. Soc. Lond. A. 2000. 358, 367-386). Here, we discuss an emerging frontier-enzyme photobiocatalysis. Photoenzymes are rarely found in nature. This limits 'difficult-to-achieve' reactions in biology that are generally accessible to chemical photocatalysts. We discuss here the emergence of photobiocatalysis as a new frontier. We review knowledge of natural photoenzymes and identify challenges and limitations in their use as photobiocatalysts. We consider emerging reports on repurposing natural enzymes as photobiocatalysts. We also discuss prospects for de novo design of photobiocatalysts which as a general concept would transform catalysis science.This article is part of the theme issue 'Science into the next millennium: 25 years on'.

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.

Decreased Transition-State Analogue Affinity in Isotopically Heavy MTAN with Increased Catalysis

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase from Helicobacter pylori (HpMTAN) demonstrated faster chemistry when expressed as an isotopically heavy protein, with 2H, 13C, and 15N replacing the bulk of normal isotopes. The inverse heavy enzyme isotope effect has been attributed to improved enzyme-reactant interactions causing more frequent transition-state formation ( Proc. Natl. Acad. Sci. U.S.A. 2021, 118, e2109118118). Transition-state analogues stabilize the transient dynamic geometry of the transition state and inform on transition-state dynamics. Here, a slow-onset, tight-binding transition-state analogue of HpMTAN is characterized with heavy and light enzymes. Dissociation constants for the initial encounter complex (Ki) and for the tightly bound complex after slow-onset inhibition (Ki*) with hexylthio-DADMe-Immucillin-A (HTDIA) gave Ki values for light and heavy HpMTAN = 52 ± 10 and 85 ± 13 pM and Ki* values = 5.9 ± 0.3 and 10.0 ± 1.2 pM, respectively. HTDIA dissociates from heavy HpMTAN at 0.063 ± 0.002 min-1, faster than that from light HpMTAN at 0.032 ± 0.004 min-1. These values are consistent with transition-state formation by an improved catalytic site dynamic search and inconsistent with catalytic efficiency proportional to tight binding of the transition state.

Photoinduced Electron Transfer from a 1,4,5,6-Tetrahydro Nicotinamide Adenine Dinucleotide (Phosphate) Analogue to Oxidized Flavin in an Ene-Reductase Flavoenzyme

Recent reports have described the use of ene-reductase flavoenzymes to catalyze non-natural photochemical reactions. These studies have focused on using reduced flavoenzyme, yet oxidized flavins have superior light harvesting properties. In a binary complex of the oxidized ene-reductase pentaerythritol tetranitrate reductase with the nonreactive nicotinamide coenzyme analogs 1,4,5,6-tetrahydro NAD(P)H, visible photoexcitation of the flavin mononucleotide (FMN) leads to one-electron transfer from the NAD(P)H₄ to FMN, generating a NAD(P)H₄ cation radical and anionic FMN semiquinone. This electron transfer occurs in ∼1 ps and appears to kinetically outcompete reductive quenching from aromatic residues in the active site. Time-resolved infrared measurements show that relaxation processes appear to be largely localized on the FMN and the charge-separated state is short-lived, with relaxation, presumably via back electron transfer, occurring over ∼3-30 ps. While this demonstrates the potential for non-natural photoactivity, useful photocatalysis will likely require longer-lived excited states, which may be accessible by enzyme engineering and/or a judicious choice of substrate.

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.

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.

Covalent Attachment of Active Enzymes to Upconversion Phosphors Allows Ratiometric Detection of Substrates

Upconverting phosphors (UCPs) convert multiple low energy photons into higher energy emission via the process of photon upconversion and offer an attractive alternative to organic fluorophores for use as luminescent probes. Here, UCPs were capped with functionalized silica in order to provide a surface to covalently conjugate proteins with surface-accessible cysteines. Variants of green fluorescent protein (GFP) and the flavoenzyme pentaerythritol tetranitrate reductase (PETNR) were then attached via maleimide-thiol coupling in order to allow energy transfer from the UCP to the GFP or flavin cofactor of PETNR, respectively. PETNR retains its activity when coupled to the UCPs, which allows reversible detection of enzyme substrates via ratiometric sensing of the enzyme redox state.

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.

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.

Temperature, Dynamics, and Enzyme-Catalyzed Reaction Rates

We review the adaptations of enzyme activity to different temperatures. Psychrophilic (cold-adapted) enzymes show significantly different activation parameters (lower activation enthalpies and entropies) from their mesophilic counterparts. Furthermore, there is increasing evidence that the temperature dependence of many enzyme-catalyzed reactions is more complex than is widely believed. Many enzymes show curvature in plots of activity versus temperature that is not accounted for by denaturation or unfolding. This is explained by macromolecular rate theory: A negative activation heat capacity for the rate-limiting chemical step leads directly to predictions of temperature optima; both entropy and enthalpy are temperature dependent. Fluctuations in the transition state ensemble are reduced compared to the ground state. We show how investigations combining experiment with molecular simulation are revealing fundamental details of enzyme thermoadaptation that are relevant for understanding aspects of enzyme evolution. Simulations can calculate relevant thermodynamic properties (such as activation enthalpies, entropies, and heat capacities) and reveal the molecular mechanisms underlying experimentally observed behavior.

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.

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.

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.

Exposing the Interplay Between Enzyme Turnover, Protein Dynamics, and the Membrane Environment in Monoamine Oxidase B

There is an increasing realization that structure-based drug design may show improved success by understanding the ensemble of conformations accessible to an enzyme and how the environment affects this ensemble. Human monoamine oxidase B (MAO-B) catalyzes the oxidation of amines and is inhibited for the treatment of both Parkinson's disease and depression. Despite its clinical importance, its catalytic mechanism remains unclear, and routes to drugging this target would be valuable. Evidence of a radical in either the transition state or the resting state of MAO-B is present throughout the literature and is suggested to be a flavin semiquinone, a tyrosyl radical, or both. Here we see evidence of a resting-state flavin semiquinone, via absorption redox studies and electron paramagnetic resonance, suggesting that the anionic semiquinone is biologically relevant. On the basis of enzyme kinetic studies, enzyme variants, and molecular dynamics simulations, we find evidence for the importance of the membrane environment in mediating the activity of MAO-B and that this mediation is related to the protein dynamics of MAO-B. Further, our MD simulations identify a hitherto undescribed entrance for substrate binding, membrane modulated substrate access, and indications for half-site reactivity: only one active site is accessible to binding at a time. Our study combines both experimental and computational evidence to illustrate the subtle interplay between enzyme activity and protein dynamics and the immediate membrane environment. Understanding key biomedical enzymes to this level of detail will be crucial to inform strategies (and binding sites) for rational drug design for these targets.

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.

Inverse enzyme isotope effects in human purine nucleoside phosphorylase with heavy asparagine labels

Transition path-sampling calculations with several enzymes have indicated that local catalytic site femtosecond motions are linked to transition state barrier crossing. Experimentally, femtosecond motions can be perturbed by labeling the protein with amino acids containing ¹³C, ¹⁵N, and nonexchangeable ²H. A slowed chemical step at the catalytic site with variable effects on steady-state kinetics is usually observed for heavy enzymes. Heavy human purine nucleoside phosphorylase (PNP) is slowed significantly (k_chemlight/k_chemheavy = 1.36). An asparagine (Asn243) at the catalytic site is involved in purine leaving-group activation in the PNP catalytic mechanism. In a PNP produced with isotopically heavy asparagines, the chemical step is faster (k_chemlight/k_chemheavy = 0.78). When all amino acids in PNP are heavy except for the asparagines, the chemical step is also faster (k_chemlight/k_chemheavy = 0.71). Substrate-trapping experiments provided independent confirmation of improved catalysis in these constructs. Transition path-sampling analysis of these partially labeled PNPs indicate altered femtosecond catalytic site motions with improved Asn243 interactions to the purine leaving group. Altered transition state barrier recrossing has been proposed as an explanation for heavy-PNP isotope effects but is incompatible with these isotope effects. Rate-limiting product release governs steady-state kinetics in this enzyme, and kinetic constants were unaffected in the labeled PNPs. The study suggests that mass-constrained femtosecond motions at the catalytic site of PNP can improve transition state barrier crossing by more frequent sampling of essential catalytic site contacts.

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.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

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.

Extracting Kinetic Isotope Effects From a Global Analysis of Reaction Progress Curves

Enzyme reaction progress curves, or time course datasets, are often rich in information, yet their analysis typically reduces their information content to a single parameter, the initial velocity. An alternative approach is described here, where the time course is described by a model constructed from rate equations. In combination with global nonlinear regression, intrinsic rate and/or equilibrium constants can be directly obtained by fitting these data. This method can be greatly enhanced when combined with the measurement of (usually deuterium) isotope effects, which selectively perturb individual step(s) within the reaction, allowing better separation of fitted parameters and more robust model testing. This chapter focuses on practical considerations when using analytical and/or numerically integrated rate equations to model enzyme reactions. The emphasis is on the underlying methodology, which is demonstrated with specific examples alongside recommendations of suitable software.

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.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.