Role of dynamics in enzyme catalysis: substantial versus semantic controversies

A role for conformational changes in enzyme catalysis

The role played by conformational changes in enzyme catalysis is controversial. In addition to examining specific enzymes, studying formal models can help identify the conditions under which conformational changes promote catalysis. Here, we present a model demonstrating how conformational changes can break a generic trade-off due to the conflicting requirements of successive steps in catalytic cycles, namely high specificity for the transition state to accelerate the chemical transformation and low affinity for the products to favor their release. The mechanism by which the trade-off is broken is a transition between conformations with different affinities for the substrate. The role of the effector that induces the transition is played by a substrate "handle," a part of the substrate that is not chemically transformed but whose interaction with the enzyme is nevertheless essential to rapidly complete the catalytic cycle. A key element of the model is the formalization of the constraints causing the trade-off that the presence of multiple states breaks, which we attribute to the strong chemical similarity between successive reaction states-substrates, transition states, and products. For the sake of clarity, we present our model for irreversible one-step unimolecular reactions. In this context, we demonstrate how the different forms that chemical similarities between reaction states can take impose limits on the overall catalytic turnover. We first analyze catalysts without internal degrees of freedom and then show how two-state catalysts can overcome their limitations. Our results recapitulate previous proposals concerning the role of conformational changes and substrate handles in a formalism that makes explicit the constraints that elicit these features. In addition, our approach establishes links with studies in the field of heterogeneous catalysis, where the same trade-offs are observed and where overcoming them is a well-recognized challenge.

Comparison of methods for the NMR measurement of motionally averaged dipolar couplings

Motionally averaged dipolar couplings are an important tool for understanding the complex dynamics of catalysts, polymers, and biomolecules. While there is a plethora of solid-state NMR pulse sequences available for their measurement, in can be difficult to gauge the methods' strengths and weaknesses. In particular, there has not been a comprehensive comparison of their performance in natural abundance samples, where 1H homonuclear dipolar couplings are important and the use of large MAS rotors may be required for sensitivity reasons. In this work, we directly compared some of the more common methods for measuring C-H dipolar couplings in natural abundance samples using L-alanine (L-Ala) and the N-formyl-L-methionyl-L-leucyl-L-phenylalanine (fMLF) tripeptide as model systems. We evaluated their performance in terms of accuracy, resolution, sensitivity, and ease of implementation. We found that, despite the presence of 1H homonuclear dipolar interactions, all methods, with the exception of REDOR, were able to yield the reasonable dipolar coupling strengths for both mobile and static moieties. Of these methods, PDLF provides the most convenient workflow and precision at the expense of low sensitivity. In low-sensitivity cases, MAS-PISEMA and DIPSHIFT appear to be the better options.

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

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

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

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

Activation and friction in enzymatic loop opening and closing dynamics

Protein loop dynamics have recently been recognized as central to enzymatic activity, specificity and stability. However, the factors controlling loop opening and closing kinetics have remained elusive. Here, we combine molecular dynamics simulations with string-method determination of complex reaction coordinates to elucidate the molecular mechanism and rate-limiting step for WPD-loop dynamics in the PTP1B enzyme. While protein conformational dynamics is often represented as diffusive motion hindered by solvent viscosity and internal friction, we demonstrate that loop opening and closing is activated. It is governed by torsional rearrangement around a single loop peptide group and by significant friction caused by backbone adjustments, which can dynamically trap the loop. Considering both torsional barrier and time-dependent friction, our calculated rate constants exhibit very good agreement with experimental measurements, reproducing the change in loop opening kinetics between proteins. Furthermore, we demonstrate the applicability of our results to other enzymatic loops, including the M20 DHFR loop, thereby offering prospects for loop engineering potentially leading to enhanced designs.

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.

Structural Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions in Acetonitrile: Charge-Transfer Complex Tightness Is a Key

It has recently frequently been found that the kinetic isotope effect (KIE) is independent of temperature (T) in H-tunneling reactions in enzymes but becomes dependent on T in their mutants. Many enzymologists found that the trend is related to different donor-acceptor distances (DADs) at tunneling-ready states (TRSs), which could be sampled by protein dynamics. That is, a more rigid system of densely populated short DADs gives rise to a weaker T dependence of KIEs. Theoreticians have attempted to develop H-tunneling theories to explain the observations, but none have been universally accepted. It is reasonable to assume that the DAD sampling concept, if it exists, applies to the H-transfer reactions in solution, as well. In this work, we designed NADH/NAD+ model reactions to investigate their structural effects on the T dependence of hydride KIEs in acetonitrile. Hammett correlations together with N-CH3/CD3 secondary KIEs were used to provide the electronic structure of the TRSs and thus the rigidity of their charge-transfer complexation vibrations. In all three pairs of reactions, a weaker T dependence of KIEs always corresponds to a steeper Hammett slope on the substituted hydride acceptors. It was found that a tighter/rigid charge-transfer complexation system corresponds with a weaker T dependence of KIEs, consistent with the observations in enzymes.

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.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Ratcheting synthesis

Synthetic chemistry has traditionally relied on reactions between reactants of high chemical potential and transformations that proceed energetically downhill to either a global or local minimum (thermodynamic or kinetic control). Catalysts can be used to manipulate kinetic control, lowering activation energies to influence reaction outcomes. However, such chemistry is still constrained by the shape of one-dimensional reaction coordinates. Coupling synthesis to an orthogonal energy input can allow ratcheting of chemical reaction outcomes, reminiscent of the ways that molecular machines ratchet random thermal motion to bias conformational dynamics. This fundamentally distinct approach to synthesis allows multi-dimensional potential energy surfaces to be navigated, enabling reaction outcomes that cannot be achieved under conventional kinetic or thermodynamic control. In this Review, we discuss how ratcheted synthesis is ubiquitous throughout biology and consider how chemists might harness ratchet mechanisms to accelerate catalysis, drive chemical reactions uphill and programme complex reaction sequences.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

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.

Dynamic self-assembly of supramolecular catalysts from precision macromolecules

We show the emergence of strong catalytic activity at low concentrations in dynamic libraries of complementary sequence-defined oligomeric chains comprising pendant functional catalytic groups and terminal recognition units. In solution, the dynamic constitutional library created from pairs of such complementary oligomers comprises free oligomers, self-assembled di(oligomeric) macrocycles, and a virtually infinite collection of linear poly(oligomeric) chains. We demonstrate, on an exemplary catalytic system requiring the cooperation of no less than five chemical groups, that supramolecular di(oligomeric) macrocycles exhibit a catalytic turnover frequency ca. 20 times larger than the whole collection of linear poly(oligomers) and free chains. Molecular dynamics simulations and network analysis indicate that self-assembled supramolecular di(oligomeric) macrocycles are stabilized by different interactions, among which chain end pairing. We mathematically model the catalytic properties of such complex dynamic libraries with a small set of physically relevant parameters, which provides guidelines for the synthesis of oligomers capable to self-assemble into functionally-active supramolecular macrocycles over a larger range of concentrations.

How Flexibility Can Enhance Catalysis

Conformational changes are observed in many enzymes, but their role in catalysis is highly controversial. Here we present a theoretical model that illustrates how rigid catalysts can be fundamentally limited and how a conformational change induced by substrate binding can overcome this limitation, ultimately enabling barrier-free catalysis. The model is deliberately minimal, but the principle it illustrates is general and consistent with unique features of proteins as well as with previous informal proposals to explain the superiority of enzymes over other classes of catalysts. Implementing the discriminative switch suggested by the model could help overcome limitations currently encountered in the design of artificial catalysts.

Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces

The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.

Electron Tunneling in Biology: When Does it Matter?

Electrons can tunnel between cofactor molecules positioned along biological electron transport chains up to a distance of ≃ 20 Å on the millisecond time scale of enzymatic turnover. This tunneling range determines the design of biological energy chains facilitating the cross-membrane transport of electrons. Tunneling distance and cofactors' redox potentials become the main physical parameters affecting the rate of electron transport. In addition, universal charge-transport properties are assigned to all proteins, making protein identity, flexibility, and dynamics insignificant. This paradigm is challenged by dynamical models of electron transfer, showing that the electron hopping rate is constant within the crossover distance R* ≃ 12 Å, followed with an exponential falloff at longer distances. If this hypothesis is fully confirmed, natural and man-made energy chains for electron transport should be best designed by placing redox cofactors near the crossover distance R*. Protein flexibility and dynamics affect the magnitude of the maximum hopping rate within the crossover distance. Changes in protein flexibility between forward and backward transitions contribute to vectorial charge transport. For biological energy chains, charge transport through proteins is not defined by universal parameters, and protein identity matters.

Correlation of temperature dependence of hydride kinetic isotope effects with donor-acceptor distances in two solvents of different polarities

Recently observed nearly temperature (T)-independent kinetic isotope effects (KIEs) in wild-type enzymes and T-dependent KIEs in variants were used to suggest that H-tunneling in enzymes is assisted by the fast protein vibrations that help sample short donor-acceptor distances (DADs). This supports the recently proposed role of protein vibrations in DAD sampling catalysis. However, use of T-dependence of KIEs to suggest DAD sampling associated with protein vibrations is debated. We have formulated a hypothesis regarding the correlation and designed experiments in solution to investigate it. The hypothesis is, a more rigid system with shorter DADTRS's at the tunneling ready states (TRSs) gives rise to a weaker T-dependence of KIEs, i.e., a smaller ΔEa (= EaD - EaH). In a former work, the solvent effects of acetonitrile versus chloroform on the ΔEa of NADH/NAD+ model reactions were determined, and the DADPRC's of the productive reactant complexes (PRCs) were computed to substitute the DADTRS for the DADTRS-ΔEa correlation study. A smaller ΔEa was found in the more polar acetonitrile where the positively charged PRC is better solvated and has a shorter DADPRC, indirectly supporting the hypothesis. In this work, the TRS structures of different DADTRS's for the hydride tunneling reaction from 1,3-dimethyl-2-phenylimidazoline to 10-methylacridinium were computed. The N-CH3/CD3 secondary KIEs on both reactants were calculated and fitted to the observed values to find the DADTRS order in both solutions. It was found that the equilibrium DADTRS is shorter in acetonitrile than in chloroform. Results directly support the DADTRS-ΔEa correlation hypothesis as well as the explanation that links T-dependence of KIEs to DAD sampling catalysis in enzymes.

The Conformational Transitions and Dynamics of Burkholderia cepacia Lipase Regulated by Water-Oil Interfaces

Structural dynamics and conformational transitions are crucial for the activities of enzymes. As one of the most widely used industrial biocatalysts, lipase could be activated by the water-oil interfaces. The interface activations were believed to be dominated by the close-to-open transitions of the lid subdomains. However, the detailed mechanism and the roles of structure transitions are still under debate. In this study, the dynamic structures and conformational transitions of Burkholderia cepacia lipase (LipA) were investigated by combining all-atom molecular dynamics simulations, enhanced sampling simulation, and spectrophotometric assay experiments. The conformational transitions between the lid-open and lid-closed states of LipA in aqueous solution are directly observed by the computational simulation methods. The interactions between the hydrophobic residues on the two lid-subdomains are the driven forces for the LipA closing. Meanwhile, the hydrophobic environment provided by the oil interfaces would separate the interactions between the lid-subdomains and promote the structure opening of LipA. Moreover, our studies demonstrate the opening of the lids structure is insufficient to initiate the interfacial activation, providing explanations for the inability of interfacial activation of many lipases with lid structures.

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.

Profiling Enzyme Activity of l-Asparaginase II by NMR-Based Methyl Fingerprinting at Natural Abundance

l-asparaginase II (MW 135 kDa) from E. coli is an FDA-approved protein drug used for the treatment of childhood leukemia. Despite its long history as a chemotherapeutic, the structural basis of enzyme action, in solution, remains widely contested. In this work, methyl-based 2D [¹H-¹³C]-heteronuclear single-quantum correlation (HSQC) NMR, at natural abundance, has been used to profile the enzymatic activity of the commercially available enzyme drug. The [¹H-¹³C]-HSQC NMR spectra of the protein reveal the role of a flexible loop segment in the activity of the enzyme, in solution. Addition of asparagine to the protein results in distinct conformational changes of the loop that could be signatures of intermediates formed in the catalytic reaction. To this end, an isothermal titration calorimetry (ITC)-based assay has been developed to measure the enzymatic reaction enthalpy, as a marker for its activity. Combining both ITC and NMR, it was shown that the disruption of the protein conformation can result in the loss of function. The scope, robustness, and validity of the loop fingerprints in relation to enzyme activity have been tested under different solution conditions. Overall, our results indicate that 2D NMR can be used reliably to gauge the structure-function of this enzyme, bypassing the need to label the protein. Such natural abundant NMR methods can be potentially extended to probe the structure-function aspects of high-molecular-weight protein therapeutics (glycosylated protein drugs, enzymes, therapeutic monoclonal antibodies, antibody-drug conjugates, and Fc-fusion proteins), where (a) flexible loops are required for their function and (b) isotope labeling may not be straightforward.

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.

Enzyme Grafting with a Cofactor-Decorated Metal-Organic Capsule for Solar-to-Chemical Conversion

Semi-artificial approaches to solar-to-chemical conversion can achieve chemical transformations that are beyond the capability of natural enzymes, but face marked challenges to facilitate in vivo cascades, due to their inevitable need for cofactor shuttling and regeneration. Here, we report on an enzyme grafting strategy to build a metal-organic capsule-docking artificial enzyme (metal-organic-enzyme, MOE) that comprised the self-assembly of a cofactor-decorated capsule and the supramolecular enzyme-recognition features between the enzyme scaffold and the capsule to bypass cofactor shuttling and regeneration. The incorporated NADH mimics within the metal-organic capsule interacted with the imine intermediate that formed from the condensation of the amines and the dehydrogenation of alcohol substrates in the microenvironment to form complexes within the capsule and subsequently served as an in situ-generated photoresponsive cofactor. Upon illumination, the photoresponsive cofactor facilitates efficient proton/electron transport between the inner space (supramolecular hydrogenation) and outer space (enzymatic dehydrogenation) of the capsule to dehydrogenize the alcohols and hydrogenize the imine intermediates, respectively, circumventing the conventionally complex multistep cofactor shuttling and regeneration. The semi-artificial enzyme endows the conversion of diverse types of alcohol to amine products in both aqueous/organic solutions and Escherichia coli with high efficiency, offering a wide range of opportunities for sustainable and environmentally friendly biomanufacturing of commodity and fine chemicals.

1H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?

Dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular ²H quadrupolar coupling and ¹H-X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of ¹H NMR, ¹H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out ¹H CSA and ¹H-¹³C dipolar coupling measurements in a series of model supported complexes to understand how ¹H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, ¹H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.

An engineered variant of MECR reductase reveals indispensability of long-chain acyl-ACPs for mitochondrial respiration

Mitochondrial fatty acid synthesis (mtFAS) is essential for respiratory function. MtFAS generates the octanoic acid precursor for lipoic acid synthesis, but the role of longer fatty acid products has remained unclear. The structurally well-characterized component of mtFAS, human 2E-enoyl-ACP reductase (MECR) rescues respiratory growth and lipoylation defects of a Saccharomyces cerevisiae Δetr1 strain lacking native mtFAS enoyl reductase. To address the role of longer products of mtFAS, we employed in silico molecular simulations to design a MECR variant with a shortened substrate binding cavity. Our in vitro and in vivo analyses indicate that the MECR G165Q variant allows synthesis of octanoyl groups but not long chain fatty acids, confirming the validity of our computational approach to engineer substrate length specificity. Furthermore, our data imply that restoring lipoylation in mtFAS deficient yeast strains is not sufficient to support respiration and that long chain acyl-ACPs generated by mtFAS are required for mitochondrial function.

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.

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

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

Dependence of crystallographic atomic displacement parameters on temperature (25-150 K) for complexes of horse liver alcohol dehydrogenase

Enzymes catalyze reactions by binding and orienting substrates with dynamic interactions. Horse liver alcohol dehydrogenase catalyzes hydrogen transfer with quantum-mechanical tunneling that involves fast motions in the active site. The structures and B factors of ternary complexes of the enzyme with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol or NAD+ and 2,2,2-trifluoroethanol were determined to 1.1-1.3 Å resolution below the `glassy transition' in order to extract information about the temperature-dependent harmonic motions, which are reflected in the crystallographic B factors. The refinement statistics and structures are essentially the same for each structure at all temperatures. The B factors were corrected for a small amount of radiation decay. The overall B factors for the complexes are similar (13-16 Å2) over the range 25-100 K, but increase somewhat at 150 K. Applying TLS refinement to remove the contribution of pseudo-rigid-body displacements of coenzyme binding and catalytic domains provided residual B factors of 7-10 Å2 for the overall complexes and of 5-10 Å2 for C4N of NAD+ and the methylene carbon of the alcohols. These residual B factors have a very small dependence on temperature and include local harmonic motions and apparently contributions from other sources. Structures at 100 K show complexes that are poised for hydrogen transfer, which involves atomic displacements of ∼0.3 Å and is compatible with the motions estimated from the residual B factors and molecular-dynamics simulations. At 298 K local conformational changes are also involved in catalysis, as enzymes with substitutions of amino acids in the substrate-binding site have similar positions of NAD+ and pentafluorobenzyl alcohol and similar residual B factors, but differ by tenfold in the rate constants for hydride transfer.

Conformational dynamics modulating electron transfer

Diffusional dynamics of the donor-acceptor distance are responsible for the appearance of a new time scale of diffusion over the distance of electronic tunneling in electron-transfer reactions. The distance dynamics compete with the medium polarization dynamics in the dynamics-controlled electron-transfer kinetics. The pre-exponential factor of the electron-transfer rate constant switches, at the crossover distance, between a distance-independent, dynamics-controlled plateau and exponential distance decay. The crossover between two regimes is controlled by an effective relaxation time slowed down by a factor exponentially depending on the variance of the donor-acceptor displacement. Flexible donor-acceptor complexes must show a greater tendency for dynamics-controlled electron transfer. Energy chains based on electron transport are best designed by placing the redox cofactors near the crossover distance.

Molecular Insights into the Intrinsic Dynamics and Their Roles During Catalysis in Pin1 Peptidyl-prolyl Isomerase

Proteins are intrinsically dynamic and change conformations over a wide range of time scales. While the conformational dynamics have been realized to be important for protein functions, e.g., in activity-stability trade-offs, how they play a role during enzyme catalysis has been of debate over decades. By studying Pin1 peptidyl-prolyl isomerase using extensive molecular dynamics simulations, here we discuss how the slow intrinsic dynamics of Pin1 observed in the NMR relaxation dispersion experiment occur and couple to isomerization reactions in molecular detail. In particular, we analyze the angular correlation functions of the backbone N-H bonds and find that slow conformational transitions occur at about the 3₁₀ helix in the apo state. These events at the helical region further affect the residues at about the ligand binding site. Unfolding of this helix leads to a tight hydrogen bond between the helical region and the ligand binding loop, thus forming a stable coiled structure. The helical and coiled structures are found to be characteristic of the Pin1-ligand complex with the ligand in the trans and cis states, respectively. These results indicate that the changes in the slow dynamics of Pin1 by the isomerization reaction occur via the shift in populations of the helical and coiled states, where the balance is dependent on the ligand isomerization states.

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.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

Thylakoid Membrane-Inspired Capsules with Fortified Cofactor Shuttling for Enzyme-Photocoupled Catalysis

Enzyme-photocoupled catalytic systems (EPCSs), combining the natural enzyme with a library of semiconductor photocatalysts, may break the constraint of natural evolution, realizing sustainable solar-to-chemical conversion and non-natural reactivity of the enzyme. The overall efficiency of EPCSs strongly relies on the shuttling of energy-carrying molecules, e.g., NAD⁺/NADH cofactor, between active centers of enzyme and photocatalyst. However, few efforts have been devoted to NAD⁺/NADH shuttling. Herein, we propose a strategy of constructing a thylakoid membrane-inspired capsule (TMC) with fortified and tunable NAD⁺/NADH shuttling to boost the enzyme-photocoupled catalytic process. The apparent shuttling number (ASN) of NAD⁺/NADH for TMC could reach 17.1, ∼8 times as high as that of non-integrated EPCS. Accordingly, our TMC exhibits a turnover frequency (TOF) of 38 000 ± 365 h^-1 with a solar-to-chemical efficiency (STC) of 0.69 ± 0.12%, ∼6 times higher than that of non-integrated EPCS.

Mapping Enzyme Landscapes by Time-Resolved Crystallography with Synchrotron and X-Ray Free Electron Laser Light

Directly observing enzyme catalysis in real time at the molecular level has been a long-standing goal of structural enzymology. Time-resolved serial crystallography methods at synchrotron and X-ray free electron laser (XFEL) sources have enabled researchers to follow enzyme catalysis and other nonequilibrium events at ambient conditions with unprecedented time resolution. X-ray crystallography provides detailed information about conformational heterogeneity and protein dynamics, which is enhanced when time-resolved approaches are used. This review outlines the ways in which information about the underlying energy landscape of a protein can be extracted from X-ray crystallographic data, with an emphasis on new developments in XFEL and synchrotron time-resolved crystallography. The emerging view of enzyme catalysis afforded by these techniques can be interpreted as enzymes moving on a time-dependent energy landscape. Some consequences of this view are discussed, including the proposal that irreversible enzymes or enzymes that use covalent catalytic mechanisms may commonly exhibit catalysis-activated motions. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

High throughput and quantitative enzymology in the genomic era

Accurate predictions from models based on physical principles are the ultimate metric of our biophysical understanding. Although there has been stunning progress toward structure prediction, quantitative prediction of enzyme function has remained challenging. Realizing this goal will require large numbers of quantitative measurements of rate and binding constants and the use of these ground-truth data sets to guide the development and testing of these quantitative models. Ground truth data more closely linked to the underlying physical forces are also desired. Here, we describe technological advances that enable both types of ground truth measurements. These advances allow classic models to be tested, provide novel mechanistic insights, and place us on the path toward a predictive understanding of enzyme structure and function.

Dynamics and Entropy of Cyclohexane Rings Control pH-Responsive Reactivity

Activation entropy (ΔS ^‡) is not normally considered the main factor in determining the reactivity of unimolecular reactions. Here, we report that the intramolecular degradation of six-membered ring compounds is mainly determined by the ΔS ^‡, which is strongly influenced by the ring-flipping motion and substituent geometry. Starting from the unique difference between the pH-dependent degradation kinetics of geometric isomers of 1,2-cyclohexanecarboxylic acid amide (1,2-CHCAA), where only the cis isomer can readily degrade under weakly acidic conditions (pH < 5.5), we found that the difference originated from the large difference in ΔS ^‡ of 16.02 cal·mol^-1·K^-1. While cis-1,2-CHCAA maintains a preference for the classical chair cyclohexane conformation, trans-1,2-CHCAA shows dynamic interconversion between the chair and twisted boat conformations, which was supported by both MD simulations and VT-NMR analysis. Steric repulsion between the bulky 1,2-substituents of the trans isomer is one of the main reasons for the reduced energy barrier between ring conformations that facilitates dynamic ring inversion motions. Consequently, the more dynamic trans isomer exhibits much a larger loss in entropy during the activation process due to the prepositioning of the reactant than the cis isomer, and the pH-dependent degradation of the trans isomer is effectively suppressed. When the ring inversion motion is inhibited by an additional methyl substituent on the cyclohexane ring, the pH degradability can be dramatically enhanced for even the trans isomer. This study shows a unique example in which spatial arrangement and dynamic properties can strongly influence molecular reactivity in unimolecular reactions, and it will be helpful for the future design of a reactive structure depending on dynamic conformational changes.

Probing Ligand Effects on the Ultrafast Dynamics of Copper Complexes via Midinfrared Pump-Probe and 2DIR Spectroscopies

The effects of ligand structural variation on the ultrafast dynamics of a series of copper coordination complexes were investigated using polarization-dependent mid-IR pump-probe spectroscopy and two-dimensional infrared (2DIR) spectroscopy. The series consists of three copper complexes [(^R3P₃tren)Cu^IIN₃]BAr₄^F (1^PR3, ^R3P₃tren = tris[2-(phosphiniminato)ethyl]amine, BAr₄^F = tetrakis(pentafluorophenyl)borate) where the number of methyl and phenyl groups in the PR₃ ligand are systematically varied across the series (PR₃ = PMe₃, PMe₂Ph, PMePh₂). The asymmetric stretching mode of azide in the 1^PR3 series is used as a vibrational probe of the small-molecule binding site. The results of the pump-probe measurements indicate that the vibrational energy of azide dissipates through intramolecular pathways and that the bulkier phenyl groups lead to an increase in the spatial restriction of the diffusive reorientation of bound azide. From 2DIR experiments, we characterize the spectral diffusion of the azide group and find that an increase in the number of phenyl groups maps to a broader inhomogeneous frequency distribution (Δ₂). This indicates that an increase in the steric bulk of the secondary coordination sphere acts to create more distinct configurations in the local environment that are accessible to the azide group. This work demonstrates how ligand structural variation affects the ultrafast dynamics of a small molecular group bound to the metal center, which could provide insight into the structure-function relationship of the copper coordination complexes and transition-metal coordination complexes in general.

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.

Molecular dynamics evidence for nonthermal effects of electric fields on pectin methylesterase activity

Experimental studies relevant to the nonthermal effects of electric fields on biological systems are emerging. However, these effects are poorly understood at the molecular level. The present study investigates pectin methylesterase, a cell wall modifying enzyme in plants, exposed to various electric field strengths. Molecular dynamics (MD) of the enzyme were studied with and without (thermal-only) electric field applications. The measurements were interpreted on the basis of equivalent energy input to gain insights into the effect of electric field treatment time at a constant temperature (50 °C). Results reveal that electric fields exert nonthermal effects on both local and global protein structure. In 1 μs simulations, the results show significant (P ≤ 0.05) shrinkage of the catalytic domain and shortening of enzyme-water hydrogen bond lifetime by a 50 V cm-1 electric field. Unwinding of the helical segments, altered intra- and intermolecular hydrogen bond patterns, and increased hydration are also caused by the 50 V cm-1 electric field. This study serves to understand the electric field influence on the functional role of proteins.

The structure of an archaeal oligosaccharyltransferase provides insight into the strict exclusion of proline from the N-glycosylation sequon

Oligosaccharyltransferase (OST) catalyzes oligosaccharide transfer to the Asn residue in the N-glycosylation sequon, Asn-X-Ser/Thr, where Pro is strictly excluded at position X. Considering the unique structural properties of proline, this exclusion may not be surprising, but the structural basis for the rejection of Pro residues should be explained explicitly. Here we determined the crystal structure of an archaeal OST in a complex with a sequon-containing peptide and dolichol-phosphate to a 2.7 Å resolution. The sequon part in the peptide forms two inter-chain hydrogen bonds with a conserved amino acid motif, TIXE. We confirmed the essential role of the TIXE motif and the adjacent regions by extensive alanine-scanning of the external loop 5. A Ramachandran plot revealed that the ring structure of the Pro side chain is incompatible with the ϕ backbone dihedral angle around -150° in the rigid sequon-TIXE structure. The present structure clearly provides the structural basis for the exclusion of Pro residues from the N-glycosylation sequon.

Correlated Motions in Structural Biology

Correlated motions in proteins arising from the collective movements of residues have long been proposed to be fundamentally important to key properties of proteins, from allostery and catalysis to evolvability. Recent breakthroughs in structural biology have made it possible to capture proteins undergoing complex conformational changes, yet intrinsic correlated motions within a conformation remain one of the least understood facets of protein structure. For many decades, the analysis of total X-ray scattering held the promise of animating crystal structures with correlated motions. With recent advances in both X-ray detectors and data interpretation methods, this long-held promise can now be met. In this Perspective, we will introduce how correlated motions are captured in total scattering and provide guidelines for the collection, interpretation, and validation of data. As structural biology continues to push the boundaries, we see an opportunity to gain atomistic insight into correlated motions using total scattering as a bridge between theory and experiment.

Direct Correlation between Donor-Acceptor Distance and Temperature Dependence of Kinetic Isotope Effects in Hydride-Tunneling Reactions of NADH/NAD+ Analogues

Recent study of structural effects on primary kinetic isotope effects (1° KIEs) of H-transfer reactions in enzymes and solution revealed that a more rigid reaction system gave rise to a weaker temperature dependence of 1° KIEs, i.e., a smaller isotopic activation energy difference (ΔE_a = E_aD - E_aH). This has been explained within the contemporary vibrationally assisted activated H-tunneling (VA-AHT) model in which rigidity is defined according to the density of donor-acceptor distance (DAD_TRS) populations at the tunneling ready state (TRS) sampled by heavy atom motions. To test the relationship between DAD_TRS and ΔE_a in the model, we developed a computational method to obtain the TRS structures for H-transfer reactions. The method was applied to three hydride transfer reactions of NADH/NAD⁺ analogues for which the ΔE_a's as well as secondary (2°) KIEs have been reported. The 2° KIEs computed from each TRS structure were fitted to the observed values to obtain the optimal TRSs/DAD_TRS's. It was found that a shorter DAD_TRS does correspond with a smaller ΔE_a. This appears to support the VA-AHT model. Moreover, an analysis of hybridizations at the bent TRS structures shows that rehybridizations at the donor-acceptor centers are much more advanced than predicted from the classical mechanism. This implies that more orbital preparations are required for the nonclassical H-tunneling to take place.

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

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

Phase-sensitive γ-encoded recoupling of heteronuclear dipolar interactions and 1H chemical shift anisotropy

γ-encoded recoupling sequences are known to produce strong amplitude modulations that lead to sharp doublets when Fourier transformed. These doublets depend very little on the recoupled tensor asymmetry and thus enable for the straightforward determination of dynamic order parameters. It can, however, be difficult to measure small anisotropies, or small order parameters, using such sequences; the resonances from the doublet may overlap with each other, or with the zero-frequency glitch. This limitation has prevented the widespread use of ¹H chemical shift anisotropy (CSA) for the measurement of dynamics, particularly for CH protons which typically have CSAs of only a few ppm when immobile. Here, we introduce a simple modification to the traditional ¹H CSA and proton-detected local field pulse sequences that enables the acquisition of a hypercomplex dataset and the removal of the uncorrelated magnetization that results in the zero-frequency glitch. These new sequences then yield a frequency shift in the indirect dimension, rather than a splitting, which is easily identifiable even in cases of weak interactions.

Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts

As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.