1
|
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
|
2
|
Lu Y, Wilhelm S, Bai M, Maness P, Ma L. 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? Biochemistry 2019; 58:4035-4046. [PMID: 31478638 DOI: 10.1021/acs.biochem.9b00574] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The change from the temperature independence of the primary (1°) H/D kinetic isotope effects (KIEs) in wild-type enzyme-catalyzed H-transfer reactions (ΔEa = EaD - EaH ∼ 0) to a strong temperature dependence with the mutated enzymes (ΔEa ≫ 0) has recently been frequently observed. This has prompted some enzymologists to develop new H-tunneling models to correlate ΔEa with the donor-acceptor distance (DAD) at the tunneling-ready state (TRS) as well as the protein thermal motions/dynamics that sample the short DADTRS'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 ΔEa'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 ΔEa'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, DADPRC's, and ΔEa's were discussed. The more rigid system with more suppressed 2° C-H vibrations at the TRS and more narrowly distributed DADPRC's in PRCs gave a smaller ΔEa. The results replicated the trend observed in enzymes versus mutated enzymes and appeared to support the concepts of different thermally activated DADTRS sampling processes in response to the rigid versus flexible donor-acceptor centers.
Collapse
Affiliation(s)
- Yun Lu
- Department of Chemistry , Southern Illinois University Edwardsville , Edwardsville , Illinois 62026 , United States
| | - Samantha Wilhelm
- Department of Chemistry , Southern Illinois University Edwardsville , Edwardsville , Illinois 62026 , United States
| | - Mingxuan Bai
- Department of Chemistry , Southern Illinois University Edwardsville , Edwardsville , Illinois 62026 , United States
| | - Peter Maness
- Department of Chemistry , Southern Illinois University Edwardsville , Edwardsville , Illinois 62026 , United States
| | - Li Ma
- Department of Chemistry , Southern Illinois University Edwardsville , Edwardsville , Illinois 62026 , United States
| |
Collapse
|
3
|
Mhashal AR, Pshetitsky Y, Cheatum CM, Kohen A, Major DT. Evolutionary Effects on Bound Substrate pKa in Dihydrofolate Reductase. J Am Chem Soc 2018; 140:16650-16660. [DOI: 10.1021/jacs.8b09089] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Anil R. Mhashal
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Yaron Pshetitsky
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | | | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - Dan Thomas Major
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
4
|
Mhashal AR, Pshetitsky Y, Eitan R, Cheatum CM, Kohen A, Major DT. Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility. J Phys Chem B 2018; 122:8006-8017. [PMID: 30040418 DOI: 10.1021/acs.jpcb.8b05556] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.
Collapse
Affiliation(s)
- Anil R Mhashal
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - 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
| | - Christopher M Cheatum
- Department of Chemistry , University of Iowa , Iowa City , Iowa 52242 , United States
| | - 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
|
5
|
Boehr DD, D'Amico RN, O'Rourke KF. Engineered control of enzyme structural dynamics and function. Protein Sci 2018; 27:825-838. [PMID: 29380452 DOI: 10.1002/pro.3379] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2017] [Revised: 01/20/2018] [Accepted: 01/24/2018] [Indexed: 12/20/2022]
Abstract
Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.
Collapse
Affiliation(s)
- David D Boehr
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania, 16802, USA
| | - Rebecca N D'Amico
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania, 16802, USA
| | - Kathleen F O'Rourke
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania, 16802, USA
| |
Collapse
|
6
|
Mhashal AR, Vardi-Kilshtain A, Kohen A, Major DT. The role of the Met 20 loop in the hydride transfer in Escherichia coli dihydrofolate reductase. J Biol Chem 2017; 292:14229-14239. [PMID: 28620051 DOI: 10.1074/jbc.m117.777136] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Revised: 05/24/2017] [Indexed: 11/06/2022] Open
Abstract
A key question concerning the catalytic cycle of Escherichia coli dihydrofolate reductase (ecDHFR) is whether the Met20 loop is dynamically coupled to the chemical step during catalysis. A more basic, yet unanswered question is whether the Met20 loop adopts a closed conformation during the chemical hydride transfer step. To examine the most likely conformation of the Met20 loop during the chemical step, we studied the hydride transfer in wild type (WT) ecDHFR using hybrid quantum mechanics-molecular mechanics free energy simulations with the Met20 loop in a closed and disordered conformation. Additionally, we investigated three mutant forms (I14X; X = Val, Ala, Gly) of the enzyme that have increased active site flexibility and donor-acceptor distance dynamics in closed and disordered Met20 loop states. We found that the conformation of the Met20 loop has a dramatic effect on the ordering of active site hydration, although the Met20 loop conformation only has a moderate effect on the hydride transfer rate and donor-acceptor distance dynamics. Finally, we evaluated the pKa of the substrate N5 position in closed and disordered Met20 loop states and found a strong correlation between N5 basicity and the conformation of the Met20 loop.
Collapse
Affiliation(s)
- Anil R Mhashal
- From the Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Alexandra Vardi-Kilshtain
- From the 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
| | - Dan Thomas Major
- From the Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel.
| |
Collapse
|
7
|
Effects of isoleucine 135 side chain length on the cofactor donor-acceptor distance within F 420H 2:NADP + oxidoreductase: A kinetic analysis. Biochem Biophys Rep 2016; 9:114-120. [PMID: 28955995 PMCID: PMC5614548 DOI: 10.1016/j.bbrep.2016.11.012] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2016] [Accepted: 11/22/2016] [Indexed: 10/29/2022] Open
Abstract
F420H2:NADP+ Oxidoreductase (Fno) catalyzes the reversible reduction of NADP+ to NADPH by transferring a hydride from the reduced F420 cofactor. Here, we have employed binding studies, steady-state and pre steady-state kinetic methods upon wtFno and isoleucine 135 (I135) Fno variants in order to study the effects of side chain length on the donor-acceptor distance between NADP+ and the F420 precursor, FO. The conserved I135 residue of Fno was converted to a valine, alanine and glycine, thereby shortening the side chain length. The steady-state kinetic analysis of wtFno and the variants showed classic Michaelis-Menten kinetics with varying FO concentrations. The data revealed a decreased kcat as side chain length decreased, with varying FO concentrations. The steady-state plots revealed non-Michaelis-Menten kinetic behavior when NADPH was varied. The double reciprocal plot of the varying NADPH concentrations displays a downward concave shape, while the NADPH binding curves gave Hill coefficients of less than 1. These data suggest that negative cooperativity occurs between the two identical monomers. The pre steady-state Abs420 versus time trace revealed biphasic kinetics, with a fast phase (hydride transfer) and a slow phase. The fast phase displayed an increased rate constant as side chain length decreased. The rate constant for the second phase, remained ~2 s-1 for each variant. Our data suggest that I135 plays a key role in sustaining the donor-acceptor distance between the two cofactors, thereby regulating the rate at which the hydride is transferred from FOH2 to NADP+. Therefore, Fno is a dynamic enzyme that regulates NADPH production.
Collapse
Key Words
- Dissociation constants
- E. coli,, Escherichia coli
- F420 cofactor
- F420 cofactor, 7,8-didemethyl-8-hydroxy-5-deazariboflavin-5′-phosphoryllactyl(glutamyl)nglutamate, A. fulgidus, Archaeoglobus fulgidus
- F420H2: NADP+ oxidoreductase
- FO, precursor of F420 cofactor
- Fno, F420H2:NADP+, oxidoreductase
- Half-site reactivity
- I135, Isoleucine 135
- IPTG, isopropyl β-D-1-thiogalactopyranoside
- Kd,, dissociation constant
- Km, Michaelis-Menten constant
- LB, Luria Bertani broth
- NADP
- NADP+, nicotinamide adenine dinucleotide phosphate
- Negative cooperativity
- PEI, Polyethyleneimine
- Pre steady-state kinetics
- Steady-state kinetics
- k, rate constant
- kcat, catalytic rate constant (turnover number), kcat /Km, catalytic efficiency
- wtFno, wild-type Fno
Collapse
|
8
|
Joseph E, Le CQ, Nguyen T, Oyugi M, Hossain MS, Foss FW, Johnson-Winters K. Evidence of Negative Cooperativity and Half-Site Reactivity within an F420-Dependent Enzyme: Kinetic Analysis of F420H2:NADP(+) Oxidoreductase. Biochemistry 2016; 55:1082-90. [PMID: 26811861 DOI: 10.1021/acs.biochem.5b00762] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Here, we report the very first example of half-site reactivity and negative cooperativity involving an important F420 cofactor-dependent enzyme. F420H2:NADP(+) oxidoreductase (Fno) is an F420 cofactor-dependent enzyme that catalyzes the reversible reduction of NADP(+) through the transfer of a hydride from the reduced F420 cofactor. These catalytic processes are of major significance in numerous biochemical processes. While the steady-state kinetic analysis showed classic Michaelis-Menten kinetics with varying concentrations of the F420 redox moiety, FO, such plots revealed non-Michaelis-Menten kinetic behavior when NADPH was varied. The double reciprocal plot of the varying concentrations of NADPH displays a downward concave shape, suggesting that negative cooperativity occurs between the two identical monomers. The transient state kinetic data show a burst prior to entering steady-state turnover. The burst suggests that product release is rate-limiting, and the amplitude of the burst phase corresponds to production of product in only one of the active sites of the functional dimer. These results suggest either half-site reactivity or an alternate sites model wherein the reduction of the cofactor, FO occurs at one active site at a time followed by reduction at the second active site. Thus, the data imply that Fno may be a functional regulatory enzyme.
Collapse
Affiliation(s)
- Ebenezer Joseph
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Cuong Quang Le
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Toan Nguyen
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Mercy Oyugi
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Mohammad Shawkat Hossain
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Frank W Foss
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| | - Kayunta Johnson-Winters
- Department of Chemistry and Biochemistry, University of Texas at Arlington , Arlington, Texas 76019-0065, United States
| |
Collapse
|
9
|
Luk LYP, Ruiz-Pernía JJ, Adesina AS, Loveridge EJ, Tuñón I, Moliner V, Allemann RK. Chemical Ligation and Isotope Labeling to Locate Dynamic Effects during Catalysis by Dihydrofolate Reductase. Angew Chem Int Ed Engl 2015; 54:9016-20. [PMID: 26079622 PMCID: PMC4985705 DOI: 10.1002/anie.201503968] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2015] [Indexed: 11/21/2022]
Abstract
Chemical ligation has been used to alter motions in specific regions of dihydrofolate reductase from E. coli and to investigate the effects of localized motional changes on enzyme catalysis. Two isotopic hybrids were prepared; one with the mobile N-terminal segment containing heavy isotopes ((2) H, (13) C, (15) N) and the remainder of the protein with natural isotopic abundance, and the other one with only the C-terminal segment isotopically labeled. Kinetic investigations indicated that isotopic substitution of the N-terminal segment affected only a physical step of catalysis, whereas the enzyme chemistry was affected by protein motions from the C-terminal segment. QM/MM studies support the idea that dynamic effects on catalysis mostly originate from the C-terminal segment. The use of isotope hybrids provides insights into the microscopic mechanism of dynamic coupling, which is difficult to obtain with other studies, and helps define the dynamic networks of intramolecular interactions central to enzyme catalysis.
Collapse
Affiliation(s)
- Louis Y P Luk
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT (UK)
| | - J Javier Ruiz-Pernía
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló (Spain)
| | | | - E Joel Loveridge
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT (UK)
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot (Spain).
| | - Vincent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló (Spain).
| | - Rudolf K Allemann
- School of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT (UK).
| |
Collapse
|
10
|
Luk LYP, Ruiz-Pernía JJ, Adesina AS, Loveridge EJ, Tuñón I, Moliner V, Allemann RK. Chemical Ligation and Isotope Labeling to Locate Dynamic Effects during Catalysis by Dihydrofolate Reductase. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201503968] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
11
|
Maharjan B, Raghibi Boroujeni M, Lefton J, White OR, Razzaghi M, Hammann BA, Derakhshani-Molayousefi M, Eilers JE, Lu Y. Steric Effects on the Primary Isotope Dependence of Secondary Kinetic Isotope Effects in Hydride Transfer Reactions in Solution: Caused by the Isotopically Different Tunneling Ready State Conformations? J Am Chem Soc 2015; 137:6653-61. [DOI: 10.1021/jacs.5b03085] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Binita Maharjan
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Mahdi Raghibi Boroujeni
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Jonathan Lefton
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Ormacinda R. White
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Mortezaali Razzaghi
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Blake A. Hammann
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | | | - James E. Eilers
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| | - Yun Lu
- Department of Chemistry, Southern Illinois University Edwardsville, Edwardsville, Illinois 62026, United States
| |
Collapse
|
12
|
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
|
13
|
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
|
14
|
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
|
15
|
Francis K, Kohen A. Protein motions and the activation of the CH bond catalyzed by dihydrofolate reductase. Curr Opin Chem Biol 2014; 21:19-24. [PMID: 24742825 PMCID: PMC4149937 DOI: 10.1016/j.cbpa.2014.03.009] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2014] [Revised: 03/03/2014] [Accepted: 03/05/2014] [Indexed: 11/28/2022]
Abstract
The role of protein motions in enzymatic CH→C transfer is an area of great contemporary debate. An effective tool in probing such a role is the temperature dependence of the intrinsic kinetic isotope effects for the enzyme-catalyzed reaction. The outcome of those experiments is interpreted within the context of phenomenological Marcus-like models of hydrogen tunneling. The current review focuses on recent studies of dihydrofolate reductase (DHFR) and how the role of protein motions in the catalyzed reaction has been demonstrated. The motions in DHFR are controlled by local effects of active site residues, global effects involving remote residues across the enzyme and appear to be preserved during the evolution of the enzyme from bacteria to human.
Collapse
Affiliation(s)
- Kevin Francis
- Department of Chemistry, The University of Iowa, Iowa City, IA 52242, United States
| | - Amnon Kohen
- Department of Chemistry, The University of Iowa, Iowa City, IA 52242, United States.
| |
Collapse
|
16
|
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
|
17
|
Francis K, Kohen A. Standards for the reporting of kinetic isotope effects in enzymology. ACTA ACUST UNITED AC 2014. [DOI: 10.1016/j.pisc.2014.02.009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
18
|
Guo J, Loveridge EJ, Luk LYP, Allemann RK. Effect of Dimerization on Dihydrofolate Reductase Catalysis. Biochemistry 2013; 52:3881-7. [DOI: 10.1021/bi4005073] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Jiannan Guo
- School of Chemistry and
Cardiff Catalysis Institute, Cardiff University, Main Building, Park Place, Cardiff
CF10 3AT, United Kingdom
| | - E. Joel Loveridge
- School of Chemistry and
Cardiff Catalysis Institute, Cardiff University, Main Building, Park Place, Cardiff
CF10 3AT, United Kingdom
| | - Louis Y. P. Luk
- School of Chemistry and
Cardiff Catalysis Institute, Cardiff University, Main Building, Park Place, Cardiff
CF10 3AT, United Kingdom
| | - Rudolf K. Allemann
- School of Chemistry and
Cardiff Catalysis Institute, Cardiff University, Main Building, Park Place, Cardiff
CF10 3AT, United Kingdom
| |
Collapse
|
19
|
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
|
20
|
Roston D, Cheatum CM, Kohen A. Hydrogen donor-acceptor fluctuations from kinetic isotope effects: a phenomenological model. Biochemistry 2012; 51:6860-70. [PMID: 22857146 PMCID: PMC3448806 DOI: 10.1021/bi300613e] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Kinetic isotope effects (KIEs) and their temperature dependence can probe the structural and dynamic nature of enzyme-catalyzed proton or hydride transfers. The molecular interpretation of their temperature dependence requires expensive and specialized quantum mechanics/molecular mechanics (QM/MM) calculations to provide a quantitative molecular understanding. Currently available phenomenological models use a nonadiabatic assumption that is not appropriate for most hydride and proton-transfer reactions, while others require more parameters than the experimental data justify. Here we propose a phenomenological interpretation of KIEs based on a simple method to quantitatively link the size and temperature dependence of KIEs to a conformational distribution of the catalyzed reaction. This model assumes adiabatic hydrogen tunneling, and by fitting experimental KIE data, the model yields a population distribution for fluctuations of the distance between donor and acceptor atoms. Fits to data from a variety of proton and hydride transfers catalyzed by enzymes and their mutants, as well as nonenzymatic reactions, reveal that steeply temperature-dependent KIEs indicate the presence of at least two distinct conformational populations, each with different kinetic behaviors. We present the results of these calculations for several published cases and discuss how the predictions of the calculations might be experimentally tested. This analysis does not replace molecular QM/MM investigations, but it provides a fast and accessible way to quantitatively interpret KIEs in the context of a Marcus-like model.
Collapse
Affiliation(s)
- Daniel Roston
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| | | | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA 52242
| |
Collapse
|
21
|
Doron D, Kohen A, Major DT. Collective Reaction Coordinate for Hybrid Quantum and Molecular Mechanics Simulations: A Case Study of the Hydride Transfer in Dihydrofolate Reductase. J Chem Theory Comput 2012; 8:2484-96. [PMID: 26588977 DOI: 10.1021/ct300235k] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The optimal description of the reaction coordinate in chemical systems is of great importance in simulating condensed phase reactions. In the current work, we present a collective reaction coordinate which is composed of several geometric coordinates which represent structural progress during the course of a hydride transfer reaction: the antisymmetric reactive stretch coordinate, the donor-acceptor distance (DAD) coordinate, and an orbital rehybridization coordinate. In this approach, the former coordinate serves as a distinguished reaction coordinate, while the latter two serve as environmental, Marcus-type inner-sphere reorganization coordinates. The classical free energy surface is obtained from multidimensional quantum mechanics-molecular mechanics (QM/MM) potential of mean force (PMF) simulations in conjunction with a general and efficient multidimensional weighted histogram method implementation. The minimum free energy path, or the collective reaction coordinate, connecting the dividing hypersurface to reactants and products, is obtained using an iterative scheme. In this approach, the string method is used to find the minimum free energy path. This path guides the multidimensional sampling, while the path is adaptively refined until convergence is achieved. As a model system, we choose the hydride transfer reaction in Escherichia coli dihydrofolate reductase (ecDHFR) using a recently developed accurate semiempirical potential energy surface. To estimate the advantages of the collective reaction coordinate, we perform activated dynamics simulations to obtain the reaction transmission coefficient. The results show that the combination of a distinguished reaction coordinate and an inner-sphere reorganization coordinate considerably reduces the dividing surface recrossing.
Collapse
Affiliation(s)
- Dvir Doron
- 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
|
22
|
Engel H, Doron D, Kohen A, Major DT. Momentum Distribution as a Fingerprint of Quantum Delocalization in Enzymatic Reactions: Open-Chain Path-Integral Simulations of Model Systems and the Hydride Transfer in Dihydrofolate Reductase. J Chem Theory Comput 2012; 8:1223-34. [DOI: 10.1021/ct200874q] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hamutal Engel
- Department of Chemistry and
the Lise Meitner−Minerva Center of Computational Quantum Chemistry,
Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dvir Doron
- 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
|
23
|
Stojković V, Perissinotti LL, Willmer D, Benkovic SJ, Kohen A. Effects of the donor-acceptor distance and dynamics on hydride tunneling in the dihydrofolate reductase catalyzed reaction. J Am Chem Soc 2012; 134:1738-45. [PMID: 22171795 DOI: 10.1021/ja209425w] [Citation(s) in RCA: 92] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A significant contemporary question in enzymology involves the role of protein dynamics and hydrogen tunneling in enhancing enzyme catalyzed reactions. Here, we report a correlation between the donor-acceptor distance (DAD) distribution and intrinsic kinetic isotope effects (KIEs) for the dihydrofolate reductase (DHFR) catalyzed reaction. This study compares the nature of the hydride-transfer step for a series of active-site mutants, where the size of a side chain that modulates the DAD (I14 in E. coli DHFR) is systematically reduced (I14V, I14A, and I14G). The contributions of the DAD and its dynamics to the hydride-transfer step were examined by the temperature dependence of intrinsic KIEs, hydride-transfer rates, activation parameters, and classical molecular dynamics (MD) simulations. Results are interpreted within the framework of the Marcus-like model where the increase in the temperature dependence of KIEs arises as a direct consequence of the deviation of the DAD from its distribution in the wild type enzyme. Classical MD simulations suggest new populations with larger average DADs, as well as broader distributions, and a reduction in the population of the reactive conformers correlated with the decrease in the size of the hydrophobic residue. The more flexible active site in the mutants required more substantial thermally activated motions for effective H-tunneling, consistent with the hypothesis that the role of the hydrophobic side chain of I14 is to restrict the distribution and dynamics of the DAD and thus assist the hydride-transfer. These studies establish relationships between the distribution of DADs, the hydride-transfer rates, and the DAD's rearrangement toward tunneling-ready states. This structure-function correlation shall assist in the interpretation of the temperature dependence of KIEs caused by mutants far from the active site in this and other enzymes, and may apply generally to C-H→C transfer reactions.
Collapse
Affiliation(s)
- Vanja Stojković
- Department of Chemistry, The University of Iowa, Iowa City, Iowa 52242, USA
| | | | | | | | | |
Collapse
|
24
|
Lee J, Goodey NM. Catalytic contributions from remote regions of enzyme structure. Chem Rev 2011; 111:7595-624. [PMID: 21923192 DOI: 10.1021/cr100042n] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Jeeyeon Lee
- Department of Chemistry, 413 Wartik Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802, USA.
| | | |
Collapse
|
25
|
Brut M, Estève A, Landa G, Renvez G, Djafari Rouhani M, Vaisset M. Atomic Scale Determination of Enzyme Flexibility and Active Site Stability through Static Modes: Case of Dihydrofolate Reductase. J Phys Chem B 2011; 115:1616-22. [DOI: 10.1021/jp109874z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marie Brut
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California, United States
| | - Alain Estève
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Georges Landa
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Guillaume Renvez
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Mehdi Djafari Rouhani
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| | - Marc Vaisset
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France
- Université de Toulouse, UPS, INSA, INP, ISAE, LAAS, F-31077 Toulouse, France
| |
Collapse
|