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

Structures of lactaldehyde reductase, FucO, link enzyme activity to hydrogen bond networks and conformational dynamics

A group-III iron containing 1,2-propanediol oxidoreductase, FucO, (also known as lactaldehyde reductase) from Escherichia coli was examined regarding its structure-dynamics-function relationships in the catalysis of the NADH-dependent reduction of (2S)-lactaldehyde. Crystal structures of FucO variants in the presence or absence of cofactors have been determined, illustrating large domain movements between the apo and holo enzyme structures. Different structures of FucO variants co-crystallized with NAD⁺ or NADH together with substrate further suggest dynamic properties of the nicotinamide moiety of the coenzyme that are important for the reaction mechanism. Modelling of the native substrate (2S)-lactaldehyde into the active site can explain the stereoselectivity exhibited by the enzyme, with a critical hydrogen bond interaction between the (2S)-hydroxyl and the side-chain of N151, as well as the previously experimentally demonstrated pro-(R) selectivity in hydride transfer from NADH to the aldehydic carbon. Furthermore, the deuterium kinetic isotope effect of hydride transfer suggests that reduction chemistry is the main rate-limiting step for turnover which is not the case in FucO catalysed alcohol oxidation. We further propose that a water molecule in the active site - hydrogen bonded to a conserved histidine (H267) and the 2'-hydroxyl of the coenzyme ribose - functions as a catalytic proton donor in the protonation of the product alcohol. A hydrogen bond network of water molecules and the side-chains of amino acid residues D360 and H267 links bulk solvent to this proposed catalytic water molecule.

Raman and Quantum Yield Studies of Trp48-d5 in Azurin: Closed-Shell and Neutral Radical Species

Isotopologues are valuable vibrational probes that shift features in a vibrational spectrum while preserving the electronic structure of the molecule. We report the vibrational and electronic spectra of perdeuterated tryptophan in solution (l-Trp-d₅), as Trp48-d₅ in azurin, and as the photogenerated neutral tryptophan radical, Trp48-d₅^•, in azurin. The UV resonance Raman bands of the perdeuterated closed-shell tryptophan in solution and in azurin are lower in frequency relative to the protiated counterpart. The observed decrease in frequencies of l-Trp-d₅ bands relative to l-Trp-h₅ enables the analysis of vibrational markers of other amino acids, e.g., phenylalanine, that overlap with some modes of l-Trp-h₅. The Raman intensities vary between l-Trp-d₅ and l-Trp-h₅; these differences likely reflect modifications in normal mode composition upon perdeuteration. Analysis of the W3, W6, and W17 modes suggests that the W3 mode retains its utility as a conformational marker; however, the H-bond markers W6 and W17 appear to be less sensitive upon perdeuteration. The neutral tryptophan radical, Trp48-d₅^•, was generated in azurin with a slightly lower radical quantum yield than for Trp48-h₅^•. The visible resonance Raman spectrum of Trp48-d₅^• is different from that of Trp48-h₅^•, especially in terms of relative intensities, and all assignable peaks decreased in frequency upon perdeuteration. The absorption and emission spectra of the perdeuterated closed-shell and radical species exhibited hypsochromic shifts of less than 1 nm relative to the protiated species. The data presented here indicate that l-Trp-d₅ is a valuable probe of vibrational structure, with minimal modification of photoreactivity and photophysics compared to l-Trp-h₅.

C-H Bond Cleavage Is Rate-Limiting for Oxidative C-P Bond Cleavage by the Mixed Valence Diiron-Dependent Oxygenase PhnZ

PhnZ utilizes a mixed valence diiron(II/III) cofactor and O₂ to oxidatively cleave the carbon-phosphorus bond of (R)-2-amino-1-hydroxyethylphosphonic acid to form glycine and orthophosphate. The active site residues Y24 and E27 are proposed to mediate induced-fit recognition of the substrate and access of O₂ to one of the active site Fe ions. H62 is proposed to deprotonate the C1-hydroxyl of the substrate during catalysis. Kinetic isotope effects (KIEs), pH-rate dependence, and site-directed mutagenesis were used to probe the rate-determining transition state and the roles of these three active site residues. Primary deuterium KIE values of 5.5 ± 0.3 for ^D(V) and 2.2 ± 0.4 for ^D(V/K) were measured with (R)-2-amino[1-²H₁]-1-hydroxyethylphosphonic acid, indicating that cleavage of the C1-H bond of the substrate is rate-limiting. This step is also rate-limiting for PhnZ Y24F, as shown by a significant deuterium KIE value of 2.3 ± 0.1 for ^D(V). In contrast, a different reaction step appears to be rate-limiting for the PhnZ E27A and H62A variants, which exhibited ^D(V) values near unity. A solvent KIE of 2.2 ± 0.3 for ^D₂O(V) is observed for PhnZ. Significant solvent KIE values are also observed for the PhnZ Y24F and E27A variants. In contrast, the PhnZ H62A variant does not show a significant solvent KIE, suggesting that H62 is mediating proton transfer in the transition state. A proton inventory study with PhnZ indicates that 1.5 ± 0.6 protons are in flight in the rate-determining step. Overall, the rate-determining transition state for oxidative C-P bond cleavage by PhnZ is proposed to involve C-H bond cleavage that is coupled to deprotonation of the substrate C1-hydroxyl by H62.

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.

Engineered control of enzyme structural dynamics and function

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

Protein dynamics promote hydride tunnelling in substrate oxidation by aryl-alcohol oxidase

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

Decarboxylation involving a ferryl, propionate, and a tyrosyl group in a radical relay yields heme b

The H₂O₂-dependent oxidative decarboxylation of coproheme III is the final step in the biosynthesis of heme b in many microbes. However, the coproheme decarboxylase reaction mechanism is unclear. The structure of the decarboxylase in complex with coproheme III suggested that the substrate iron, reactive propionates, and an active-site tyrosine convey a net 2e^-/2H⁺ from each propionate to an activated form of H₂O₂ Time-resolved EPR spectroscopy revealed that Tyr-145 formed a radical species within 30 s of the reaction of the enzyme-coproheme complex with H₂O₂ This radical disappeared over the next 270 s, consistent with a catalytic intermediate. Use of the harderoheme III intermediate as substrate or substitutions of redox-active side chains (W198F, W157F, or Y113S) did not strongly affect the appearance or intensity of the radical spectrum measured 30 s after initiating the reaction with H₂O₂, nor did it change the ∼270 s required for the radical signal to recede to ≤10% of its initial intensity. These results suggested Tyr-145 as the site of a catalytic radical involved in decarboxylating both propionates. Tyr-145^• was accompanied by partial loss of the initially present Fe(III) EPR signal intensity, consistent with the possible formation of Fe(IV)=O. Site-specifically deuterated coproheme gave rise to a kinetic isotope effect of ∼2 on the decarboxylation rate constant, indicating that cleavage of the propionate Cβ-H bond was partly rate-limiting. The inferred mechanism requires two consecutive hydrogen atom transfers, first from Tyr-145 to the substrate Fe/H₂O₂ intermediate and then from the propionate Cβ-H to Tyr-145^•.

Theory of single-molecule controlled rotation experiments, predictions, tests, and comparison with stalling experiments in F1-ATPase

A recently proposed chemomechanical group transfer theory of rotary biomolecular motors is applied to treat single-molecule controlled rotation experiments. In these experiments, single-molecule fluorescence is used to measure the binding and release rate constants of nucleotides by monitoring the occupancy of binding sites. It is shown how missed events of nucleotide binding and release in these experiments can be corrected using theory, with F₁-ATP synthase as an example. The missed events are significant when the reverse rate is very fast. Using the theory the actual rate constants in the controlled rotation experiments and the corrections are predicted from independent data, including other single-molecule rotation and ensemble biochemical experiments. The effective torsional elastic constant is found to depend on the binding/releasing nucleotide, and it is smaller for ADP than for ATP. There is a good agreement, with no adjustable parameters, between the theoretical and experimental results of controlled rotation experiments and stalling experiments, for the range of angles where the data overlap. This agreement is perhaps all the more surprising because it occurs even though the binding and release of fluorescent nucleotides is monitored at single-site occupancy concentrations, whereas the stalling and free rotation experiments have multiple-site occupancy.

The quantum mitochondrion and optimal health

A sufficiently complex set of molecules, if subject to perturbation, will self-organize and show emergent behaviour. If such a system can take on information it will become subject to natural selection. This could explain how self-replicating molecules evolved into life and how intelligence arose. A pivotal step in this evolutionary process was of course the emergence of the eukaryote and the advent of the mitochondrion, which both enhanced energy production per cell and increased the ability to process, store and utilize information. Recent research suggest that from its inception life embraced quantum effects such as 'tunnelling' and 'coherence' while competition and stressful conditions provided a constant driver for natural selection. We believe that the biphasic adaptive response to stress described by hormesis-a process that captures information to enable adaptability, is central to this whole process. Critically, hormesis could improve mitochondrial quantum efficiency, improving the ATP/ROS ratio, whereas inflammation, which is tightly associated with the aging process, might do the opposite. This all suggests that to achieve optimal health and healthy aging, one has to sufficiently stress the system to ensure peak mitochondrial function, which itself could reflect selection of optimum efficiency at the quantum level.

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.

Structural and Kinetic Studies of Formate Dehydrogenase from Candida boidinii

The structure of formate dehydrogenase from Candida boidinii (CbFDH) is of both academic and practical interests. First, this enzyme represents a unique model system for studies on the role of protein dynamics in catalysis, but so far these studies have been limited by the availability of structural information. Second, CbFDH and its mutants can be used in various industrial applications (e.g., CO2 fixation or nicotinamide recycling systems), and the lack of structural information has been a limiting factor in commercial development. Here, we report the crystallization and structural determination of both holo- and apo-CbFDH. The free-energy barrier for the catalyzed reaction was computed and indicates that this structure indeed represents a catalytically competent form of the enzyme. Complementing kinetic examinations demonstrate that the recombinant CbFDH has a well-organized reactive state. Finally, a fortuitous observation has been made: the apoenzyme crystal was obtained under cocrystallization conditions with a saturating concentration of both the cofactor (NAD(+)) and inhibitor (azide), which has a nanomolar dissociation constant. It was found that the fraction of the apoenzyme present in the solution is less than 1.7 × 10(-7) (i.e., the solution is 99.9999% holoenzyme). This is an extreme case where the crystal structure represents an insignificant fraction of the enzyme in solution, and a mechanism rationalizing this phenomenon is presented.

Enzymatic Kinetic Isotope Effects from First-Principles Path Sampling Calculations

In this study, we develop and test a method to determine the rate of particle transfer and kinetic isotope effects in enzymatic reactions, specifically yeast alcohol dehydrogenase (YADH), from first-principles. Transition path sampling (TPS) and normal mode centroid dynamics (CMD) are used to simulate these enzymatic reactions without knowledge of their reaction coordinates and with the inclusion of quantum effects, such as zero-point energy and tunneling, on the transferring particle. Though previous studies have used TPS to calculate reaction rate constants in various model and real systems, it has not been applied to a system as large as YADH. The calculated primary H/D kinetic isotope effect agrees with previously reported experimental results, within experimental error. The kinetic isotope effects calculated with this method correspond to the kinetic isotope effect of the transfer event itself. The results reported here show that the kinetic isotope effects calculated from first-principles, purely for barrier passage, can be used to predict experimental kinetic isotope effects in enzymatic systems.

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.

Theory for rates, equilibrium constants, and Brønsted slopes in F1-ATPase single molecule imaging experiments

A theoretical model of elastically coupled reactions is proposed for single molecule imaging and rotor manipulation experiments on F1-ATPase. Stalling experiments are considered in which rates of individual ligand binding, ligand release, and chemical reaction steps have an exponential dependence on rotor angle. These data are treated in terms of the effect of thermodynamic driving forces on reaction rates, and lead to equations relating rate constants and free energies to the stalling angle. These relations, in turn, are modeled using a formalism originally developed to treat electron and other transfer reactions. During stalling the free energy profile of the enzymatic steps is altered by a work term due to elastic structural twisting. Using biochemical and single molecule data, the dependence of the rate constant and equilibrium constant on the stall angle, as well as the Børnsted slope are predicted and compared with experiment. Reasonable agreement is found with stalling experiments for ATP and GTP binding. The model can be applied to other torque-generating steps of reversible ligand binding, such as ADP and Pi release, when sufficient data become available.

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.

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.