Dynamically achieved active site precision in enzyme catalysis

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Enzyme dynamics-a brief review

Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.

Atomistic description of the relationship between protein dynamics and catalysis with transition path sampling

Despite initial resistance, it has been increasingly accepted that protein dynamics plays a role in enzymatic catalysis. There have been two lines of research. Some works study slow conformational motions that are not coupled to the reaction coordinate, but guide the system towards catalytically competent conformations. Understanding at the atomistic level how this is accomplished has remained elusive except for a few systems. In this review we focus on fast sub-picosecond motions that are coupled to the reaction coordinate. The use of Transition Path Sampling has allowed us an atomistic description of how these rate-promoting vibrational motions are incorporated in the reaction mechanism. We will also show how we used insights from rate-promoting motions in protein design.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Non-active Site Residue in Loop L4 Alters Substrate Capture and Product Release in d-Arginine Dehydrogenase

Numerous studies demonstrate that enzymes undergo multiple conformational changes during catalysis. The malleability of enzymes forms the basis for allosteric regulation: residues located far from the active site can exert long-range dynamical effects on the active site residues to modulate catalysis. The structure of Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) shows four loops (L1, L2, L3, and L4) that span the substrate and the FAD-binding domains. Loop L4 comprises residues 329-336, spanning over the flavin cofactor. The I335 residue on loop L4 is ∼10 Å away from the active site and ∼3.8 Å from N(1)-C(2)═O atoms of the flavin. In this study, we used molecular dynamics and biochemical techniques to investigate the effect of the mutation of I335 to histidine on the catalytic function of PaDADH. Molecular dynamics showed that the conformational dynamics of PaDADH are shifted to a more closed conformation in the I335H variant. In agreement with an enzyme that samples more in a closed conformation, the kinetic data of the I335H variant showed a 40-fold decrease in the rate constant of substrate association (k₁), a 340-fold reduction in the rate constant of substrate dissociation from the enzyme-substrate complex (k₂), and a 24-fold decrease in the rate constant of product release (k₅), compared to that of the wild-type. Surprisingly, the kinetic data are consistent with the mutation having a negligible effect on the reactivity of the flavin. Altogether, the data indicate that the residue at position 335 has a long-range dynamical effect on the catalytic function in PaDADH.

Obtaining anomalous and ensemble information from protein crystals from 220 K up to physiological temperatures

X-ray crystallography has been invaluable in delivering structural information about proteins. Previously, an approach has been developed that allows high-quality X-ray diffraction data to be obtained from protein crystals at and above room temperature. Here, this previous work is built on and extended by showing that high-quality anomalous signal can be obtained from single protein crystals using diffraction data collected at 220 K up to physiological temperatures. The anomalous signal can be used to directly determine the structure of a protein, i.e. to phase the data, as is routinely performed under cryoconditions. This ability is demonstrated by obtaining diffraction data from model lysozyme, thaumatin and proteinase K crystals, the anomalous signal from which allowed their structures to be solved experimentally at 7.1 keV X-ray energy and at room temperature with relatively low data redundancy. It is also demonstrated that the anomalous signal from diffraction data obtained at 310 K (37°C) can be used to solve the structure of proteinase K and to identify ordered ions. The method provides useful anomalous signal at temperatures down to 220 K, resulting in an extended crystal lifetime and increased data redundancy. Finally, we show that useful anomalous signal can be obtained at room temperature using X-rays of 12 keV energy as typically used for routine data collection, allowing this type of experiment to be carried out at widely accessible synchrotron beamline energies and enabling the simultaneous extraction of high-resolution data and anomalous signal. With the recent emphasis on obtaining conformational ensemble information for proteins, the high resolution of the data allows such ensembles to be built, while the anomalous signal allows the structure to be experimentally solved, ions to be identified, and water molecules and ions to be differentiated. Because bound metal-, phosphorus- and sulfur-containing ions all have anomalous signal, obtaining anomalous signal across temperatures and up to physiological temperatures will provide a more complete description of protein conformational ensembles, function and energetics.

Enzymatische C4-Epimerisierung von UDP-Glucuronsäure: präzise gesteuerte Rotation eines transienten 4-Ketointermediats für eine invertierende Reaktion ohne Decarboxylierung

UDP‐Glucuronsäure(UDP‐GlcA)‐4‐Epimerase repräsentiert eine wichtige Fragestellung in der Enzymkatalyse: die Balance zwischen konformativer Flexibilität und genauer Positionierung. Das Enzym koordiniert die C4‐Oxidation des Substrats durch NAD+ mit der Rotation eines leicht decarboxylierbaren β‐Ketosäure‐Intermediats im aktiven Zentrum zur Ermöglichung der stereoinvertierenden Reduktion der Ketogruppe durch NADH. Wir zeigen hier die nur schwer erfassbare Rotationskoordinate des 4‐Ketointermediats. Distorsion des Zuckerrings in eine Boot‐Konformation erzeugt torsionale Mobilität in der Bindungstasche des Enzyms. Die Endpunkte der Rotation zeigen den 4‐Ketozucker in einer unverformten 4C1‐Sesselkonformation. Die äquatorial positionierte Carboxylatgruppe ist ungünstig für die 4‐Ketozucker‐Decarboxylierung. Varianten der Epimerase zeigen Decarboxylierung, wenn sie die Bindung mit der Carboxylatgruppe im entgegengesetzten Rotationsisomer des Substrats entfernen. R185A/D‐Substitutionen wandeln die Epimerase in UDP‐Xylose‐Synthasen um, welche UDP‐GlcA in stereospezifischen, konfigurationserhaltenden Reaktionen decarboxylieren.

Enzymatic C4-Epimerization of UDP-Glucuronic Acid: Precisely Steered Rotation of a Transient 4-Keto Intermediate for an Inverted Reaction without Decarboxylation

UDP-glucuronic acid (UDP-GlcA) 4-epimerase illustrates an important problem regarding enzyme catalysis: balancing conformational flexibility with precise positioning. The enzyme coordinates the C4-oxidation of the substrate by NAD⁺ and rotation of a decarboxylation-prone β-keto acid intermediate in the active site, enabling stereoinverting reduction of the keto group by NADH. We reveal the elusive rotational landscape of the 4-keto intermediate. Distortion of the sugar ring into boat conformations induces torsional mobility in the enzyme's binding pocket. The rotational endpoints show that the 4-keto sugar has an undistorted ⁴ C₁ chair conformation. The equatorially placed carboxylate group disfavors decarboxylation of the 4-keto sugar. Epimerase variants lead to decarboxylation upon removal of the binding interactions with the carboxylate group in the opposite rotational isomer of the substrate. Substitutions R185A/D convert the epimerase into UDP-xylose synthases that decarboxylate UDP-GlcA in stereospecific, configuration-retaining reactions.

Graph-Theoretic Molecular Fragmentation for Potential Surfaces Leads Naturally to a Tensor Network Form and Allows Accurate and Efficient Quantum Nuclear Dynamics

Molecular fragmentation methods have revolutionized quantum chemistry. Here, we use a graph-theoretically generated molecular fragmentation method, to obtain accurate and efficient representations for multidimensional potential energy surfaces and the quantum time-evolution operator, which plays a critical role in quantum chemical dynamics. In doing so, we find that the graph-theoretic fragmentation approach naturally reduces the potential portion of the time-evolution operator into a tensor network that contains a stream of coupled lower-dimensional propagation steps to potentially achieve quantum dynamics with reduced complexity. Furthermore, the fragmentation approach used here has previously been shown to allow accurate and efficient computation of post-Hartree-Fock electronic potential energy surfaces, which in many cases has been shown to be at density functional theory cost. Thus, by combining the advantages of molecular fragmentation with the tensor network formalism, the approach yields an on-the-fly quantum dynamics scheme where both the electronic potential calculation and nuclear propagation portion are enormously simplified through a single stroke. The method is demonstrated by computing approximations to the propagator and to potential surfaces for a set of coupled nuclear dimensions within a protonated water wire problem exhibiting the Grotthuss mechanism of proton transport. In all cases, our approach has been shown to reduce the complexity of representing the quantum propagator, and by extension action of the propagator on an initial wavepacket, by several orders, with minimal loss in accuracy.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Bypassing evolutionary dead ends and switching the rate-limiting step of a human immunotherapeutic enzyme

The Trp metabolite kynurenine (KYN) accumulates in numerous solid tumours and mediates potent immunosuppression. Bacterial kynureninases (KYNases), which preferentially degrade kynurenine, can relieve immunosuppression in multiple cancer models, but immunogenicity concerns preclude their clinical use, while the human enzyme (HsKYNase) has very low activity for kynurenine and shows no therapeutic effect. Using fitness selections, we evolved a HsKYNase variant with 27-fold higher activity, beyond which exploration of >30 evolutionary trajectories involving the interrogation of >109 variants led to no further improvements. Introduction of two amino acid substitutions conserved in bacterial KYNases reduced enzyme fitness but potentiated rapid evolution of variants with ~500-fold improved activity and reversed substrate specificity, resulting in an enzyme capable of mediating strong anti-tumour effects in mice. Pre-steady-state kinetics revealed a switch in rate-determining step attributable to changes in both enzyme structure and conformational dynamics. Apart from its clinical significance, our work highlights how rationally designed substitutions can potentiate trajectories that overcome barriers in protein evolution.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate

Transient oxidation–reduction through hydride transfer with tightly bound NAD coenzyme is used by a large class of sugar nucleotide epimerases to promote configurational inversion of carbon stereocenters in carbohydrate substrates. A requirement for the epimerases to coordinate hydride abstraction and re-addition with substrate rotation in the binding pocket poses a challenge for dynamical protein conformational selection linked to enzyme catalysis. Here, we studied the thermophilic C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose substrate (k_cat ≈ 1.0 min^–1; 60 °C) to explore the sensitivity of the enzymatic hydride transfer toward environmental fluctuations affected by temperature (20–80 °C). We determined noncompetitive primary kinetic isotope effects (KIE) due to ²H at the glucose C2 and showed that a normal KIE on the k_cat (^Dk_cat) reflects isotope sensitivity of the hydrogen abstraction to enzyme-NAD⁺ in a rate-limiting transient oxidation. The ^Dk_cat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C) and 3.3 at high temperature (80 °C). The temperature profiles for k_cat with the ¹H and ²H substrate showed a decrease in the rate below a dynamically important breakpoint (∼40 °C), suggesting an equilibrium shift to an impaired conformational landscape relevant for catalysis in the low-temperature region. Full Marcus-like model fits of the rate and KIE profiles provided evidence for a high-temperature reaction via low-frequency conformational sampling associated with a broad distribution of hydride donor–acceptor distances (long-distance population centered at 3.31 ± 0.02 Å), only poorly suitable for quantum mechanical tunneling. Collectively, dynamical characteristics of TaCPa2E-catalyzed hydride transfer during transient oxidation of CDP-glucose reveal important analogies to mechanistically simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase. A loose-fit substrate (in TaCPa2E) resembles structural variants of these enzymes by extensive dynamical sampling to balance conformational flexibility and catalytic efficiency.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Temperature-Dependent Low-Frequency Modes in the Active Site of Bovine Carbonic Anhydrase II Probed by 2D-IR Spectroscopy

Enzyme catalysis achieves tremendous rate accelerations. Enzyme reaction centers provide a constraint geometry that preferentially binds an activated form of the substrate and thus lowers the energy barrier. However, this transition state picture neglects the flexibility of proteins and its role in enzymatic catalysis. Especially for proton transfer reactions, it has been suggested that motions of the protein modulate the donor-acceptor distance and prepare a tunneling-ready state. We report the detection of frequency fluctuations of an azide anion (N₃^-) bound in the active site of the protein carbonic anhydrase II, where a low-frequency mode of the protein has been proposed to facilitate proton transfer over two water molecules during the catalyzed reaction. 2D-IR spectroscopy resolves an underdamped low-frequency mode at about 1 THz (30 cm^-1). We find its frequency to be viscosity- and temperature-dependent and to decrease by 6 cm^-1 between 230 and 320 K, reporting the softening of the mode's potential.

X-ray free-electron laser studies reveal correlated motion during isopenicillin N synthase catalysis

Isopenicillin N synthase (IPNS) catalyzes the unique reaction of l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) with dioxygen giving isopenicillin N (IPN), the precursor of all natural penicillins and cephalosporins. X-ray free-electron laser studies including time-resolved crystallography and emission spectroscopy reveal how reaction of IPNS:Fe(II):ACV with dioxygen to yield an Fe(III) superoxide causes differences in active site volume and unexpected conformational changes that propagate to structurally remote regions. Combined with solution studies, the results reveal the importance of protein dynamics in regulating intermediate conformations during conversion of ACV to IPN. The results have implications for catalysis by multiple IPNS-related oxygenases, including those involved in the human hypoxic response, and highlight the power of serial femtosecond crystallography to provide insight into long-range enzyme dynamics during reactions presently impossible for nonprotein catalysts.

Non-Coordinative Binding of O2 at the Active Center of a Copper-Dependent Enzyme

Molecular oxygen (O₂ ) is a sustainable oxidation reagent. O₂ is strongly oxidizing but kinetically stable and its final reaction product is water. For these reasons learning how to activate O₂ and how to steer its reactivity along desired reaction pathways is a longstanding challenge in chemical research.^[1] Activation of ground-state diradical O₂ can occur either via conversion to singlet oxygen or by one-electron reduction to superoxide. Many enzymes facilitate activation of O₂ by direct fomation of a metal-oxygen coordination complex concomitant with inner sphere electron transfer. The formylglycine generating enzyme (FGE) is an unusual mononuclear copper enzyme that appears to follow a different strategy. Atomic-resolution crystal structures of the precatalytic complex of FGE demonstrate that this enzyme binds O₂ juxtaposed, but not coordinated to the catalytic Cu^I . Isostructural complexes that contain Ag^I instead of Cu^I or nitric oxide instead of O₂ confirm that formation of the initial oxygenated complex of FGE does not depend on redox activity. A stepwise mechanism that decouples binding and activation of O₂ is unprecedented for metal-dependent oxidases, but is reminiscent of flavin-dependent enzymes.

Identification of Thermal Conduits That Link the Protein-Water Interface to the Active Site Loop and Catalytic Base in Enolase

We report here on the salient role of protein mobility in accessing conformational landscapes that enable efficient enzyme catalysis. We are focused on yeast enolase, a highly conserved lyase with a TIM barrel domain and catalytic loop, as part of a larger study of the relationship of site selective protein motions to chemical reactivity within superfamilies. Enthalpically hindered variants were developed by replacement of a conserved hydrophobic side chain (Leu 343) with smaller side chains. Leu343 is proximal to the active site base in enolase, and comparative pH rate profiles for the valine and alanine variants indicate a role for side chain hydrophobicity in tuning the pK_a of the catalytic base. However, the magnitude of a substrate deuterium isotope effect is almost identical for wild-type (WT) and Leu343Ala, supporting an unchanged rate-determining proton abstraction step. The introduced hydrophobic side chains at position 343 lead to a discontinuous break in both activity and activation energy as a function of side chain volume. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments were performed as a function of time and temperature for WT and Leu343Ala, and provide a spatially resolved map of changes in protein flexibility following mutation. Impacts on protein flexibility are localized to specific networks that arise at the protein-solvent interface and terminate in a loop that has been shown by X-ray crystallography to close over the active site. These interrelated effects are discussed in the context of long-range, solvent-accessible and thermally activated networks that play key roles in tuning the precise distances and interactions among reactants.

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.

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.

Hydrogen-Deuterium Exchange within Adenosine Deaminase, a TIM Barrel Hydrolase, Identifies Networks for Thermal Activation of Catalysis

Proteins are intrinsically flexible macromolecules that undergo internal motions with time scales spanning femtoseconds to milliseconds. These fluctuations are implicated in the optimization of reaction barriers for enzyme catalyzed reactions. Time, temperature, and mutation dependent hydrogen-deuterium exchange coupled to mass spectrometry (HDX-MS) has been previously employed to identify spatially resolved, catalysis-linked dynamical regions of enzymes. We now extend this technique to pursue the correlation of protein flexibility and chemical reactivity within the diverse and widespread TIM barrel proteins, targeting murine adenosine deaminase (mADA) that catalyzes the irreversible deamination of adenosine to inosine and ammonia. Following a structure-function analysis of rate and activation energy for a series of mutations at a second sphere phenylalanine positioned in proximity to the bound substrate, the catalytically impaired Phe61Ala with an elevated activation energy (Ea = 7.5 kcal/mol) and the wild type (WT) mADA (Ea = 5.0 kcal/mol) were selected for HDX-MS experiments. The rate constants and activation energies of HDX for peptide segments are quantified and used to assess mutation-dependent changes in local and distal motions. Analyses reveal that approximately 50% of the protein sequence of Phe61Ala displays significant changes in the temperature dependence of HDX behaviors, with the dominant change being an increase in protein flexibility. Utilizing Phe61Ile, which displays the same activation energy for k_cat as WT, as a control, we were able to further refine the HDX analysis, highlighting the regions of mADA that are altered in a functionally relevant manner. A map is constructed that illustrates the regions of protein that are proposed to be essential for the thermal optimization of active site configurations that dominate reaction barrier crossings in the native enzyme.

The dynamic nature of the Mre11-Rad50 DNA break repair complex

The Mre11-Rad50-Nbs1/Xrs2 protein complex plays a pivotal role in the detection and repair of DNA double strand breaks. Through traditional and emerging structural biology techniques, various functional structural states of this complex have been visualized; however, relatively little is known about the transitions between these states. Indeed, it is these structural transitions that are important for Mre11-Rad50-mediated DNA unwinding at a break and the activation of downstream repair signaling events. Here, we present a brief overview of the current understanding of the structure of the core Mre11-Rad50 complex. We then highlight our recent studies emphasizing the contributions of solution state NMR spectroscopy and other biophysical techniques in providing insight into the structures and dynamics associated with Mre11-Rad50 functions.

Biological concepts for catalysis and reactivity: empowering bioinspiration

Biological systems provide attractive reactivity blueprints for the design of challenging chemical transformations. Emulating the operating mode of natural systems may however not be so easy and direct translation of structural observations does not always afford the anticipated efficiency. Metalloenzymes rely on earth-abundant metals to perform an incredibly wide range of chemical transformations. To do so, enzymes in general have evolved tools and tricks to enable control of such reactivity. The underlying concepts related to these tools are usually well-known to enzymologists and bio(inorganic) chemists but may be a little less familiar to organometallic chemists. So far, the field of bioinspired catalysis has greatly focused on the coordination sphere and electronic effects for the design of functional enzyme models but might benefit from a paradigm shift related to recent findings in biological systems. The goal of this review is to bring these fields closer together as this could likely result in the development of a new generation of highly efficient bioinspired systems. This contribution covers the fields of redox-active ligands, entatic state reactivity, energy conservation through electron bifurcation, and quantum tunneling for C-H activation.

High-Resolution XFEL Structure of the Soluble Methane Monooxygenase Hydroxylase Complex with its Regulatory Component at Ambient Temperature in Two Oxidation States

Soluble methane monooxygenase (sMMO) is a multicomponent metalloenzyme that catalyzes the conversion of methane to methanol at ambient temperature using a nonheme, oxygen-bridged dinuclear iron cluster in the active site. Structural changes in the hydroxylase component (sMMOH) containing the diiron cluster caused by complex formation with a regulatory component (MMOB) and by iron reduction are important for the regulation of O2 activation and substrate hydroxylation. Structural studies of metalloenzymes using traditional synchrotron-based X-ray crystallography are often complicated by partial X-ray-induced photoreduction of the metal center, thereby obviating determination of the structure of the enzyme in pure oxidation states. Here, microcrystals of the sMMOH:MMOB complex from Methylosinus trichosporium OB3b were serially exposed to X-ray free electron laser (XFEL) pulses, where the ≤35 fs duration of exposure of an individual crystal yields diffraction data before photoreduction-induced structural changes can manifest. Merging diffraction patterns obtained from thousands of crystals generates radiation damage-free, 1.95 Å resolution crystal structures for the fully oxidized and fully reduced states of the sMMOH:MMOB complex for the first time. The results provide new insight into the manner by which the diiron cluster and the active site environment are reorganized by the regulatory protein component in order to enhance the steps of oxygen activation and methane oxidation. This study also emphasizes the value of XFEL and serial femtosecond crystallography (SFX) methods for investigating the structures of metalloenzymes with radiation sensitive metal active sites.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Enzyme-Substrate-Cofactor Dynamical Networks Revealed by High-Resolution Field Cycling Relaxometry

The remarkable power and specificity of enzyme catalysis rely on the dynamic alignment of the enzyme, substrates, and cofactors, yet the role of dynamics has usually been approached from the perspective of the protein. We have been using an underappreciated NMR technique, subtesla high-resolution field cycling 31P NMR relaxometry, to investigate the dynamics of the enzyme-bound substrates and cofactor on guanosine-5'-monophosphate reductase (GMPR). GMPR forms two dead end, yet catalytically competent, complexes that mimic distinct steps in the catalytic cycle: E·IMP·NADP+ undergoes a partial hydride transfer reaction, while E·GMP·NADP+ undergoes a partial deamination reaction. A different cofactor conformation is required for each partial reaction. Here we report the effects of mutations designed to perturb cofactor conformation and ammonia binding with the goal of identifying the structural features that contribute to the distinct dynamic signatures of the hydride transfer and deamination complexes. These experiments suggest that Asp129 is a central cog in a dynamic network required for both hydride transfer and deamination. In contrast, Lys77 modulates the conformation and mobility of substrates and cofactors in a reaction-specific manner. Thr105 and Tyr318 are part of a deamination-specific dynamic network that includes the 2'-OH of GMP. These residues have comparatively little effect on the dynamic properties of the hydride transfer complex. These results further illustrate the potential of high-resolution field cycling NMR relaxometry for the investigation of ligand dynamics. In addition, exchange experiments indicate that NH3/NH4+ has a high affinity for the deamination complex but a low affinity for the hydride transfer complex, suggesting that the movement of ammonia may gate the cofactor conformational change. Collectively, these experiments reinforce the view that the enzyme, substrates, and cofactor are linked in intricate, reaction-specific, dynamic networks and demonstrate that distal portions of the substrates and cofactors are critical features in these networks.

Substitutions of Amino Acid Residues in the Substrate Binding Site of Horse Liver Alcohol Dehydrogenase Have Small Effects on the Structures but Significantly Affect Catalysis of Hydrogen Transfer

Previous studies showed that the L57F and F93W alcohol dehydrogenases catalyze the oxidation of benzyl alcohol with some quantum mechanical hydrogen tunneling, whereas the V203A enzyme has diminished tunneling. Here, steady-state kinetics for the L57F and F93W enzymes were studied, and microscopic rate constants for the ordered bi-bi mechanism were estimated from simulations of transient kinetics for the S48T, F93A, S48T/F93A, F93W, and L57F enzymes. Catalytic efficiencies for benzyl alcohol oxidation (V₁/E_tK_b) vary over a range of ∼100-fold for the less active enzymes up to the L57F enzyme and are mostly associated with the binding of alcohol rather than the rate constants for hydride transfer. In contrast, catalytic efficiencies for benzaldehyde reduction (V₂/E_tK_p) are ∼500-fold higher for the L57F enzyme than for the less active enzymes and are mostly associated with the rate constants for hydride transfer. Atomic-resolution structures (1.1 Å) for the F93W and L57F enzymes complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol are almost identical to previous structures for the wild-type, S48T, and V203A enzymes. Least-squares refinement with SHELXL shows that the nicotinamide ring is slightly strained in all complexes and that the apparent donor-acceptor distances from C4N of NAD to C7 of pentafluorobenzyl alcohol range from 3.28 to 3.49 Å (±0.02 Å) and are not correlated with the rate constants for hydride transfer or hydrogen tunneling. How the substitutions affect the dynamics of reorganization during hydrogen transfer and the extent of tunneling remain to be determined.

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.

Detecting and Characterizing the Kinetic Activation of Thermal Networks in Proteins: Thermal Transfer from a Distal, Solvent-Exposed Loop to the Active Site in Soybean Lipoxygenase

The rate-limiting chemical reaction catalyzed by soybean lipoxygenase (SLO) involves quantum mechanical tunneling of a hydrogen atom from substrate to its active site ferric-hydroxide cofactor. SLO has emerged as a prototypical system for linking the thermal activation of a protein scaffold to the efficiency of active site chemistry. Significantly, hydrogen-deuterium exchange-mass spectrometry (HDX-MS) experiments on wild type and mutant forms of SLO have uncovered trends in the enthalpic barriers for HDX within a solvent-exposed loop (positions 317-334) that correlate well with trends in the corresponding enthalpic barriers for kcat. A model for this behavior posits that collisions between water and loop 317-334 initiate thermal activation at the protein surface that is then propagated 15-34 Å inward toward the reactive carbon of substrate in proximity to the iron catalyst. In this study, we have prepared protein samples containing cysteine residues either at the tip of the loop 317-334 (Q322C) or on a control loop, 586-603 (S596C). Chemical modification of cysteines with the fluorophore 6-bromoacetyl-2-dimethylaminonaphthalene (Badan, BD) provides site-specific probes for the measurement of fluorescence relaxation lifetimes and Stokes shift decays as a function of temperature. Computational studies indicate that surface water structure is likely to be largely preserved in each sample. While both loops exhibit temperature-independent fluorescence relaxation lifetimes as do the Stokes shifts for S596C-BD, the activation enthalpy for the nanosecond solvent reorganization at Q322C-BD (Ea(ksolv) = 2.8(0.9) kcal/mol)) approximates the enthalpy of activation for catalytic C-H activation (Ea(kcat) = 2.3(0.4) kcal/mol). This study establishes and validates the methodology for measuring rates of rapid local motions at the protein/solvent interface of SLO. These new findings, when combined with previously published correlations between protein motions and the rate-limiting hydride transfer in a thermophilic alcohol dehydrogenase, provide experimental evidence for thermally induced "protein quakes" as the origin of enthalpic barriers in catalysis.

High-resolution crystal structure of human asparagine synthetase enables analysis of inhibitor binding and selectivity

Expression of human asparagine synthetase (ASNS) promotes metastatic progression and tumor cell invasiveness in colorectal and breast cancer, presumably by altering cellular levels of L-asparagine. Human ASNS is therefore emerging as a bona fide drug target for cancer therapy. Here we show that a slow-onset, tight binding inhibitor, which exhibits nanomolar affinity for human ASNS in vitro, exhibits excellent selectivity at 10 μM concentration in HCT-116 cell lysates with almost no off-target binding. The high-resolution (1.85 Å) crystal structure of human ASNS has enabled us to identify a cluster of negatively charged side chains in the synthetase domain that plays a key role in inhibitor binding. Comparing this structure with those of evolutionarily related AMP-forming enzymes provides insights into intermolecular interactions that give rise to the observed binding selectivity. Our findings demonstrate the feasibility of developing second generation human ASNS inhibitors as lead compounds for the discovery of drugs against metastasis.

Wen Zhu et al. report the crystal structure of human asparagine synthetase at a 1.85 Å resolution, enabling computational analysis of inhibitor binding. They also find new insights into the intermolecular interactions contributing to binding specificity of inhibitors.

The Structural Dynamics of Engineered β-Lactamases Vary Broadly on Three Timescales yet Sustain Native Function

Understanding the principles of protein dynamics will help guide engineering of protein function: altering protein motions may be a barrier to success or may be an enabling tool for protein engineering. The impact of dynamics on protein function is typically reported over a fraction of the full scope of motional timescales. If motional patterns vary significantly at different timescales, then only by monitoring motions broadly will we understand the impact of protein dynamics on engineering functional proteins. Using an integrative approach combining experimental and in silico methodologies, we elucidate protein dynamics over the entire span of fast to slow timescales (ps to ms) for a laboratory-engineered system composed of five interrelated β-lactamases: two natural homologs and three laboratory-recombined variants. Fast (ps-ns) and intermediate (ns-µs) dynamics were mostly conserved. However, slow motions (µs-ms) were few and conserved in the natural homologs yet were numerous and widely dispersed in their recombinants. Nonetheless, modified slow dynamics were functionally tolerated. Crystallographic B-factors from high-resolution X-ray structures were partly predictive of the conserved motions but not of the new slow motions captured in our solution studies. Our inspection of protein dynamics over a continuous range of timescales vividly illustrates the complexity of dynamic impacts of protein engineering as well as the functional tolerance of an engineered enzyme system to new slow motions.

On the use of noncompetitive kinetic isotope effects to investigate flavoenzyme mechanism

This account describes the application of kinetic isotope effects (KIEs) to investigate the mechanistic properties of flavin dependent enzymes. Assays can be conducted during steady-state catalytic turnover of the flavoenzyme with its substrate or by using rapid-kinetic techniques to measure either the reductive or oxidative half-reactions of the enzyme. Great care should be taken to ensure that the observed effects are due to isotopic substitution and not other factors such as pH effects or changes in the solvent viscosity of the reaction mixture. Different types of KIEs are described along with a physical description of their origins and the unique information each can provide about the mechanism of an enzyme. Detailed experimental techniques are outlined with special emphasis on the proper controls and data analysis that must be carried out to avoid erroneous conclusions. Examples are provided for each type of KIE measurement from references in the literature. It is our hope that this article will clarify any confusion concerning the utility of KIEs in the study of flavoprotein mechanism and encourage their use by the community.

Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example

Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.

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.

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 dynamics promote hydride tunnelling in substrate oxidation by aryl-alcohol oxidase

The temperature dependence of hydride transfer from the substrate to the N5 of the FAD cofactor during the reductive half-reaction of Pleurotus eryngii aryl-alcohol oxidase (AAO) is assessed here. Kinetic isotope effects on both the pre-steady state reduction of the enzyme and its steady-state kinetics, with differently deuterated substrates, suggest an environmentally-coupled quantum-mechanical tunnelling process. Moreover, those kinetic data, along with the crystallographic structure of the enzyme in complex with a substrate analogue, indicate that AAO shows a pre-organized active site that would only require the approaching of the hydride donor and acceptor for the tunnelled transfer to take place. Modification of the enzyme's active-site architecture by replacement of Tyr92, a residue establishing hydrophobic interactions with the substrate analogue in the crystal structure, in the Y92F, Y92L and Y92W variants resulted in different temperature dependence patterns that indicated a role of this residue in modulating the transfer reaction.

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.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

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.

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.