Understanding Biological Hydrogen Transfer Through the Lens of Temperature Dependent Kinetic Isotope Effects

Thermal Rates and High-Temperature Tunneling from Surface Reaction Dynamics and First-Principles

Studying dynamics of the dissociative adsorption and recombinative desorption of hydrogen on copper surfaces has shaped our atomic-scale understanding of surface chemistry, yet experimentally determining the thermal rates for these processes, which dictate the outcome of catalytic reactions, has been impossible so far. In this work, we determine the thermal rate constants for dissociative adsorption and recombinative desorption of hydrogen on Cu(111) between 200 and 1000 K using data from reaction dynamics experiments. Contrary to current understanding, our findings demonstrate the predominant role of quantum tunneling, even at temperatures as high as 400 K. We also provide precise values for the reaction barrier (0.619 ± 0.020 eV) and adsorption energy (0.348 ± 0.026 eV) for H2 on Cu(111). Remarkably, the thermal rate constants are in excellent agreement with a first-principles quantum rate theory based on a new implementation of ring polymer molecular dynamics for reactions on surfaces, paving the way to discovering better catalysts using reliable and efficient computational methods.

Environment- and Conformation-Induced Frequency Shifts of C-D Vibrational Stark Probes in NAD(P)H Cofactors

NAD(P)H cofactors are found in all forms of life and are essential for electron and hydrogen atom transfer. The linear response of a carbon-deuterium (C-D) vibration based on the vibrational Stark effect can facilitate measurements of electric fields for critical biological reactions including cofactor-mediated hydride transfer. We find both inter- and intramolecular electric fields influence the C-D frequency in NAD(P)H and nicotinamide-like models where the reactive C4-hydrogen has been deuterated. Hence, the C-D frequency can report both environmental electrostatics and conformational changes of the nicotinamide ring. Conformation-dependent effects are mediated through space as electrostatic effects, rather than through-bond. A Stark tuning rate of ∼0.57 cm-1/(MV/cm) was determined using both experimental and computational approaches, including vibrational solvatochromism, molecular dynamics simulations, and in silico Stark calculations. The vibrational probe's Stark tuning rate is shown to be robust and suitable for measuring fields along hydride transfer reaction coordinates in enzymes.

Heterogeneous CuS QDs/BiVO4@Y2O2S Nanoreactor for Monitorable Photocatalysis

Exploration of multifunctional integrated catalysts is of great significance for photocatalysis toward practical application. Herein, a 1D confined nanoreactor with a heterogeneous core-shell structure is designed for synergies of efficient catalysis and temperature monitoring by custom encapsulation of Z-scheme heterojunction CuS quantum dots/BiVO4 (CuS QDs/BiVO4) and Y2O2S-Er, Yb. The dispersed active sites created by the QDs with high surface energy improve the mass transfer efficiency, and the efficient electron transport channels at the heterogeneous interface extend the carrier lifetime, which endows the nanoreactor with excellent catalytic performance. Meanwhile, real-time temperature monitoring is realized based on the thermally coupled levels 2H11/2/4S3/2→4I15/2 of Er3+ using fluorescence intensity ratio, which enables the monitorable photocatalysis. Furthermore, the nanoreactor with a multidimensional structure increases effective intermolecular collisions to facilitate the catalytic process by restricting the reaction within distinct enclosed spaces and circumvents potential unknown interaction effects. The design of multi-space nanoconfined reactors opens up a new avenue to modulate catalyst function, providing a unique perspective for photocatalytic applications in the mineralization of organic pollutants, hydrogen production, and nitrogen fixation.

Melanin: Nature's 4th bioorganic polymer

The pigments known as the melanins are widely recognized for their responsibility in the coloration of human skin, eyes, hair, and minimising the harmful effects of solar ultraviolet radiation. But specialists are aware that the melanins are present in all living kingdoms, barring viruses, and have functionality that extends beyond neutralizing ionising radiation. The ubiquitous presence of melanin in almost all human organs, recognized in recent years, as well as the presence of melanin in organisms that are evolutionarily distant from each other, indicate the fundamental importance of this class of material for all life forms. In this review, we argue for the need to accept melanins as the fourth primordial class of biological polymers, along with nucleic acids, proteins and polysaccharides. We consistently compare the properties of these canonical biological polymers with the properties of melanin and highlight key features that fundamentally distinguish melanins, their function and its mysteries.

Hydrogen Tunneling Exhibiting Unexpectedly Small Primary Kinetic Isotope Effects

Probing quantum mechanical tunneling (QMT) in chemical reactions is crucial to understanding and developing new transformations. Primary H/D kinetic isotopic effects (KIEs) beyond the semiclassical maximum values of 7-10 (room temperature) are commonly used to assess substantial QMT contributions in one-step hydrogen transfer reactions, because of the much greater QMT probability of protium vs. deuterium. Nevertheless, we report here the discovery of a reaction model occurring exclusively by H-atom QMT with residual primary H/D KIEs. 2-Hydroxyphenylnitrene, generated in N2 matrix, was found to isomerize to an imino-ketone via sequential (domino) QMT involving anti to syn OH-rotamerization (rate determining step) and [1,4]-H shift reactions. These sequential QMT transformations were also observed in the OD-deuterated sample, and unexpected primary H/D KIEs between 3 and 4 were measured at 3 to 20 K. Analogous residual primary H/D KIEs were found in the anti to syn OH-rotamerization QMT of 2-cyanophenol in a N2 matrix. Evidence strongly indicates that these intriguing isotope-insensitive QMT reactivities arise due to the solvation effects of the N2 matrix medium, putatively through coupling with the moving H/D tunneling particle. Should a similar scenario be extrapolated to conventional solution conditions, then QMT may have been overlooked in many chemical reactions.

Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions

Olive oil phenols are recognized as molecules with numerous positive health effects, many of which rely on their antioxidative activity, i.e., the ability to transfer hydrogen to radicals. Proton-coupled electron transfer reactions and hydrogen tunneling are ubiquitous in biological systems. Reactions of olive oil phenols, hydroxytyrosol, tyrosol, oleuropein, oleacein, oleocanthal, homovanillyl alcohol, vanillin, and a few phenolic acids with a DPPH• (2,2-diphenyl-1-picrylhydrazyl) radical in a 1,4-dioxane:water = 95:5 or 99:1 v/v solvent mixture were studied through an experimental kinetic analysis and computational chemistry calculations. The highest rate constants corresponding to the highest antioxidative activity are obtained for the ortho-diphenols hydroxytyrosol, oleuropein, and oleacein. The experimentally determined kinetic isotope effects (KIEs) for hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions are 16.0, 15.4, and 16.7, respectively. Based on these KIEs, thermodynamic activation parameters, and an intrinsic bond orbital (IBO) analysis along the IRC path calculations, we propose a proton-coupled electron transfer mechanism. The average local ionization energy and electron donor Fukui function obtained for the phenolic compounds show that the most reactive electron-donating sites are associated with π electrons above and below the aromatic ring, in support of the IBO analysis and proposed PCET reaction mechanism. Large KIEs and isotopic values of Arrhenius pre-exponential factor AH/AD determined for the hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions of 0.6, 1.3, and 0.3, respectively, reveal the involvement of hydrogen tunneling in the process.

Impact of N-Glycosylation on Protein Structure and Dynamics Linked to Enzymatic C-H Activation in the M. oryzae Lipoxygenase

Lipoxygenases (LOXs) from pathogenic fungi are potential therapeutic targets for defense against plant and select human diseases. In contrast to the canonical LOXs in plants and animals, fungal LOXs are unique in having appended N-linked glycans. Such important post-translational modifications (PTMs) endow proteins with altered structure, stability, and/or function. In this study, we present the structural and functional outcomes of removing or altering these surface carbohydrates on the LOX from the devastating rice blast fungus, M. oryzae, MoLOX. Alteration of the PTMs did notinfluence the active site enzyme-substrate ground state structures as visualized by electron-nuclear double resonance (ENDOR) spectroscopy. However, removal of the eight N-linked glycans by asparagine-to-glutamine mutagenesis nonetheless led to a change in substrate selectivity and an elevated activation energy for the reaction with substrate linoleic acid, as determined by kinetic measurements. Comparative hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis of wild-type and Asn-to-Gln MoLOX variants revealed a regionally defined impact on the dynamics of the arched helix that covers the active site. Guided by these HDX results, a single glycan sequon knockout was generated at position 72, and its comparative substrate selectivity from kinetics nearly matched that of the Asn-to-Gln variant. The cumulative data from model glyco-enzyme MoLOX showcase how the presence, alteration, or removal of even a single N-linked glycan can influence the structural integrity and dynamics of the protein that are linked to an enzyme's catalytic proficiency, while indicating that extensive glycosylation protects the enzyme during pathogenesis by protecting it from protease degradation.

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.

What triggers the coupling of proton transfer and electron transfer at the active site of nitrogenase?

In converting N2 to NH3 the enzyme nitrogenase utilises 8 electrons and 8 protons in the complete catalytic cycle. The source of the electrons is an Fe4S4 reductase protein (Fe-protein) which temporarily docks with the MoFe-protein that contains the catalytic active cofactor, FeMo-co, and an electron transfer cluster called the P cluster. The overall mechanism involves 8 repetitions of a cycle in which reduced Fe-protein docks with the MoFe-protein, one electron transfers to the P-cluster, and then to FeMo-co, followed by dissociation of the two proteins and re-reduction of the Fe-protein. Protons are supplied serially to FeMo-co by a Grotthuss proton translocation mechanism from the protein surface along a conserved chain of water molecules (a proton wire) that terminates near S atoms of the FeMo-co cluster [CFe7S9Mo(homocitrate)] where the multiple steps of the chemical conversions are effected. It is assumed that the chemical mechanisms use proton-coupled electron-transfer (PCET) and that H atoms (e- + H+) are involved in each of the hydrogenation steps. However there is neither evidence for, or mechanism proposed, for this coupling. Here I report calculations of cluster charge distribution upon electron addition, revealing that the added negative charge is on the S atoms of FeMo-co, which thereby become more basic, and able to trigger proton transfer from H3O+ waiting at the near end of the proton wire. This mechanism is supported by calculations of the dynamics of the proton transfer step, in which the barrier is reduced by ca. 3.5 kcal mol-1 and the product stabilised by ca. 7 kcal mol-1 upon electron addition. H tunneling is probable in this step. In nitrogenase it is electron transfer that triggers proton transfer.

C-H Bond Oxidation by MnIV-Oxo Complexes: Hydrogen-Atom Tunneling and Multistate Reactivity

The reactivity of six MnIV-oxo complexes in C-H bond oxidation has been examined using a combination of kinetic experiments and computational methods. Variable-temperature studies of the oxidation of 9,10-dihydroanthracene (DHA) and ethylbenzene by these MnIV-oxo complexes yielded activation parameters suitable for evaluating electronic structure computations. Complementary kinetic experiments of the oxidation of deuterated DHA provided evidence for hydrogen-atom tunneling in C-H bond oxidation for all MnIV-oxo complexes. These results are in accordance with the Bell model, where tunneling occurs near the top of the transition-state barrier. Density functional theory (DFT) and DLPNO-CCSD(T1) computations were performed for three of the six MnIV-oxo complexes to probe a previously predicted multistate reactivity model. The DFT computations predicted a thermal crossing from the 4B1 ground state to a 4E state along the C-H bond oxidation reaction coordinate. DLPNO-CCSD(T1) calculations further confirm that the 4E transition state offers a lower energy barrier, reinforcing the multistate reactivity model for these complexes. We discuss how this multistate model can be reconciled with recent computations that revealed that the kinetics of C-H bond oxidation by this set of MnIV-oxo complexes can be well-predicted on the basis of the thermodynamic driving force for these reactions.

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

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.

C-H Bond Activation by Iridium(III) and Iridium(IV) Oxo Complexes

Oxidation of an iridium(III) oxo precursor enabled the structural, spectroscopic, and quantum-chemical characterization of the first well-defined iridium(IV) oxo complex. Side-by-side examination of the proton-coupled electron transfer thermochemistry revealed similar driving forces for the isostructural oxo complexes in two redox states due to compensating contributions from H+ and e- transfer. However, C-H activation of dihydroanthracene revealed significant hydrogen tunneling for the distinctly more basic iridium(III) oxo complex. Our findings complement the growing body of data that relate tunneling to ground state properties as predictors for the selectivity of C-H bond activation.

Robust fragment-based method of calculating hydrogen atom transfer activation barrier in complex molecules

Dynamic processes driven by non-covalent interactions (NCI), such as conformational exchange, molecular binding, and solvation, can strongly influence the rate constants of reactions with low activation barriers, especially at low temperatures. Examples of this may include hydrogen-atom-transfer (HAT) reactions involved in the oxidative stress of an active pharmaceutical ingredient (API). Here, we develop an automated workflow to generate HAT transition-state (TS) geometries for complex and flexible APIs and then systematically evaluate the influences of NCI on the free activation energies, based on the multi-conformational transition-state theory (MC-TST) within the framework of a multi-step reaction path. The two APIs studied: fesoterodine and imipramine, display considerable conformational complexity and have multiple ways of forming hydrogen bonds with the abstracting radical-a hydroxymethyl peroxyl radical. Our results underscore the significance of considering conformational exchange and multiple activation pathways in activation calculations. We also show that structural elements and NCIs outside the reaction site minimally influence TS core geometry and covalent activation barrier, although they more strongly affect reactant binding and consequently the overall activation barrier. We further propose a robust and economical fragment-based method to obtain overall activation barriers, by combining the covalent activation barrier calculated for a small molecular fragment with the binding free energy calculated for the whole molecule.

Harnessing Radicals in Confined Supramolecular Environments Made Possible by MOFs

Researching and utilizing radical intermediates in organic synthetic chemistry have innovated discoveries in methodology and theory. Reactions concerning free radical species opened new pathways beyond the frame of the two-electron mechanism while commonly characterized as rampant processes lacking selectivity. As a result, research in this field has always focused on the controllable generation of radical species and determining factors of selectivity. Metal-organic frameworks (MOFs) have emerged as compelling candidates as catalysts in radical chemistry. From a catalytic point of view, the porous nature of MOFs entails an inner phase for the reaction that could offer possibilities for the regulation of reactivity and selectivity. From a material science perspecti ve, MOFs are organic-inorganic hybrid materials that integrate functional units in organic compounds and complex forms in the tunable long-ranged periodic structure. In this account, we summarized our progress in the application of MOFs in radical chemistry in three parts: (1) The generation of radical species; (2) The weak interactions and site selectivity; (3) Regio- and stereo-selectivity. The unique role of MOFs play in these paradigms is demonstrated in a supramolecular narrative through the analyses of the multi-constituent collaboration within the MOF and the interactions between MOFs and the intermediates during the reactions.

Effect of solvent viscosity on the activation barrier of hydrogen tunneling in the lipoxygenase reaction

Hydrogen tunneling in enzyme reactions has played an important role in linking protein thermal motions to the chemical steps of catalysis. Lipoxygenases (LOXs) have served as model systems for such reactions, showcasing deep hydrogen tunneling mechanisms associated with enzymatic C-H bond cleavage from polyunsaturated fatty acids. Here, we examined the effect of solvent viscosity on the protein thermal motions associated with LOX catalysis using trehalose and glucose as viscogens. Kinetic analysis of the reaction of the paradigm plant orthologue, soybean lipoxygenase (SLO), with linoleic acid revealed no effect on the first-order rate constants, kcat, or activation energy, Ea. Further studies of SLO active site mutants displaying varying Eas, which have been used to probe catalytically relevant motions, likewise provided no evidence for viscogen-dependent motions. Kinetic analyses were extended to a representative fungal LOX from M. oryzae, MoLOX, and a human LOX, 15-LOX-2. While MoLOX behaved similarly to SLO, we show that viscogens inhibit 15-LOX-2 activity. The latter implicates viscogen sensitive, conformational motions in animal LOX reactions. The data provide insight into the role of water hydration layers in facilitating hydrogen (quantum) tunneling in LOX.

Impact of Catalyst Deuteration on the Reactivity of Chiral Phase-Transfer Organocatalysts

Site-specifically deuterated organocatalysts were prepared and found to show improved reactivity over the non-deuterated analogs. Two privileged C2 -symmetric chiral binaphthyl-modified tetraalkylammonium salts were selected for this study. The stability of these phase-transfer catalysts was generally improved by site-specific deuteration, though the degree of improvement was structure dependent. In particular, a large secondary kinetic isotope effect was observed for the tetradeuterated phase-transfer catalyst. The performance of these deuterated catalysts in the asymmetric catalytic alkylation of amino acid derivatives was better than that of non-deuterated analogs at low catalyst loadings. The results suggest that catalyst deuteration is a promising strategy for enhancing the stability and performance of organocatalysts.

Investigating membrane-binding properties of lipoxygenases using surface plasmon resonance

Lipoxygenases (LOXs) catalyze the oxidation of polyunsaturated fatty acids and synthesize oxylipin products that drive important cellular signaling processes in plants and animals. While there has been indirect evidence presented for the interaction of mammalian LOXs with membranes, a quantitative study of the molecular details of LOX-membrane interactions is lacking. Here, we mimicked biological membranes using surface plasmon resonance (SPR) sensor chips derivatized with 2-D planar lipophilic anchors (2D LP) to capture liposomes of varying phospholipid compositions that self-assemble into lipid bilayers on the SPR chip. The sensor chip surfaces were then used to investigate the membrane-binding properties of model LOX enzymes. SPR binding assays displayed reproducible and stable liposome capture to the sensor chip surface that allowed for the detailed characterization of LOX-membrane interactions. Our studies demonstrate a calcium-dependence for the membrane binding activities of coral 8R-LOX and human 15-LOX-2. Furthermore, our data confirm the importance of key membrane insertion loop residues in each of these LOX enzymes for membrane binding activity. Experiments utilizing model plant and human LOXs reveal differences in membrane-binding specificities. Our study establishes and validates a robust SPR-based platform using 2D LP sensor chips that allows for the detailed study of LOX-membrane interactions under different experimental conditions, including altered membrane compositions. Collectively, this investigation improves our overall understanding of LOX-membrane interaction properties, and our SPR-based approach holds potential for future use in the development of LOX-based therapeutics.

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.

13C Electron Nuclear Double Resonance Spectroscopy-Guided Molecular Dynamics Computations Reveal the Structure of the Enzyme-Substrate Complex of an Active, N-Linked Glycosylated Lipoxygenase

Lipoxygenase (LOX) enzymes produce important cell-signaling mediators, yet attempts to capture and characterize LOX-substrate complexes by X-ray co-crystallography are commonly unsuccessful, requiring development of alternative structural methods. We previously reported the structure of the complex of soybean lipoxygenase, SLO, with substrate linoleic acid (LA), as visualized through the integration of 13C/1H electron nuclear double resonance (ENDOR) spectroscopy and molecular dynamics (MD) computations. However, this required substitution of the catalytic mononuclear, nonheme iron by the structurally faithful, yet inactive Mn2+ ion as a spin probe. Unlike canonical Fe-LOXs from plants and animals, LOXs from pathogenic fungi contain active mononuclear Mn2+ metallocenters. Here, we report the ground-state active-site structure of the native, fully glycosylated fungal LOX from rice blast pathogen Magnaporthe oryzae, MoLOX complexed with LA, as obtained through the 13C/1H ENDOR-guided MD approach. The catalytically important distance between the hydrogen donor, carbon-11 (C11), and the acceptor, Mn-bound oxygen, (donor-acceptor distance, DAD) for the MoLOX-LA complex derived in this fashion is 3.4 ± 0.1 Å. The difference of the MoLOX-LA DAD from that of the SLO-LA complex, 3.1 ± 0.1 Å, is functionally important, although is only 0.3 Å, despite the MoLOX complex having a Mn-C11 distance of 5.4 Å and a "carboxylate-out" substrate-binding orientation, whereas the SLO complex has a 4.9 Å Mn-C11 distance and a "carboxylate-in" substrate orientation. The results provide structural insights into reactivity differences across the LOX family, give a foundation for guiding development of MoLOX inhibitors, and highlight the robustness of the ENDOR-guided MD approach to describe LOX-substrate structures.

How proton transfer impacts hachimoji DNA

Hachimoji DNA is a synthetic nucleic acid extension of DNA, formed by an additional four bases, Z, P, S, and B, that can encode information and sustain Darwinian evolution. In this paper, we aim to look into the properties of hachimoji DNA and investigate the probability of proton transfer between the bases, resulting in base mismatch under replication. First, we present a proton transfer mechanism for hachimoji DNA, analogous to the one presented by Löwdin years prior. Then, we use density functional theory to calculate proton transfer rates, tunnelling factors and the kinetic isotope effect in hachimoji DNA. We determined that the reaction barriers are sufficiently low that proton transfer is likely to occur even at biological temperatures. Furthermore, the rates of proton transfer of hachimoji DNA are much faster than in Watson-Crick DNA due to the barrier for Z-P and S-B being 30% lower than in G-C and A-T. Suggesting that proton transfer occurs more frequently in hachimoji DNA than canonical DNA, potentially leading to a higher mutation rate.

Effects of the 10B/11B isotopic substitution on shear relaxation in supercooled B2O3 liquid: A validation of the elastic model of viscous flow

The effects of atomic mass in terms of its zero-point vibrational energy, on molar volume, glass transition temperature T_g, and viscosity are studied in glassy and supercooled B₂O₃ liquids using boron isotope substitutions. The molar volume decreases and T_g and isothermal viscosity increase on the substitution of lighter ¹⁰B isotopes with the heavier ¹¹B isotopes. These effects are argued to be a manifestation of the higher zero-point vibrational energy of the lighter isotope, which along with the anharmonicity of the potential well, results in a longer equilibrium inter-atomic distance and larger mean-square displacement with respect to that for the heavier isotope. The isotope effect on viscosity is increasingly enhanced as the temperature approaches T_g, which is shown to be consistent with the prediction of the elastic models of viscous flow and shear relaxation.

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

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

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.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Metallocavitins as Advanced Enzyme Mimics and Promising Chemical Catalysts

The supramolecular approach is becoming increasingly dominant in biomimetics and chemical catalysis due to the expansion of the enzyme active center idea, which now includes binding cavities (hydrophobic pockets), channels and canals for transporting substrates and products. For a long time, the mimetic strategy was mainly focused on the first coordination sphere of the metal ion. Understanding that a highly organized cavity-like enzymatic pocket plays a key role in the sophisticated functionality of enzymes and that the activity and selectivity of natural metalloenzymes are due to the effects of the second coordination sphere, created by the protein framework, opens up new perspectives in biomimetic chemistry and catalysis. There are two main goals of mimicking enzymatic catalysis: (1) scientific curiosity to gain insight into the mysterious nature of enzymes, and (2) practical tasks of mankind: to learn from nature and adopt from its many years of evolutionary experience. Understanding the chemistry within the enzyme nanocavity (confinement effect) requires the use of relatively simple model systems. The performance of the transition metal catalyst increases due to its retention in molecular nanocontainers (cavitins). Given the greater potential of chemical synthesis, it is hoped that these promising bioinspired catalysts will achieve catalytic efficiency and selectivity comparable to and even superior to the creations of nature. Now it is obvious that the cavity structure of molecular nanocontainers and the real possibility of modifying their cavities provide unlimited possibilities for simulating the active centers of metalloenzymes. This review will focus on how chemical reactivity is controlled in a well-defined cavitin nanospace. The author also intends to discuss advanced metal–cavitin catalysts related to the study of the main stages of artificial photosynthesis, including energy transfer and storage, water oxidation and proton reduction, as well as highlight the current challenges of activating small molecules, such as H2O, CO2, N2, O2, H2, and CH4.

Oxygen versus Sulfur Coordination in Cobalt Superoxo Complexes: Spectroscopic Properties, O2 Binding, and H-Atom Abstraction Reactivity

A five-coordinate, disiloxide-ligated cobalt(II) (S = 3/2) complex (1) was prepared as an oxygen-ligated analogue to the previously reported silanedithiolate-ligated CoII(Me3TACN)(S2SiMe2) (J. Am. Chem. Soc., 2019, 141, 3641-3653). The structural and spectroscopic properties of 1 were analyzed by single-crystal X-ray diffraction, electron paramagnetic resonance (EPR), and NMR spectroscopies. The reactivity of 1 with dioxygen was examined, and it was shown to bind O2 reversibly in a range of solvents at low temperatures. A cobalt(III)-superoxo complex, CoIII(O2·-)(Me3TACN)((OSi2Ph)2O) (2), was generated, and was analyzed by UV-vis, EPR, and resonance Raman spectroscopies. Unlike its sulfur-ligated analogue, complex 2 can thermally release O2 to regenerate 1. Vibrational assignments for selective 18O isotopic labeling of both O2 and disiloxide ligands in 2 are consistent with a 6-coordinate, Co(η1-O2·-)("end-on") complex. Complex 2 reacts with the O-H bond of 4-methoxy-2,2,6,6-tetramethylpiperidin-1-ol (4-MeO-TEMPOH) via H-atom abstraction with a rate of 0.58(2) M-1 s-1 at -105 °C, but it is unable to oxidize phenol substrates. This bracketed reactivity suggests that the O-H bond being formed in the putative CoIII(OOH) product has a relatively weak O-H bond strength (BDFE ∼66-74 kcal mol-1). These thermodynamic and kinetic parameters are similar to those seen for the sulfur-ligated Co(O2)(Me3TACN)(S2SiMe2), indicating that the differences in the electronic structure for O versus S ligation do not have a large impact on H-atom abstraction reactivity.

Dynamical activation of function in metalloenzymes

Formulations of hydrogen tunneling in enzyme-catalysed C-H activation reactions indicate enthalpic barriers to reaction that are independent of chemical steps and dependent on the protein scaffold. A tool to identify catalytically relevant site-specific protein thermal networks has emerged from temperature-dependent hydrogen deuterium exchange (TDHDX). Focusing on mutant enzyme forms with altered activation energies for catalysis, TDHDX provides a comparative analysis of the impact of mutation on Ea for local protein unfolding. Identified thermal networks appear unrelated to protein scaffold conservation and track to the dictates of the catalysed reaction, including sites for metal binding. The positions of thermal networks provide a framework for further understanding of time-dependent, functionally relevant protein motions. Measurement of nanosecond Stokes shifts at the surface of the thermal network in soybean lipoxygenase yields activation energies that are identical to Ea values measured for kcat . This finding identifies a rapid (> nanosecond), long-range and cooperative structural reorganization as the thermal barrier to catalysis. A model for protein dynamics is put forward that integrates broadly distributed protein conformational sampling with protein embedded thermal networks.

Perspective: Path Sampling Methods Applied to Enzymatic Catalysis

This Perspective reviews the use of Transition Path Sampling methods to study enzymatically catalyzed chemical reactions. First applied by our group to an enzymatic reaction over 15 years ago, the method has uncovered basic principles in enzymatic catalysis such as the protein promoting vibration, and it has also helped harmonize such ideas as electrostatic preorganization with dynamic views of enzyme function. It is now being used to help uncover principles of protein design necessary to artificial enzyme creation.

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.

Comparison of the Role of Protein Dynamics in Catalysis by Dihydrofolate Reductase from E. coli and H. sapiens

Dihydrofolate reductase (DHFR) is a well-studied, clinically relevant enzyme known for being highly dynamic over the course of its catalytic cycle. However, the role dynamic motions play in the explicit hydride transfer from the nicotinamide cofactor to the dihydrofolate substrate remains unclear because reaction initiation and direct spectroscopic examination on the appropriate time scale for such femtosecond to picosecond motions is challenging. Here, we employ pre-steady-state kinetics to observe the hydride transfer as directly as possible in two different species of DHFR: Escherichia coli and Homo sapiens. While the hydride transfer has been well-characterized in DHFR from E. coli, improvements in time resolution now allow for sub-millisecond dead times for stopped-flow spectroscopy, which reveals that the maximum rate is indeed faster than previously recorded. The rate in the human enzyme, previously only estimated, is also able to be directly observed using cutting-edge stopped-flow instrumentation. In addition to the pH dependence of the hydride transfer rates for both enzymes, we examine the primary H/D kinetic isotope effect to reveal a temperature dependence in the human enzyme that is absent from the E. coli counterpart. This dependence, which appears above a temperature of 15 °C is a shared feature among other hydride transfer enzymes and is also consistent with computational work suggesting the presence of a fast promoting-vibration that provides donor-acceptor compression on the time scale of catalysis to facilitate the chemistry step.

Coherent state representation of thermal correlation functions with applications to rate theory

A coherent state phase space representation of operators, based on the Husimi distribution, is used to derive an exact expression for the symmetrized version of thermal correlation functions. In addition to the time and temperature independent phase space representation of the two operators whose correlation function is of interest, the integrand includes a non-negative distribution function where only one imaginary time and one real time propagation are needed to compute it. The methodology is exemplified for the flux side correlation function used in rate theory. The coherent state representation necessitates the use of a smeared Gaussian flux operator whose coherent state phase space representation is identical to the classical flux expression. The resulting coherent state expression for the flux side correlation function has a number of advantages as compared to previous formulations. Since only one time propagation is needed, it is much easier to converge it with a semiclassical initial value representation. There is no need for forward–backward approximations, and in principle, the computation may be implemented on the fly. It also provides a route for analytic semiclassical approximations for the thermal rate, as exemplified by a computation of the transmission factor through symmetric and asymmetric Eckart barriers using a thawed Gaussian approximation for both imaginary and real time propagations. As a by-product, this example shows that one may obtain “good” tunneling rates using only above barrier classical trajectories even in the deep tunneling regime.

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.

Condensed-phase isomerization through tunnelling gateways

Quantum mechanical tunnelling describes transmission of matter waves through a barrier with height larger than the energy of the wave¹. Tunnelling becomes important when the de Broglie wavelength of the particle exceeds the barrier thickness; because wavelength increases with decreasing mass, lighter particles tunnel more efficiently than heavier ones. However, there exist examples in condensed-phase chemistry where increasing mass leads to increased tunnelling rates². In contrast to the textbook approach, which considers transitions between continuum states, condensed-phase reactions involve transitions between bound states of reactants and products. Here this conceptual distinction is highlighted by experimental measurements of isotopologue-specific tunnelling rates for CO rotational isomerization at an NaCl surface^3,4, showing nonmonotonic mass dependence. A quantum rate theory of isomerization is developed wherein transitions between sub-barrier reactant and product states occur through interaction with the environment. Tunnelling is fastest for specific pairs of states (gateways), the quantum mechanical details of which lead to enhanced cross-barrier coupling; the energies of these gateways arise nonsystematically, giving an erratic mass dependence. Gateways also accelerate ground-state isomerization, acting as leaky holes through the reaction barrier. This simple model provides a way to account for tunnelling in condensed-phase chemistry, and indicates that heavy-atom tunnelling may be more important than typically assumed.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.

Tunneling in the Hydrogen-Transfer Reaction from a Vitamin E Analog to an Inclusion Complex of 2,2-Diphenyl-1-picrylhydrazyl Radical with β-Cyclodextrin in an Aqueous Buffer Solution at Ambient Temperature

Recently, increasing attention has been paid to quantum mechanical behavior in biology. In this study, we investigated the involvement of quantum mechanical tunneling in the hydrogen-transfer reaction from Trolox, a water-soluble analog of vitamin E (α-tocopherol), to 2,2-diphenyl-1-picrylhydrazyl radical (DPPH•) in a phosphate buffer solution (0.05 M, pH 7.0). DPPH• was used as a reactivity model of reactive oxygen species and solubilized in water using β-cyclodextrin (β-CD). The second-order rate constants, kH and kD, in 0.05 M phosphate buffer solutions prepared with H2O (pH 7.0) and D2O (pD 7.0), respectively, were determined for the reaction between Trolox and DPPH•, using a stopped-flow technique at various temperatures (283-303 K). Large kinetic isotope effects (KIE, kH/kD) were observed for the hydrogen-transfer reaction from Trolox to the β-CD-solubilized DPPH• in the whole temperature range. The isotopic ratio of the Arrhenius prefactor (AH/AD = 0.003), as well as the isotopic difference in the activation energies (19 kJ mol-1), indicated that quantum mechanical tunneling plays a role in the reaction.

Sterically Stabilized End-On Superoxocopper(II) Complexes and Mechanistic Insights into Their Reactivity with O-H, N-H, and C-H Substrates

Instability of end-on superoxocopper(II) complexes, with respect to conversion to peroxo-bridged dicopper(II) complexes, has largely constrained their study to very low temperatures. This limits their kinetic capacity to oxidize substrates. In response, we have developed a series of bulky ligands, Ar₃-TMPA (Ar = tpb, dpb, dtbpb), and used them to support copper(I) complexes that react with O₂ to yield [Cu^II(η¹-O₂^•-)(Ar₃-TMPA)]⁺ species, which are stable against dimerization at all temperatures. Binding of O₂ saturates at subambient temperatures and can be reversed by warming. The onset of oxygenation for the Ar = tpb and dpb systems is observed at 25 °C, and all three [Cu^II(η¹-O₂^•-)(Ar₃-TMPA)]⁺ complexes are stable against self-decay at temperatures of ≤-20 °C. This provides a wide temperature window for study of these complexes, which was exploited by performing extensive reaction kinetics measurements for [Cu^II(η¹-O₂^•-)(tpb₃-TMPA)]⁺ using a broad range of O-H, N-H, and C-H bond substrates. This includes correlation of second order rate constants (k₂) versus oxidation potentials (E_ox) for a range of phenols, construction of Eyring plots, and temperature-dependent kinetic isotope effect (KIE) measurements. The data obtained indicate that reaction with all substrates proceeds via H atom transfer (HAT), reaction with the phenols proceeds with significant charge transfer, and full tunneling of both H and D atoms occurs in the case of 1,2-diphenylhydrazine and 4-methoxy-2,6-di-tert-butylphenol. Oxidation of C-H bonds proved to be kinetically challenging, and whereas [Cu^II(η¹-O₂^•-)(tpb₃-TMPA)]⁺ can oxidize moderately strong O-H and N-H bonds, it is only able to oxidize very weak C-H bonds.

Mimicking Elementary Reactions of Manganese Lipoxygenase Using Mn-hydroxo and Mn-alkylperoxo Complexes

Manganese lipoxygenase (MnLOX) is an enzyme that converts polyunsaturated fatty acids to alkyl hydroperoxides. In proposed mechanisms for this enzyme, the transfer of a hydrogen atom from a substrate C-H bond to an active-site Mn^III-hydroxo center initiates substrate oxidation. In some proposed mechanisms, the active-site Mn^III-hydroxo complex is regenerated by the reaction of a Mn^III-alkylperoxo intermediate with water by a ligand substitution reaction. In a recent study, we described a pair of Mn^III-hydroxo and Mn^III-alkylperoxo complexes supported by the same amide-containing pentadentate ligand (^6Medpaq). In this present work, we describe the reaction of the Mn^III-hydroxo unit in C-H and O-H bond oxidation processes, thus mimicking one of the elementary reactions of the MnLOX enzyme. An analysis of kinetic data shows that the Mn^III-hydroxo complex [Mn^III(OH)(^6Medpaq)]⁺ oxidizes TEMPOH (2,2′-6,6′-tetramethylpiperidine-1-ol) faster than the majority of previously reported Mn^III-hydroxo complexes. Using a combination of cyclic voltammetry and electronic structure computations, we demonstrate that the weak Mn^III-N(pyridine) bonds lead to a higher Mn^III/II reduction potential, increasing the driving force for substrate oxidation reactions and accounting for the faster reaction rate. In addition, we demonstrate that the Mn^III-alkylperoxo complex [Mn^III(OO^tBu)(^6Medpaq)]⁺ reacts with water to obtain the corresponding Mn^III-hydroxo species, thus mimicking the ligand substitution step proposed for MnLOX.

Transition States for Psychrophilic and Mesophilic (R)-3-Hydroxybutyrate Dehydrogenase-Catalyzed Hydride Transfer at Sub-zero Temperatures

(R)-3-Hydroxybutyrate dehydrogenase (HBDH) catalyzes the NADH-dependent reduction of 3-oxocarboxylates to (R)-3-hydroxycarboxylates. The active sites of a pair of cold- and warm-adapted HBDHs are identical except for a single residue, yet kinetics evaluated at -5, 0, and 5 °C show a much higher steady-state rate constant (k_cat) for the cold-adapted than for the warm-adapted HBDH. Intriguingly, single-turnover rate constants (k_STO) are strikingly similar between the two orthologues. Psychrophilic HBDH primary deuterium kinetic isotope effects on k_cat (^Dk_cat) and k_STO (^Dk_STO) decrease at lower temperatures, suggesting more efficient hydride transfer relative to other steps as the temperature decreases. However, mesophilic HBDH ^Dk_cat and ^Dk_STO are generally temperature-independent. The ^Dk_STO data allowed calculation of intrinsic primary deuterium kinetic isotope effects. Intrinsic isotope effects of 4.2 and 3.9 for cold- and warm-adapted HBDH, respectively, at 5 °C, supported by quantum mechanics/molecular mechanics calculations, point to a late transition state for both orthologues. Conversely, intrinsic isotope effects of 5.7 and 3.1 for cold- and warm-adapted HBDH, respectively, at -5 °C indicate the transition state becomes nearly symmetric for the psychrophilic enzyme, but more asymmetric for the mesophilic enzyme. His-to-Asn and Asn-to-His mutations in the psychrophilic and mesophilic HBDH active sites, respectively, swap the single active-site position where these orthologues diverge. At 5 °C, the His-to-Asn mutation in psychrophilic HBDH decreases ^Dk_cat to 3.1, suggesting a decrease in transition-state symmetry, while the His-to-Asn mutation in mesophilic HBDH increases ^Dk_cat to 4.4, indicating an increase in transition-state symmetry. Hence, temperature adaptation and a single divergent active-site residue may influence transition-state geometry in HBDHs.

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.

Small-Molecule Tunnels in Metalloenzymes Viewed as Extensions of the Active Site

Rigorous substrate selectivity is a hallmark of enzyme catalysis. This selectivity is generally ascribed to a thermodynamically favorable process of substrate binding to the enzyme active site based upon complementary physiochemical characteristics, which allows both acquisition and orientation. However, this chemical selectivity is more difficult to rationalize for diminutive molecules that possess too narrow a range of physical characteristics to allow either precise positioning or discrimination between a substrate and an inhibitor. Foremost among these small molecules are dissolved gases such as H₂, N₂, O₂, CO, CO₂, NO, N₂O, NH₃, and CH₄ so often encountered in metalloenzyme catalysis. Nevertheless, metalloenzymes have evolved to metabolize these small-molecule substrates with high selectivity and efficiency.The soluble methane monooxygenase enzyme (sMMO) acts upon two of these small molecules, O₂ and CH₄, to generate methanol as part of the C1 metabolic pathway of methanotrophic organisms. sMMO is capable of oxidizing many alternative hydrocarbon substrates. Remarkably, however, it will preferentially oxidize methane, the substrate with the fewest discriminating physical characteristics and the strongest C-H bond. Early studies led us to broadly attribute this specificity to the formation of a "molecular sieve" in which a methane- and oxygen-sized tunnel provides a size-selective route from bulk solvent to the completely buried sMMO active site. Indeed, recent cryogenic and serial femtosecond ambient temperature crystallographic studies have revealed such a route in sMMO. A detailed study of the sMMO tunnel considered here in the context of small-molecule tunnels identified in other metalloenzymes reveals three discrete characteristics that contribute to substrate selectivity and positioning beyond that which can be provided by the active site itself. Moreover, the dynamic nature of many tunnels allows an exquisite coordination of substrate binding and reaction phases of the catalytic cycle. Here we differentiate between the highly selective molecular tunnel, which allows only the one-dimensional transit of small molecules, and the larger, less-selective channels found in typical enzymes. Methods are described to identify and characterize tunnels as well as to differentiate them from channels. In metalloenzymes which metabolize dissolved gases, we posit that the contribution of tunnels is so great that they should be considered to be extensions of the active site itself. A full understanding of catalysis by these enzymes requires an appreciation of the roles played by tunnels. Such an understanding will also facilitate the use of the enzymes or their synthetic mimics in industrial or pharmaceutical applications.

A large kinetic isotope effect in the reaction of ascorbic acid with 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl 3-oxide (PTIO˙) in aqueous buffer solutions

A large kinetic isotope effect (KIE, kH/kD) of 12.8 was observed for the hydrogen-transfer reaction from ascorbic acid to 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl 3-oxide (PTIO˙) in a phosphate buffer solution (0.05 M, pH/pD 7.0) at 298 K. The isotopic difference in the activation energies (6.8 kJ mol-1) determined from the temperature dependence of the KIE suggests that quantum mechanical tunneling may partly play a role in the reaction, although the isotopic ratio of the Arrhenius prefactor (AH/AD = 0.86) is within the semiclassical limits.