Hydrogen tunnelling in enzyme-catalysed H-transfer reactions: flavoprotein and quinoprotein systems

It is now widely accepted that enzyme-catalysed C-H bond breakage occurs by quantum mechanical tunnelling. This paradigm shift in the conceptual framework for these reactions away from semi-classical transition state theory (TST, i.e. including zero-point energy, but with no tunnelling correction) has been driven over the recent years by experimental studies of the temperature dependence of kinetic isotope effects (KIEs) for these reactions in a range of enzymes, including the tryptophan tryptophylquinone-dependent enzymes such as methylamine dehydrogenase and aromatic amine dehydrogenase, and the flavoenzymes such as morphinone reductase and pentaerythritol tetranitrate reductase, which produced observations that are also inconsistent with the simple Bell-correction model of tunnelling. However, these data-especially, the strong temperature dependence of reaction rates and the variable temperature dependence of KIEs-are consistent with other tunnelling models (termed full tunnelling models), in which protein and/or substrate fluctuations generate a configuration compatible with tunnelling. These models accommodate substrate/protein (environment) fluctuations required to attain a configuration with degenerate nuclear quantum states and, when necessary, motion required to increase the probability of tunnelling in these states. Furthermore, tunnelling mechanisms in enzymes are supported by atomistic computational studies performed within the framework of modern TST, which incorporates quantum nuclear effects.     

Colworth Medal Lecture. Enzymes in the quantum world

Studies on those enzymes and electron-transfer proteins involved in the catabolism of 'C1' substrates in methylotrophic bacteria have provided a wealth of information concerning the transfer of electrons and hydrogen by quantum tunnelling mechanisms. With regard to H-transfer, studies with MADH have provided the first example of ground-state tunnelling of hydrogen driven by the natural, thermally activated, low-frequency motions of the enzyme molecule. Subsequent studies with related enzymes (e.g. TMADH and bacterial sarcosine oxidase) and with thermophilic alcohol dehydrogenase suggest that vibrationally assisted tunnelling of hydrogen may be more widespread than originally assumed. Our studies of electron transfer in TMADH and ETF have established a role for large-scale protein dynamics in interprotein electron transfer, and have made a contribution to the ongoing debate concerning the mechanism of amine oxidation by enzymes. Moreover, our work has identified a hitherto unknown mechanism for the control of electron density in reduced flavin that influences the rate of electron transfer between redox centres within a protein molecule. Despite this progress, however, many questions still remain to be resolved. With the development of more sophisticated experimental techniques (and also continued financial support from the funding agencies!), the mechanistic uncertainties surrounding the quantum mechanical transfer of electrons and hydrogen in biological molecules should be transmogrified into the certainties one more readily acquaints with the classical world.     

Protein dynamics and enzyme catalysis: the ghost in the machine?

One of the most controversial questions in enzymology today is whether protein dynamics are significant in enzyme catalysis. A particular issue in these debates is the unusual temperature-dependence of some kinetic isotope effects for enzyme-catalysed reactions. In the present paper, we review our recent model [Glowacki, Harvey and Mulholland (2012) Nat. Chem. 4, 169-176] that is capable of reproducing intriguing temperature-dependences of enzyme reactions involving significant quantum tunnelling. This model relies on treating multiple conformations of the enzyme-substrate complex. The results show that direct 'driving' motions of proteins are not necessary to explain experimental observations, and show that enzyme reactivity can be understood and accounted for in the framework of transition state theory.     

Unusual origins of isotope effects in enzyme-catalysed reactions

High hydrostatic pressure is a neglected tool for probing the origins of isotope effects. In chemical reactions, normal primary deuterium isotope effects (DIEs) arising solely from differences in zero point energies are unaffected by pressure; but some anomalous isotope effects in which hydrogen tunnelling is suspected are partially suppressed. In some enzymatic reactions, high pressure completely suppresses the DIE. We have now measured the effects of high pressure on the parallel 13C heavy atom isotope effect of yeast alcohol dehydrogenase and found that it is also suppressed by high pressure and, similarly, suppressed in its entirety. Moreover, the volume changes associated with the suppression of both deuterium and heavy atom isotope effects are virtually identical. The equivalent decrease in activation volumes for hydride transfer, when one mass unit is added to the carbon end of a scissile C-H bond as when one mass unit is added to the hydrogen end, suggests a common origin. Given that carbon is highly unlikely to undergo tunnelling, it follows that hydrogen is not doing so either. The origin of these isotope effects must lie elsewhere. We offer protein domain motions as a possibility.     

Binding isotope effects: boon and bane

Kinetic isotope effects are increasingly applied to investigate enzyme reactions and have been used to understand transition state structure, reaction mechanisms, quantum mechanical hydride ion tunneling and to design transition state analogue inhibitors. Binding isotope effects are an inherent part of most isotope effect measurements but are usually assumed to be negligible. More detailed studies have established surprisingly large binding isotope effects with lactate dehydrogenase, hexokinase, thymidine phosphorylase, and purine nucleoside phosphorylase. Binding reactants into catalytic sites immobilizes conformationally flexible groups, polarizes bonds, and distorts bond angle geometry, all of which generate binding isotope effects. Binding isotope effects are easily measured and provide high-resolution and detailed information on the atomic changes resulting from ligand-macromolecular interactions. Although binding isotope effects complicate kinetic isotope effect analysis, they also provide a powerful tool for finding atomic distortion in molecular interactions.     

Molecular basis for complex formation between methylamine dehydrogenase and amicyanin revealed by inverse mutagenesis of an interprotein salt bridge

Methylamine dehydrogenase (MADH) and amicyanin form a physiologic complex which is required for interprotein electron transfer. The crystal structure of this protein complex is known, and the importance of certain residues on amicyanin in its interaction with MADH has been demonstrated by site-directed mutagenesis. In this study, site-directed mutagenesis of MADH, kinetic data, and thermodynamic analysis are used to probe the molecular basis for stabilization of the protein complex by an interprotein salt bridge between Arg99 of amicyanin and Asp180 of the alpha subunit of MADH. This paper reports the first site-directed mutagenesis of MADH, as well as the construction, heterologous expression, and characterization of a six-His-tagged MADH. alpha Asp180 of MADH was converted to arginine to examine the effect on complex formation with native and mutant amicyanins. This mutation had no effect on the parameters for methylamine oxidation by MADH, but significantly affected its interaction with amicyanin. Of the native and mutant proteins that were studied, their observed order of affinity for each other was as follows: native MADH and native amicyanin > native MADH and R99D amicyanin > alpha D180R MADH and native amicyanin > alpha D180R MADH and R99D amicyanin, and alpha D180R MADH and R99L amicyanin. The alpha D180R mutation also eliminated the ionic strength dependence of the reaction of MADH with amicyanin that is observed with wild-type MADH. Interestingly, the inverse mutation pair of alpha D180R MADH and R99D amicyanin did not restore the favorable salt bridge, but instead disrupted complex formation much more severely than did either individual mutation. These results are explained using molecular modeling and thermodynamic analysis of the kinetic data to correlate the energy contributions of specific stabilizing and destabilizing interactions that are present in the wild-type and mutant complexes. A model is also proposed to describe the sequence of events that leads to stable complex formation between MADH and amicyanin.     

Active-site residues are critical for the folding and stability of methylamine dehydrogenase.

Site-directed mutagenesis was used to alter active-site residues of methylamine dehydrogenase (MADH) from Paracoccus denitrificans. Four residues of the beta subunit of MADH which are in close proximity to the tryptophan tryptophylquinone (TTQ) prosthetic group were modified. The crystal structure of MADH reveals that each of these residues participates in hydrogen bonding interactions with other active-site residues, TTQ or water. Relatively conservative mutations which removed the potentially reactive oxygens on the side chains of Thr122, Tyr119, Asp76 and Asp32 each resulted in greatly reduced or undetectable levels of MADH production. The reduction of MADH levels was determined by assays of activity and Western blots of crude extracts with antisera specific for the MADH beta subunit. No activity or cross-reactive protein was detected in extracts of cells expressing D76N, T122A and T122C MADH mutants. Very low levels of active MADH were produced by cells expressing D32N, Y119F, Y119E and Y119K MADH mutants. The Y119F and D32N mutants were purified from cell extracts and found to be significantly less stable than wild-type MADH. Only the T122S MADH mutant was produced at near wild-type levels. Possible roles for these amino acid residues in stabilizing unusual structural features of the MADH beta subunit, protein folding and TTQ biosynthesis are discussed.     

Simulating large nuclear quantum mechanical corrections in hydrogen atom transfer reactions in metalloenzymes

The role of nuclear quantum mechanical effects in enzyme catalysis has recently attracted significant interest both from theoretical and experimental points of view. From a theoretical point of view, it is undoubtedly a challenge to try to account for the observed tunneling in the protein by microscopic simulations without adjustable parameters. One of the most spectacular examples of nuclear quantum mechanical effects is the reaction of lipoxygenase, which is characterized by a very large kinetic isotope effect and, thus, provides an excellent benchmark for simulation approaches. In the present study, we report a microscopic simulation of the large kinetic isotope effect in soybean lipoxygenase and its temperature dependence. This is, to the best of our knowledge, the first time that a very large nuclear quantum mechanical contribution to the activation free energy of a hydrogen atom transfer reaction and its temperature dependence have been evaluated by microscopic simulation. The simulation reproduces quite well the experimental kinetic information and the corresponding difference between the classical and quantum mechanical activation free energies for the H and D transfer reactions.     

Linking protein structure and dynamics to catalysis: the role of hydrogen tunnelling

Early studies of enzyme-catalysed hydride transfer reactions indicated kinetic anomalies that were initially interpreted in the context of a 'tunnelling correction'. An alternate model for tunnelling emerged following studies of the hydrogen atom transfer catalysed by the enzyme soybean lipoxygenase. This invokes full tunnelling of all isotopes of hydrogen, with reaction barriers reflecting the heavy atom, environmental reorganization terms. Using the latter approach, we offer an integration of the aggregate data implicating hydrogen tunnelling in enzymes (i.e. deviations from Swain-Schaad relationships and the semi-classical temperature dependence of the hydrogen isotope effect). The impact of site-specific mutations of enzymes plays a critical role in our understanding of the factors that control tunnelling in enzyme reactions.     

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.     

Enzymatic H-transfer requires vibration-driven extreme tunneling

Enzymatic breakage of the substrate C-H bond by Methylophilus methyltrophus (sp. W3A1) methylamine dehydrogenase (MADH) has been studied by stopped-flow spectroscopy. The rate of reduction of the tryptophan tryptophylquinone (TTQ) cofactor has a large kinetic isotope effect (KIE = 16.8 +/- 0.5), and the KIE is independent of temperature. Analysis of the temperature dependence of C-H bond breakage revealed that extreme (ground state) quantum tunneling is responsible for the transfer of the hydrogen nucleus. Reaction rates are strongly dependent on temperature, indicating thermally induced, vibrational motion drives the H-transfer reaction. The data provide direct experimental evidence for enzymatic bond breakage by extreme tunneling driven by vibrational motion of the protein scaffold. The results demonstrate that classical transition state theory and its tunneling derivatives do not adequately describe this enzymatic reaction.     

Advances in studying bioinorganic reaction mechanisms: isotopic probes of activated oxygen intermediates in metalloenzymes

Metalloenzymes catalyze reactions of molecular oxygen and its reduced forms through the controlled formation of metal-bound, activated oxygen intermediates. These intermediates have been a challenge to characterize and new experimental approaches capable of relating structure to reactivity under physiologically relevant conditions are needed. The application of a competitive isotope fractionation technique has enabled changes in O-O bonding to be probed during enzyme-catalyzed reactions. The derived isotope effects provide insights into the reaction mechanisms of O2 and O2*-, which probably could not have been obtained using more conventional methods.     

Quantum catalysis in B12-dependent methylmalonyl-CoA mutase: experimental and computational insights

B12-dependent methylmalonyl-CoA mutase catalyses the interchange of a hydrogen atom and the carbonyl-CoA group on adjacent carbons of methylmalonyl-CoA to give the rearranged product, succinyl-CoA. The first step in this reaction involves the transient generation of cofactor radicals by homolytic rupture of the cobalt-carbon bond to generate the deoxyadenosyl radical and cob(II)alamin. This step exhibits a curious sensitivity to isotopic substitution in the substrate, methylmalonyl-CoA, which has been interpreted as evidence for kinetic coupling. The magnitude of the isotopic discrimination is large and a deuterium isotope effect ranging from 35.6 at 20 degrees C to 49.9 at 5 degrees C has been recorded. Arrhenius analysis of the temperature dependence of this isotope effect provides evidence for quantum tunnelling in this hydrogen transfer step. The mechanistic complexity of the observed rate constant for cobalt-carbon bond homolysis together with the spectroscopically silent nature of many of the component steps limits the insights that can be derived by experimental approaches alone. Computational studies using a newly developed geometry optimization scheme that allows determination of the transition state in the full quantum mechanical/molecular mechanical coordinate space have yielded novel insights into the strategy deployed for labilizing the cobalt-carbon bond and poising the resulting deoxyadenosyl radical for subsequent hydrogen atom abstraction.     

Effects of high pressure on enzymatic activity

Effects of high pressure on enzymatic reactions are poised to revolutionize enzyme kinetics. The reason for this is that experimental designs are at hand to separate effects on equilibria between reactant states from effects on catalytic transition states and both yield new information. The first of the former runs contrary to Pauling's hypothesis that substrates are bound more tightly in the transition state, while the latter penetrates the 'black box' of catalysis, the stabilized transition state itself, and returns a precise measure of a physical parameter, deltaV. This in turn opens the door to new forms of structure-activity relationships. The first of these has been described, the effect of pressure on isotope effects, with the surprising finding that the entire isotope effect comes from a transition state phenomenon such as quantum mechanical hydrogen tunneling.     

Incorporation of hydrostatic pressure into models of hydrogen tunneling highlights a role for pressure-modulated promoting vibrations

Hydrostatic pressure offers an alternative to temperature as an experimental probe of hydrogen-transfer reactions. H tunneling reactions have been shown to exhibit kinetic isotope effects (KIEs) that are sensitive to pressure, and environmentally coupled H tunneling reactions, those reactions in which H transfer is coupled to atomic fluctuations (a promoting vibration) along the reaction coordinate, often have quite temperature-dependent KIEs. We present here a theoretical treatment of the combined effect of temperature and pressure on environmentally coupled H tunneling reactions. We develop a generalized expression for the KIE, which can be used as a simple fitting function for combined experimental temperature- and pressure-dependent KIE data sets. With this expression, we are able to extract information about the pressure dependence of both the apparent tunneling distance and the frequency of the promoting vibration. The KIE expression is tested on two data sets {the reduction of chloranil by leuco crystal violet [Isaacs, N. S., Javaid, K., and Rannala, E. (1998) J. Chem. Soc., Perkin Trans. 2, 709-711] and the reduction of morphinone reductase by NADH [Hay, S., Sutcliffe, M. J., and Scrutton, N. S. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 507-512]} and suggests that hydrostatic pressure is a sensitive probe of nuclear quantum mechanical effects in H-transfer reactions.     

Arrhenius curves of hydrogen transfers: tunnel effects, isotope effects and effects of pre-equilibria

In this paper, the Arrhenius curves of selected hydrogen-transfer reactions for which kinetic data are available in a large temperature range are reviewed. The curves are discussed in terms of the one-dimensional Bell-Limbach tunnelling model. The main parameters of this model are the barrier heights of the isotopic reactions, barrier width of the H-reaction, tunnelling masses, pre-exponential factor and minimum energy for tunnelling to occur. The model allows one to compare different reactions in a simple way and prepare the kinetic data for more-dimensional treatments. The first type of reactions is concerned with reactions where the geometries of the reacting molecules are well established and the kinetic data of the isotopic reactions are available in a large temperature range. Here, it is possible to study the relation between kinetic isotope effects (KIEs) and chemical structure. Examples are the tautomerism of porphyrin, the porphyrin anion and related compounds exhibiting intramolecular hydrogen bonds of medium strength. We observe pre-exponential factors of the order of kT/h congruent with 10(13) s-1 corresponding to vanishing activation entropies in terms of transition state theory. This result is important for the second type of reactions discussed in this paper, referring mostly to liquid solutions. Here, the reacting molecular configurations may be involved in equilibria with non- or less-reactive forms. Several cases are discussed, where the less-reactive forms dominate at low or at high temperature, leading to unusual Arrhenius curves. These cases include examples from small molecule solution chemistry like the base-catalysed intramolecular H-transfer in diaryltriazene, 2-(2'-hydroxyphenyl)-benzoxazole, 2-hydroxy-phenoxyl radicals, as well as in the case of an enzymatic system, thermophilic alcohol dehydrogenase. In the latter case, temperature-dependent KIEs are interpreted in terms of a transition between two regimes with different temperature-independent KIEs.     

Hydrogen tunneling in quinoproteins

It is now widely accepted that substrate C-H bond breakage by quinoprotein enzymes occurs by quantum mechanical tunneling. This paradigm shift in the conceptual framework for these reactions away from semi-classical transition state theory (i.e., including zero-point energy but with no tunneling correction) has been driven over recent years by experimental studies of the temperature dependence of kinetic isotope effects for these reactions in the TTQ-dependent enzymes methylamine dehydrogenase and aromatic amine dehydrogenase, which produced observations also inconsistent with the simple Bell correction model of tunneling. However, these data-specifically, the strong temperature dependence of reaction rates and the variable temperature dependence of kinetic isotope effects-are consistent with other tunneling models (denoted full tunneling models) in which protein and/or substrate fluctuations generate a configuration compatible with tunneling. These models accommodate substrate/protein (environment) fluctuations required to attain a configuration with degenerate quantum states and, when necessary, motion required to increase the probability of tunneling in these states. Furthermore, tunneling mechanisms in quinoproteins are supported by computational studies employing variational transition state theory with multidimensional tunneling corrections; these studies are also discussed in this review. Potential pitfalls in analyzing the temperature dependence of kinetic isotope effects as probes of tunneling are also discussed with reference to PQQ-dependent methanol dehydrogenase.     

H-transfers in Photosystem II: what can we learn from recent lessons in the enzyme community?

Over the last 10 years, studies of enzyme systems have demonstrated that, in many cases, H-transfers occur by a quantum mechanical tunneling mechanism analogous to long-range electron transfer. H-transfer reactions can be described by an extension of Marcus theory and, by substituting hydrogen with deuterium (or even tritium), it is possible to explore this theory in new ways by employing kinetic isotope effects. Because hydrogen has a relatively short deBroglie wavelength, H-transfers are controlled by the width of the reaction barrier. By coupling protein dynamics to the reaction coordinate, enzymes have the potential ability to facilitate more efficient H-tunneling by modulating barrier properties. In this review, we describe recent advances in both experimental and theoretical studies of enzymatic H-transfer, in particular the role of protein dynamics or promoting motions. We then discuss possible consequences with regard to tyrosine oxidation/reduction kinetics in Photosystem II.     

Effectors of Escherichia coli aspartate transcarbamoylase differentially perturb aspartate binding rather than the T-R transition

New systematic methods developed for equilibrium isotope exchange kinetics have been used to analyze the effects of activator ATP and inhibitor CTP with Escherichia coli aspartate transcarbamoylase. This indepth approach requires (a) variation of [modifier] with fixed subsaturating levels of substrates, and (b) variation of at least three combinations of reactant-product pairs in constant ratio at equilibrium: [A,B,P,Q], [A,P], and [B,Q] with the co-substrates held constant, in the presence and absence of added modifier. Both ATP and CTP had much stronger effects on the [14C]Asp in equilibrium C-Asp exchange rate than on [32P]C-P in equilibrium Pi. The bisubstrate analog N-phosphonacetyl-L-aspartate activated, then inhibited, Asp in equilibrium C-Asp more strongly than C-P in equilibrium Pi. N-Phosphonacetyl-L-aspartate gave complete (100%) inhibition, whereas CTP inhibition of either exchange was only partial. Substrate saturation curves in the presence and absence of effectors indicate that ATP and CTP perturb the observed values of Rmax and S0.5 in different fashions without appreciably changing the observed Hill number. Computer simulations indicate that the primary site of ATP and CTP action is the association rate for Asp, not the allosteric T-R transition. This finding is substantiated by previous studies in which modified aspartate transcarbamoylase had lost cooperative Asp binding without loss of sensitivity to effectors, or in which sensitivity to one effector could be deleted selectively. The present results, with newly devised computer simulation and analysis methods, illustrate the usefulness of equilibrium isotope exchange kinetic probes for providing unique insights to enzyme regulatory mechanisms, to define exactly which steps are altered in a given kinetic mechanism.