Hydride transfer reaction catalyzed by hyperthermophilic dihydrofolate reductase is dominated by quantum mechanical tunneling and is promoted by both inter- and intramonomeric correlated motions

Simulations of hydride and deuteride transfer catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) are presented. TmDHFR was modeled with its active homodimeric quaternary structure, where each monomer has three subdomains. The potential energy function was a combined quantum mechanical and molecular mechanical potential (69 atoms were treated quantum mechanically, and 35 287, by molecular mechanics). The calculations of the rate constants by ensemble-averaged variational transition state theory with multidimensional tunneling predicted that hydride and deuteride transfer at 278 K proceeded with 81 and 80% by tunneling. These percentages decreased to 50 and 49% at 338 K. The kinetic isotope effect was dominated by contributions of bound vibrations and decreased from 3.0 to 2.2 over the temperature range. The calculated rates for hydride and deuteride transfer catalyzed by the hypothetical monomer were smaller by approximately 2 orders of magnitude. At 298 K tunneling contributed 73 and 66% to hydride and deuteride transfer in the monomer. The decreased catalytic efficiency of the monomer was therefore not the result of a decrease of the tunneling contribution but an increase in the quasi-classical activation free energy. The catalytic effect was associated in the dimer with correlated motions between domains as well as within and between subunits. The intrasubunit correlated motions were decreased in the monomer when compared to both native dimeric TmDHFR and monomeric E. coli enzyme. TmDHFR and its E. coli homologue involve similar patterns of correlated interactions that affect the free energy barrier of hydride transfer despite only 27% sequence identity and different quaternary structures.     

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 β‐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 ～1.0 and both the 1° and 2° 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 ～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.     

Hydride reduction of NAD+ analogues by isopropyl alcohol: kinetics, deuterium isotope effects and mechanism

Observed pseudo-first-order rate constants (k(obs)) of the hydride-transfer reactions from isopropyl alcohol (i-PrOH) to two NAD(+) analogues, 9-phenylxanthylium ion (PhXn(+)) and 10-methylacridinium ion (MA(+)), were determined at temperatures ranging from 49 to 82 degrees C in i-PrOH containing various amounts of AN or water. Formations of the alcohol-cation ether adducts (ROPr-i) were observed as side equilibria. The equilibrium constants for the conversion of PhXn(+) to PhXnOPr-i in i-PrOH/AN (v/v = 1) were determined, and the equilibrium isotope effect (EIE = K(i-PrOH)/K(i-PrOD)) at 62 degrees C was calculated to be 2.67. The k(H) of the hydride-transfer step for both reactions were calculated on the basis of the k(obs) and K. The corresponding deuterium kinetic isotope effects (e.g., KIE(OD)(H) = k(H)(i-PrOH)/k(H)(i-PrOD) and KIE(beta-D6)(H) = k(obs)(i-PrOH)/k(obs)((CD3)2CHOH)), as well as the activation parameters, were derived. For the reaction of PhXn(+) (62 degrees C) and MA(+) (67 degrees C), primary KIE(alpha-D)(H) (4.4 and 2.1, respectively) as well as secondary KIE(OD)(H) (1.07 and 1.18) and KIE(beta-D6)(H) (1.1 and 1.5) were observed. The observed EIE and KIE(OD)(H) were explained in terms of the fractionation factors for deuterium between OH and OH(+)(OH(delta+)) sites. The observed inverse kinetic solvent isotope effect for the reaction of PhXn(+) (k(obs)(i-PrOH)/k(obs)(i-PrOD) = 0.39) is consistent with the intermolecular hydride-transfer mechanism. The dramatic reduction of the reaction rate for MA(+), when the water or i-PrOH cosolvent was replaced by AN, suggests that the hydride-transfer T.S. is stabilized by H-bonding between O of the solvent OH and the substrate alcohol OH(delta+). This result suggests an H-bonding stabilization effect on the T.S. of the alcohol dehydrogenase reactions.     

Effects of the donor-acceptor distance and dynamics on hydride tunneling in the dihydrofolate reductase catalyzed reaction

A significant contemporary question in enzymology involves the role of protein dynamics and hydrogen tunneling in enhancing enzyme catalyzed reactions. Here, we report a correlation between the donor-acceptor distance (DAD) distribution and intrinsic kinetic isotope effects (KIEs) for the dihydrofolate reductase (DHFR) catalyzed reaction. This study compares the nature of the hydride-transfer step for a series of active-site mutants, where the size of a side chain that modulates the DAD (I14 in E. coli DHFR) is systematically reduced (I14V, I14A, and I14G). The contributions of the DAD and its dynamics to the hydride-transfer step were examined by the temperature dependence of intrinsic KIEs, hydride-transfer rates, activation parameters, and classical molecular dynamics (MD) simulations. Results are interpreted within the framework of the Marcus-like model where the increase in the temperature dependence of KIEs arises as a direct consequence of the deviation of the DAD from its distribution in the wild type enzyme. Classical MD simulations suggest new populations with larger average DADs, as well as broader distributions, and a reduction in the population of the reactive conformers correlated with the decrease in the size of the hydrophobic residue. The more flexible active site in the mutants required more substantial thermally activated motions for effective H-tunneling, consistent with the hypothesis that the role of the hydrophobic side chain of I14 is to restrict the distribution and dynamics of the DAD and thus assist the hydride-transfer. These studies establish relationships between the distribution of DADs, the hydride-transfer rates, and the DAD's rearrangement toward tunneling-ready states. This structure-function correlation shall assist in the interpretation of the temperature dependence of KIEs caused by mutants far from the active site in this and other enzymes, and may apply generally to C-H→C transfer reactions.     

A computational study of the intramolecular deprotonation of a carbon acid in aqueous solution

Proton transfer reactions are the rate-limiting steps in many biological and synthetic chemical processes, often requiring complex cofactors or catalysts to overcome the generally unfavourable thermodynamic process of carbanion intermediate formation. It has been suggested that quantum tunnelling processes enhance the kinetics of some of these reactions, which when coupled to protein motions may be an important consideration for enzyme catalysis. To obtain a better fundamental and quantitative understanding of these proton transfer mechanisms, a computational analysis of the intramolecular proton transfer from a carbon acid in the small molecule, 4-nitropentanoic acid, in aqueous solution is presented. Potential-energy surfaces from gas-phase, implicit and QM/MM (quantum mechanical/molecular mechanical) explicit solvation quantum chemistry models are compared, and the potential of mean force, for the full reaction coordinate, using umbrella-sampling molecular dynamics is analysed. Semi-classical multidimensional tunnelling corrections are also used to estimate the quantum tunnelling contributions and to understand the origin of the primary deuterium kinetic isotope effects (KIEs). The computational results are found to be in excellent agreement with the KIEs and the energetics obtained experimentally.     

Mechanism of oxidation of diazepam by Chloramine-B: A kinetic approach

The kinetics of the oxidation of diazepam (DZ) by Chloramine-B (CAB) has been studied in aqueous hydrochloric acid medium. The oxidation reaction follows the rate law: The dependence of the reaction rate on temperature is studied and activation parameters for the rate-determining step are evaluated. The dielectric constant of the medium has a small effect on the rate. Ionic strength and the reaction product benzenesulfonamide have no effect on the reaction rate. The solvent isotope effect is studied. A probable mechanism for the observed kinetic data is proposed. © 1998 John Wiley & Sons, Inc. Int J Chem Kinet: 30: 605–611, 1998     

Examining the importance of dynamics,barrier compression and hydrogen tunnelling in enzyme catalysed reactions

Nuclear quantum mechanical tunnelling is important in enzyme-catalysed H-transfer reactions. This viewpoint has arisen after a number of experimental studies have described enzymatic reactions with kinetic isotope effects that are significantly larger than the semiclassical limit. Other experimental evidence for tunnelling, and the potential role of promoting vibrations that transiently compress the reaction barrier, is more indirect, being derived from the interpretation of e.g. mutational analyses of enzyme systems and temperature perturbation studies of reaction rates/kinetic isotope effects. Computational simulations have, in some cases, determined exalted kinetic isotope effects and tunnelling contributions, and identified putative promoting vibrations. In this review, we present the available evidence – both experimental and computational – for environmentally-coupled Htunnelling in several enzyme systems, namely aromatic amine dehydrogenase and members of the Old Yellow Enzyme family. We then consider the relative importance of tunnelling contributions to these reactions. We find that the tunnelling contribution to these reactions confers a rate enhancement of ～1000-fold. Without tunnelling, a 1000-fold reduction in activity would seriously impair cellular metabolism. We therefore infer that tunnelling is crucial to host organism viability thereby emphasising the general importance of tunnelling in biology.     

Quantum dynamics of hydride transfer catalyzed by bimetallic electrophilic catalysis: synchronous motion of Mg(2+) and H(-) in xylose isomerase

Xylose isomerase exhibits a bridged-bimetallic active-site motif in which the substrate is bound to two metals connected by a glutamate bridge, and X-ray crystallographic studies suggest that metal movement is involved in the hydride transfer rate-controlling catalytic step. Here we report classical/quantal dynamical simulations of this step that provide new insight into the metal motion. The potential energy surface is calculated by treating xylose with semiempirical molecular orbital theory augmented by a simple valence bond potential and the rest of the system by molecular mechanics. The rate constant for the hydride-transfer step was calculated by ensemble-averaged dynamical simulations including both variational transition-state theory for determination of the statistically averaged dynamical bottleneck and optimized multidimensional tunneling calculations. The dynamics calculations include 25 317 atoms, with quantized vibrational free energy in 89 active-site degrees of freedom, and with 32 atoms moving through static secondary zone transition-state configurations in the quantum tunneling simulation. Our simulations show that the average Mg-Mg distance R increases monotonically as a function of the hydride-transfer progress variable z. The range of the average R along the reaction path is consistent with the X-ray structure, thus providing a dynamical demonstration of the postulated role of Mg in catalysis. We also predicted the primary deuterium kinetic isotope effect (KIE) for the chemical step. We calculated a KIE of 3.8 for xylose at 298 K, which is consistent with somewhat smaller experimentally observed KIEs for glucose substrate at higher temperatures. More than half of our KIE is due to tunneling; neglecting quantum effects on the reaction coordinate reduces the calculated KIE to 1.8.     

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.     

Proton-coupled electron transfer in soybean lipoxygenase: dynamical behavior and temperature dependence of kinetic isotope effects

The dynamical behavior and the temperature dependence of the kinetic isotope effects (KIEs) are examined for the proton-coupled electron transfer reaction catalyzed by the enzyme soybean lipoxygenase. The calculations are based on a vibronically nonadiabatic formulation that includes the quantum mechanical effects of the active electrons and the transferring proton, as well as the motions of all atoms in the complete solvated enzyme system. The rate constant is represented by the time integral of a probability flux correlation function that depends on the vibronic coupling and on time correlation functions of the energy gap and the proton donor-acceptor mode, which can be calculated from classical molecular dynamics simulations of the entire system. The dynamical behavior of the probability flux correlation function is dominated by the equilibrium protein and solvent motions and is not significantly influenced by the proton donor-acceptor motion. The magnitude of the overall rate is strongly influenced by the proton donor-acceptor frequency, the vibronic coupling, and the protein/solvent reorganization energy. The calculations reproduce the experimentally observed magnitude and temperature dependence of the KIE for the soybean lipoxygenase reaction without fitting any parameters directly to the experimental kinetic data. The temperature dependence of the KIE is determined predominantly by the proton donor-acceptor frequency and the distance dependence of the vibronic couplings for hydrogen and deuterium. The ratio of the overlaps of the hydrogen and deuterium vibrational wavefunctions strongly impacts the magnitude of the KIE but does not significantly influence its temperature dependence. For this enzyme reaction, the large magnitude of the KIE arises mainly from the dominance of tunneling between the ground vibronic states and the relatively large ratio of the overlaps between the corresponding hydrogen and deuterium vibrational wavefunctions. The weak temperature dependence of the KIE is due in part to the dominance of the local component of the proton donor-acceptor motion.     

Kinetics and mechanism of the oxidation of some vicinal and non-vicinal diols by tetrabutylammonium tribromide

Kinetics of oxidation of five vicinal and four non-vicinal diols, and two of their monoethers, by tetrabutylammonium tribromide (TBATB) has been studied. The vicinal diols yield products arising out of glycol-bond fission, while the non-vicinal diols produce the hydroxycarbonyl compounds. The reaction is first-order with respect to TBATB. Michaelis-Menten type kinetics is observed with respect to diols. The reaction fails to induce the polymerization of acrylonitrile. There is no effect of tetrabutylammonium chloride on the reaction rate. The proposed reactive oxidizing species is the tribromide ion. The effect of solvent composition indicates that the rate increases with increase in the polarity of the solvent. The oxidation of [1,1,2,2-²H₄] ethanediol shows the absence of any primary kinetic isotope effect. Values of solvent isotope effect, k(H₂O)/k(D₂O), at 288 K for the oxidation of ethanediol, propane-1,3-diol and 3-methoxybutan-1-ol are 3.41, 0.98 and 1.02 respectively. A mechanism involving a glycol-bond fission has been proposed for the oxidation of vicinal diols. Non-vicinal diols are oxidised by a hydride-transfer mechanism, as they are monohydric alcohols.     

One-step versus stepwise mechanism in protonated amino acid-promoted electron-transfer reduction of a quinone by electron donors and two-electron reduction by a dihydronicotinamide adenine dinucleotide analogue. Interplay between electron transfer and hydrogen bonding

Semiquinone radical anion of 1-(p-tolylsulfinyl)-2,5-benzoquinone (TolSQ(*-)) forms a strong hydrogen bond with protonated histidine (TolSQ(*-)/His x 2 H(+)), which was successfully detected by electron spin resonance. Strong hydrogen bonding between TolSQ(*-) and His x 2 H(+) results in acceleration of electron transfer (ET) from ferrocenes [R2Fc, R = C5H5, C5H4(n-Bu), C5H4Me] to TolSQ, when the one-electron reduction potential of TolSQ is largely shifted to the positive direction in the presence of His x 2 H(+). The rates of His x 2 H(+)-promoted ET from R2Fc to TolSQ exhibit deuterium kinetic isotope effects due to partial dissociation of the N-H bond in His x 2 H(+) at the transition state, when His x 2 H(+) is replaced by the deuterated compound (His x 2 D(+)-d6). The observed deuterium kinetic isotope effect (kH/kD) decreases continuously with increasing the driving force of ET to approach kH/kD = 1.0. On the other hand, His x 2 H(+) also promotes a hydride reduction of TolSQ by an NADH analogue, 9,10-dihydro-10-methylacridine (AcrH2). The hydride reduction proceeds via the one-step hydride-transfer pathway. In such a case, a large deuterium kinetic isotope effect is observed in the rate of the hydride transfer, when AcrH2 is replaced by the dideuterated compound (AcrD2). In sharp contrast to this, no deuterium kinetic isotope effect is observed, when His x 2 H(+) is replaced by His x 2 D(+)-d6. On the other hand, direct protonation of TolSQ and 9,10-phenanthrenequinone (PQ) also results in efficient reductions of TolSQH(+) and PQH(+) by AcrH2, respectively. In this case, however, the hydride-transfer reactions occur via the ET pathway, that is, ET from AcrH2 to TolSQH(+) and PQH(+) occurs in preference to direct hydride transfer from AcrH2 to TolSQH(+) and PQH(+), respectively. The AcrH2(*+) produced by the ET oxidation of AcrH2 by TolSQH(+) and PQH(+) was directly detected by using a stopped-flow technique.     

Primary semiclassical kinetic hydrogen isotope effects in identity carbon-to-carbon proton- and hydride-transfer reactions, an ab initio and DFT computational study

Enthalpies of activation, transition state (ts) geometries, and primary semiclassical (without tunneling) kinetic isotope effects (KIEs) have been calculated for eleven bimolecular identity proton-transfer reactions, four intramolecular proton transfers, four nonidentity proton-transfer reactions, eleven identity hydride transfers, and two 1,2-intramolecular hydride shifts at the HF/6-311+G, MP2/6-311+G, and B3LYP/6-311++G levels. We find the KIEs to be systematically smaller for hydride transfers than for proton transfers. This outcome is not the result of "bent" transition states, although extreme bending can lower the KIE. Rather, it is a consequence of generally greater total bonding in a hydride-transfer ts than in a proton-transfer ts, most prominently manifested as a reduced contribution from the zero-point vibrational energy difference between reactant and transition states (the DeltaZPVE factor) for hydride transfers relative to proton transfers. This and other differences between proton and hydride transfers are rationalized by modeling the central .C...H...C unit of a proton-transfer ts as a 4-electron, 3-center (4-e 3-c) system and the same unit of a hydride-transfer ts as a 2-e 3-c system. Inclusion of tunneling is most likely to magnify the observed differences between proton-transfer and hydride-transfer KIEs, leaving our qualitative conclusions unchanged.     

A mechanistic probe for asymmetric reactions: deuterium isotope effects at enantiotopic groups

The development of a mechanistic probe that is especially suitable for the study of asymmetric reactions is presented. Chemically innocuous enantiotopic methyl groups are utilized as probes for the distinct environments that develop at the transition state for the (-)-B-chlorodiisopinocampheylborane reduction of 4'-methylisobutyrophenone. 2H kinetic isotope effects (KIEs) are determined for both enantiotopic methyl groups using two types of competition reactions. One competition is that between the d3-methyl enantiomeric isotopomers. The other competition reaction is that between the d6-dimethyl and perprotiated isotopologues. The rate constant ratios can be converted into kinetic isotope effects upon each of the individual enantiotopic methyl groups by invoking the rule of the geometric mean. The resulting isotope effect measurements yield highly precise values and contribute further understanding to the transition structure for this stereoselective reduction. The results are discussed in the context of steric isotope effects and the origins of these effects, which arise from the impact of steric crowding upon the anharmonicity of C-H bonds in the transition structure relative to the reactant state.     

Kinetics and mechanism of oxidation of chloramphenicol by 1-chlorobenzotriazole in acidic medium

Chloramphenicol (CAP) is an antibiotic drug having a wide spectrum of activity. The kinetics of oxidation of chloramphenicol by 1-chlorobenzotriazole (CBT) in HClO₄ medium over the temperature range 293–323 K has been investigated. The reaction exhibits first-order kinetics with respect to [CBT]_o and zero-order with respect to [CAP]_o. The fractional-order dependence of rate on [H⁺] suggests complex formation between CBT and H⁺. It fails to induce polymerization of acrylonitrile under the experimental conditions employed. Activation parameters are evaluated. The observed solvent isotope effect indicates the absence of hydride transfer during oxidation. Effects of dielectric constant and ionic strength of the medium on the reaction rate have been studied. Oxidation products are identified. A suitable reaction scheme is proposed and an appropriate rate law is deduced to account for the observed kinetic data.     

Large kinetic isotope effects with unsymmetrical transition states

The results of a series of calculations for a wide, unrestricted variation in the force constants for the making and breaking bonds and their interaction constant are presented for the abstraction reactions of CH₂D₂ with Cl atoms. A wide range of asymmetrical force constants leads to a high kinetic isotope effect as has been pointed out by others for a more restricted range of consideration. These results pointedly contradict the assumed connection between a high kinetic isotope effect and a symmetrical transition state. It is found by examining the atomic displacements of the normal mode that the motion of the H or D atom in the real stretch of the transition state will often have little influence on the isotope effect because the mode can be dominated by end group motions. It is further found that a 3-center model multiplied by a constant factor to account for the contributions of the other vibrations is capable of very satisfactorily reproducing the more rigorous 6-center calculations.     

Micellar effect on hydrolysis of 4-methyl-2-nitroaniline phosphate

The effect of anionic surfactant sodium dodecyl sulfate (SDS) on the hydrolysis of a substrate (mono-4-methyl-2-nitroaniline phosphate) by HCl was studied at 303 K. The reaction followed the first-order kinetics with respect to both HCl and the substrate. SDS effectively catalyzes this reaction, which rate increases with the concentration of SDS due to an increase of dielectric constant of the medium. The kinetic data were fitted to Menger-Portnoy, Piszkiewicz and Berezin kinetic models to explain the observed micellar effects. The various activation parameters both in the presence and absence of SDS were evaluated; a reasonable mechanism was proposed. The rate constant in micellar phase, binding constant and index of cooperativity were calculated accordingly.     

Preorganization and protein dynamics in enzyme catalysis

Recently, an alternative has been offered to the concept of transition state (TS) stabilization as an explanation for rate enhancements in enzyme-catalyzed reactions. Instead, most of the rate increase has been ascribed to preorganization of the enzyme active site to bind substrates in a geometry close to that of the TS, which then transit the activation barrier impelled by motions along the reaction coordinate. The question as to how an enzyme achieves such preorganization and concomitant TS stabilization as well as potential coupled motions along the reaction coordinate leads directly to the role of protein dynamic motion. Dihydrofolate reductase (DHFR) is a paradigm in which the role of dynamics in catalysis continues to be unraveled by a wealth of kinetic, structural, and computational studies. DHFR has flexible loop regions adjacent to the active site whose motions modulate passage through the kinetically preferred pathway. The participation of residues distant from the DHFR active site in enhancing the rate of hydride transfer, however, is unanticipated and may signify the importance of long range protein motions. The general significance of protein dynamics in understanding other biological processes is briefly discussed.     

Conformational substates modulate hydride transfer in dihydrofolate reductase

In earlier studies of the hydride-transfer reaction catalyzed by dihydrofolate reductase (DHFR) we identified features of the protein correlated with variations in the reaction barrier. We extend the scope of those studies by carrying out potential of mean force (PMF) simulations to determine the hydride-transfer barrier in the wild-type protein as well as the G121V and G121S mutants. While our prior studies focused on the reactant state, our current work addresses the full reaction pathway and directly probes the reactive event. The free energy barriers and structural ensembles resulting from these PMF calculations exhibit the same trends reported in our previous work. Fluctuations present in these simulations also exhibit trends associated with differences in the hydride-transfer barrier height. Moreover, vibrational modes anticipated to promote hydride transfer exhibit larger amplitudes in simulations that generate lowered barriers. The results of our study indicate that discrete basins (substates) on a potential energy landscape of the enzyme give rise to distinct hydride-transfer barriers. We suggest that the long-range effects of mutations at position 121 within DHFR are mediated by differentially preorganized protein environments in the context of distinct substate distributions, with concomitant changes to the dynamic properties of the enzyme.