Structural and mechanistic aspects of flavoproteins: probes of hydrogen tunnelling

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 Hydride Transfer Process in NADP-dependent Methylene-tetrahydromethanopterin Dehydrogenase

NADP-dependent methylene-tetrahydromethanopterin (methylene-H₄MPT) dehydrogenase (MtdA) catalyzes the reversible dehydrogenation of methylene-H₄MPT to form methenyl-H₄MPT⁺ by using NADP⁺ as a hydride acceptor. This hydride transfer reaction is involved in the oxidative metabolism from formaldehyde to CO₂ in methylotrophic and methanotrophic bacteria. Here, we report on the crystal structures of the ternary MtdA-substrate complexes from Methylorubrum extorquens AM1 obtained in open and closed forms. Their conversion is accomplished by opening/closing the active site cleft via a 15° rotation of the NADP, relative to the pterin domain. The 1.08 Å structure of the closed and active enzyme-NADP-methylene-H₄MPT complex allows a detailed geometric analysis of the bulky substrates and a precise prediction of the hydride trajectory. Upon domain closure, the bulky substrate rings become compressed resulting in a tilt of the imidazolidine group of methylene-H₄MPT that optimizes the geometry for hydride transfer. An additional 1.5 Å structure of MtdA in complex with the nonreactive NADP⁺ and methenyl-H₄MPT⁺ revealed an extremely short distance between nicotinamide-C4 and imidazoline-C14a of 2.5 Å, which demonstrates the strong pressure imposed. The pterin-imidazolidine-phenyl butterfly angle of methylene-H₄MPT bound to MtdA is smaller than that in the enzyme-free state but is similar to that in H₂- and F₄₂₀-dependent methylene-H₄MPT dehydrogenases. The concept of compression-driven hydride transfer including quantum mechanical hydrogen tunneling effects, which are established for flavin- and NADP-dependent enzymes, can be expanded to hydride-transferring H₄MPT-dependent enzymes.

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.

On the use of noncompetitive kinetic isotope effects to investigate flavoenzyme mechanism

This account describes the application of kinetic isotope effects (KIEs) to investigate the mechanistic properties of flavin dependent enzymes. Assays can be conducted during steady-state catalytic turnover of the flavoenzyme with its substrate or by using rapid-kinetic techniques to measure either the reductive or oxidative half-reactions of the enzyme. Great care should be taken to ensure that the observed effects are due to isotopic substitution and not other factors such as pH effects or changes in the solvent viscosity of the reaction mixture. Different types of KIEs are described along with a physical description of their origins and the unique information each can provide about the mechanism of an enzyme. Detailed experimental techniques are outlined with special emphasis on the proper controls and data analysis that must be carried out to avoid erroneous conclusions. Examples are provided for each type of KIE measurement from references in the literature. It is our hope that this article will clarify any confusion concerning the utility of KIEs in the study of flavoprotein mechanism and encourage their use by the community.

Isotopically labeled flavoenzymes and their uses in probing reaction mechanisms

The incorporation of stable isotopes into proteins is beneficial or essential for a range of experiments, including NMR, neutron scattering and reflectometry, proteomic mass spectrometry, vibrational spectroscopy and "heavy" enzyme kinetic isotope effect (KIE) measurements. Here, we present detailed protocols for the stable isotopic labeling of pentaerythritol tetranitrate reductase (PETNR) via recombinant expression in E. coli. PETNR is an ene-reductase belonging to the Old Yellow Enzyme (OYE) family of flavoenzymes, and is regarded as a model system for studying hydride transfer reactions. Included is a discussion of how efficient back-exchange of amide protons in the protein core can be achieved and how the intrinsic flavin mononucleotide (FMN) cofactor can be exchanged, allowing the production of isotopologues with differentially labeled protein and cofactor. In addition to a thorough description of labeling strategies, we briefly exemplify how data analysis and interpretation of "heavy" enzyme KIEs can be performed.

Nonequivalence of Second Sphere "Noncatalytic" Residues in Pentaerythritol Tetranitrate Reductase in Relation to Local Dynamics Linked to H-Transfer in Reactions with NADH and NADPH Coenzymes

Many enzymes that catalyze hydride transfer reactions work via a mechanism dominated by quantum mechanical tunneling. The involvement of fast vibrational modes of the reactive complex is often inferred in these reactions, as in the case of the NAD(P)H-dependent pentaerythritol tetranitrate reductase (PETNR). Herein, we interrogated the H-transfer mechanism in PETNR by designing conservative (L25I and I107L) and side chain shortening (L25A and I107A) PETNR variants and using a combination of experimental approaches (stopped-flow rapid kinetics, X-ray crystallography, isotope/temperature dependence studies of H-transfer and NMR spectroscopy). X-ray data show subtle changes in the local environment of the targeted side chains but no major structural perturbation caused by mutagenesis of these two second sphere active site residues. However, temperature dependence studies of H-transfer revealed a coenzyme-specific and complex thermodynamic equilibrium between different reactive configurations in PETNR–coenzyme complexes. We find that mutagenesis of these second sphere “noncatalytic” residues affects differently the reactivity of PETNR with NADPH and NADH coenzymes. We attribute this to subtle, dynamic structural changes in the PETNR active site, the effects of which impact differently in the nonequivalent reactive geometries of PETNR−NADH and PETNR−NADPH complexes. This inference is confirmed through changes observed in the NMR chemical shift data for PETNR complexes with unreactive 1,4,5,6-tetrahydro-NAD(P) analogues. We show that H-transfer rates can (to some extent) be buffered through entropy–enthalpy compensation, but that use of integrated experimental tools reveals hidden complexities that implicate a role for dynamics in this relatively simple H-transfer reaction. Similar approaches are likely to be informative in other enzymes to understand the relative importance of (distal) hydrophobic side chains and dynamics in controlling the rates of enzymatic H-transfer.

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.

A qualitative quantum rate model for hydrogen transfer in soybean lipoxygenase

The hydrogen transfer reaction catalysed by soybean lipoxygenase (SLO) has been the focus of intense study following observations of a high kinetic isotope effect (KIE). Today high KIEs are generally thought to indicate departure from classical rate theory and are seen as a strong signature of tunnelling of the transferring particle, hydrogen or one of its isotopes, through the reaction energy barrier. In this paper, we build a qualitative quantum rate model with few free parameters that describes the dynamics of the transferring particle when it is exposed to energetic potentials exerted by the donor and the acceptor. The enzyme's impact on the dynamics is modelled by an additional energetic term, an oscillatory contribution known as "gating." By varying two key parameters, the gating frequency and the mean donor-acceptor separation, the model is able to reproduce well the KIE data for SLO wild-type and a variety of SLO mutants over the experimentally accessible temperature range. While SLO-specific constants have been considered here, it is possible to adapt these for other enzymes.

Nuclear quantum tunnelling in enzymatic reactions – an enzymologist's perspective

The roles of nuclear quantum tunnelling and dynamics in enzyme reactions are discussed in this perspective on H-transfer reactions.

Practical aspects on the use of kinetic isotope effects as probes of flavoprotein enzyme mechanisms

The measurement of kinetic isotope effects (KIEs) has proved useful in many mechanistic studies of enzyme activity, most notably in enzyme-catalyzed hydrogen-transfer reactions. Primary KIEs (1° KIE) greater than unity indicate that transfer of the hydrogen species of interest is partially or fully rate limiting, and studies of the magnitude of the temperature and pressure dependence of these KIEs can inform on the mechanism of transfer. For example, KIE measurements have proved crucial in understanding the role of quantum mechanical tunneling in enzyme systems. The measurement of secondary KIEs (2° KIEs) is also informative and can be used to infer a significant tunneling contribution and details of transition state geometry. Here the deuterium label is introduced next to that of the transferred hydrogen. Measurements of 1° and 2° KIEs are being used increasingly in studies of H-transfer in flavoprotein enzymes and this requires the preparation of high purity and stereospecific labeled isotopologues. Strategies for the synthesis of labeled substrates are dependent on the enzyme system being studied. However, the nicotinamide coenzymes are often used in studies of flavoprotein enzyme mechanisms. Here, we provide practical details for the enzymatic synthesis of high purity deuterated isotopologues of the common biological coenzymes NADH and NADPH as well as the corresponding nonreactive mimics, tetrahydroNAD(P)H. Both forms of the coenzyme have proven useful in the study of mechanisms, particularly in relation to the involvement of quantum mechanical tunneling and dynamics in enzymatic H-transfer chemistry. The focus here is on practical considerations in the synthesis of these compounds. We also provide an abbreviated description of how measurements of KIEs can inform on flavoprotein mechanisms. The aim of this contribution is not to give a detailed description of the underlying theory (which has been reviewed extensively in the literature), but to provide a basic introduction and practical considerations for nonexpert readers who wish to incorporate such measurements in studies of enzyme mechanisms.

Excited state dynamics can be used to probe donor-acceptor distances for H-tunneling reactions catalyzed by flavoproteins

In enzyme systems where fast motions are thought to contribute to H-transfer efficiency, the distance between hydrogen donor and acceptor is a very important factor. Sub-ångstrom changes in donor-acceptor distance can have a large effect on the rate of reaction, so a sensitive probe of these changes is a vital tool in our understanding of enzyme function. In this study we use ultrafast transient absorption spectroscopy to investigate the photoinduced electron transfer rates, which are also very sensitive to small changes in distance, between coenzyme analog, NAD(P)H4, and the isoalloxazine center in the model flavoenzymes morphinone reductase (wild-type and selected variants) and pentaerythritol tetranitrate reductase (wild-type). It is shown that upon addition of coenzyme to the protein the rate of photoinduced electron transfer is increased. By comparing the magnitude of this increase with existing values for NAD(P)H4-FMN distances, based on charge-transfer complex absorbance and experimental kinetic isotope effect reaction data, we show that this method can be used as a sensitive probe of donor-acceptor distance in a range of enzyme systems.

Increased dynamic effects in a catalytically compromised variant of Escherichia coli dihydrofolate reductase

Isotopic substitution ((15)N, (13)C, (2)H) of a catalytically compromised variant of Escherichia coli dihydrofolate reductase, EcDHFR-N23PP/S148A, has been used to investigate the effect of these mutations on catalysis. The reduction of the rate constant of the chemical step in the EcDHFR-N23PP/S148A catalyzed reaction is essentially a consequence of an increase of the quasi-classical free energy barrier and to a minor extent of an increased number of recrossing trajectories on the transition state dividing surface. Since the variant enzyme is less well set up to catalyze the reaction, a higher degree of active site reorganization is needed to reach the TS. Although millisecond active site motions are lost in the variant, there is greater flexibility on the femtosecond time scale. The "dynamic knockout" EcDHFR-N23PP/S148A is therefore a "dynamic knock-in" at the level of the chemical step, and the increased dynamic coupling to the chemical coordinate is in fact detrimental to catalysis. This finding is most likely applicable not just to hydrogen transfer in EcDHFR but also to other enzymatic systems.

Hydrogen transfer reactions in viscous media — Potential and free energy surfaces in solvent–solute coordinates and their kinetic implications

Two-dimensional potential energy and free energy surfaces are obtained using quantum mechanical and molecular dynamics calculations for four hydrogen transfer reactions in n-hexane solvent: the methyl–methane, n-propyl–n-propane, n-pentyl–n-pentane, and t-butyl–isobutane systems. The resultant surfaces have similar landscapes despite the fact the equilibrated solvent cavities for these systems are notably different in size and shape. The kinetic implications of these landscapes are discussed. The Arrhenius and tunneling kinetics of hydrogen transfer in nonpolar hydrocarbon solute–solvent systems are not expected to show any significant viscosity dependence.

Evidence that a 'dynamic knockout' in Escherichia coli dihydrofolate reductase does not affect the chemical step of catalysis

The question of whether protein motions play a role in the chemical step of enzymatic catalysis has generated much controversy in recent years. Debate has recently reignited over possible dynamic contributions to catalysis in dihydrofolate reductase, following conflicting conclusions from studies of the N23PP/S148A variant of the Escherichia coli enzyme. By investigating the temperature dependence of kinetic isotope effects, we present evidence that the reduction in the hydride transfer rate constants in this variant is not a direct result of impairment of conformational fluctuations. Instead, the conformational state of the enzyme immediately before hydride transfer, which determines the electrostatic environment of the active site, affects the rate constant for the reaction. Although protein motions are clearly important for binding and release of substrates and products, there appears to be no detectable dynamic coupling of protein motions to the hydride transfer step itself.

Good vibrations in enzyme-catalysed reactions

Fast motions (femtosecond to picosecond) and their potential involvement during enzyme-catalysed reactions have ignited considerable interest in recent years. Their influence on reaction chemistry has been inferred indirectly from studies of the anomalous temperature dependence of kinetic isotope effects and computational simulations. But can such motion reduce the width and height of energy barriers along the reaction coordinate, and contribute to quantum mechanical and/or classical nuclear-transfer chemistry? Here we discuss contemporary ideas for enzymatic reactions invoking a role for fast 'promoting' (or 'compressive') motions that, in principle, can aid hydrogen-transfer reactions. Of key importance is the direct demonstration of a role for compressive motions and the ability to understand in atomic detail the structural origin of these fast motions, but so far this has not been achieved. Here we discuss both indirect experimental evidence that supports a role for compressive motion and the additional insight gained from computational simulations.

Effect of pH on hydride transfer by Escherichia coli dihydrofolate reductase

The kinetic isotope effect (KIE) on hydride transfer in the reaction catalysed by dihydrofolate reductase from Escherichia coli (EcDHFR) is known to be temperature dependent at pH 7, but essentially independent of temperature at elevated pH. Here, we show that the transition from the temperature-dependent regime to the temperature-independent regime occurs sharply between pH 7.5 and 8. The activation energy for hydride transfer is independent of pH. The mechanism leading to the change in behaviour of the KIEs is not clear, but probably involves a conformational change in the enzyme brought about by deprotonation of a key residue (or residues) at high pH. The KIE on hydride transfer at low pH suggests that the rate constant for the reaction is not limited by a conformational change to the enzyme under these conditions. The effect of pH on the temperature dependence of the rate constants and KIEs for hydride transfer catalysed by EcDHFR suggests that enzyme motions and conformational changes do not directly influence the chemistry, but that the reaction conditions affect the conformational ensemble of the enzyme prior to reaction and control the reaction though this route.

Solvent effects on catalysis by Escherichia coli dihydrofolate reductase

Hydride transfer catalyzed by dihydrofolate reductase (DHFR) has been described previously within an environmentally coupled model of hydrogen tunneling, where protein motions control binding of substrate and cofactor to generate a tunneling ready conformation and modulate the width of the activation barrier and hence the reaction rate. Changes to the composition of the reaction medium are known to perturb protein motions. We have measured kinetic parameters of the reaction catalyzed by DHFR from Escherichia coli in the presence of various cosolvents and cosolutes and show that the dielectric constant, but not the viscosity, of the reaction medium affects the rate of reaction. Neither the primary kinetic isotope effect on the reaction nor its temperature dependence were affected by changes to the bulk solvent properties. These results are in agreement with our previous report on the effect of solvent composition on catalysis by DHFR from the hyperthermophile Thermotoga maritima. However, the effect of solvent on the temperature dependence of the kinetic isotope effect on hydride transfer catalyzed by E. coli DHFR is difficult to explain within a model, in which long-range motions couple to the chemical step of the reaction, but may indicate the existence of a short-range promoting vibration or the presence of multiple nearly isoenergetic conformational substates of enzymes with similar but distinct catalytic properties.

Probing active site geometry using high pressure and secondary isotope effects in an enzyme-catalysed 'deep' H-tunnelling reaction

We report the first study of the effects of hydrostatic pressure on α-2° KIEs for an enzyme-catalysed H-transfer reaction that occurs by 'deep' tunnelling. High pressure causes a significant decrease in the observed α-2° KIE on the pre-steady-state hydride transfer from NADH to FMN in the flavoprotein morphinone reductase. We have recently shown that high pressure causes a reduction in macroscopic reaction barrier width for this reaction. Using DFT vibrational analysis of a simple active site model, we posit that the decrease in α-2° KIE with pressure may arise due to a decrease in the vibrational coupling between the NADH primary (transferred) and secondary hydrogens in the 'tunnelling ready configuration', which more closely resembles the reactant state than the transition state.

Evidence to support the hypothesis that promoting vibrations enhance the rate of an enzyme catalyzed H-tunneling reaction

In recent years there has been a shift away from transition state theory models for H-transfer reactions. Models that incorporate tunneling as the mechanism of H-transfer are now recognized as a better description of such reactions. Central to many models of H-tunneling is the notion that specific vibrational modes of the protein and/or substrate can increase the probability of a H-tunneling reaction, modes that are termed promoting vibrations. Thus far there has been limited evidence that promoting vibrations can increase the rate of H-transfer. In the present communication we examine the single hydride transfer from both NADPH and NADH to FMN in the reductive half-reaction of pentaerythritol tetranitrate reductase (PETNR). We find that there is a significant promoting vibration with NADPH but not with NADH and that the observed rate of hydride transfer is significantly (approximately 15x) faster with NADPH. We rule out differences in rate due to variation in driving force and the donor-acceptor distance, suggesting it is the promoting vibration with NADPH that is the origin of the increased observed rate. This study therefore provides direct evidence that promoting vibrations can lead to an increase in rate.