1
|
Eberhart ME, Wilson TR, Johnston NW, Alexandrova AN. Geometry of Charge Density as a Reporter on the Role of the Protein Scaffold in Enzymatic Catalysis: Electrostatic Preorganization and Beyond. J Chem Theory Comput 2023; 19:694-704. [PMID: 36562645 DOI: 10.1021/acs.jctc.2c01060] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Enzymes host active sites inside protein macromolecules, which have diverse, often incredibly complex, and atom-expensive structures. It is an outstanding question what the role of these expensive scaffolds might be in enzymatic catalysis. Answering this question is essential to both enzymology and the design of artificial enzymes with proficiencies that will match those of the best natural enzymes. Protein rigidifying the active site, contrasted with the dynamics and vibrational motion promoting the reaction, as well as long-range electrostatics (also known as electrostatic preorganization) were all proposed as central contributions of the scaffold to the catalysis. Here, we show that all these effects inevitably produce changes in the quantum mechanical electron density in the active site, which in turn defines the reactivity. The phenomena are therefore fundamentally inseparable. The geometry of the electron density-a scalar field characterized by a number of mathematical features such as critical points-is a rigorous and convenient descriptor of enzymatic catalysis and a reporter on the role of the protein. We show how this geometry can be analyzed, linked to the reaction barriers, and report in particular on intramolecular electric fields in enzymes. We illustrate these tools on the studies of electrostatic preorganization in several representative enzyme classes, both natural and artificial. We highlight the forward-looking aspects of the approach.
Collapse
Affiliation(s)
- Mark E Eberhart
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80401, United States
| | - Timothy R Wilson
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80401, United States
| | - Nathaniel W Johnston
- Department of Chemistry, and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095, United States
| | - Anastassia N Alexandrova
- Department of Chemistry, and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095, United States
| |
Collapse
|
2
|
Abstract
This Perspective reviews the use of Transition Path Sampling methods to study enzymatically catalyzed chemical reactions. First applied by our group to an enzymatic reaction over 15 years ago, the method has uncovered basic principles in enzymatic catalysis such as the protein promoting vibration, and it has also helped harmonize such ideas as electrostatic preorganization with dynamic views of enzyme function. It is now being used to help uncover principles of protein design necessary to artificial enzyme creation.
Collapse
Affiliation(s)
- Steven D Schwartz
- Department of Chemistry and Biochemistry University of Arizona Tucson, Arizona 85721, United States
| |
Collapse
|
3
|
Andrews BA, Dyer RB. Comparison of the Role of Protein Dynamics in Catalysis by Dihydrofolate Reductase from E. coli and H. sapiens. J Phys Chem B 2022; 126:7126-7134. [PMID: 36069763 DOI: 10.1021/acs.jpcb.2c05112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dihydrofolate reductase (DHFR) is a well-studied, clinically relevant enzyme known for being highly dynamic over the course of its catalytic cycle. However, the role dynamic motions play in the explicit hydride transfer from the nicotinamide cofactor to the dihydrofolate substrate remains unclear because reaction initiation and direct spectroscopic examination on the appropriate time scale for such femtosecond to picosecond motions is challenging. Here, we employ pre-steady-state kinetics to observe the hydride transfer as directly as possible in two different species of DHFR: Escherichia coli and Homo sapiens. While the hydride transfer has been well-characterized in DHFR from E. coli, improvements in time resolution now allow for sub-millisecond dead times for stopped-flow spectroscopy, which reveals that the maximum rate is indeed faster than previously recorded. The rate in the human enzyme, previously only estimated, is also able to be directly observed using cutting-edge stopped-flow instrumentation. In addition to the pH dependence of the hydride transfer rates for both enzymes, we examine the primary H/D kinetic isotope effect to reveal a temperature dependence in the human enzyme that is absent from the E. coli counterpart. This dependence, which appears above a temperature of 15 °C is a shared feature among other hydride transfer enzymes and is also consistent with computational work suggesting the presence of a fast promoting-vibration that provides donor-acceptor compression on the time scale of catalysis to facilitate the chemistry step.
Collapse
Affiliation(s)
- Brooke A Andrews
- Chemistry Department, Emory University, Atlanta, Georgia 30322, United States
| | - R Brian Dyer
- Chemistry Department, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
4
|
Mhashal AR, Major DT. Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations. J Phys Chem B 2021; 125:1369-1377. [PMID: 33522797 PMCID: PMC7883348 DOI: 10.1021/acs.jpcb.0c10318] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 01/05/2021] [Indexed: 11/28/2022]
Abstract
Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.
Collapse
Affiliation(s)
- Anil R. Mhashal
- Department of Chemistry and Institute
for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry and Institute
for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
5
|
Chen X, Schwartz SD. Multiple Reaction Pathways in the Morphinone Reductase-Catalyzed Hydride Transfer Reaction. ACS OMEGA 2020; 5:23468-23480. [PMID: 32954200 PMCID: PMC7496013 DOI: 10.1021/acsomega.0c03472] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 08/20/2020] [Indexed: 06/11/2023]
Abstract
Morphinone reductase (MR) is an important model system for studying the contribution of protein motions to H-transfer reactions. In this research, we used quantum mechanical/molecular mechanics (QM/MM) simulation together with transition path sampling (TPS) simulation to study two important topics of current research on MR: the existence of multiple catalytic reaction pathways and the involvement of fast protein motions in the catalytic process. We have discovered two reaction pathways for the wild type and three reaction pathways for the N189A mutant. With the committor distribution analysis method, we found reaction coordinates for all five reaction pathways. Only one wild-type reaction pathway has a rate-promoting vibration from His186, while all of the other four pathways do not involve any protein motions in their catalytic process through the transition state. The rate-promoting vibration in the wild-type MR, which comes from a direction perpendicular to the donor-acceptor axis, functions to decrease the donor-acceptor distance by causing a subtle "out-of-plane" motion of a donor atom. By comparing reaction pathways between the two enzymes, we concluded that the major effect of the N189A point mutation is to increase the active site volume by altering the active site backbone and eliminating the Asn189 side chain. This effect causes a different NADH geometry at the reactant state, which very well explains the different reaction mechanisms between the two enzymes, as well as the disappearance of the His186 rate-promoting vibrations in the N189A mutant. The unfavorable geometry of the NADH pyridine ring induced by the N189A point mutation is the potential cause of multiple reaction pathways in N189A mutants.
Collapse
|
6
|
Li J, Fortunato G, Lin J, Agarwal PK, Kohen A, Singh P, Cheatum CM. Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase. Biochemistry 2019; 58:3861-3868. [PMID: 31423766 PMCID: PMC7296831 DOI: 10.1021/acs.biochem.9b00460] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- Jiayue Li
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| | | | - Jennifer Lin
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| | - Pratul K. Agarwal
- Department of Biochemistry & Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996
| | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| | - Priyanka Singh
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| | | |
Collapse
|
7
|
Schramm VL, Schwartz SD. Promoting Vibrations and the Function of Enzymes. Emerging Theoretical and Experimental Convergence. Biochemistry 2018; 57:3299-3308. [PMID: 29608286 DOI: 10.1021/acs.biochem.8b00201] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A complete understanding of enzyme catalysis requires knowledge of both transition state features and the detailed motions of atoms that cause reactant molecules to form and traverse the transition state. The seeming intractability of the problem arises from the femtosecond lifetime of chemical transition states, preventing most experimental access. Computational chemistry is admirably suited to short time scale analysis but can be misled by inappropriate starting points or by biased assumptions. Kinetic isotope effects provide an experimental approach to transition state structure and a method for obtaining transition state analogues but, alone, do not inform how that transition state is reached. Enzyme structures with transition state analogues provide computational starting points near the transition state geometry. These well-conditioned starting points, combined with the unbiased computational method of transition path sampling, provide realistic atomistic motions involved in transition state formation and passage. In many, but not all, enzymatic systems, femtosecond local protein motions near the catalytic site are linked to transition state formation. These motions are not inherently revealed by most approaches of transition state theory, because transition state theory replaces dynamics with the statistics of the transition state. Experimental and theoretical convergence of the link between local catalytic site vibrational modes and catalysis comes from heavy atom ("Born-Oppenheimer") enzymes. Fully labeled and catalytic site local heavy atom labels perturb the probability of finding enzymatic transition states in ways that can be analyzed and predicted by transition path sampling. Recent applications of these experimental and computational approaches reveal how subpicosecond local catalytic site protein modes play important roles in creating the transition state.
Collapse
Affiliation(s)
- Vern L Schramm
- Department of Biochemistry , Albert Einstein College of Medicine , Bronx , New York 10461 , United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry , University of Arizona , Tucson , Arizona 85721 , United States
| |
Collapse
|
8
|
Chen X, Schwartz SD. Directed Evolution as a Probe of Rate Promoting Vibrations Introduced via Mutational Change. Biochemistry 2018; 57:3289-3298. [PMID: 29553716 DOI: 10.1021/acs.biochem.8b00185] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this article, we study with transition path sampling and reaction coordinate analysis how directed evolution in the Kemp eliminase family of artificial enzymes makes differential use of rapid rate promoting vibrations as a component of their chemical mechanism. Even though this family was initially created by placing the expected active site in a fixed protein matrix, we find a shift from largely static to more dynamic active sites that make use of donor-acceptor compression as the evolutionary process proceeds. We see that this introduction of dynamics significantly shifts the order of processes in the reaction. We also suggest that the lack of "design for dynamics" may help explain the relatively low proficiency of such designed enzymes.
Collapse
Affiliation(s)
- Xi Chen
- Department of Chemistry and Biochemistry , University of Arizona , 1306 East University Boulevard , Tucson , Arizona 85721 , United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry , University of Arizona , 1306 East University Boulevard , Tucson , Arizona 85721 , United States
| |
Collapse
|
9
|
Shanmuganatham KK, Wallace RS, Ting-I Lee A, Plapp BV. Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis. Protein Sci 2018; 27:750-768. [PMID: 29271062 DOI: 10.1002/pro.3370] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Revised: 12/18/2017] [Accepted: 12/20/2017] [Indexed: 01/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Karthik K Shanmuganatham
- Department of Biochemistry, The University of Iowa, Iowa City, IA, 52242-1109.,Diagnostic Virology Laboratory, USDA, Ames, IA, 50010
| | - Rachel S Wallace
- Department of Biochemistry, The University of Iowa, Iowa City, IA, 52242-1109.,Department of Physiology, School of Biomedical Sciences, University of Otago, Dunedin, 9054, New Zealand
| | - Ann Ting-I Lee
- Department of Biochemistry, The University of Iowa, Iowa City, IA, 52242-1109.,No 92, Jing Mao 1st Rd., Taichung, Taiwan, 406, Republic of China
| | - Bryce V Plapp
- Department of Biochemistry, The University of Iowa, Iowa City, IA, 52242-1109
| |
Collapse
|
10
|
Zinovjev K, Tuñón I. Reaction coordinates and transition states in enzymatic catalysis. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1329] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Kirill Zinovjev
- Departament de Química FísicaUniversitat de València Valencia Spain
| | - Iñaki Tuñón
- Departament de Química FísicaUniversitat de València Valencia Spain
| |
Collapse
|
11
|
Abstract
Enzyme isotope effects, or the kinetic effects of "heavy" enzymes, refer to the effect of isotopically labeled protein residues on the enzyme's activity or physical properties. These effects are increasingly employed in the examination of the possible contributions of protein dynamics to enzyme catalysis. One hypothesis assumed that isotopic substitution of all 12C, 14N, and nonexchangeable 1H by 13C, 15N, and 2H, would slow down protein picosecond to femtosecond dynamics without any effect on the system's electrostatics following the Born-Oppenheimer approximation. It was suggested that reduced reaction rates reported for several "heavy" enzymes accords with that hypothesis. However, numerous deviations from the predictions of that hypothesis were also reported. Current studies also attempt to test the role of individual residues by site-specific labeling or by labeling a pattern of residues on activity. It appears that in several systems the protein's fast dynamics are indeed reduced in "heavy" enzymes in a way that reduces the probability of barrier crossing of its chemical step. Other observations, however, indicated that slower protein dynamics are electrostatically altered in isotopically labeled enzymes. Interestingly, these effects appear to be system dependent, thus it might be premature to suggest a general role of "heavy" enzymes' effect on catalysis.
Collapse
|
12
|
Catalytic-site design for inverse heavy-enzyme isotope effects in human purine nucleoside phosphorylase. Proc Natl Acad Sci U S A 2017; 114:6456-6461. [PMID: 28584087 DOI: 10.1073/pnas.1704786114] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Heavy-enzyme isotope effects (15N-, 13C-, and 2H-labeled protein) explore mass-dependent vibrational modes linked to catalysis. Transition path-sampling (TPS) calculations have predicted femtosecond dynamic coupling at the catalytic site of human purine nucleoside phosphorylase (PNP). Coupling is observed in heavy PNPs, where slowed barrier crossing caused a normal heavy-enzyme isotope effect (kchemlight/kchemheavy > 1.0). We used TPS to design mutant F159Y PNP, predicted to improve barrier crossing for heavy F159Y PNP, an attempt to generate a rare inverse heavy-enzyme isotope effect (kchemlight/kchemheavy < 1.0). Steady-state kinetic comparison of light and heavy native PNPs to light and heavy F159Y PNPs revealed similar kinetic properties. Pre-steady-state chemistry was slowed 32-fold in F159Y PNP. Pre-steady-state chemistry compared heavy and light native and F159Y PNPs and found a normal heavy-enzyme isotope effect of 1.31 for native PNP and an inverse effect of 0.75 for F159Y PNP. Increased isotopic mass in F159Y PNP causes more efficient transition state formation. Independent validation of the inverse isotope effect for heavy F159Y PNP came from commitment to catalysis experiments. Most heavy enzymes demonstrate normal heavy-enzyme isotope effects, and F159Y PNP is a rare example of an inverse effect. Crystal structures and TPS dynamics of native and F159Y PNPs explore the catalytic-site geometry associated with these catalytic changes. Experimental validation of TPS predictions for barrier crossing establishes the connection of rapid protein dynamics and vibrational coupling to enzymatic transition state passage.
Collapse
|
13
|
Varga MJ, Dzierlenga MW, Schwartz SD. Structurally Linked Dynamics in Lactate Dehydrogenases of Evolutionarily Distinct Species. Biochemistry 2017; 56:2488-2496. [PMID: 28445027 DOI: 10.1021/acs.biochem.7b00245] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present new findings about how primary and secondary structure affects the role of fast protein motions in the reaction coordinates of enzymatic reactions. Using transition path sampling and committor distribution analysis, we examined the difference in the role of these fast protein motions in the reaction coordinate of lactate dehydrogenases (LDHs) of Apicomplexa organisms Plasmodium falciparum and Cryptosporidium parvum. Having evolved separately from a common malate dehydrogenase ancestor, the two enzymes exhibit several important structural differences, notably a five-amino acid insertion in the active site loop of P. falciparum LDH. We find that these active site differences between the two organisms' LDHs likely cause a decrease in the contribution of the previously determined LDH rate-promoting vibration to the reaction coordinate of P. falciparum LDH compared to that of C. parvum LDH, specifically in the coupling of the rate-promoting vibration and the hydride transfer. This effect, while subtle, directly shows how changes in structure near the active site of LDH alter catalytically important motions. Insights provided by studying these alterations would prove to be useful in identifying LDH inhibitors that specifically target the isozymes of these parasitic organisms.
Collapse
Affiliation(s)
- Matthew J Varga
- Department of Chemistry and Biochemistry, University of Arizona , 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Michael W Dzierlenga
- Department of Chemistry and Biochemistry, University of Arizona , 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry, University of Arizona , 1306 East University Boulevard, Tucson, Arizona 85721, United States
| |
Collapse
|
14
|
Wang Z, Chang EP, Schramm VL. Triple Isotope Effects Support Concerted Hydride and Proton Transfer and Promoting Vibrations in Human Heart Lactate Dehydrogenase. J Am Chem Soc 2016; 138:15004-15010. [PMID: 27766841 DOI: 10.1021/jacs.6b09049] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Transition path sampling simulations have proposed that human heart lactate dehydrogenase (LDH) employs protein promoting vibrations (PPVs) on the femtosecond (fs) to picosecond (ps) time scale to promote crossing of the chemical barrier. This chemical barrier involves both hydride and proton transfers to pyruvate to form l-lactate, using reduced nicotinamide adenine dinucleotide (NADH) as the cofactor. Here we report experimental evidence from three types of isotope effect experiments that support coupling of the promoting vibrations to barrier crossing and the coincidence of hydride and proton transfer. We prepared the native (light) LDH and a heavy LDH labeled with 13C, 15N, and nonexchangeable 2H (D) to perturb the predicted PPVs. Heavy LDH has slowed chemistry in single turnover experiments, supporting a contribution of PPVs to transition state formation. Both the [4-2H]NADH (NADD) kinetic isotope effect and the D2O solvent isotope effect were increased in dual-label experiments combining both NADD and D2O, a pattern maintained with both light and heavy LDHs. These isotope effects support concerted hydride and proton transfer for both light and heavy LDHs. Although the transition state barrier-crossing probability is reduced in heavy LDH, the concerted mechanism of the hydride-proton transfer reaction is not altered. This study takes advantage of triple isotope effects to resolve the chemical mechanism of LDH and establish the coupling of fs-ps protein dynamics to barrier crossing.
Collapse
Affiliation(s)
- Zhen Wang
- Department of Biochemistry, Albert Einstein College of Medicine , Bronx, New York 10461, United States
| | - Eric P Chang
- Department of Biochemistry, Albert Einstein College of Medicine , Bronx, New York 10461, United States
| | - Vern L Schramm
- Department of Biochemistry, Albert Einstein College of Medicine , Bronx, New York 10461, United States
| |
Collapse
|
15
|
Luk LYP, Loveridge EJ, Allemann RK. Protein motions and dynamic effects in enzyme catalysis. Phys Chem Chem Phys 2016; 17:30817-27. [PMID: 25854702 DOI: 10.1039/c5cp00794a] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.
Collapse
Affiliation(s)
- Louis Y P Luk
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, UK.
| | - E Joel Loveridge
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, UK.
| | - Rudolf K Allemann
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, UK.
| |
Collapse
|
16
|
Pshetitsky Y, Eitan R, Verner G, Kohen A, Major DT. Improved Sugar Puckering Profiles for Nicotinamide Ribonucleoside for Hybrid QM/MM Simulations. J Chem Theory Comput 2016; 12:5179-5189. [PMID: 27490188 DOI: 10.1021/acs.jctc.6b00401] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The coenzyme nicotinamide adenine dinucleotide (NAD+) and its reduced form (NADH) play ubiquitous roles as oxidizing and reducing agents in nature. The binding, and possibly the chemical redox step, of NAD+/NADH may be influenced by the cofactor conformational distribution and, in particular, by the ribose puckering of its nicotinamide-ribonucleoside (NR) moiety. In many hybrid quantum mechanics-molecular mechanics (QM/MM) studies of NAD+/NADH dependent enzymes, the QM region is treated by semiempirical (SE) methods. Recent work suggests that SE methods do not adequately describe the ring puckering in sugar molecules. In the present work we adopt an efficient and practical strategy to correct for this deficiency for NAD+/NADH. We have implemented a cost-effective correction to a SE Hamiltonian by adding a correction potential, which is defined as the difference between an accurate benchmark density functional theory (DFT) potential energy surface (PES) and the SE PES. In practice, this is implemented via a B-spline interpolation scheme for the grid-based potential energy difference surface. We find that the puckering population distributions obtained from free energy QM(SE)/MM simulations are in good agreement with DFT and in fair accord with experimental results. The corrected PES should facilitate a more accurate description of the ribose puckering in the NAD+/NADH cofactor in simulations of biological systems.
Collapse
Affiliation(s)
- Yaron Pshetitsky
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Reuven Eitan
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Gilit Verner
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| | - Amnon Kohen
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 52900, Israel
| |
Collapse
|
17
|
Abstract
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.
Collapse
Affiliation(s)
- P Singh
- University of Iowa, Iowa City, IA, United States
| | - Z Islam
- University of Iowa, Iowa City, IA, United States
| | - A Kohen
- University of Iowa, Iowa City, IA, United States.
| |
Collapse
|
18
|
Zoi I, Suarez J, Antoniou D, Cameron SA, Schramm VL, Schwartz SD. Modulating Enzyme Catalysis through Mutations Designed to Alter Rapid Protein Dynamics. J Am Chem Soc 2016; 138:3403-9. [PMID: 26927977 PMCID: PMC4794390 DOI: 10.1021/jacs.5b12551] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The relevance of sub-picosecond protein motions to the catalytic event remains a topic of debate. Heavy enzymes (isotopically substituted) provide an experimental tool for bond-vibrational links to enzyme catalysis. A recent transition path sampling study with heavy purine nucleoside phosphorylase (PNP) characterized the experimentally observed mass-dependent slowing of barrier crossing (Antoniou, D.; Ge, X.; Schramm, V. L.; Schwartz, S. D. J. Phys. Chem. Lett. 2012, 3, 3538). Here we computationally identify second-sphere amino acid residues predicted to influence the freedom of the catalytic site vibrational modes linked to heavy enzyme effects in PNP. We mutated heavy and light PNPs to increase the catalytic site vibrational freedom. Enzymatic barrier-crossing rates were converted from mass-dependent to mass-independent as a result of the mutations. The mutagenic uncoupling of femtosecond motions between catalytic site groups and reactants decreased transition state barrier crossing by 2 orders of magnitude, an indication of the femtosecond dynamic contributions to catalysis.
Collapse
Affiliation(s)
- Ioanna Zoi
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| | - Javier Suarez
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Dimitri Antoniou
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| | - Scott A. Cameron
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Vern L. Schramm
- Department of Biochemistry, Albert Einstein College of Medicine, 1300 Morris Park Ave, Bronx, New York 10461, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd, Tucson, Arizona 85721, United States
| |
Collapse
|
19
|
Varga MJ, Schwartz SD. Enzymatic Kinetic Isotope Effects from First-Principles Path Sampling Calculations. J Chem Theory Comput 2016; 12:2047-54. [PMID: 26949835 DOI: 10.1021/acs.jctc.5b01169] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this study, we develop and test a method to determine the rate of particle transfer and kinetic isotope effects in enzymatic reactions, specifically yeast alcohol dehydrogenase (YADH), from first-principles. Transition path sampling (TPS) and normal mode centroid dynamics (CMD) are used to simulate these enzymatic reactions without knowledge of their reaction coordinates and with the inclusion of quantum effects, such as zero-point energy and tunneling, on the transferring particle. Though previous studies have used TPS to calculate reaction rate constants in various model and real systems, it has not been applied to a system as large as YADH. The calculated primary H/D kinetic isotope effect agrees with previously reported experimental results, within experimental error. The kinetic isotope effects calculated with this method correspond to the kinetic isotope effect of the transfer event itself. The results reported here show that the kinetic isotope effects calculated from first-principles, purely for barrier passage, can be used to predict experimental kinetic isotope effects in enzymatic systems.
Collapse
Affiliation(s)
- Matthew J Varga
- Department of Chemistry and Biochemistry, University of Arizona , Tucson, Arizona 85721, United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry, University of Arizona , Tucson, Arizona 85721, United States
| |
Collapse
|
20
|
Francis K, Sapienza PJ, Lee AL, Kohen A. The Effect of Protein Mass Modulation on Human Dihydrofolate Reductase. Biochemistry 2016; 55:1100-6. [PMID: 26813442 DOI: 10.1021/acs.biochem.5b00945] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dihydrofolate reductase (DHFR) from Escherichia coli has long served as a model enzyme with which to elucidate possible links between protein dynamics and the catalyzed reaction. Such physical properties of its human counterpart have not been rigorously studied so far, but recent computer-based simulations suggest that these two DHFRs differ significantly in how closely coupled the protein dynamics and the catalyzed C-H → C hydride transfer step are. To test this prediction, two contemporary probes for studying the effect of protein dynamics on catalysis were combined here: temperature dependence of intrinsic kinetic isotope effects (KIEs), which are sensitive to the physical nature of the chemical step, and protein mass modulation, which slows down fast dynamics (femto- to picosecond time scale) throughout the protein. The intrinsic H/T KIEs of human DHFR, like those of E. coli DHFR, are shown to be temperature-independent in the range from 5 to 45 °C, indicating fast sampling of donor and acceptor distances (DADs) at the reaction's transition state (or tunneling ready state, TRS). Mass modulation of these enzymes through isotopic labeling with (13)C, (15)N, and (2)H at nonexchangeable hydrogens yields an 11% heavier enzyme. The additional mass has no effect on the intrinsic KIEs of the human enzyme. This finding indicates that the mass modulation of the human DHFR affects neither DAD distribution nor the DAD's conformational sampling dynamics. Furthermore, reduction in the enzymatic turnover number and the dissociation rate constant for the product indicate that the isotopic substitution affects kinetic steps that are not the catalyzed C-H → C hydride transfer. The findings are discussed in terms of fast dynamics and their role in catalysis, the comparison of calculations and experiments, and the interpretation of isotopically modulated heavy enzymes in general.
Collapse
Affiliation(s)
- Kevin Francis
- The Department of Chemistry, The University of Iowa , Iowa City, Iowa 52242, United States
| | - Paul J Sapienza
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina , Chapel Hill, North Carolina 27599, United States
| | - Andrew L Lee
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina , Chapel Hill, North Carolina 27599, United States
| | - Amnon Kohen
- The Department of Chemistry, The University of Iowa , Iowa City, Iowa 52242, United States
| |
Collapse
|
21
|
Antoniou D, Schwartz SD. Phase Space Bottlenecks in Enzymatic Reactions. J Phys Chem B 2016; 120:433-9. [PMID: 26756622 PMCID: PMC4734068 DOI: 10.1021/acs.jpcb.5b11157] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2015] [Revised: 01/11/2016] [Indexed: 01/09/2023]
Abstract
The definition of a transition state on an individual reactive trajectory is made via a committor analysis. In the past, the bottleneck definition has often been applied in configuration space. This is an approximation, and in order to expand this definition, we are revisiting an enzyme in which we had identified a fast subpicosecond motion that makes the reaction possible. First we used a time-series analysis method to identify the exact time when this motion initiates donor-acceptor compression. Then we modified the standard committor analysis of transition path sampling to identify events in phase space and found that there is a dividing surface in phase space significantly earlier than the configurationally defined transition-state crossing.
Collapse
Affiliation(s)
- Dimitri Antoniou
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| |
Collapse
|
22
|
Wang Z, Antoniou D, Schwartz SD, Schramm VL. Hydride Transfer in DHFR by Transition Path Sampling, Kinetic Isotope Effects, and Heavy Enzyme Studies. Biochemistry 2015; 55:157-66. [PMID: 26652185 DOI: 10.1021/acs.biochem.5b01241] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Escherichia coli dihydrofolate reductase (ecDHFR) is used to study fundamental principles of enzyme catalysis. It remains controversial whether fast protein motions are coupled to the hydride transfer catalyzed by ecDHFR. Previous studies with heavy ecDHFR proteins labeled with (13)C, (15)N, and nonexchangeable (2)H reported enzyme mass-dependent hydride transfer kinetics for ecDHFR. Here, we report refined experimental and computational studies to establish that hydride transfer is independent of protein mass. Instead, we found the rate constant for substrate dissociation to be faster for heavy DHFR. Previously reported kinetic differences between light and heavy DHFRs likely arise from kinetic steps other than the chemical step. This study confirms that fast (femtosecond to picosecond) protein motions in ecDHFR are not coupled to hydride transfer and provides an integrative computational and experimental approach to resolve fast dynamics coupled to chemical steps in enzyme catalysis.
Collapse
Affiliation(s)
- Zhen Wang
- Department of Biochemistry, Albert Einstein College of Medicine , 1300 Morris Park Avenue, Bronx, New York 10461, United States
| | - Dimitri Antoniou
- Department of Chemistry and Biochemistry, University of Arizona , P.O. Box 210041, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Steven D Schwartz
- Department of Chemistry and Biochemistry, University of Arizona , P.O. Box 210041, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Vern L Schramm
- Department of Biochemistry, Albert Einstein College of Medicine , 1300 Morris Park Avenue, Bronx, New York 10461, United States
| |
Collapse
|
23
|
Abstract
Dihydrofolate reductase from
Escherichia coli (ecDHFR) serves as a model system for investigating the role of protein dynamics in enzyme catalysis. We discuss calculations predicting a network of dynamic motions that is coupled to the chemical step catalyzed by this enzyme. Kinetic studies testing these predictions are presented, and their potential use in better understanding the role of these dynamics in enzyme catalysis is considered. The cumulative results implicate motions across the entire protein in catalysis.
Collapse
Affiliation(s)
- Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA, USA
| |
Collapse
|
24
|
Change in heat capacity accurately predicts vibrational coupling in enzyme catalyzed reactions. FEBS Lett 2015; 589:2200-6. [DOI: 10.1016/j.febslet.2015.06.042] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2015] [Revised: 06/24/2015] [Accepted: 06/30/2015] [Indexed: 11/23/2022]
|
25
|
Singh P, Francis K, Kohen A. Network of remote and local protein dynamics in dihydrofolate reductase catalysis. ACS Catal 2015; 5:3067-3073. [PMID: 27182453 DOI: 10.1021/acscatal.5b00331] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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 (ΔEa(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔEa(T-H)G121V +ΔEa(T-H) I14A < ΔEa(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.
Collapse
Affiliation(s)
- Priyanka Singh
- The Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| | - Kevin Francis
- The Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| | - Amnon Kohen
- The Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
26
|
Abstract
Enzymatic transition states have lifetimes of a few femtoseconds (fs). Computational analysis of enzyme motions leading to transition state formation suggests that local catalytic site motions on the fs time scale provide the mechanism to locate transition states. An experimental test of protein fs motion and its relation to transition state formation can be provided by isotopically heavy proteins. Heavy enzymes have predictable mass-altered bond vibration states without altered electrostatic properties, according to the Born-Oppenheimer approximation. On-enzyme chemistry is slowed in most heavy proteins, consistent with altered protein bond frequencies slowing the search for the transition state. In other heavy enzymes, structural changes involved in reactant binding and release are also influenced. Slow protein motions associated with substrate binding and catalytic site preorganization are essential to allow the subsequent fs motions to locate the transition state and to facilitate the efficient release of products. In the catalytically competent geometry, local groups move in stochastic atomic motion on the fs time scale, within transition state-accessible conformations created by slower protein motions. The fs time scale for the transition state motions does not permit thermodynamic equilibrium between the transition state and stable enzyme states. Isotopically heavy enzymes provide a diagnostic tool for fast coupled protein motions to transition state formation and mass-dependent conformational changes. The binding of transition state analogue inhibitors is the opposite in catalytic time scale to formation of the transition state but is related by similar geometries of the enzyme-transition state and enzyme-inhibitor interactions. While enzymatic transition states have lifetimes as short as 10(-15) s, transition state analogues can bind tightly to enzymes with release rates greater than 10(3) s. Tight-binding transition state analogues stabilize the rare but evolved enzymatic geometry to form the transition state. Evolution to efficient catalysis optimized this geometry and its stabilization by a transition state mimic results in tight binding. Release rates of transition state analogues are orders of magnitude slower than product release in normal catalytic function. During catalysis, product release is facilitated by altered chemistry. Compared to the weak associations found in Michaelis complexes, transition state analogues involve strong interactions related to those in the transition state. Optimum binding of transition state analogues occurs when the complex retains the system motions intrinsic to transition state formation. Conserved dynamic motion retains the entropic components of inhibitor complexes, improving the thermodynamics of analogue binding.
Collapse
Affiliation(s)
- Vern L. Schramm
- Department
of Biochemistry, Albert Einstein College of Medicine, 1300 Morris
Park Avenue, Bronx, New York 10461, United States
| |
Collapse
|
27
|
Bhabha G, Biel JT, Fraser JS. Keep on moving: discovering and perturbing the conformational dynamics of enzymes. Acc Chem Res 2015; 48:423-30. [PMID: 25539415 PMCID: PMC4334266 DOI: 10.1021/ar5003158] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
CONSPECTUS: Because living organisms are in constant motion, the word "dynamics" can hold many meanings to biologists. Here we focus specifically on the conformational changes that occur in proteins and how studying these protein dynamics may provide insights into enzymatic catalysis. Advances in integrating techniques such as X-ray crystallography, nuclear magnetic resonance, and electron cryomicroscopy (cryo EM) allow us to model the dominant structures and exchange rates for many proteins and protein complexes. For proteins amenable to atomic resolution techniques, the major questions shift from simply describing the motions to discovering their role in function. Concurrently, there is an increasing need for using perturbations to test predictive models of dynamics-function relationships. Examples are the catalytic cycles of dihydrofolate reductase (DHFR) and cyclophilin A (CypA). In DHFR, mutations that alter the ability of the active site to sample productive higher energy states on the millisecond time scale reduce the rate of hydride transfer significantly. Recently identified rescue mutations restore function, but the mechanism by which they do so remains unclear. The exact role of any changes in the dynamics remains an open question. For CypA, a network of side chains that exchange between two conformations is critical for catalysis. Mutations that lock the network in one state also reduce catalytic activity. A further understanding of enzyme dynamics of well-studied enzymes such as dihydrofolate reductase and cyclophilin A will lead to improvement in ability to modulate the functions of computationally designed enzymes and large macromolecular machines. In designed enzymes, directed evolution experiments increase catalytic rates. Detailed X-ray studies suggest that these mutations likely limit the conformational space explored by residues in the active site. For proteins where atomic resolution information is currently inaccessible, other techniques such as cryo-EM and high-resolution single molecule microscopy continue to advance. Understanding the conformational dynamics of larger systems such as protein machines will likely become more accessible and provide new opportunities to rationally modulate protein function.
Collapse
Affiliation(s)
- Gira Bhabha
- Department
of Cellular and Molecular Pharmacology and the Howard Hughes Medical
Institute, University of California, San Francisco, California 94158, United States
| | - Justin T. Biel
- Biophysics Graduate Group, University of California, San Francisco, California 94158, United States
| | - James S. Fraser
- Department of Bioengineering and Therapeutic
Science and California Institute for Quantitative Biology, University of California, San Francisco, California 94158, United States
| |
Collapse
|
28
|
Singh P, Abeysinghe T, Kohen A. Linking protein motion to enzyme catalysis. Molecules 2015; 20:1192-209. [PMID: 25591120 PMCID: PMC4341894 DOI: 10.3390/molecules20011192] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2014] [Accepted: 01/07/2015] [Indexed: 12/01/2022] Open
Abstract
Enzyme motions on a broad range of time scales can play an important role in various intra- and intermolecular events, including substrate binding, catalysis of the chemical conversion, and product release. The relationship between protein motions and catalytic activity is of contemporary interest in enzymology. To understand the factors influencing the rates of enzyme-catalyzed reactions, the dynamics of the protein-solvent-ligand complex must be considered. The current review presents two case studies of enzymes—dihydrofolate reductase (DHFR) and thymidylate synthase (TSase)—and discusses the role of protein motions in their catalyzed reactions. Specifically, we will discuss the utility of kinetic isotope effects (KIEs) and their temperature dependence as tools in probing such phenomena.
Collapse
Affiliation(s)
- Priyanka Singh
- Department of Chemistry, University of Iowa, Iowa City, IA 52242, USA.
| | - Thelma Abeysinghe
- Department of Chemistry, University of Iowa, Iowa City, IA 52242, USA.
| | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA 52242, USA.
| |
Collapse
|
29
|
Doron D, Stojković V, Gakhar L, Vardi-Kilshtain A, Kohen A, Major DT. Free energy simulations of active-site mutants of dihydrofolate reductase. J Phys Chem B 2014; 119:906-16. [PMID: 25382260 DOI: 10.1021/jp5059963] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
This study employs hybrid quantum mechanics-molecular mechanics (QM/MM) simulations to investigate the effect of mutations of the active-site residue I14 of E. coli dihydrofolate reductase (DHFR) on the hydride transfer. Recent kinetic measurements of the I14X mutants (X = V, A, and G) indicated slower hydride transfer rates and increasingly temperature-dependent kinetic isotope effects (KIEs) with systematic reduction of the I14 side chain. The QM/MM simulations show that when the original isoleucine residue is substituted in silico by valine, alanine, or glycine (I14V, I14A, and I14G DHFR, respectively), the free energy barrier height of the hydride transfer reaction increases relative to the wild-type enzyme. These trends are in line with the single-turnover rate measurements reported for these systems. In addition, extended dynamics simulations of the reactive Michaelis complex reveal enhanced flexibility in the mutants, and in particular for the I14G mutant, including considerable fluctuations of the donor-acceptor distance (DAD) and the active-site hydrogen bonding network compared with those detected in the native enzyme. These observations suggest that the perturbations induced by the mutations partly impair the active-site environment in the reactant state. On the other hand, the average DADs at the transition state of all DHFR variants are similar. Crystal structures of I14 mutants (V, A, and G) confirmed the trend of increased flexibility of the M20 and other loops.
Collapse
Affiliation(s)
- Dvir Doron
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 5290002, Israel
| | | | | | | | | | | |
Collapse
|
30
|
Masterson JE, Schwartz SD. Evolution alters the enzymatic reaction coordinate of dihydrofolate reductase. J Phys Chem B 2014; 119:989-96. [PMID: 25369552 PMCID: PMC4306496 DOI: 10.1021/jp506373q] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
How evolution has affected enzyme function is a topic of great interest in the field of biophysical chemistry. Evolutionary changes from Escherichia coli dihydrofolate reductase (ecDHFR) to human dihydrofolate reductase (hsDHFR) have resulted in increased catalytic efficiency and an altered dynamic landscape in the human enzyme. Here, we show that a subpicosecond protein motion is dynamically coupled to hydride transfer catalyzed by hsDHFR but not ecDHFR. This motion propagates through residues that correspond to mutational events along the evolutionary path from ecDHFR to hsDHFR. We observe an increase in the variability of the transition states, reactive conformations, and times of barrier crossing in the human system. In the hsDHFR active site, we detect structural changes that have enabled the coupling of fast protein dynamics to the reaction coordinate. These results indicate a shift in the DHFR family to a form of catalysis that incorporates rapid protein dynamics and a concomitant shift to a more flexible path through reactive phase space.
Collapse
Affiliation(s)
- Jean E Masterson
- Department of Chemistry and Biochemistry, University of Arizona , 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | | |
Collapse
|
31
|
Wang Z, Singh P, Czekster CM, Kohen A, Schramm VL. Protein mass-modulated effects in the catalytic mechanism of dihydrofolate reductase: beyond promoting vibrations. J Am Chem Soc 2014; 136:8333-41. [PMID: 24820793 PMCID: PMC4063187 DOI: 10.1021/ja501936d] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The role of fast protein dynamics
in enzyme catalysis has been
of great interest in the past decade. Recent “heavy enzyme”
studies demonstrate that protein mass-modulated vibrations are linked
to the energy barrier for the chemical step of catalyzed reactions.
However, the role of fast dynamics in the overall catalytic mechanism
of an enzyme has not been addressed. Protein mass-modulated effects
in the catalytic mechanism of Escherichia coli dihydrofolate
reductase (ecDHFR) are explored by isotopic substitution (13C, 15N, and non-exchangeable 2H) of the wild-type
ecDHFR (l-DHFR) to generate a vibrationally perturbed
“heavy ecDHFR” (h-DHFR). Steady-state,
pre-steady-state, and ligand binding kinetics, intrinsic kinetic isotope
effects (KIEint) on the chemical step, and thermal unfolding
experiments of both l- and h-DHFR
show that the altered protein mass affects the conformational ensembles
and protein–ligand interactions, but does not affect the hydride
transfer at physiological temperatures (25–45 °C). Below
25 °C, h-DHFR shows altered transition state
(TS) structure and increased barrier-crossing probability of the chemical
step compared with l-DHFR, indicating temperature-dependent
protein vibrational coupling to the chemical step. Protein mass-modulated
vibrations in ecDHFR are involved in TS interactions at cold temperatures
and are linked to dynamic motions involved in ligand binding at physiological
temperatures. Thus, mass effects can affect enzymatic catalysis beyond
alterations in promoting vibrations linked to chemistry.
Collapse
Affiliation(s)
- Zhen Wang
- Department of Biochemistry, Albert Einstein College of Medicine , Bronx, New York 10461, United States
| | | | | | | | | |
Collapse
|
32
|
Doron D, Kohen A, Nam K, Major DT. How Accurate Are Transition States from Simulations of Enzymatic Reactions? J Chem Theory Comput 2014; 10:1863-1871. [PMID: 24860275 PMCID: PMC4025581 DOI: 10.1021/ct5000742] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2014] [Indexed: 11/30/2022]
Abstract
The rate expression of traditional transition state theory (TST) assumes no recrossing of the transition state (TS) and thermal quasi-equilibrium between the ground state and the TS. Currently, it is not well understood to what extent these assumptions influence the nature of the activated complex obtained in traditional TST-based simulations of processes in the condensed phase in general and in enzymes in particular. Here we scrutinize these assumptions by characterizing the TSs for hydride transfer catalyzed by the enzyme Escherichia coli dihydrofolate reductase obtained using various simulation approaches. Specifically, we compare the TSs obtained with common TST-based methods and a dynamics-based method. Using a recently developed accurate hybrid quantum mechanics/molecular mechanics potential, we find that the TST-based and dynamics-based methods give considerably different TS ensembles. This discrepancy, which could be due equilibrium solvation effects and the nature of the reaction coordinate employed and its motion, raises major questions about how to interpret the TSs determined by common simulation methods. We conclude that further investigation is needed to characterize the impact of various TST assumptions on the TS phase-space ensemble and on the reaction kinetics.
Collapse
Affiliation(s)
- Dvir Doron
- Department
of Chemistry and the Lise Meitner-Minerva Center of Computational
Quantum Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Amnon Kohen
- Department
of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - Kwangho Nam
- Department
of Chemistry and Computational Life Science Cluster (CLiC), Umeå University, 901 87 Umeå, Sweden
| | - Dan Thomas Major
- Department
of Chemistry and the Lise Meitner-Minerva Center of Computational
Quantum Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
33
|
Roston D, Kohen A, Doron D, Major DT. Simulations of remote mutants of dihydrofolate reductase reveal the nature of a network of residues coupled to hydride transfer. J Comput Chem 2014; 35:1411-7. [PMID: 24798860 DOI: 10.1002/jcc.23629] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2014] [Revised: 03/11/2014] [Accepted: 04/06/2014] [Indexed: 11/07/2022]
Abstract
Recent experimental and theoretical studies have proposed that enzymes involve networks of coupled residues throughout the protein that participate in motions accompanying chemical barrier crossing. Here, we have examined portions of a proposed network in dihydrofolate reductase (DHFR) using quantum mechanics/molecular mechanics simulations. The simulations use a hybrid quantum mechanics-molecular mechanics approach with a recently developed semiempirical AM1-SRP Hamiltonian that provides accurate results for this reaction. The simulations reproduce experimentally determined catalytic rates for the wild type and distant mutants of E. coli DHFR, underscoring the accuracy of the simulation protocol. Additionally, the simulations provide detailed insight into how residues remote from the active site affect the catalyzed chemistry, through changes in the thermally averaged properties along the reaction coordinate. The mutations do not greatly affect the structure of the transition state near the bond activation, but we observe differences somewhat removed from the point of C-H cleavage that affect the rate. The mutations have global effects on the thermally averaged structure that propagate throughout the enzyme and the current simulations highlight several interactions that appear to be particularly important.
Collapse
Affiliation(s)
- Daniel Roston
- Department of Chemistry, University of Iowa, Iowa City, Iowa, 52242
| | | | | | | |
Collapse
|
34
|
Francis K, Stojković V, Kohen A. Preservation of protein dynamics in dihydrofolate reductase evolution. J Biol Chem 2013; 288:35961-8. [PMID: 24158440 PMCID: PMC3861645 DOI: 10.1074/jbc.m113.507632] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2013] [Revised: 10/09/2013] [Indexed: 11/06/2022] Open
Abstract
The hydride transfer reaction catalyzed by dihydrofolate reductase (DHFR) is a model for examining how protein dynamics contribute to enzymatic function. The relationship between functional motions and enzyme evolution has attracted significant attention. Recent studies on N23PP Escherichia coli DHFR (ecDHFR) mutant, designed to resemble parts of the human enzyme, indicated a reduced single turnover rate. NMR relaxation dispersion experiments with that enzyme showed rigidification of millisecond Met-20 loop motions (Bhabha, G., Lee, J., Ekiert, D. C., Gam, J., Wilson, I. A., Dyson, H. J., Benkovic, S. J., and Wright, P. E. (2011) Science 332, 234-238). A more recent study of this mutant, however, indicated that fast motions along the reaction coordinate are actually more dispersed than for wild-type ecDHFR (WT). Furthermore, a double mutant (N23PP/G51PEKN) that better mimics the human enzyme seems to restore both the single turnover rates and narrow distribution of fast dynamics (Liu, C. T., Hanoian, P., French, T. H., Hammes-Schiffer, S., and Benkovic, S. J. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 10159-11064). Here, we measured intrinsic kinetic isotope effects for both N23PP and N23PP/G51PEKN double mutant DHFRs over a temperature range. The findings indicate that although the C-H→C transfer and dynamics along the reaction coordinate are impaired in the altered N23PP mutant, both seem to be restored in the N23PP/G51PEKN double mutant. This indicates that the evolution of G51PEKN, although remote from the Met-20 loop, alleviated the loop rigidification that would have been caused by N23PP, enabling WT-like H-tunneling. The correlation between the calculated dynamics, the nature of C-H→C transfer, and a phylogenetic analysis of DHFR sequences are consistent with evolutionary preservation of the protein dynamics to enable H-tunneling from well reorganized active sites.
Collapse
Affiliation(s)
- Kevin Francis
- From the Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242
| | - Vanja Stojković
- From the Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242
| | - Amnon Kohen
- From the Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242
| |
Collapse
|
35
|
Ruiz-Pernia JJ, Luk LYP, García-Meseguer R, Martí S, Loveridge EJ, Tuñón I, Moliner V, Allemann RK. Increased dynamic effects in a catalytically compromised variant of Escherichia coli dihydrofolate reductase. J Am Chem Soc 2013; 135:18689-96. [PMID: 24252106 PMCID: PMC3949409 DOI: 10.1021/ja410519h] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2013] [Indexed: 01/19/2023]
Abstract
Isotopic substitution ((15)N, (13)C, (2)H) of a catalytically compromised variant of Escherichia coli dihydrofolate reductase, EcDHFR-N23PP/S148A, has been used to investigate the effect of these mutations on catalysis. The reduction of the rate constant of the chemical step in the EcDHFR-N23PP/S148A catalyzed reaction is essentially a consequence of an increase of the quasi-classical free energy barrier and to a minor extent of an increased number of recrossing trajectories on the transition state dividing surface. Since the variant enzyme is less well set up to catalyze the reaction, a higher degree of active site reorganization is needed to reach the TS. Although millisecond active site motions are lost in the variant, there is greater flexibility on the femtosecond time scale. The "dynamic knockout" EcDHFR-N23PP/S148A is therefore a "dynamic knock-in" at the level of the chemical step, and the increased dynamic coupling to the chemical coordinate is in fact detrimental to catalysis. This finding is most likely applicable not just to hydrogen transfer in EcDHFR but also to other enzymatic systems.
Collapse
Affiliation(s)
- J. Javier Ruiz-Pernia
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castello, Spain
| | - Louis Y. P. Luk
- School
of Chemistry & Cardiff Catalysis Institute, Cardiff University, Park Place, Cardiff, CF10
3AT, U.K.
| | | | - Sergio Martí
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castello, Spain
| | - E. Joel Loveridge
- School
of Chemistry & Cardiff Catalysis Institute, Cardiff University, Park Place, Cardiff, CF10
3AT, U.K.
| | - Iñaki Tuñón
- Departament
de Química Física, Universitat
de València, 46100 Burjassot, Spain
| | - Vicent Moliner
- Departament
de Química Física i Analítica, Universitat Jaume I, 12071 Castello, Spain
| | - Rudolf K. Allemann
- School
of Chemistry & Cardiff Catalysis Institute, Cardiff University, Park Place, Cardiff, CF10
3AT, U.K.
| |
Collapse
|
36
|
Roston D, Islam Z, Kohen A. Isotope effects as probes for enzyme catalyzed hydrogen-transfer reactions. Molecules 2013; 18:5543-67. [PMID: 23673528 PMCID: PMC4342783 DOI: 10.3390/molecules18055543] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2013] [Revised: 04/30/2013] [Accepted: 05/03/2013] [Indexed: 11/16/2022] Open
Abstract
Kinetic Isotope effects (KIEs) have long served as a probe for the mechanisms of both enzymatic and solution reactions. Here, we discuss various models for the physical sources of KIEs, how experimentalists can use those models to interpret their data, and how the focus of traditional models has grown to a model that includes motion of the enzyme and quantum mechanical nuclear tunneling. We then present two case studies of enzymes, thymidylate synthase and alcohol dehydrogenase, and discuss how KIEs have shed light on the C-H bond cleavages those enzymes catalyze. We will show how the combination of both experimental and computational studies has changed our notion of how these enzymes exert their catalytic powers.
Collapse
Affiliation(s)
| | | | - Amnon Kohen
- Department of Chemistry, The University of Iowa, Iowa City, IA 52242, USA
| |
Collapse
|
37
|
Masterson JE, Schwartz SD. Changes in protein architecture and subpicosecond protein dynamics impact the reaction catalyzed by lactate dehydrogenase. J Phys Chem A 2013; 117:7107-13. [PMID: 23441954 DOI: 10.1021/jp400376h] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have previously established the importance of a promoting vibration, a subpicosecond protein motion that propagates through a specific axis of residues, in the reaction coordinate of lactate dehydrogenase (LDH). To test the effect that perturbation of this motion would have on the enzymatic reaction, we employ transition path sampling to obtain transition path ensembles for four independent LDH enzymatic systems: the wild type enzyme, a version of the enzyme expressing heavy isotopic substitution, and two enzymes with mutations in the promoting vibration axis. We show that even slight changes in the promoting vibration of LDH result in dramatic changes in enzymatic chemistry. In the "heavy" version of the enzyme, we find that the dampening of the subpicosecond dynamics from heavy isotopic substitution leads to a drastic increase in the time of barrier crossing. Furthermore, we see that mutation of the promoting vibration axis causes a decrease in the variability of transition paths available to the enzymatic reaction. The combined results reveal the importance of the protein architecture of LDH in enzymatic catalysis by establishing how the promoting vibration is finely tuned to facilitate chemistry.
Collapse
Affiliation(s)
- Jean E Masterson
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Blvd., Tucson, Arizona 85721, USA
| | | |
Collapse
|
38
|
Toney MD, Castro JN, Addington TA. Heavy-enzyme kinetic isotope effects on proton transfer in alanine racemase. J Am Chem Soc 2013; 135:2509-11. [PMID: 23373756 DOI: 10.1021/ja3101243] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The catalytic effects of perdeuterating the pyridoxal phosphate-dependent enzyme alanine racemase from Geobacillus stearothermophilus are reported. The mass of the heavy perdeuterated form is ~5.5% greater than that of the protiated form, causing kinetic isotope effects (KIEs) of ~1.3 on k(cat) and k(cat)/K(M) for both L- and D-alanine. These values increase when Cα-deuterated alanine is used as the substrate. The heavy-enzyme KIEs of ~3 on k(cat)/K(M) with deuterated substrates are greater than the product of the individual heavy-enzyme and primary substrate KIEs. This breakdown of the rule of the geometric mean is likely due to coupled motion between the protein and the proton-transfer reaction coordinate in the rate-limiting step. These data implicate a direct role for protein vibrational motions in barrier crossing for proton-transfer steps in alanine racemase.
Collapse
Affiliation(s)
- Michael D Toney
- Department of Chemistry, University of California-Davis, One Shields Avenue, Davis, California 95616, USA.
| | | | | |
Collapse
|
39
|
Fan Y, Cembran A, Ma S, Gao J. Connecting protein conformational dynamics with catalytic function as illustrated in dihydrofolate reductase. Biochemistry 2013; 52:2036-49. [PMID: 23297871 DOI: 10.1021/bi301559q] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Combined quantum mechanics/molecular mechanics molecular dynamics simulations reveal that the M20 loop conformational dynamics of dihydrofolate reductase (DHFR) is severely restricted at the transition state of the hydride transfer as a result of the M42W/G121V double mutation. Consequently, the double-mutant enzyme has a reduced entropy of activation, i.e., increased entropic barrier, and altered temperature dependence of kinetic isotope effects in comparison with those of wild-type DHFR. Interestingly, in both wild-type DHFR and the double mutant, the average donor-acceptor distances are essentially the same in the Michaelis complex state (~3.5 Å) and the transition state (2.7 Å). It was found that an additional hydrogen bond is formed to stabilize the M20 loop in the closed conformation in the M42W/G121V double mutant. The computational results reflect a similar aim designed to knock out precisely the dynamic flexibility of the M20 loop in a different double mutant, N23PP/S148A.
Collapse
Affiliation(s)
- Yao Fan
- Department of Chemistry, Digital Technology Center, and Supercomputing Institute, University of Minnesota , 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | | | | | | |
Collapse
|
40
|
Abstract
The relationship between protein dynamics and function is a subject of considerable contemporary interest. Although protein motions are frequently observed during ligand binding and release steps, the contribution of protein motions to the catalysis of bond making/breaking processes is more difficult to probe and verify. Here, we show how the quantum mechanical hydrogen tunneling associated with enzymatic C-H bond cleavage provides a unique window into the necessity of protein dynamics for achieving optimal catalysis. Experimental findings support a hierarchy of thermodynamically equilibrated motions that control the H-donor and -acceptor distance and active-site electrostatics, creating an ensemble of conformations suitable for H-tunneling. A possible extension of this view to methyl transfer and other catalyzed reactions is also presented. The impact of understanding these dynamics on the conceptual framework for enzyme activity, inhibitor/drug design, and biomimetic catalyst design is likely to be substantial.
Collapse
Affiliation(s)
- Judith P. Klinman
- Department of Chemistry, Department of Molecular and Cell Biology, and the California Institute for Quantitative Sciences, University of California, Berkeley, California 94720;
| | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1294;
| |
Collapse
|
41
|
Abstract
This chapter discusses progress over the past 15 years in understanding the role of protein dynamics in enzymatically catalyzed chemical reactions. Research has shown that protein motion on all timescales from femtoseconds to milliseconds can contribute to function, and in particular in some enzymes there are sub-picosecond motions, on the same timescale as barrier passage, the couple directly to chemical transformation, and are thus part of the reaction coordinate. Approaches such as transition path sampling and committor analysis have greatly enhanced our understanding of these processes.
Collapse
Affiliation(s)
- Steven D Schwartz
- Department of Chemistry, University of Arizona, 1306 East University Blvd., Tucson, AZ, 85721, USA,
| |
Collapse
|
42
|
Pollak E. Quantum variational transition state theory for hydrogen tunneling in enzyme catalysis. J Phys Chem B 2012; 116:12966-71. [PMID: 22992044 DOI: 10.1021/jp307556j] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Experiments in recent years have shown that there is a large kinetic isotope effect in the rate of transfer of hydrogen or deuterium in enzymatic reactions of soybean lipoxygenase-1. The kinetic isotope effect (KIE) is only weakly temperature dependent but varies significantly in the presence of mutants whose functional groups are located rather far from the reaction center. In this paper we suggest that variational transition state theory as applied to dissipative systems, above the crossover temperature between deep tunneling and thermal activation, may be used as a paradigm for understanding the dynamics of these reactions. We find that the theory fits the experimental data rather well. The effects of different mutants are readily interpreted in terms of the friction they exert on the reaction center. Increasing the distal functional group increases the friction and thus lowers the kinetic isotope effect.
Collapse
Affiliation(s)
- Eli Pollak
- Chemical Physics Department, Weizmann Institute of Science, 76100 Rehovoth, Israel.
| |
Collapse
|