Are environmentally coupled enzymatic hydrogen tunneling reactions influenced by changes in solution viscosity?

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

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

UbiD domain dynamics underpins aromatic decarboxylation

The widespread UbiD enzyme family utilises the prFMN cofactor to achieve reversible decarboxylation of acrylic and (hetero)aromatic compounds. The reaction with acrylic compounds based on reversible 1,3-dipolar cycloaddition between substrate and prFMN occurs within the confines of the active site. In contrast, during aromatic acid decarboxylation, substantial rearrangement of the substrate aromatic moiety associated with covalent catalysis presents a molecular dynamic challenge. Here we determine the crystal structures of the multi-subunit vanillic acid decarboxylase VdcCD. We demonstrate that the small VdcD subunit acts as an allosteric activator of the UbiD-like VdcC. Comparison of distinct VdcCD structures reveals domain motion of the prFMN-binding domain directly affects active site architecture. Docking of substrate and prFMN-adduct species reveals active site reorganisation coupled to domain motion supports rearrangement of the substrate aromatic moiety. Together with kinetic solvent viscosity effects, this establishes prFMN covalent catalysis of aromatic (de)carboxylation is afforded by UbiD dynamics.

Influence of local microenvironment on the double hydrogen transfer in porphycene

We performed time-resolved transient absorption and fluorescence anisotropy measurements in order to study tautomerization of porphycene in rigid polymer matrices at cryogenic temperatures. Studies were carried out in poly(methyl methacrylate) (PMMA), poly(vinyl butyral) (PVB), and poly(vinyl alcohol) (PVA). The results prove that in all studied media hydrogen tunnelling plays a significant role in the double hydrogen transfer which becomes very sensitive to properties of the environment below approx. 150 K. We also demonstrate that there exist two populations of porphycene molecules in rigid media: "hydrogen-transferring" molecules, in which tautomerization occurs on time scales below 1 ns and "frozen" molecules in which double hydrogen transfer is too slow to be monitored with nanosecond techniques. The number of "frozen" molecules increases when the sample is cooled. We explain this effect by interactions of guest molecules with a rigid host matrix which disturbs symmetry of porphycene and hinders tunnelling. Temperature dependence of the number of hydrogen-transferring molecules suggests that the factor which restores the symmetry of the double-minimum potential well in porphycene are intermolecular vibrations localized in separated regions of the amorphous polymer.

Structures of carboxylic acid reductase reveal domain dynamics underlying catalysis

Carboxylic acid reductase (CAR) catalyzes the ATP- and NADPH-dependent reduction of carboxylic acids to the corresponding aldehydes. The enzyme is related to the nonribosomal peptide synthetases, consisting of an adenylation domain fused via a peptidyl carrier protein (PCP) to a reductase termination domain. Crystal structures of the CAR adenylation-PCP didomain demonstrate that large-scale domain motions occur between the adenylation and thiolation states. Crystal structures of the PCP-reductase didomain reveal that phosphopantetheine binding alters the orientation of a key Asp, resulting in a productive orientation of the bound nicotinamide. This ensures that further reduction of the aldehyde product does not occur. Combining crystallography with small-angle X-ray scattering (SAXS), we propose that molecular interactions between initiation and termination domains are limited to competing PCP docking sites. This theory is supported by the fact that (R)-pantetheine can support CAR activity for mixtures of the isolated domains. Our model suggests directions for further development of CAR as a biocatalyst.

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.

Evidence for Dominant Role of Tunneling in Condensed Phases and at High Temperatures: Double Hydrogen Transfer in Porphycenes

Investigation of the double hydrogen transfer in porphycene, its 2,7,12,17-tetra-tert-butyl derivative, and their N-deuterated isotopologues revealed the dominant role of tunneling, even at room temperature in condensed phase. Ultrafast optical spectroscopy with polarized light employed in a wide range of temperatures allowed the identification and evaluation of contributions of two tunneling modes: vibrational ground-state tunneling, occurring from the zero vibrational level, and vibrationally activated, via a large amplitude, low-frequency mode. Good correspondence was found between the rates of incoherent tunneling occurring in condensed phase and the values estimated on the basis of tunneling splittings observed in molecules isolated in supersonic jets or helium nanodroplets. The results provide solid experimental insight into widely proposed quantum facets of ubiquitous hydrogen-transfer phenomena.

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.

Multiple intermediates, diverse conformations, and cooperative conformational changes underlie the catalytic hydride transfer reaction of dihydrofolate reductase

It has become increasingly clear that protein motions play an essential role in enzyme catalysis. However, exactly how these motions are related to an enzyme's chemical step is still intensely debated. This chapter examines the possible role of protein motions that display a hierarchy of timescales in enzyme catalysis. The linkage between protein motions and catalysis is investigated in the context of a model enzyme, E. coli dihydrofolate reductase (DHFR), that catalyzes the hydride transfer reaction in the conversion of dihydrofolate to tetrahydrofolate. The results of extensive computer simulations probing the protein motions that are manifest during different steps along the turnover cycle of DHFR are summarized. Evidence is presented that the protein motions modulate the catalytic efficacy of DHFR by generating a conformational ensemble conducive to the hydride transfer. The alteration of the equilibrium conformational ensemble rather than any protein dynamical effects is found to be sufficient to explain the rate-diminishing effects of mutation on the kinetics of the enzyme. These data support the view that the protein motions facilitate catalysis by establishing reaction competent conformations of the enzyme, but they do not directly couple to the chemical reaction itself. These findings have broad implications for our understanding of enzyme mechanisms and the design of novel protein catalysts.

Relationship of femtosecond-picosecond dynamics to enzyme-catalyzed H-transfer

At physiological temperatures, enzymes exhibit a broad spectrum of conformations, which interchange via thermally activated dynamics. These conformations are sampled differently in different complexes of the protein and its ligands, and the dynamics of exchange between these conformers depends on the mass of the group that is moving and the length scale of the motion, as well as restrictions imposed by the globular fold of the enzymatic complex. Many of these motions have been examined and their role in the enzyme function illuminated, yet most experimental tools applied so far have identified dynamics at time scales of seconds to nanoseconds, which are much slower than the time scale for H-transfer between two heavy atoms. This chemical conversion and other processes involving cleavage of covalent bonds occur on picosecond to femtosecond time scales, where slower processes mask both the kinetics and dynamics. Here we present a combination of kinetic and spectroscopic methods that may enable closer examination of the relationship between enzymatic C-H → C transfer and the dynamics of the active site environment at the chemically relevant time scale. These methods include kinetic isotope effects and their temperature dependence, which are used to study the kinetic nature of the H-transfer, and 2D IR spectroscopy, which is used to study the dynamics of transition-state- and ground-state-analog complexes. The combination of these tools is likely to provide a new approach to examine the protein dynamics that directly influence the chemical conversion catalyzed by enzymes.

The role of large-scale motions in catalysis by dihydrofolate reductase

Dihydrofolate reductase has long been used as a model system to study the coupling of protein motions to enzymatic hydride transfer. By studying environmental effects on hydride transfer in dihydrofolate reductase (DHFR) from the cold-adapted bacterium Moritella profunda (MpDHFR) and comparing the flexibility of this enzyme to that of DHFR from Escherichia coli (EcDHFR), we demonstrate that factors that affect large-scale (i.e., long-range, but not necessarily large amplitude) protein motions have no effect on the kinetic isotope effect on hydride transfer or its temperature dependence, although the rates of the catalyzed reaction are affected. Hydrogen/deuterium exchange studies by NMR-spectroscopy show that MpDHFR is a more flexible enzyme than EcDHFR. NMR experiments with EcDHFR in the presence of cosolvents suggest differences in the conformational ensemble of the enzyme. The fact that enzymes from different environmental niches and with different flexibilities display the same behavior of the kinetic isotope effect on hydride transfer strongly suggests that, while protein motions are important to generate the reaction ready conformation, an optimal conformation with the correct electrostatics and geometry for the reaction to occur, they do not influence the nature of the chemical step itself; large-scale motions do not couple directly to hydride transfer proper in DHFR.

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.

Mutagenesis alters the catalytic mechanism of the light-driven enzyme protochlorophyllide oxidoreductase

The light-activated enzyme protochlorophyllide oxidoreductase (POR) catalyzes an essential step in the synthesis of the most abundant pigment on Earth, chlorophyll. This unique reaction involves the sequential addition of a hydride and proton across the C17=C18 double bond of protochlorophyllide (Pchlide) by dynamically coupled quantum tunneling and is an important model system for studying the mechanism of hydrogen transfer reactions. In the present work, we have combined site-directed mutagenesis studies with a variety of sensitive spectroscopic and kinetic measurements to provide new insights into the mechanistic role of three universally conserved Cys residues in POR. We show that mutation of Cys-226 dramatically alters the catalytic mechanism of the enzyme. In contrast to wild-type POR, the characteristic charge-transfer intermediate, formed upon hydride transfer from NADPH to the C17 position of Pchlide, is absent in C226S variant enzymes. This suggests a concerted hydrogen transfer mechanism where proton transfer only is rate-limiting. Moreover, Pchlide reduction does not require the network of solvent-coupled conformational changes that play a key role in the proton transfer step of wild-type POR. We conclude that this globally important enzyme is finely tuned to facilitate efficient photochemistry, and the removal of a key interaction with Pchlide in the C226S variants significantly affects the local active site structure in POR, resulting in a shorter donor-acceptor distance for proton transfer.

Correction of pre-steady-state KIEs for isotopic impurities and the consequences of kinetic isotope fractionation

We show, both experimentally and by kinetic modeling, that enzymatic single-turnover (pre-steady-state) H-transfer reactions can be significantly complicated by kinetic isotope fractionation. This fractionation results in the formation of more protiated than deuterated product and is a unique problem for pre-steady-state reactions. When observed rate constants are measured using rapid-mixing (e.g., stopped flow) methodologies, kinetic isotope fractionation can lead to a large underestimation of both the magnitude and temperature dependence of kinetic isotope effects (KIEs). This fractionation is related to the isotopic purity of the substrates used and highlights a major problem with experimental studies which measure KIEs with substrates that are not isotopically pure. As it is not always possible to prepare isotopically pure substrates, we describe two general methods for the correction, for known isotope impurities, of KIEs calculated from pre-steady-state measurements.

Solvent-slaved protein motions accompany proton but not hydride tunneling in light-activated protochlorophyllide oxidoreductase

H(+) but not H(-): The reduction reaction of protochlorophyllide catalyzed by protochlorophyllide oxidoreductase features solvent-slaved motions that control the proton- but not the hydride-tunneling mechanism. These motions imply a long-range dynamic network from the solvent to the enzyme active site that facilitate proton transfer (see picture, left). Motions for hydride transfer are more localized and are not slaved by the solvent (see picture, right).

Barrier compression enhances an enzymatic hydrogen-transfer reaction

Putting the squeeze on: Hydrostatic pressure causes a shortening of the charge-transfer bond in the binary complex of morphinone reductase and NADH(4) (see diagram). Molecular dynamics simulations suggest that pressure reduces the average reaction barrier width by restricting the conformational space available to the flavin mononucleotide and NADH within the active site. The apparent rate of catalysis increases with pressure.

Probing coupled motions in enzymatic hydrogen tunnelling reactions

Much work has gone into understanding the physical basis of the enormous catalytic power of enzymes over the last 50 years or so. Nevertheless, the detailed mechanism used by Nature's catalysts to speed chemical transformations remains elusive. DHFR (dihydrofolate reductase) has served as a paradigm to study the relationship between the structure, function and dynamics of enzymatic transformations. A complex reaction cascade, which involves rearrangements and movements of loops and domains of the enzyme, is used to orientate cofactor and substrate in a reactive configuration from which hydride is transferred by quantum mechanical tunnelling. In the present paper, we review results from experiments that probe the influence of protein dynamics on the chemical step of the reaction catalysed by TmDHFR (DHFR from Thermotoga maritima). This enzyme appears to have evolved an optimal structure that can maintain a catalytically competent conformation under extreme conditions.

Solvent effects on environmentally coupled hydrogen tunnelling during catalysis by dihydrofolate reductase from Thermotoga maritima

Protein motions may be perturbed by altering the properties of the reaction medium. Here we show that dielectric constant, but not viscosity, affects the rate of the hydride-transfer reaction catalysed by dihydrofolate reductase from Thermotoga maritima (TmDHFR), in which quantum-mechanical tunnelling has previously been shown to be driven by protein motions. Neither dielectric constant nor viscosity directly alters the kinetic isotope effect of the reaction or the mechanism of coupling of protein motions to tunnelling. Glycerol and sucrose cause a significant increase in the rate of hydride transfer, but lead to a reduction in the magnitude of the kinetic isotope effect as well as an extension of the temperature range over which "passive" protein dynamics (rather than "active" gating motions) dominate the reaction. Our results are in agreement with the proposal that non-equilibrium dynamical processes (promoting motions) drive the hydride-transfer reaction in TmDHFR.

Solvent as a probe of active site motion and chemistry during the hydrogen tunnelling reaction in morphinone reductase

The reductive half-reaction of morphinone reductase involves a hydride transfer from enzyme-bound beta-nicotinamide adenine dinucleotide (NADH) to a flavin mononucleotide (FMN). We have previously demonstrated that this step proceeds via a quantum mechanical tunnelling mechanism. Herein, we probe the effect of the solvent on the active site chemistry. The pK(a) of the reduced FMN N1 is 7.4+/-0.7, based on the pH-dependence of the FMN midpoint potential. We rule out that protonation of the reduced FMN N1 is coupled to the preceding H-transfer as both the rate and temperature-dependence of the reaction are insensitive to changes in solution pH above and below this pK(a). Further, the solvent kinetic isotope effect is approximately 1.0 and both the 1 degrees and 2 degrees KIEs are insensitive to solution pH. The effect of the solvent's dielectric constant is investigated and the rate of H-transfer is found to be unaffected by changes in the dielectric constant between approximately 60 and 80. We suggest that, while there is crystallographic evidence for some water in the active site, the putative promoting motion involved in the H-tunnelling reaction is insensitive to such changes.

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.