Update 1 of: Tunneling and dynamics in enzymatic hydride transfer

Hydrogen Tunneling Exhibiting Unexpectedly Small Primary Kinetic Isotope Effects

Probing quantum mechanical tunneling (QMT) in chemical reactions is crucial to understanding and developing new transformations. Primary H/D kinetic isotopic effects (KIEs) beyond the semiclassical maximum values of 7-10 (room temperature) are commonly used to assess substantial QMT contributions in one-step hydrogen transfer reactions, because of the much greater QMT probability of protium vs. deuterium. Nevertheless, we report here the discovery of a reaction model occurring exclusively by H-atom QMT with residual primary H/D KIEs. 2-Hydroxyphenylnitrene, generated in N2 matrix, was found to isomerize to an imino-ketone via sequential (domino) QMT involving anti to syn OH-rotamerization (rate determining step) and [1,4]-H shift reactions. These sequential QMT transformations were also observed in the OD-deuterated sample, and unexpected primary H/D KIEs between 3 and 4 were measured at 3 to 20 K. Analogous residual primary H/D KIEs were found in the anti to syn OH-rotamerization QMT of 2-cyanophenol in a N2 matrix. Evidence strongly indicates that these intriguing isotope-insensitive QMT reactivities arise due to the solvation effects of the N2 matrix medium, putatively through coupling with the moving H/D tunneling particle. Should a similar scenario be extrapolated to conventional solution conditions, then QMT may have been overlooked in many chemical reactions.

Study of the Effects of Remote Heavy Group Vibrations on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions

It has recently been observed that the temperature(T)-dependence of KIEs in H-tunneling reactions, characterized by isotopic activation energy difference (ΔEa = EaD - EaH), is correlated to the rigidity of the tunneling ready states (TRSs) in enzymes. A more rigid system with narrowly distributed H-donor-acceptor distances (DADs) at the TRSs gives rise to a weaker T-dependence of KIEs (i.e., a smaller ΔEa). Theoreticians have attempted to develop new H-tunneling models to explain this, but none has been universally accepted. In order to further understand the observations in enzymes and provide useful data to build new theoretical models, we have studied the electronic and solvent effects on ΔEa's for the hydride-tunneling reactions of NADH/NAD+ analogues. We found that a tighter charge-transfer (CT) complex system gives rises to a smaller ΔEa, consistent with the enzyme observations. In this paper, we use the remote heavy group (R) vibrational effects to mediate the system rigidity to study the rigidity-ΔEa relationship. The specific hypothesis is that slower vibrations of a heavier remote group would broaden the DAD distributions and increase the ΔEa value. Four NADH/NAD+ systems were studied in acetonitrile but most of such heavy group vibrations do not appear to significantly increase the ΔEa. The remote heavy group vibrations in these systems may have not affected the CT complexation rigidity to a degree that can significantly increase the DADs, and further, the ΔEa values.

Structural Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions in Acetonitrile: Charge-Transfer Complex Tightness Is a Key

It has recently frequently been found that the kinetic isotope effect (KIE) is independent of temperature (T) in H-tunneling reactions in enzymes but becomes dependent on T in their mutants. Many enzymologists found that the trend is related to different donor-acceptor distances (DADs) at tunneling-ready states (TRSs), which could be sampled by protein dynamics. That is, a more rigid system of densely populated short DADs gives rise to a weaker T dependence of KIEs. Theoreticians have attempted to develop H-tunneling theories to explain the observations, but none have been universally accepted. It is reasonable to assume that the DAD sampling concept, if it exists, applies to the H-transfer reactions in solution, as well. In this work, we designed NADH/NAD+ model reactions to investigate their structural effects on the T dependence of hydride KIEs in acetonitrile. Hammett correlations together with N-CH3/CD3 secondary KIEs were used to provide the electronic structure of the TRSs and thus the rigidity of their charge-transfer complexation vibrations. In all three pairs of reactions, a weaker T dependence of KIEs always corresponds to a steeper Hammett slope on the substituted hydride acceptors. It was found that a tighter/rigid charge-transfer complexation system corresponds with a weaker T dependence of KIEs, consistent with the observations in enzymes.

Correlation of temperature dependence of hydride kinetic isotope effects with donor-acceptor distances in two solvents of different polarities

Recently observed nearly temperature (T)-independent kinetic isotope effects (KIEs) in wild-type enzymes and T-dependent KIEs in variants were used to suggest that H-tunneling in enzymes is assisted by the fast protein vibrations that help sample short donor-acceptor distances (DADs). This supports the recently proposed role of protein vibrations in DAD sampling catalysis. However, use of T-dependence of KIEs to suggest DAD sampling associated with protein vibrations is debated. We have formulated a hypothesis regarding the correlation and designed experiments in solution to investigate it. The hypothesis is, a more rigid system with shorter DADTRS's at the tunneling ready states (TRSs) gives rise to a weaker T-dependence of KIEs, i.e., a smaller ΔEa (= EaD - EaH). In a former work, the solvent effects of acetonitrile versus chloroform on the ΔEa of NADH/NAD+ model reactions were determined, and the DADPRC's of the productive reactant complexes (PRCs) were computed to substitute the DADTRS for the DADTRS-ΔEa correlation study. A smaller ΔEa was found in the more polar acetonitrile where the positively charged PRC is better solvated and has a shorter DADPRC, indirectly supporting the hypothesis. In this work, the TRS structures of different DADTRS's for the hydride tunneling reaction from 1,3-dimethyl-2-phenylimidazoline to 10-methylacridinium were computed. The N-CH3/CD3 secondary KIEs on both reactants were calculated and fitted to the observed values to find the DADTRS order in both solutions. It was found that the equilibrium DADTRS is shorter in acetonitrile than in chloroform. Results directly support the DADTRS-ΔEa correlation hypothesis as well as the explanation that links T-dependence of KIEs to DAD sampling catalysis in enzymes.

Solvent Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects: Correlation to the Donor-Acceptor Distances

Protein structural effects on the temperature (T) dependence of kinetic isotope effects (KIEs) in H-tunneling reactions have recently been used to discuss about the role of enzyme thermal motions in catalysis. Frequently observed nearly T-independent KIEs in the wild-type enzymes and T-dependent KIEs in variants suggest that H-tunneling in the former is assisted by the naturally evolved protein constructive vibrations that help sample short donor-acceptor distances (DADs) needed. This explanation that correlates the T-dependence of KIEs with DAD sampling has been highly debated as simulations following other H-tunneling models sometimes gave alternative explanations. In this paper, solvent effects on the T-dependence of KIEs of two hydride tunneling reactions of NADH/NAD⁺ analogues (represented by ΔE_a = E_aD - E_aH) were determined in attempts to replicate the observations in enzymes and test the protein vibration-assisted DAD sampling concept. Effects of selected aprotic solvents on the DAD_PRC's of the productive reactant complexes (PRCs) and the DAD_TRS's of the activated tunneling ready states (TRSs) were obtained through computations and analyses of the kinetic data, including 2° KIEs, respectively. A weaker T-dependence of KIEs (i.e., smaller ΔE_a) was found in a more polar aprotic solvent in which the system has a shorter average DAD_PRC and DAD_TRS. Further results show that a charge-transfer (CT) complexation made of a stronger donor/acceptor gives rise to a smaller ΔE_a. Overall, the shorter and less broadly distributed DADs resulting from the stronger CT complexation vibrations give rise to a smaller ΔE_a. Our results appear to support the explanation that links the T-dependence of KIEs to the donor-acceptor rigidity in enzymes.

Dependence of crystallographic atomic displacement parameters on temperature (25-150 K) for complexes of horse liver alcohol dehydrogenase

Enzymes catalyze reactions by binding and orienting substrates with dynamic interactions. Horse liver alcohol dehydrogenase catalyzes hydrogen transfer with quantum-mechanical tunneling that involves fast motions in the active site. The structures and B factors of ternary complexes of the enzyme with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol or NAD+ and 2,2,2-trifluoroethanol were determined to 1.1-1.3 Å resolution below the `glassy transition' in order to extract information about the temperature-dependent harmonic motions, which are reflected in the crystallographic B factors. The refinement statistics and structures are essentially the same for each structure at all temperatures. The B factors were corrected for a small amount of radiation decay. The overall B factors for the complexes are similar (13-16 Å2) over the range 25-100 K, but increase somewhat at 150 K. Applying TLS refinement to remove the contribution of pseudo-rigid-body displacements of coenzyme binding and catalytic domains provided residual B factors of 7-10 Å2 for the overall complexes and of 5-10 Å2 for C4N of NAD+ and the methylene carbon of the alcohols. These residual B factors have a very small dependence on temperature and include local harmonic motions and apparently contributions from other sources. Structures at 100 K show complexes that are poised for hydrogen transfer, which involves atomic displacements of ∼0.3 Å and is compatible with the motions estimated from the residual B factors and molecular-dynamics simulations. At 298 K local conformational changes are also involved in catalysis, as enzymes with substitutions of amino acids in the substrate-binding site have similar positions of NAD+ and pentafluorobenzyl alcohol and similar residual B factors, but differ by tenfold in the rate constants for hydride transfer.

Allosteric Communication in the Multifunctional and Redox NQO1 Protein Studied by Cavity-Making Mutations

Allosterism is a common phenomenon in protein biochemistry that allows rapid regulation of protein stability; dynamics and function. However, the mechanisms by which allosterism occurs (by mutations or post-translational modifications (PTMs)) may be complex, particularly due to long-range propagation of the perturbation across protein structures. In this work, we have investigated allosteric communication in the multifunctional, cancer-related and antioxidant protein NQO1 by mutating several fully buried leucine residues (L7, L10 and L30) to smaller residues (V, A and G) at sites in the N-terminal domain. In almost all cases, mutated residues were not close to the FAD or the active site. Mutations L→G strongly compromised conformational stability and solubility, and L30A and L30V also notably decreased solubility. The mutation L10A, closer to the FAD binding site, severely decreased FAD binding affinity (≈20 fold vs. WT) through long-range and context-dependent effects. Using a combination of experimental and computational analyses, we show that most of the effects are found in the apo state of the protein, in contrast to other common polymorphisms and PTMs previously characterized in NQO1. The integrated study presented here is a first step towards a detailed structural–functional mapping of the mutational landscape of NQO1, a multifunctional and redox signaling protein of high biomedical relevance.

Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate

Transient oxidation–reduction through hydride transfer with tightly bound NAD coenzyme is used by a large class of sugar nucleotide epimerases to promote configurational inversion of carbon stereocenters in carbohydrate substrates. A requirement for the epimerases to coordinate hydride abstraction and re-addition with substrate rotation in the binding pocket poses a challenge for dynamical protein conformational selection linked to enzyme catalysis. Here, we studied the thermophilic C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose substrate (k_cat ≈ 1.0 min^–1; 60 °C) to explore the sensitivity of the enzymatic hydride transfer toward environmental fluctuations affected by temperature (20–80 °C). We determined noncompetitive primary kinetic isotope effects (KIE) due to ²H at the glucose C2 and showed that a normal KIE on the k_cat (^Dk_cat) reflects isotope sensitivity of the hydrogen abstraction to enzyme-NAD⁺ in a rate-limiting transient oxidation. The ^Dk_cat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C) and 3.3 at high temperature (80 °C). The temperature profiles for k_cat with the ¹H and ²H substrate showed a decrease in the rate below a dynamically important breakpoint (∼40 °C), suggesting an equilibrium shift to an impaired conformational landscape relevant for catalysis in the low-temperature region. Full Marcus-like model fits of the rate and KIE profiles provided evidence for a high-temperature reaction via low-frequency conformational sampling associated with a broad distribution of hydride donor–acceptor distances (long-distance population centered at 3.31 ± 0.02 Å), only poorly suitable for quantum mechanical tunneling. Collectively, dynamical characteristics of TaCPa2E-catalyzed hydride transfer during transient oxidation of CDP-glucose reveal important analogies to mechanistically simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase. A loose-fit substrate (in TaCPa2E) resembles structural variants of these enzymes by extensive dynamical sampling to balance conformational flexibility and catalytic efficiency.

Pairing structural reconstruction with catalytic competence to evaluate the mechanisms of key enzymes in the folate-mediated one-carbon pathway

Mammalian metabolism comprises a series of interlinking pathways that include two major cycles: the folate and methionine cycles. The folate-mediated metabolic cycle uses several oxidation states of tetrahydrofolate to carry activated one-carbon units to be readily used and interconverted within the cell. They are required for nucleotide synthesis, methylation and metabolism, and particularly for proliferation of cancer cells. Based on the latest progress in genome-wide CRISPR loss-of-function viability screening of 789 cell lines, we focus on the most cancer-dependent enzymes in this pathway, especially those that are hyperactivated in cancer, to provide new insight into the chemical basis for cancer drug development. Since the complete 3D structure of several of these enzymes of the one-carbon pathway in their active form are not available, we used homology modelling integrated with the interpretation of the reaction mechanism. In addition, have reconstructed the most likely scenario for the reactions taking place paired with their catalytic competence that provides a testable framework for this pathway.

Thermodynamics of the elementary steps of organic hydride chemistry determined in acetonitrile and their applications

This review focuses on the thermodynamics of the elementary step of 421 organic hydrides and unsaturated compounds releasing or accepting hydride or hydrogen determined in acetonitrile as well as their potential applications.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.

Large Isotope Effects in Organometallic Chemistry

The kinetic isotope effect (KIE) is key to understanding reaction mechanisms in many areas of chemistry and chemical biology, including organometallic chemistry. This ratio of rate constants, k_H /k_D , typically falls between 1-7. However, KIEs up to 105 have been reported, and can even be so large that reactivity with deuterium is unobserved. We collect here examples of large KIEs across organometallic chemistry, in catalytic and stoichiometric reactions, along with their mechanistic interpretations. Large KIEs occur in proton transfer reactions such as protonation of organometallic complexes and clusters, protonolysis of metal-carbon bonds, and dihydrogen reactivity. C-H activation reactions with large KIEs occur with late and early transition metals, photogenerated intermediates, and abstraction by metal-oxo complexes. We categorize the mechanistic interpretations of large KIEs into the following three types: (a) proton tunneling, (b) compound effects from multiple steps, and (c) semi-classical effects on a single step. This comprehensive collection of large KIEs in organometallics provides context for future mechanistic interpretation.

Structural basis of the pleiotropic and specific phenotypic consequences of missense mutations in the multifunctional NAD(P)H:quinone oxidoreductase 1 and their pharmacological rescue

The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer-associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by long-range structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models.

Direct Correlation between Donor-Acceptor Distance and Temperature Dependence of Kinetic Isotope Effects in Hydride-Tunneling Reactions of NADH/NAD+ Analogues

Recent study of structural effects on primary kinetic isotope effects (1° KIEs) of H-transfer reactions in enzymes and solution revealed that a more rigid reaction system gave rise to a weaker temperature dependence of 1° KIEs, i.e., a smaller isotopic activation energy difference (ΔE_a = E_aD - E_aH). This has been explained within the contemporary vibrationally assisted activated H-tunneling (VA-AHT) model in which rigidity is defined according to the density of donor-acceptor distance (DAD_TRS) populations at the tunneling ready state (TRS) sampled by heavy atom motions. To test the relationship between DAD_TRS and ΔE_a in the model, we developed a computational method to obtain the TRS structures for H-transfer reactions. The method was applied to three hydride transfer reactions of NADH/NAD⁺ analogues for which the ΔE_a's as well as secondary (2°) KIEs have been reported. The 2° KIEs computed from each TRS structure were fitted to the observed values to obtain the optimal TRSs/DAD_TRS's. It was found that a shorter DAD_TRS does correspond with a smaller ΔE_a. This appears to support the VA-AHT model. Moreover, an analysis of hybridizations at the bent TRS structures shows that rehybridizations at the donor-acceptor centers are much more advanced than predicted from the classical mechanism. This implies that more orbital preparations are required for the nonclassical H-tunneling to take place.

Protein Mass-Modulated Effects in Alkaline Phosphatase

Recent experimental studies engaging isotopically substituted protein (heavy protein) have revealed that many, but not all, enzymatic systems exhibit altered chemical steps in response to an altered mass. The results have been interpreted as femtosecond protein dynamics at the active site being linked (or not) to transition-state barrier crossing. An altered enzyme mass can influence several kinetic parameters (kcat, Km, and kchem) in amounts of ≤30% relative to light enzymes. An early report on deuterium-labeled Escherichia coli alkaline phosphatase (AP) showed an unusually large enzyme kinetic isotope effect on kcat. We examined steady-state and chemical step properties of native AP, [2H]AP, and [2H,13C,15N]AP to characterize the role of heavy enzyme protein dynamics in reactions catalyzed by AP. Both [2H]- and [2H,13C,15N]APs showed unaltered steady-state and single-turnover rate constants. These findings characterize AP as one of the enzymes in which mass-dependent catalytic site dynamics is dominated by reactant-linked atomic motions. Two catalytic site zinc ions activate the oxygen nucleophiles in the catalytic site of AP. The mass of the zinc ions is unchanged in light and heavy APs. They are essentially linked to catalysis and provide a possible explanation for the loss of linkage between catalysis and protein mass in these enzymes.

The Catalytic Cycle of the Antioxidant and Cancer-Associated Human NQO1 Enzyme: Hydride Transfer, Conformational Dynamics and Functional Cooperativity

Human NQO1 [NAD(H):quinone oxidoreductase 1] is a multi-functional and stress-inducible dimeric protein involved in the antioxidant defense, the activation of cancer prodrugs and the stabilization of oncosuppressors. Despite its roles in human diseases, such as cancer and neurological disorders, a detailed characterization of its enzymatic cycle is still lacking. In this work, we provide a comprehensive analysis of the NQO1 catalytic cycle using rapid mixing techniques, including multiwavelength and spectral deconvolution studies, kinetic modeling and temperature-dependent kinetic isotope effects (KIEs). Our results systematically support the existence of two pathways for hydride transfer throughout the NQO1 catalytic cycle, likely reflecting that the two active sites in the dimer catalyze two-electron reduction with different rates, consistent with the cooperative binding of inhibitors such as dicoumarol. This negative cooperativity in NQO1 redox activity represents a sort of half-of-sites activity. Analysis of KIEs and their temperature dependence also show significantly different contributions from quantum tunneling, structural dynamics and reorganizations to catalysis at the two active sites. Our work will improve our understanding of the effects of cancer-associated single amino acid variants and post-translational modifications in this protein of high relevance in cancer progression and treatment.

Substituent Effects on Temperature Dependence of Kinetic Isotope Effects in Hydride-Transfer Reactions of NADH/NAD+ Analogues in Solution: Reaction Center Rigidity Is the Key

Substituent effects on the temperature dependence of primary kinetic isotope effects, characterized by ΔE_a = E_aD - E_aH, for two series of the title reactions in acetonitrile were studied. The change from ΔE_a ≈ 0 for a highly rigid system to ΔE_a > 0 for systems with reduced rigidities was observed. The rigidities were controlled by the electronic and steric effects. This work replicates the observations in enzymes and opens a new research direction that studies structure-ΔE_a relationship.

Kinetic Isotope Effects and Hydrogen Tunnelling in PCET Oxidations of Ascorbate: New Insights into Aqueous Chemistry?

Recent experimental studies of kinetic isotope effects (KIE-s) and hydrogen tunnelling comprising three proton-coupled electron transfer (PCET) oxidations of ascorbate monoanion, (a) in aqueous reaction solutions, (b) in the mixed water-organic cosolvent systems, (c) in aqueous solutions of various salts and (d) in fairly diluted aqueous solutions of the various partial hydrophobes are reviewed. A number of new insights into the wealth of the kinetic isotope phenomena in the PCET reactions have been obtained. The modulation of KIE-s and hydrogen tunnelling observed when partially hydrophobic solutes are added into water reaction solution, in the case of fairly diluted solutions is revealed as the strong linear correlation of the isotopic ratios of the Arrhenius prefactors Ah/Ad and the isotopic differences in activation energies ΔE_a (D,H). The observation has been proposed to be a signature of the involvement of the collective intermolecular excitonic vibrational dynamics of water in activation processes and aqueous chemistry.

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.

Substitutions of Amino Acid Residues in the Substrate Binding Site of Horse Liver Alcohol Dehydrogenase Have Small Effects on the Structures but Significantly Affect Catalysis of Hydrogen Transfer

Previous studies showed that the L57F and F93W alcohol dehydrogenases catalyze the oxidation of benzyl alcohol with some quantum mechanical hydrogen tunneling, whereas the V203A enzyme has diminished tunneling. Here, steady-state kinetics for the L57F and F93W enzymes were studied, and microscopic rate constants for the ordered bi-bi mechanism were estimated from simulations of transient kinetics for the S48T, F93A, S48T/F93A, F93W, and L57F enzymes. Catalytic efficiencies for benzyl alcohol oxidation (V₁/E_tK_b) vary over a range of ∼100-fold for the less active enzymes up to the L57F enzyme and are mostly associated with the binding of alcohol rather than the rate constants for hydride transfer. In contrast, catalytic efficiencies for benzaldehyde reduction (V₂/E_tK_p) are ∼500-fold higher for the L57F enzyme than for the less active enzymes and are mostly associated with the rate constants for hydride transfer. Atomic-resolution structures (1.1 Å) for the F93W and L57F enzymes complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol are almost identical to previous structures for the wild-type, S48T, and V203A enzymes. Least-squares refinement with SHELXL shows that the nicotinamide ring is slightly strained in all complexes and that the apparent donor-acceptor distances from C4N of NAD to C7 of pentafluorobenzyl alcohol range from 3.28 to 3.49 Å (±0.02 Å) and are not correlated with the rate constants for hydride transfer or hydrogen tunneling. How the substitutions affect the dynamics of reorganization during hydrogen transfer and the extent of tunneling remain to be determined.

Odd Loop Regions of XenA and XenB Enzymes Modulate Their Interaction with Nitro-explosives Compounds and Provide Structural Support for Their Regioselectivity

The nitro-explosive compounds 2,4,6-trinitrotoluene, 2,4,6-trinitrophenol, and 1,2,3-trinitroglycerol are persistent environmental contaminants. The presence of different functional groups in these molecules represents a great challenge to enzymatic catalysis. The chemical variety of these three substrates is such that they do not bind and interact with catalytic residues within an enzyme with the same affinity. In this context, two Xenobiotic Reductase enzymes produced by the bacteria Pseudomonas putida can catalyze the reduction of these compounds with different affinities and regioselectivity. The structural bases that support this substrate promiscuity and catalytic preferences are unknown. Therefore, through molecular dynamics simulations and free energy calculations, we explored the structural properties driving the specific interactions of these enzymes with their substrates and cofactor. Models of Xenobiotic Reductase A and B enzymes in complex with 2,4,6-trinitrotoluene, 2,4,6-trinitrophenol, or 1,2,3-trinitroglycerol were built, and the ligand enzyme interaction was simulated by molecular dynamics. The structural analysis of the molecular dynamics simulations shows that loops 3, 5, 7, 9, 11, and 13 of Xenobiotic Reductase B, and loops 4, 5, 7, 11, 13, and 15 Xenobiotic Reductase A, are in contact with the ligands during the first stages of the molecular recognition. These loops are the most flexible regions for both enzymes; however, Xenobiotic Reductase B presents a greater range of movement and a higher number of residues interacting with the ligands. Finally, the distance between the cofactor and the different reactive groups in the substrate reflects the regioselectivity of the enzymes, and the free energy calculations are consistent with the substrate specificity of both enzymes studied. The simulation shows a stable interaction between the aromatic ring of the substrates and Xenobiotic Reductase B. In contrast, a less stable interaction with the different nitro groups of the aromatic ligands was observed. In the case of 1,2,3-trinitroglycerol, Xenobiotic Reductase B interacts more closely with the nitro groups of carbon 1 or 3, while Xenobiotic Reductase A is more selective by nitro groups of carbon 2. The obtained results suggest that the flexibility of the loops in Xenobiotic Reductase B and the presence of polar and aromatic residues present in loops 5 and 7 are fundamental to determine the affinity of the enzyme with the different substrates, and they also contribute to the proper orientation of the ligands that directs the catalytic reaction.

Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations as well as bioinformatic analyses. A few specific residues (M42, G121, F125, and I14) were shown to function synergistically with measurements of single-turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEsint) of site-directed mutants. This study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEsint and computer simulations based on the empirical valence bond method. We study mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

Replication of the Enzymatic Temperature Dependency of the Primary Hydride Kinetic Isotope Effects in Solution: Caused by the Protein-Controlled Rigidity of the Donor-Acceptor Centers?

The change from the temperature independence of the primary (1°) H/D kinetic isotope effects (KIEs) in wild-type enzyme-catalyzed H-transfer reactions (ΔE_a = E_aD - E_aH ∼ 0) to a strong temperature dependence with the mutated enzymes (ΔE_a ≫ 0) has recently been frequently observed. This has prompted some enzymologists to develop new H-tunneling models to correlate ΔE_a with the donor-acceptor distance (DAD) at the tunneling-ready state (TRS) as well as the protein thermal motions/dynamics that sample the short DAD_TRS's for H-tunneling to occur. While extensive evidence supporting or disproving the thermally activated DAD sampling concept has emerged, a comparable study of the simpler bimolecular H-tunneling reactions in solution has not been carried out. In particular, small ΔE_a's (∼0) have not been found. In this paper, we report a study of the hydride-transfer reactions from four NADH models to the same hydride acceptor in acetonitrile. The ΔE_a's were determined: 0.37 (small), 0.60, 0.99, and 1.53 kcal/mol (large). The α-secondary (2°) KIEs on the acceptor that serve as a ruler for the rigidity of reaction centers were previously reported or determined. All possible productive reactant complex (PRC) configurations were computed to provide insight into the structures of the TRS's. Relationships among structures, 2° KIEs, DAD_PRC's, and ΔE_a's were discussed. The more rigid system with more suppressed 2° C-H vibrations at the TRS and more narrowly distributed DAD_PRC's in PRCs gave a smaller ΔE_a. The results replicated the trend observed in enzymes versus mutated enzymes and appeared to support the concepts of different thermally activated DAD_TRS sampling processes in response to the rigid versus flexible donor-acceptor centers.

Molecular modeling of conformational dynamics and its role in enzyme evolution

With increasing computational power, biomolecular simulations have become an invaluable tool for understanding enzyme mechanisms and the origins of enzyme catalysis. More recently, computational studies have started to focus on understanding how enzyme activity itself evolves, both in terms of enhancing the native or new activities on existing enzyme scaffolds, or completely de novo on previously non-catalytic scaffolds. In this context, both experiment and molecular modeling provided strong evidence for an important role of conformational dynamics in the evolution of enzyme functions. This contribution will present a brief overview of the current state of the art for computationally exploring enzyme conformational dynamics in enzyme evolution, and, using several showcase studies, illustrate the ways molecular modeling can be used to shed light on how enzyme function evolves, at the most fundamental molecular level.

Substitutions of a buried glutamate residue hinder the conformational change in horse liver alcohol dehydrogenase and yield a surprising complex with endogenous 3'-Dephosphocoenzyme A

Glu-267 is highly conserved in alcohol dehydrogenases and buried as a negatively-charged residue in a loop of the NAD coenzyme binding domain. Glu-267 might have a structural role and contribute to a rate-promoting vibration that facilitates catalysis. Substitutions of Glu-267 with histidine or asparagine residues increase the dissociation constants for the coenzymes (NAD⁺ by ∼40-fold, NADH by ∼200-fold) and significantly decrease catalytic efficiencies by 16-1200-fold various substrates and substituted enzymes. The turnover numbers modestly change with the substitutions, but hydride transfer is at least partially rate-limiting for turnover for alcohol oxidation. X-ray structures of the E267H and E267 N enzymes are similar to the apoenzyme (open) conformation of the wild-type enzyme, and the substitutions are accommodated by local changes in the structure. Surprisingly, the E267H and E267 N enzymes have endogenous (from the expression in E. coli) 3'-dephosphocoenzyme A bound in the active site with the ADP moiety in the NAD binding site and the pantethiene sulfhydryl bound to the catalytic zinc. The kinetics and crystallography show that the substitutions of Glu-267 hinder the conformational change, which occurs when wild-type enzyme binds coenzymes, and affect productive binding of substrates.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Bacterial versus human thymidylate synthase: Kinetics and functionality

Thymidylate Synthase (TSase) is a highly conserved enzyme that catalyzes the production of the DNA building block thymidylate. Structurally, functionally and mechanistically, bacterial and mammalian TSases share remarkable similarities. Because of this closeness, bacterial enzymes have long been used as model systems for human TSase. Furthermore, while TSase inhibitors have long served as chemotherapeutic drugs, no TSase inhibitor serves as an antibiotic. Despite their high resemblance, the mammalian TSases are distinct in a few known aspects, such as having a N-terminal tail and two insertions in the primary sequence and active/inactive conformations. Here, we aim to comprehensively characterize human (hs) TSase and delineate its contrasts and the similarities to the well-studied Escherichia coli (ec) TSase. We found that, in contrast to ecTSase, Mg²⁺ does not enhance reaction rates for hsTSase. The temperature dependence of intrinsic kinetic isotope effects (KIEs), on the other hand, suggests that Mg²⁺ has little or no impact on the transition state of hydride transfer in either enzyme, and that the transition state for the hydride transfer in hsTSase is looser than in ecTSase. Additionally, the substrates’ binding order is strictly ordered for ecTSase but slightly less ordered for hsTSase. The observed kinetic and functional differences between bacterial and human enzymes may aid in the development of antibiotic drugs with reduced toxicity.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Protein Mass Effects on Formate Dehydrogenase

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Separating Proton and Electron Transfer Effects in Three-Component Concerted Proton-Coupled Electron Transfer Reactions

Multiple-site concerted proton-electron transfer (MS-CPET) reactions were studied in a three-component system. 1-Hydroxy-2,2,6,6-tetramethylpiperidine (TEMPOH) was oxidized to the stable radical TEMPO by electron transfer to ferrocenium oxidants coupled to proton transfer to various pyridine bases. These MS-CPET reactions contrast with the usual reactivity of TEMPOH by hydrogen atom transfer (HAT) to a single e^-/H⁺ acceptor. The three-component reactions proceed by pre-equilibrium formation of a hydrogen-bonded adduct between TEMPOH and the pyridine base, and the adduct is then oxidized by the ferrocenium in a bimolecular MS-CPET step. The second-order rate constants, measured using stopped-flow kinetic techniques, spanned 4 orders of magnitude. An advantage of this system is that the MS-CPET driving force could be independently varied by changing either the pK_a of the base or the reduction potential (E°) of the oxidant. Changes in ΔG°_MS-CPET from either source had the same effect on the MS-CPET rate constants, and a combined Brønsted plot of ln(k_MS-CPET) vs ln(K_eq) was linear with a slope of 0.46. These results imply a synchronous concerted mechanism, in which the proton and electron transfer components of the CPET process make equal contributions to the rate constants. The only outliers to the Brønsted correlation are the reactions with sterically hindered pyridines, which apparently hinder the close approach of proton donor and acceptor that facilitates MS-CPET. These three-component reactions are compared with a related HAT reaction of TEMPOH, with the 2,4,6-tri-tert-butylphenoxyl radical. The MS-CPET and HAT oxidations of TEMPOH at the same driving force occurred with similar rate constants. While this is an imperfect comparison, the data suggest that the separation of the proton and electron to different reagents does not significantly inhibit the proton-coupled electron transfer process.

Mechanistic Studies of an Amine Oxidase Derived from d-Amino Acid Oxidase

The flavoprotein d-amino acid oxidase has long served as a paradigm for understanding the mechanism of oxidation of amino acids by flavoproteins. Recently, a mutant d-amino acid oxidase (Y228L/R283G) that catalyzed the oxidation of amines rather than amino acids was described [Yasukawa, K., et al. (2014) Angew. Chem., Int. Ed. 53, 4428-4431]. We describe here the use of pH and kinetic isotope effects with (R)-α-methylbenzylamine as a substrate to determine whether the mutant enzyme utilizes the same catalytic mechanism as the wild-type enzyme. The effects of pH on the steady-state and rapid-reaction kinetics establish that the neutral amine is the substrate, while an active-site residue, likely Tyr224, must be uncharged for productive binding. There is no solvent isotope effect on the k_cat/K_m value for the amine, consistent with the neutral amine being the substrate. The deuterium isotope effect on the k_cat/K_m value is pH-independent, with an average value of 5.3, similar to values found with amino acids as substrates for the wild-type enzyme and establishing that there is no commitment to catalysis with this substrate. The k_cat/K_O2 value is similar to that seen with amino acids as the substrate, consistent with the oxidative half-reaction being unperturbed by the mutation and with flavin oxidation preceding product release. All of the data are consistent with the mutant enzyme utilizing the same mechanism as the wild-type enzyme, transfer of hydride from the neutral amine to the flavin.

Examinations of the Chemical Step in Enzyme Catalysis

Advances in computational and experimental methods in enzymology have aided comprehension of enzyme-catalyzed chemical reactions. The main difficulty in comparing computational findings to rate measurements is that the first examines a single energy barrier, while the second frequently reflects a combination of many microscopic barriers. We present here intrinsic kinetic isotope effects and their temperature dependence as a useful experimental probe of a single chemical step in a complex kinetic cascade. Computational predictions are tested by this method for two model enzymes: dihydrofolate reductase and thymidylate synthase. The description highlights the significance of collaboration between experimentalists and theoreticians to develop a better understanding of enzyme-catalyzed chemical conversions.

Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

A reaction's transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.

Investigating Information Dynamics in Living Systems through the Structure and Function of Enzymes

Enzymes are proteins that accelerate intracellular chemical reactions often by factors of 10⁵−10¹²s⁻¹. We propose the structure and function of enzymes represent the thermodynamic expression of heritable information encoded in DNA with post-translational modifications that reflect intra- and extra-cellular environmental inputs. The 3 dimensional shape of the protein, determined by the genetically-specified amino acid sequence and post translational modifications, permits geometric interactions with substrate molecules traditionally described by the key-lock best fit model. Here we apply Kullback-Leibler (K-L) divergence as metric of this geometric “fit” and the information content of the interactions. When the K-L ‘distance’ between interspersed substrate p_n and enzyme r_n positions is minimized, the information state, reaction probability, and reaction rate are maximized. The latter obeys the Arrhenius equation, which we show can be derived from the geometrical principle of minimum K-L distance. The derivation is first limited to optimum substrate positions for fixed sets of enzyme positions. However, maximally improving the key/lock fit, called ‘induced fit,’ requires both sets of positions to be varied optimally. We demonstrate this permits and is maximally efficient if the key and lock particles p_n,r_n are quantum entangled because the level of entanglement obeys the same minimized value of the Kullback-Leibler distance that occurs when all p_n ≈ r_n. This implies interchanges p_n ⇄ br_n randomly taking place during a reaction successively improves key/lock fits, reducing the activation energy E_a and increasing the reaction rate k. Our results demonstrate the summation of heritable and environmental information that determines the enzyme spatial configuration, by decreasing the K-L divergence, is converted to thermodynamic work by reducing E_a and increasing k of intracellular reactions. Macroscopically, enzyme information increases the order in living systems, similar to the Maxwell demon gedanken, by selectively accelerating specific reaction thus generating both spatial and temporal concentration gradients.

The general base in the thymidylate synthase catalyzed proton abstraction

The enzyme thymidylate synthase (TSase), an important chemotherapeutic drug target, catalyzes the formation of 2'-deoxythymidine-5'-monophosphate (dTMP), a precursor of one of the DNA building blocks. TSase catalyzes a multi-step mechanism that includes the abstraction of a proton from the C5 of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). Previous studies on ecTSase proposed that an active-site residue, Y94 serves the role of the general base abstracting this proton. However, since Y94 is neither very basic, nor connected to basic residues, nor located close enough to the pyrimidine proton to be abstracted, the actual identity of this base remains enigmatic. Based on crystal structures, an alternative hypothesis is that the nearest potential proton-acceptor of C5 of dUMP is a water molecule that is part of a hydrogen bond (H-bond) network comprised of several water molecules and several protein residues including H147, E58, N177, and Y94. Here, we examine the role of the residue Y94 in the proton abstraction step by removing its hydroxyl group (Y94F mutant). We investigated the effect of the mutation on the temperature dependence of intrinsic kinetic isotope effects (KIEs) and found that these KIEs are more temperature dependent than those of the wild-type enzyme (WT). These results suggest that the phenolic -OH of Y94 is a component of the transition state for the proton abstraction step. The findings further support the hypothesis that no single functional group is the general base, but a network of bases and hydroxyls (from water molecules and tyrosine) sharing H-bonds across the active site can serve the role of the general base to remove the pyrimidine proton.

Activation of Two Sequential H-transfers in the Thymidylate Synthase Catalyzed Reaction

Thymidylate synthase (TSase) catalyzes the de novo biosynthesis of thymidylate, a precursor for DNA, and is thus an important target for chemotherapeutics and antibiotics. Two sequential C-H bond cleavages catalyzed by TSase are of particular interest: a reversible proton abstraction from the 2'-deoxy-uridylate substrate, followed by an irreversible hydride transfer forming the thymidylate product. QM/MM calculations of the former predicted a mechanism where the abstraction of the proton leads to formation of a novel nucleotide-folate intermediate that is not covalently bound to the enzyme (Wang, Z.; Ferrer, S.; Moliner, V.; Kohen, A. Biochemistry2013, 52, 2348-2358). Existence of such intermediate would hold promise as a target for a new class of drugs. Calculations of the subsequent hydride transfer predicted a concerted H-transfer and elimination of the enzymatic cysteine (Kanaan, N.; Ferrer, S.; Marti, S.; Garcia-Viloca, M.; Kohen, A.; Moliner, V. J. Am. Chem. Soc.2011, 133, 6692-6702). A key to both C-H activations is a highly conserved arginine (R166) that stabilizes the transition state of both H-transfers. Here we test these predictions by studying the R166 to lysine mutant of E. coli TSase (R166K) using intrinsic kinetic isotope effects (KIEs) and their temperature dependence to assess effects of the mutation on both chemical steps. The findings confirmed the predictions made by the QM/MM calculations, implicate R166 as an integral component of both reaction coordinates, and thus provide critical support to the nucleotide-folate intermediate as a new target for rational drug design.

The Importance of the MM Environment and the Selection of the QM Method in QM/MM Calculations: Applications to Enzymatic Reactions

In this chapter, we discuss the influence of an anisotropic protein environment on the reaction mechanisms of saccharopine reductase and uroporphyrinogen decarboxylase, respectively, via the use of a quantum mechanical and molecular mechanical (QM/MM) approach. In addition, we discuss the importance of selecting a suitable DFT functional to be used in a QM/MM study of a key intermediate in the mechanism of 8R-lipoxygenase, a nonheme iron enzyme. In the case of saccharopine reductase, while the enzyme utilizes a substrate-assisted catalytic pathway, it was found that only through treating the polarizing effect of the active site, via the use of an electronic embedding formalism, was agreement with experimental kinetic data obtained. Similarly, in the case of uroporphyrinogen decarboxylase, the effect of the protein environment on the catalytic mechanism was found to be such that the calculated rate-limiting barrier is in good agreement with related experimentally determined values for the first decarboxylation of the substrate. For 8R-lipoxygenase, it was found that the geometries and energies of the multicentered open-shell intermediate complexes formed during the mechanism are quite sensitive to the choice of the density functional theory method. Thus, while density functional theory has become the method of choice in QM/MM studies, care must be taken in the selection of a particular high-level method.

Network of remote and local protein dynamics in dihydrofolate reductase catalysis

Molecular dynamics calculations and bionformatic studies of dihydrofolate reductase (DHFR) have suggested a network of coupled motions across the whole protein that is correlated to the reaction coordinate. Experimental studies demonstrated that distal residues G121, M42 and F125 in E. coli DHFR participate in that network. The missing link in our understanding of DHFR catalysis is the lack of a mechanism by which such remote residues can affect the catalyzed chemistry at the active site. Here, we present a study of the temperature dependence of intrinsic kinetic isotope effects (KIEs) that indicates synergism between a remote residue in that dynamic network, G121, and the active site's residue I14. The intrinsic KIEs for the I14A-G121V double mutant showed steeper temperature dependence (ΔE_a(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔE_a(T-H)G121V +ΔE_{a(T-H) I14A} < ΔE_{a(T-H) double mutant}, which indicates a synergism between the two residues. This finding links the remote residues in the network under investigation to the enzyme's active site, providing a mechanism by which these residues can be coupled to the catalyzed chemistry. This experimental evidence validates calculations proposing that both remote and active site residues constitute a network of coupled promoting motions correlated to the bond activation step (C-H→C hydride transfer in this case). Additionally, the effect of I14A and G121V mutations on single turnover rates was additive rather than synergistic. Although single turnover rate measurements are more readily available and thus more popular than assessing intrinsic kinetic isotope effects, the current finding demonstrates that for these rates, which in DHFR reflect several microscopic rate constants, can fall short of revealing the nature of the C-H bond activation per se.

Another Look at the Mechanisms of Hydride Transfer Enzymes with Quantum and Classical Transition Path Sampling

The mechanisms involved in enzymatic hydride transfer have been studied for years, but questions remain due, in part, to the difficulty of probing the effects of protein motion and hydrogen tunneling. In this study, we use transition path sampling (TPS) with normal mode centroid molecular dynamics (CMD) to calculate the barrier to hydride transfer in yeast alcohol dehydrogenase (YADH) and human heart lactate dehydrogenase (LDH). Calculation of the work applied to the hydride allowed for observation of the change in barrier height upon inclusion of quantum dynamics. Similar calculations were performed using deuterium as the transferring particle in order to approximate kinetic isotope effects (KIEs). The change in barrier height in YADH is indicative of a zero-point energy (ZPE) contribution and is evidence that catalysis occurs via a protein compression that mediates a near-barrierless hydride transfer. Calculation of the KIE using the difference in barrier height between the hydride and deuteride agreed well with experimental results.

Role of long-range protein dynamics in different thymidylate synthase catalyzed reactions

Recent studies of Escherichia coli thymidylate synthase (ecTSase) showed that a highly conserved residue, Y209, that is located 8 Å away from the reaction site, plays a key role in the protein’s dynamics. Those crystallographic studies indicated that Y209W mutant is a structurally identical but dynamically altered relative to the wild type (WT) enzyme, and that its turnover catalytic rate governed by a slow hydride-transfer has been affected. The most challenging test of an examination of a fast chemical conversion that precedes the rate-limiting step has been achieved here. The physical nature of both fast and slow C-H bond activations have been compared between the WT and mutant by means of observed and intrinsic kinetic isotope effects (KIEs) and their temperature dependence. The findings indicate that the proton abstraction step has not been altered as much as the hydride transfer step. Additionally, the comparison indicated that other kinetic steps in the TSase catalyzed reaction were substantially affected, including the order of the substrate binding. Enigmatically, although Y209 is H-bonded to 3'-OH of 2'-deoxyuridine-5'-mono­phosphate (dUMP), its altered dynamics is more pronounced on the binding of the remote cofactor, (6R)-N⁵,N¹⁰-methylene-5,6,7,8-tetrahydrofolate (CH₂H₄folate), revealing the importance of long-range dynamics of the enzymatic complex and its catalytic function.