Isotope effects as probes for enzyme catalyzed hydrogen-transfer reactions

Solvent isotope and mutagenesis studies on the proton relay system in yeast alcohol dehydrogenase 1

Alcohol dehydrogenase catalyzes the reversible transfer of a hydride directly from an alcohol to the nicotinamide ring of NAD+ to form an aldehyde and NADH, and the proton from the alcohol probably is transferred through a hydrogen-bonded system to the imidazole of His-48. Studies of the pH dependencies, and solvent and substrate isotope effects on the wild-type and the enzyme with His-48 substituted with Gln-48 were used to demonstrate a role for the proton relay system. The H48Q substitution increases affinities for NAD+ and NADH by ∼2-fold, suggesting that the overall protein structure is maintained. In contrast, catalytic efficiencies (V/Km) on ethanol and acetaldehyde and affinity for 2,2,2-trifluoroethanol are decreased by about 10-fold. The pH dependencies for catalytic efficiencies on ethanol and acetaldehyde (log V/Km versus pH), show pK values of about 7.5 for wild-type enzyme, but ethanol oxidation by H48Q ADH is essentially linear over the pH range from 5.5 to 9.2 with a slope of 0.47. Steady-state kinetics and substrate isotope effects suggest that the kinetic mechanism of H48Q ADH has become partly random for oxidation of ethanol. Both wild-type and H48Q ADHs have pH-independent isotope effects for oxidation (V1/Kb) of 1-butanol/1-butanol-d9 of 4, suggesting that hydride transfer is a major rate-limiting step. The pH dependence for butanol oxidation by wild type ADH shows a wavy profile over the pH range from pH 6 to 10, with a ∼2.3-fold larger V1/Kb in D2O than in H2O, an "inverse" isotope effect. The substrate isotope effect of 4 is not altered by the solvent isotope effect, suggesting concerted proton/hydride transfer. The solvent isotope effect can be explained by a ground state with a water bound to the catalytic zinc in the enzyme-NAD+ complex, and a transition state that resembles a complex with NADH and aldehyde. In contrast, the H48Q enzyme has a diminished inverse solvent isotope effect of ∼1.3 and an essentially linear pH dependence with a slope of log V1/Kb against pH of 0.49 for oxidation of 1-butanol, which together are consistent with a transition state where hydroxide ion directly accepts a proton from the 2'-hydroxyl group of the nicotinamide ribose in the proton relay system in the enzyme-NAD+-alcohol complex. The results support a catalytic role for His-48 in the proton relay system.

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.

SABRE hyperpolarized anticancer agents for use in 1 H MRI

Purpose

Enabling drug tracking (distribution/specific pathways) with magnetic resonance spectroscopy requires manipulation (via hyperpolarization) of spin state populations and targets with sufficiently long magnetic lifetimes to give the largest possible window of observation. Here, we demonstrate how the proton resonances of a group of thienopyridazines (with known anticancer properties), can be amplified using the para‐hydrogen (p‐H₂) based signal amplification by reversible exchange (SABRE) hyperpolarization technique.

Methods

Thienopyridazine isomers, including a ²H version, were synthesized in house. Iridium‐based catalysts dissolved in a methanol‐d₄ solvent facilitated polarization transfer from p‐H₂ gas to the target thienopyridazines. Subsequent SABRE ¹H responses of hyperpolarized thienopyridazines were completed (400 MHz NMR). Pseudo‐singlet state approaches were deployed to extend magnetic state lifetimes. Proof of principle spectral‐spatial images were acquired across a range of field strengths (7T‐9.4T MRI).

Results

¹H‐NMR signal enhancements of −10,130‐fold at 9.4T (~33% polarization) were achieved on thieno[2,3‐d]pyridazine (T[2,3‐d]P), using SABRE under optimal mixing/field transfer conditions. ¹H T₁ lifetimes for the thienopyridazines were ~18‐50 s. Long‐lived state approaches extended the magnetic lifetime of target proton sites in T[2,3‐d]P from an average of 25‐40 seconds. Enhanced in vitro imaging (spatial and chemical shift based) of target T[2,3‐d]P was demonstrated.

Conclusion

Here, we demonstrate the power of SABRE to deliver a fast and cost‐effective route to hyperpolarization of important chemical motifs of anticancer agents. The SABRE approach outlined here lays the foundations for realizing continuous flow, hyperpolarized tracking of drug delivery/pathways.

A guide to the Michaelis-Menten equation: steady state and beyond

The modern definition of enzymology is synonymous with the Michaelis-Menten equation instituted by Leonor Michaelis and Maud Menten. Most textbooks, or chapters within, discussing enzymology start with the derivation of the equation under the assumption of rapid equilibrium (as done by Michaelis-Menten) or steady state (as modified by Briggs and Haldane) conditions to highlight the importance of this equation as the bedrock on which interpretation of enzyme kinetic results is dependent. However, few textbooks or monographs take the effort of placing the equation within its right historical context and discuss the assumptions that have gone into its institution. This guide will dwell on these in substantial detail. Further, this guide will attempt to instil a sense of appreciation for the mathematical curve rectangular hyperbola, its unique attributes and how ubiquitous the curve is in biological systems. To conclude, this guide will discuss the limitations of the equation, and the method it embodies, and trace the journey of how investigators are attempting to move beyond the steady-state approach and the Michaelis-Menten equation into full progress curve, pre-steady state and single-turnover kinetic analysis to obtain greater insights into enzyme kinetics and catalysis.

Synthesis and Biological Evaluation of [18F]FECNT-d4 as a Novel PET Agent for Dopamine Transporter Imaging

PURPOSE

The dopamine transporter (DAT) is a marker of the occurrence and development of Parkinson's disease (PD) and other diseases with nigrostriatal degeneration. 2β-Carbomethoxy-3β-(4-chlorophenyl)-8-(2-[¹⁸F]-fluoroethyl)nortropane ([¹⁸F]FECNT), an ¹⁸F-labelled tropane derivative, was reported to be a useful positron-emitting probe for DAT. However, the rapid formation of brain-penetrating radioactive metabolites is an impediment to the proper quantitation of DAT in PET studies with [¹⁸F]FECNT. Deuterium-substituted analogues have presented better in vivo stability to reduce metabolites. This study aimed to synthesize a deuterium-substituted DAT radiotracer, [¹⁸F]FECNT-d₄, and to make a preliminary investigation of its properties as a DAT tracer in vivo.

PROCEDURES

The ligand [¹⁸F]FECNT-d₄ was obtained by one-step radiolabelling reaction. The lipophilicity was measured by the shake-flask method. Binding properties of [¹⁸F]FECNT-d₄ were estimated by in vitro binding assay, biodistribution, and microPET imaging in rats. In vivo stability of [¹⁸F]FECNT-d₄ was estimated by radio-HPLC.

RESULTS

[¹⁸F]FECNT-d₄ was synthesized at an average activity yield of 46 ± 17 % (n = 15) and the molar activity was 67 ± 12 GBq/μmol. The deuterated tracer showed suitable lipophilicity and the ability to penetrate the blood-brain barrier (brain uptake of 1.72 % ID at 5 min). [¹⁸F]FECNT-d₄ displayed a high binding affinity for DAT comparable to that of [¹⁸F]FECNT in rat striatum homogenates. Biodistribution results in normal rats showed that [¹⁸F]FECNT-d₄ exhibited a higher ratio of the target to non-target (striatum/cerebellum) at 15 min post administration (5.00 ± 0.44 vs 3.84 ± 0.24 for [¹⁸F]FECNT-d₄ vs [¹⁸F]FECNT). MicroPET imaging studies of [¹⁸F]FECNT-d₄ in normal rats showed that the ligand selectively localized to DAT-rich striatal regions and the accumulation could be blocked with DAT inhibitor. Furthermore, in the unilateral PD model rat, a significant reduction of the signal was found in the lesioned side relative to the unlesioned side. Striatal standardized uptake value of [¹⁸F]FECNT-d₄ remained ~4.02 in the striatum between 5 and 20 min, whereas that of [¹⁸F]FECNT fell rapidly from 4.11 to 2.95. Radio-HPLC analysis of the plasma demonstrated better in vivo stability of [¹⁸F]FECNT-d₄ than [¹⁸F]FECNT.

CONCLUSION

The deuterated compound [¹⁸F]FECNT-d₄ may serve as a promising PET imaging agent to assess DAT-related disorders.

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.

The effect of deuteration on the keto-enol equilibrium and photostability of the sunscreen agent avobenzone

The remarkable properties of deuterium have led to many exciting and favourable results in enhancing material properties, for applications in the physical, medical, and biological sciences. Deuterated isotopologues of avobenzone, a sunscreen active ingredient, were synthesised to examine for any changes to the equilibrium between the diketone and enol isomers, as well as their UV photostability and photoprotective properties. Prior to UV irradiation, deuteration of the diketone methylene/enol moiety (i.e. avobenzone-d2) led to an increase in the % diketone compared to non-deuterated, determined by 1H NMR experiments in CDCl3 and C6D12. This can be rationalised from two angles; mechanistically by a deuterium kinetic isotope effect for the CH vs. CD abstraction step during tautomerisation from the diketone to the enol, and a weaker chelating hydrogen bond for the enol when deuterated allowing increased equilibration to the diketone. Avobenzone-d2 was further examined by solid state 13C NMR. The higher % diketone for avobenzone-d2 was postulated to favour increased photodegradation by a non-reversible pathway. This was investigated by UV irradiation of the avobenzone isotopologues in C6D12, both in real time in situ within the NMR by fibre optic cable as well as ex situ using sunlight. An increase in the relative amount of photoproducts for avobenzone-d2 compared to non-deuterated was observed by 1H NMR upon UV irradiation ex situ. Overall, the study demonstrates that deuteration can be applied to alter complex equilibria, and has potential to be manifested as changes to the properties and behaviour of materials.

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

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

Isotopic Labeling of Formate Dehydrogenase Perturbs the Protein Dynamics

Isotope substitution of enzymes has become a means of addressing the participation of protein motions in enzyme-catalyzed reactions. The idea is that only the enzyme mass will be altered and not the electrostatics, so that the protein dynamics are essentially the same but at lower frequencies because of the mass change. In this study, we variably label all carbon atoms in formate dehydrogenase (FDH) with ¹³C, all nitrogen atoms with ¹⁵N, and all nonexchangeable hydrogen atoms with deuterium and investigate the impact that isotopic substitution has on the dynamics at the active site by two-dimensional infrared spectroscopy and compare with the measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs). We show that ¹⁵N labeling of FDH has the largest effect and makes the active site more heterogeneous, whereas the addition of nonexchangeable deuterium appears to have the opposite effect of ¹⁵N on active-site dynamics, resulting in a behavior similar to that of native FDH. Nevertheless, the temperature dependence of the KIEs shows a monotonic trend with protein mass that does not correspond with the changes in dynamics. These results suggest that isotope labeling has more than just a mass effect on enzyme dynamics and may influence electrostatics in ways that complicate the interpretation of the protein isotope effect.

Temperature-Independent Kinetic Isotope Effects as Evidence for a Marcus-like Model of Hydride Tunneling in Phosphite Dehydrogenase

Phosphite dehydrogenase catalyzes the transfer of a hydride from phosphite to NAD⁺, producing phosphate and NADH. We have evaluated the role of hydride tunneling in a thermostable variant of this enzyme (17X-PTDH) by measuring the temperature dependence of the primary ²H kinetic isotope effects (KIEs) between 5 and 45 °C. Pre-steady-state kinetic measurements were used to demonstrate that the hydride transfer is rate-determining across this temperature range and that the observed KIEs are equal to the intrinsic isotope effect on the chemical step. The KIEs on the pre-exponential factor (A_H/A_D) and the activation energy (ΔE_a) were 1.6 ± 0.1 and 0.21 ± 0.05 kcal/mol, respectively, suggesting that 17X-PTDH facilitates extensive tunneling of both isotopes via a Marcus-like model. Site-directed mutagenesis was used to evaluate the role of an active site threonine (Thr104) found on the back face of the nicotinamide in promoting the close packing of the substrates. In mutants with reduced steric bulk at this position, values of A_H/A_D and ΔE_a fall within the range describing semiclassical "over the barrier" reactivity, suggesting that Thr104 acts as a steric backstop to promote tunneling in 17X-PTDH. Whereas hydrogen tunneling is now a widely appreciated feature of C-H activating enzymes, these observations with a P-H activating system are consistent with the proposal that tunneling is likely to be a common feature on all enzymes that catalyze hydrogen transfers.

Quantum Dynamics and Kinetics of the F + H2 and F + D2 Reactions at Low and Ultra-Low Temperatures

Integral cross sections and rate constants for the prototypical chemical reactions of the fluorine atom with molecular hydrogen and deuterium have been calculated over a wide interval of collision energy and temperature ranging from the sub-thermal (50 K) down to the ultra-cold regimes (0.5 mK). Rigorous close coupling time-independent quantum reactive scattering calculations have been carried out on two potential energy surfaces, differing only at long-range in the reactants' channel. The results show that tunnel, resonance and virtual state effects enhance under-barrier reactivity giving rise to pronounced deviations from the Arrhenius law as temperature is lowered. Within the ultra-cold domain (below 1 mK), the reactivity is governed by virtual state effects and by tunneling through the reaction barrier; in the cold regime (1 mK–1 K), the shape resonances in the entrance channel of the potential energy surface make the quantum tunneling contribution larger so enhancing cross sections and rate constants by about one order of magnitude; at higher temperatures (above 10 K), the tunneling pathway enhanced by the constructive interference between two Feshbach resonances trapped in the reaction exit channel competes with the thermally activated mechanism, as the energy gets closer to the reaction barrier height. The results show that at low temperatures cross sections and rate constants are extremely sensitive to small changes in the long-range intermolecular interaction in the entrance channel of the potential energy surface, as well as to isotopic substitution.

Theoretical Studies on Mechanism of Inactivation of Kanamycin A by 4'-O-Nucleotidyltransferase

This work is focused on mechanistic studies of the transfer of an adenylyl group (Adenoside-5'-monophosfate) from adenosine 5'-triphosphate (ATP) to a OH-4' hydroxyl group of an antibiotic. Using hybrid quantum mechanics/molecular mechanics (QM/MM) techniques, we study the substrate and base-assisted mechanisms of the inactivation process of kanamycin A (KAN) catalyzed by 4'-O-Nucleotidyltransferase [ANT(4')], an active enzyme against almost all aminoglycoside antibiotics. Free energy surfaces, obtained with Free Energy Perturbation methods at the M06-2X/MM level of theory, show that the most favorable reaction path presents a barrier of 12.2 kcal·mol-1 that corresponds to the concerted activation of O4' from KAN by Glu145. In addition, the primary and secondary 18O kinetic isotope effects (KIEs) have been computed for bridge O3α, and non-bridge O1α, O2α, and O5' atoms of ATP. The observed normal 1°-KIE of 1.2% and 2°-KIE of 0.07% for the Glu145-assisted mechanism are in very good agreement with experimentally measured data. Additionally, based on the obtained results, the role of electrostatic and compression effects in enzymatic catalysis is discussed.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Bacterial versus human thymidylate synthase: Kinetics and functionality

Thymidylate Synthase (TSase) is a highly conserved enzyme that catalyzes the production of the DNA building block thymidylate. Structurally, functionally and mechanistically, bacterial and mammalian TSases share remarkable similarities. Because of this closeness, bacterial enzymes have long been used as model systems for human TSase. Furthermore, while TSase inhibitors have long served as chemotherapeutic drugs, no TSase inhibitor serves as an antibiotic. Despite their high resemblance, the mammalian TSases are distinct in a few known aspects, such as having a N-terminal tail and two insertions in the primary sequence and active/inactive conformations. Here, we aim to comprehensively characterize human (hs) TSase and delineate its contrasts and the similarities to the well-studied Escherichia coli (ec) TSase. We found that, in contrast to ecTSase, Mg²⁺ does not enhance reaction rates for hsTSase. The temperature dependence of intrinsic kinetic isotope effects (KIEs), on the other hand, suggests that Mg²⁺ has little or no impact on the transition state of hydride transfer in either enzyme, and that the transition state for the hydride transfer in hsTSase is looser than in ecTSase. Additionally, the substrates’ binding order is strictly ordered for ecTSase but slightly less ordered for hsTSase. The observed kinetic and functional differences between bacterial and human enzymes may aid in the development of antibiotic drugs with reduced toxicity.

Deuterium- and Tritium-Labelled Compounds: Applications in the Life Sciences

Hydrogen isotopes are unique tools for identifying and understanding biological and chemical processes. Hydrogen isotope labelling allows for the traceless and direct incorporation of an additional mass or radioactive tag into an organic molecule with almost no changes in its chemical structure, physical properties, or biological activity. Using deuterium-labelled isotopologues to study the unique mass-spectrometric patterns generated from mixtures of biologically relevant molecules drastically simplifies analysis. Such methods are now providing unprecedented levels of insight in a wide and continuously growing range of applications in the life sciences and beyond. Tritium (³ H), in particular, has seen an increase in utilization, especially in pharmaceutical drug discovery. The efforts and costs associated with the synthesis of labelled compounds are more than compensated for by the enhanced molecular sensitivity during analysis and the high reliability of the data obtained. In this Review, advances in the application of hydrogen isotopes in the life sciences are described.

Protein Mass Effects on Formate Dehydrogenase

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Mechanism-Informed Refinement Reveals Altered Substrate-Binding Mode for Catalytically Competent Nitroreductase

Nitroreductase (NR) from Enterobacter cloacae reduces diverse nitroaromatics including herbicides, explosives, and prodrugs, and holds promise for bioremediation, prodrug activation, and enzyme-assisted synthesis. We solved crystal structures of NR complexes with bound substrate or analog for each of its two half-reactions. We complemented these with kinetic isotope effect (KIE) measurements elucidating H-transfer steps essential to each half-reaction. KIEs indicate hydride transfer from NADH to the flavin consistent with our structure of NR with the NADH analog nicotinic acid adenine dinucleotide (NAAD). The KIE on reduction of p-nitrobenzoic acid (p-NBA) also indicates hydride transfer, and requires revision of prior computational mechanisms. Our mechanistic information provided a structural restraint for the orientation of bound substrate, placing the nitro group closer to the flavin N5 in the pocket that binds the amide of NADH. KIEs show that solvent provides a proton, enabling accommodation of different nitro group placements, consistent with the broad repertoire of NR.

Hydrogen Tunneling above Room Temperature Evidenced by Infrared Ion Spectroscopy

While hydrogen tunneling at elevated temperatures has, for instance, often been postulated in biochemical processes, spectroscopic proof is thus far limited to cryogenic conditions, under which thermal reactivity is negligible. We report spectroscopic evidence for H-tunneling in the gas phase at temperatures around 320-350 K observed in the isomerization reaction of a hydroxycarbene into an aldehyde. The charge-tagged carbene was generated in situ in a tandem mass spectrometer by decarboxylation of oxo[4-(trimethylammonio)phenyl]acetic acid upon collision induced dissociation. All ion structures involved are characterized by infrared ion spectroscopy and quantum chemical calculations. The charge-tagged phenylhydroxycarbene undergoes a 1,2-H-shift to the corresponding aldehyde with an half-life of about 10 s, evidenced by isomer-selective two-color (IR-IR) spectroscopy. In contrast, the deuterated (OD) carbene analogue showed much reduced 1,2-D-shift reactivity with an estimated half-life of at least 200 s under the experimental conditions, and provides clear evidence for hydrogen atom tunneling in the H-isotopologue. This is the first spectroscopic confirmation of hydrogen atom tunneling governing 1,2-H-shift reactions at noncryogenic temperatures, which is of broad significance for a range of (bio)chemical processes, including enzymatic transformations and organocatalysis.

Simultaneous Measurement of Glucose-6-phosphate 3-Dehydrogenase (NtdC) Catalysis and the Nonenzymatic Reaction of Its Product: Kinetics and Isotope Effects on the First Step in Kanosamine Biosynthesis

Glucose-6-phosphate 3-dehydrogenase (NtdC) is an NAD-dependent oxidoreductase encoded in the NTD operon of Bacillus subtilis. The oxidation of glucose 6-phosphate by NtdC is the first step in kanosamine biosynthesis. The product, 3-oxo-d-glucose 6-phosphate (3oG6P), has never been synthesized or isolated. The NtdC-catalyzed reaction is very slow at low and neutral pH, and its rate increases to a maximum near pH 9.5. However, under alkaline conditions, the product is not stable because of ring opening followed by deprotonation of the 1,3-dicarbonyl compound. The absorbance band due to this enolate at 310 nm overlaps with that of the other enzymatic product, NADH, complicating kinetic measurements. We report the deconvolution of the resulting spectra of the reaction to determine the rate constants and likely kinetic mechanism. In doing so, we were able to determine the extinction coefficient of the enolate of 3oG6P (23000 M^-1 cm^-1), which allowed the measurement of the first-order rate constant (5.51 × 10^-3 s^-1) and activation energy (93 kJ mol^-1) of nonenzymatic enolate formation. Using deuterium-labeled substrates, we show that hydride transfer from carbon 3 is partially rate-limiting in the enzymatic reaction, and deuterium substitution on carbon 2 has no significant effect on the enzymatic reaction but lowers the rate of deprotonation of 3oG6P 4-fold. These experiments clearly establish the regiochemistry of the reactions. Coupling of the NtdC reaction with the subsequent step in the pathway, NtdA-catalyzed glutamate-dependent amino transfer, has a small but significant effect on the rate of NAD reduction, consistent with these enzymes working together to process the unstable metabolite.

Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects

Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Examinations of the Chemical Step in Enzyme Catalysis

Advances in computational and experimental methods in enzymology have aided comprehension of enzyme-catalyzed chemical reactions. The main difficulty in comparing computational findings to rate measurements is that the first examines a single energy barrier, while the second frequently reflects a combination of many microscopic barriers. We present here intrinsic kinetic isotope effects and their temperature dependence as a useful experimental probe of a single chemical step in a complex kinetic cascade. Computational predictions are tested by this method for two model enzymes: dihydrofolate reductase and thymidylate synthase. The description highlights the significance of collaboration between experimentalists and theoreticians to develop a better understanding of enzyme-catalyzed chemical conversions.

Wide-dynamic-range kinetic investigations of deep proton tunnelling in proteins

Directional proton transport along 'wires' that feed biochemical reactions in proteins is poorly understood. Amino-acid residues with high pKa are seldom considered as active transport elements in such wires because of their large classical barrier for proton dissociation. Here, we use the light-triggered proton wire of the green fluorescent protein to study its ground-electronic-state proton-transport kinetics, revealing a large temperature-dependent kinetic isotope effect. We show that 'deep' proton tunnelling between hydrogen-bonded oxygen atoms with a typical donor-acceptor distance of 2.7-2.8 Å fully accounts for the rates at all temperatures, including the unexpectedly large value (2.5 × 10(9) s(-1)) found at room temperature. The rate-limiting step in green fluorescent protein is assigned to tunnelling of the ionization-resistant serine hydroxyl proton. This suggests how high-pKa residues within a proton wire can act as a 'tunnel diode' to kinetically trap protons and control the direction of proton flow.

Possible interstellar formation of glycine through a concerted mechanism: a computational study on the reaction of CH2NH, CO2 and H2

Computational studies on the reaction of CH₂NH, CO₂ and H₂ show the possible interstellar formation of glycine in both hot-cores and cold interstellar clouds.

Preserved hydride transfer mechanism in evolutionarily divergent thymidylate synthases

Thymidylate synthase (TSase) catalyzes a hydride transfer in the last step of the de novo biosynthesis of the DNA nucleotide thymine. We compared two isozymes, namely, TSase from Escherichia coli (ecTSase) and TSase from Bacillus subtilis (bsTSase) that represent a case of divergent evolution. Interestingly, a highly conserved histidine (H147 of ecTSase) was proposed to serve a critical role in catalysis, but in bsTSase it is naturally substituted by valine (Val). Yet, bsTSase is more active than ecTSase, and the intrinsic kinetic isotope effects (KIEs) of both are temperature-independent, suggesting a similarly well-organized transition state (TS) for the catalyzed hydride transfer. To examine the role of that histidine (His) in TSase catalysis, we examined the kinetics of H147V ecTSase, which "bridges" between these two TSases. In contrast to both wild-type TSases, the single mutation results in deficient catalysis. The mutation leads to intrinsic KIEs that are temperature-dependent, indicating a substantial imperfection in its TS. The findings reveal two important features: a direct role of H147 in the hydride transfer step catalyzed by the ecTSase and the evolutionary compensation for its deficiency in bsTSase via extensive polymorphism across the protein. Very different active site residues are observed for these evolutionarily divergent isozymes, which result in a well-organized TS for both. It is suggested that evolutionary pressure compensated for the H to V substitution at the active site of bsTSase by polymorphism leading to a well-organized TS in both enzymes.

Catalytic Mechanism of the Maltose Transporter Hydrolyzing ATP

We use quantum mechanical and molecular mechanical (QM/MM) simulations to study ATP hydrolysis catalyzed by the maltose transporter. This protein is a prototypical member of a large family that consists of ATP-binding cassette (ABC) transporters. The ABC proteins catalyze ATP hydrolysis to perform a variety of biological functions. Despite extensive research efforts, the precise molecular mechanism of ATP hydrolysis catalyzed by the ABC enzymes remains elusive. In this work, the reaction pathway for ATP hydrolysis in the maltose transporter is evaluated using a QM/MM implementation of the nudged elastic band method without presuming reaction coordinates. The potential of mean force along the reaction pathway is obtained with an activation free energy of 19.2 kcal/mol in agreement with experiments. The results demonstrate that the reaction proceeds via a dissociative-like pathway with a trigonal bipyramidal transition state in which the cleavage of the γ-phosphate P-O bond occurs and the O-H bond of the lytic water molecule is not yet broken. Our calculations clearly show that the Walker B glutamate as well as the switch histidine stabilizes the transition state via electrostatic interactions rather than serving as a catalytic base. The results are consistent with biochemical and structural experiments, providing novel insight into the molecular mechanism of ATP hydrolysis in the ABC proteins.

The general base in the thymidylate synthase catalyzed proton abstraction

The enzyme thymidylate synthase (TSase), an important chemotherapeutic drug target, catalyzes the formation of 2'-deoxythymidine-5'-monophosphate (dTMP), a precursor of one of the DNA building blocks. TSase catalyzes a multi-step mechanism that includes the abstraction of a proton from the C5 of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). Previous studies on ecTSase proposed that an active-site residue, Y94 serves the role of the general base abstracting this proton. However, since Y94 is neither very basic, nor connected to basic residues, nor located close enough to the pyrimidine proton to be abstracted, the actual identity of this base remains enigmatic. Based on crystal structures, an alternative hypothesis is that the nearest potential proton-acceptor of C5 of dUMP is a water molecule that is part of a hydrogen bond (H-bond) network comprised of several water molecules and several protein residues including H147, E58, N177, and Y94. Here, we examine the role of the residue Y94 in the proton abstraction step by removing its hydroxyl group (Y94F mutant). We investigated the effect of the mutation on the temperature dependence of intrinsic kinetic isotope effects (KIEs) and found that these KIEs are more temperature dependent than those of the wild-type enzyme (WT). These results suggest that the phenolic -OH of Y94 is a component of the transition state for the proton abstraction step. The findings further support the hypothesis that no single functional group is the general base, but a network of bases and hydroxyls (from water molecules and tyrosine) sharing H-bonds across the active site can serve the role of the general base to remove the pyrimidine proton.

Activation of Two Sequential H-transfers in the Thymidylate Synthase Catalyzed Reaction

Thymidylate synthase (TSase) catalyzes the de novo biosynthesis of thymidylate, a precursor for DNA, and is thus an important target for chemotherapeutics and antibiotics. Two sequential C-H bond cleavages catalyzed by TSase are of particular interest: a reversible proton abstraction from the 2'-deoxy-uridylate substrate, followed by an irreversible hydride transfer forming the thymidylate product. QM/MM calculations of the former predicted a mechanism where the abstraction of the proton leads to formation of a novel nucleotide-folate intermediate that is not covalently bound to the enzyme (Wang, Z.; Ferrer, S.; Moliner, V.; Kohen, A. Biochemistry2013, 52, 2348-2358). Existence of such intermediate would hold promise as a target for a new class of drugs. Calculations of the subsequent hydride transfer predicted a concerted H-transfer and elimination of the enzymatic cysteine (Kanaan, N.; Ferrer, S.; Marti, S.; Garcia-Viloca, M.; Kohen, A.; Moliner, V. J. Am. Chem. Soc.2011, 133, 6692-6702). A key to both C-H activations is a highly conserved arginine (R166) that stabilizes the transition state of both H-transfers. Here we test these predictions by studying the R166 to lysine mutant of E. coli TSase (R166K) using intrinsic kinetic isotope effects (KIEs) and their temperature dependence to assess effects of the mutation on both chemical steps. The findings confirmed the predictions made by the QM/MM calculations, implicate R166 as an integral component of both reaction coordinates, and thus provide critical support to the nucleotide-folate intermediate as a new target for rational drug design.

Methods to Increase the Metabolic Stability of (18)F-Radiotracers

The majority of pharmaceuticals and other organic compounds incorporating radiotracers that are considered foreign to the body undergo metabolic changes in vivo. Metabolic degradation of these drugs is commonly caused by a system of enzymes of low substrate specificity requirement, which is present mainly in the liver, but drug metabolism may also take place in the kidneys or other organs. Thus, radiotracers and all other pharmaceuticals are faced with enormous challenges to maintain their stability in vivo highlighting the importance of their structure. Often in practice, such biologically active molecules exhibit these properties in vitro, but fail during in vivo studies due to obtaining an increased metabolism within minutes. Many pharmacologically and biologically interesting compounds never see application due to their lack of stability. One of the most important issues of radiotracers development based on fluorine-18 is the stability in vitro and in vivo. Sometimes, the metabolism of ¹⁸F-radiotracers goes along with the cleavage of the C-F bond and with the rejection of [¹⁸F]fluoride mostly combined with high background and accumulation in the skeleton. This review deals with the impact of radiodefluorination and with approaches to stabilize the C-F bond to avoid the cleavage between fluorine and carbon.

Network of remote and local protein dynamics in dihydrofolate reductase catalysis

Molecular dynamics calculations and bionformatic studies of dihydrofolate reductase (DHFR) have suggested a network of coupled motions across the whole protein that is correlated to the reaction coordinate. Experimental studies demonstrated that distal residues G121, M42 and F125 in E. coli DHFR participate in that network. The missing link in our understanding of DHFR catalysis is the lack of a mechanism by which such remote residues can affect the catalyzed chemistry at the active site. Here, we present a study of the temperature dependence of intrinsic kinetic isotope effects (KIEs) that indicates synergism between a remote residue in that dynamic network, G121, and the active site's residue I14. The intrinsic KIEs for the I14A-G121V double mutant showed steeper temperature dependence (ΔE_a(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔE_a(T-H)G121V +ΔE_{a(T-H) I14A} < ΔE_{a(T-H) double mutant}, which indicates a synergism between the two residues. This finding links the remote residues in the network under investigation to the enzyme's active site, providing a mechanism by which these residues can be coupled to the catalyzed chemistry. This experimental evidence validates calculations proposing that both remote and active site residues constitute a network of coupled promoting motions correlated to the bond activation step (C-H→C hydride transfer in this case). Additionally, the effect of I14A and G121V mutations on single turnover rates was additive rather than synergistic. Although single turnover rate measurements are more readily available and thus more popular than assessing intrinsic kinetic isotope effects, the current finding demonstrates that for these rates, which in DHFR reflect several microscopic rate constants, can fall short of revealing the nature of the C-H bond activation per se.

Role of long-range protein dynamics in different thymidylate synthase catalyzed reactions

Recent studies of Escherichia coli thymidylate synthase (ecTSase) showed that a highly conserved residue, Y209, that is located 8 Å away from the reaction site, plays a key role in the protein’s dynamics. Those crystallographic studies indicated that Y209W mutant is a structurally identical but dynamically altered relative to the wild type (WT) enzyme, and that its turnover catalytic rate governed by a slow hydride-transfer has been affected. The most challenging test of an examination of a fast chemical conversion that precedes the rate-limiting step has been achieved here. The physical nature of both fast and slow C-H bond activations have been compared between the WT and mutant by means of observed and intrinsic kinetic isotope effects (KIEs) and their temperature dependence. The findings indicate that the proton abstraction step has not been altered as much as the hydride transfer step. Additionally, the comparison indicated that other kinetic steps in the TSase catalyzed reaction were substantially affected, including the order of the substrate binding. Enigmatically, although Y209 is H-bonded to 3'-OH of 2'-deoxyuridine-5'-mono­phosphate (dUMP), its altered dynamics is more pronounced on the binding of the remote cofactor, (6R)-N⁵,N¹⁰-methylene-5,6,7,8-tetrahydrofolate (CH₂H₄folate), revealing the importance of long-range dynamics of the enzymatic complex and its catalytic function.

Compound-specific isotope analysis as a tool to characterize biodegradation of ethylbenzene

This study applied one- and two-dimensional compound-specific isotope analysis (CSIA) for the elements carbon and hydrogen to assess different means of microbial ethylbenzene activation. Cultures incubated under nitrate-reducing conditions showed significant carbon and highly pronounced hydrogen isotope fractionation of comparable magnitudes, leading to nearly identical slopes in dual-isotope plots. The results imply that Georgfuchsia toluolica G5G6 and an enrichment culture dominated by an Azoarcus species activate ethylbenzene by anaerobic hydroxylation catalyzed by ethylbenzene dehydrogenase, similar to Aromatoleum aromaticum EbN1. The isotope enrichment pattern in dual plots from two strictly anaerobic enrichment cultures differed considerably from those for benzylic hydroxylation, indicating an alternative anaerobic activation step, most likely fumarate addition. Large hydrogen fractionation was quantified using a recently developed Rayleigh-based approach considering hydrogen atoms at reactive sites. Data from nine investigated microbial cultures clearly suggest that two-dimensional CSIA in combination with the magnitude of hydrogen isotope fractionation is a valuable tool to distinguish ethylbenzene degradation and may be of practical use for monitoring natural or technological remediation processes at field sites.

Concerted versus stepwise mechanism in thymidylate synthase

Thymidylate synthase (TSase) catalyzes the intracellular de novo formation of thymidylate (a DNA building block) in most living organisms, making it a common target for chemotherapeutic and antibiotic drugs. Two mechanisms have been proposed for the rate-limiting hydride transfer step in TSase catalysis: a stepwise mechanism in which the hydride transfer precedes the cleavage of the covalent bond between the enzymatic cysteine and the product and a mechanism where both happen concertedly. Striking similarities between the enzyme-bound enolate intermediates formed in the initial and final step of the reaction supported the first mechanism, while QM/MM calculations favored the concerted mechanism. Here, we experimentally test these two possibilities using secondary kinetic isotope effect (KIE), mutagenesis study, and primary KIEs. The findings support the concerted mechanism and demonstrate the critical role of an active site arginine in substrate binding, activation of enzymatic nucleophile, and the hydride transfer studied here. The elucidation of this reduction/substitution sheds light on the critical catalytic step in TSase and may aid future drug or biomimetic catalyst design.

Kinetic isotope effects as a probe of hydrogen transfers to and from common enzymatic cofactors

Enzymes use a number of common cofactors as sources of hydrogen to drive biological processes, but the physics of the hydrogen transfers to and from these cofactors is not fully understood. Researchers study the mechanistically important contributions from quantum tunneling and enzyme dynamics and connect those processes to the catalytic power of enzymes that use these cofactors. Here we describe some progress that has been made in studying these reactions, particularly through the use of kinetic isotope effects (KIEs). We first discuss the general theoretical framework necessary to interpret experimental KIEs, and then describe practical uses for KIEs in the context of two case studies. The first example is alcohol dehydrogenase, which uses a nicotinamide cofactor to catalyze a hydride transfer, and the second example is thymidylate synthase, which uses a folate cofactor to catalyze both a hydride and a proton transfer.

Extension and limits of the network of coupled motions correlated to hydride transfer in dihydrofolate reductase

Enzyme catalysis has been studied extensively, but the role of enzyme dynamics in the catalyzed chemical conversion is still an enigma. The enzyme dihydrofolate reductase (DHFR) is often used as a model system to assess a network of coupled motions across the protein that may affect the catalyzed chemical transformation. Molecular dynamics simulations, quantum mechanical/molecular mechanical studies, and bioinformatics studies have suggested the presence of a “global dynamic network” of residues in DHFR. Earlier studies of two DHFR distal mutants, G121V and M42W, indicated that these residues affect the chemical step synergistically. While this finding was in accordance with the concept of a network of functional motions across the protein, two residues do not constitute a network. To better define the extent and limits of the proposed network, the current work studied two remote residues predicted to be part of the same network: W133 and F125. The effect of mutations in these residues on the nature of the chemical step was examined via measurements of the temperature-dependence of the intrinsic kinetic isotope effects (KIEs) and other kinetic parameters, and double mutants were used to tie the findings to G121 and M42. The findings indicate that residue F125, which was implicated by both calculations and bioinformatic methods, is a part of the same global dynamic network as G121 and M42, while W133, implicated only by bioinformatics, is not. These findings extend our understanding of the proposed network and the relations between functional and genomic couplings. Delineating that network illuminates the need to consider remote residues and protein structural dynamics in the rational design of drugs and of biomimetic catalysts.

Effects of cavities at the nicotinamide binding site of liver alcohol dehydrogenase on structure, dynamics and catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.

Stereospecific multiple isotopic labeling of benzyl alcohol

Isotopically labeled enzymatic substrates and biological metabolites are useful for many mechanistic analyses, particularly the study of kinetic and equilibrium isotope effects, determining the stereospecificity of enzymes, and resolving metabolic pathways. Here, we present the one-pot synthesis, purification, and kinetic analysis of 7R-[(2) H]-phenyl-[(14) C]-benzyl alcohol. The procedure involves a chemoenzymatic synthesis that couples formate dehydrogenase to alcohol dehydrogenase with a catalytic amount of nicotinamide cofactor. The reaction goes to completion overnight, and the measurement of a competitive kinetic isotope effect on the enzymatic oxidation of the purified product identified no (1) H contamination. This measurement is very sensitive to such isotopic contamination and verified the high level of isotopic and enantiomeric purity yielded by the new synthetic procedure.

A critical test of the "tunneling and coupled motion" concept in enzymatic alcohol oxidation

The physical mechanism of C-H bond activation by enzymes is the subject of intense study, and we have tested the predictions of two competing models for C-H activation in the context of alcohol dehydrogenase. The kinetic isotope effects (KIEs) in this enzyme have previously suggested a model of quantum mechanical tunneling and coupled motion of primary (1°) and secondary (2°) hydrogens. Here we measure the 2° H/T KIEs with both H and D at the 1° position and find that the 2° KIE is significantly deflated with D-transfer, consistent with the predictions of recent Marcus-like models of H-transfer. The results suggest that the fast dynamics of H-tunneling result in a 1° isotope effect on the structure of the tunneling ready state: the trajectory of D-transfer goes through a shorter donor-acceptor distance than that of H-transfer.

Quantum Tunnelling to the Origin and Evolution of Life

Quantum tunnelling is a phenomenon which becomes relevant at the nanoscale and below. It is a paradox from the classical point of view as it enables elementary particles and atoms to permeate an energetic barrier without the need for sufficient energy to overcome it. Tunnelling might seem to be an exotic process only important for special physical effects and applications such as the Tunnel Diode, Scanning Tunnelling Microscopy (electron tunnelling) or Near-field Optical Microscopy operating in photon tunnelling mode. However, this review demonstrates that tunnelling can do far more, being of vital importance for life: physical and chemical processes which are crucial in theories about the origin and evolution of life can be traced directly back to the effects of quantum tunnelling. These processes include the chemical evolution in stellar interiors and within the cold interstellar medium, prebiotic chemistry in the atmosphere and subsurface of planetary bodies, planetary habitability via insolation and geothermal heat as well as the function of biomolecular nanomachines. This review shows that quantum tunnelling has many highly important implications to the field of molecular and biological evolution, prebiotic chemistry and astrobiology.