1
|
Cristobal J, Nagorski RW, Richard JP. Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces. Biochemistry 2023; 62:2314-2324. [PMID: 37463347 PMCID: PMC10399567 DOI: 10.1021/acs.biochem.3c00290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 06/29/2023] [Indexed: 07/20/2023]
Abstract
The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.
Collapse
Affiliation(s)
- Judith
R. Cristobal
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| | - Richard W. Nagorski
- Department
of Chemistry, Illinois State University, Normal, Illinois 61790-4160, United
States
| | - John P. Richard
- Department
of Chemistry, University at Buffalo, SUNY, Buffalo, New York 14260-3000, United
States
| |
Collapse
|
2
|
Kannath S, Adamczyk P, Ferro-Costas D, Fernández-Ramos A, Major DT, Dybala-Defratyka A. Role of Microsolvation and Quantum Effects in the Accurate Prediction of Kinetic Isotope Effects: The Case of Hydrogen Atom Abstraction in Ethanol by Atomic Hydrogen in Aqueous Solution. J Chem Theory Comput 2020; 16:847-859. [PMID: 31904954 PMCID: PMC7588029 DOI: 10.1021/acs.jctc.9b00774] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Hydrogen abstraction from ethanol
by atomic hydrogen in aqueous
solution is studied using two theoretical approaches: the multipath
variational transition state theory (MP-VTST) and a path-integral
formalism in combination with free-energy perturbation and umbrella
sampling (PI-FEP/UM). The performance of the models is compared to
experimental values of H kinetic isotope effects (KIE). Solvation
models used in this study ranged from purely implicit, via mixed–microsolvation
treated quantum mechanically via the density functional theory (DFT)
to fully explicit representation of the solvent, which was incorporated
using a combined quantum mechanical-molecular mechanical (QM/MM) potential.
The effects of the transition state conformation and the position
of microsolvating water molecules interacting with the solute on the
KIE are discussed. The KIEs are in good agreement with experiment
when MP-VTST is used together with a model that includes microsolvation
of the polar part of ethanol by five or six water molecules, emphasizing
the importance of explicit solvation in KIE calculations. Both, MP-VTST
and PI-FEP/UM enable detailed characterization of nuclear quantum
effects accompanying the hydrogen atom transfer reaction in aqueous
solution.
Collapse
Affiliation(s)
- Suraj Kannath
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| | - Paweł Adamczyk
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| | - David Ferro-Costas
- LAQV@REQUIMTE, Department of Chemistry & Biochemistry, Faculty of Sciences , University of Porto , Rua do Campo Alegre , 4169-007 Porto , Portugal.,Center for Research in Biological Chemistry and Molecular Materials (CIQUS) , Universidade de Santiago de Compostela , 15782 Santiago de Compostela , Spain
| | - Antonio Fernández-Ramos
- Center for Research in Biological Chemistry and Molecular Materials (CIQUS) , Universidade de Santiago de Compostela , 15782 Santiago de Compostela , Spain
| | - Dan Thomas Major
- Department of Chemistry and Institute for Nanotechnology & Advanced Materials , Bar-Ilan University , Ramat-Gan 52900 , Israel
| | - Agnieszka Dybala-Defratyka
- Institute of Applied Radiation Chemistry, Faculty of Chemistry , Lodz University of Technology , Zeromskiego 116 , 90-924 Lodz , Poland
| |
Collapse
|
3
|
Prah A, Ogrin P, Mavri J, Stare J. Nuclear quantum effects in enzymatic reactions: simulation of the kinetic isotope effect of phenylethylamine oxidation catalyzed by monoamine oxidase A. Phys Chem Chem Phys 2020; 22:6838-6847. [DOI: 10.1039/d0cp00131g] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
By using computational techniques for quantizing nuclear motion one can accurately reproduce kinetic isotope effect of enzymatic reactions, as demonstrated for phenylethylamine oxidation catalyzed by the monoamine oxidase A enzyme.
Collapse
Affiliation(s)
- Alja Prah
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
- University of Ljubljana
| | - Peter Ogrin
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
- University of Ljubljana
| | - Janez Mavri
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
| | - Jernej Stare
- Theory Department
- National Institute of Chemistry
- Ljubljana
- Slovenia
| |
Collapse
|
4
|
Examination of the performance of semiempirical methods in QM/MM studies of the SN2-like reaction of an adenylyl group transfer catalysed by ANT4′. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2507-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
5
|
Tuñón I, Williams IH. The transition state and cognate concepts. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.09.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
6
|
Roca M, Ruiz-Pernía JJ, Castillo R, Oliva M, Moliner V. Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase. Phys Chem Chem Phys 2018; 20:25722-25737. [PMID: 30280169 DOI: 10.1039/c8cp04244f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.
Collapse
Affiliation(s)
- Maite Roca
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | | | - Raquel Castillo
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | - Mónica Oliva
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain.
| |
Collapse
|
7
|
Świderek K, Tuñón I, Williams IH, Moliner V. Insights on the Origin of Catalysis on Glycine N-Methyltransferase from Computational Modeling. J Am Chem Soc 2018; 140:4327-4334. [PMID: 29460630 PMCID: PMC6613375 DOI: 10.1021/jacs.7b13655] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor-acceptor distance (DAD): Is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The computed activation free energies and KIEs are in good agreement with experimental data, but the mutations do not meaningfully affect the DAD: Compression cannot explain the experimental variations on KIEs. On the contrary, electrostatic properties in the active site correlate with the catalytic activity of wild type and mutants. The plasticity of the enzyme moderates the effects of the mutations, explaining the rather small degree of variation in KIEs and reactivities.
Collapse
Affiliation(s)
- Katarzyna Świderek
- Departament de Química Física i Analítica; Universitat Jaume I, 12071 Castellón (Spain)
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjasot (Spain)
| | - Ian H. Williams
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| | - Vicent Moliner
- Departament de Química Física i Analítica; Universitat Jaume I, 12071 Castellón (Spain)
- Department of Chemistry, University of Bath, Bath BA2 7AY (United Kingdom)
| |
Collapse
|
8
|
Das S, Nam K, Major DT. Rapid Convergence of Energy and Free Energy Profiles with Quantum Mechanical Size in Quantum Mechanical–Molecular Mechanical Simulations of Proton Transfer in DNA. J Chem Theory Comput 2018; 14:1695-1705. [DOI: 10.1021/acs.jctc.7b00964] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Susanta Das
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Kwangho Nam
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019-0065, United States
| | - Dan Thomas Major
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
9
|
Ranasinghe C, Guo Q, Sapienza PJ, Lee AL, Quinn DM, Cheatum CM, Kohen A. Protein Mass Effects on Formate Dehydrogenase. J Am Chem Soc 2017; 139:17405-17413. [PMID: 29083897 PMCID: PMC5800309 DOI: 10.1021/jacs.7b08359] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.
Collapse
Affiliation(s)
- Chethya Ranasinghe
- Department of Chemistry, University of Iowa, Iowa City, IA 52242-1727, USA
| | - Qi Guo
- Department of Chemistry, University of Iowa, Iowa City, IA 52242-1727, USA
| | - Paul J. Sapienza
- Division of Chemical Biology and Medicinal Chemistry, UNC Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Andrew L. Lee
- Division of Chemical Biology and Medicinal Chemistry, UNC Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Daniel M. Quinn
- Department of Chemistry, University of Iowa, Iowa City, IA 52242-1727, USA
| | | | - Amnon Kohen
- Department of Chemistry, University of Iowa, Iowa City, IA 52242-1727, USA
| |
Collapse
|
10
|
Zinovjev K, Tuñón I. Quantifying the limits of transition state theory in enzymatic catalysis. Proc Natl Acad Sci U S A 2017; 114:12390-12395. [PMID: 29101125 PMCID: PMC5703300 DOI: 10.1073/pnas.1710820114] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
While being one of the most popular reaction rate theories, the applicability of transition state theory to the study of enzymatic reactions has been often challenged. The complex dynamic nature of the protein environment raised the question about the validity of the nonrecrossing hypothesis, a cornerstone in this theory. We present a computational strategy to quantify the error associated to transition state theory from the number of recrossings observed at the equicommittor, which is the best possible dividing surface. Application of a direct multidimensional transition state optimization to the hydride transfer step in human dihydrofolate reductase shows that both the participation of the protein degrees of freedom in the reaction coordinate and the error associated to the nonrecrossing hypothesis are small. Thus, the use of transition state theory, even with simplified reaction coordinates, provides a good theoretical framework for the study of enzymatic catalysis.
Collapse
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
11
|
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
|
12
|
Guo Q, Gakhar L, Wickersham K, Francis K, Vardi-Kilshtain A, Major DT, Cheatum CM, Kohen A. Structural and Kinetic Studies of Formate Dehydrogenase from Candida boidinii. Biochemistry 2016; 55:2760-71. [PMID: 27100912 DOI: 10.1021/acs.biochem.6b00181] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The structure of formate dehydrogenase from Candida boidinii (CbFDH) is of both academic and practical interests. First, this enzyme represents a unique model system for studies on the role of protein dynamics in catalysis, but so far these studies have been limited by the availability of structural information. Second, CbFDH and its mutants can be used in various industrial applications (e.g., CO2 fixation or nicotinamide recycling systems), and the lack of structural information has been a limiting factor in commercial development. Here, we report the crystallization and structural determination of both holo- and apo-CbFDH. The free-energy barrier for the catalyzed reaction was computed and indicates that this structure indeed represents a catalytically competent form of the enzyme. Complementing kinetic examinations demonstrate that the recombinant CbFDH has a well-organized reactive state. Finally, a fortuitous observation has been made: the apoenzyme crystal was obtained under cocrystallization conditions with a saturating concentration of both the cofactor (NAD(+)) and inhibitor (azide), which has a nanomolar dissociation constant. It was found that the fraction of the apoenzyme present in the solution is less than 1.7 × 10(-7) (i.e., the solution is 99.9999% holoenzyme). This is an extreme case where the crystal structure represents an insignificant fraction of the enzyme in solution, and a mechanism rationalizing this phenomenon is presented.
Collapse
Affiliation(s)
- Qi Guo
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Lokesh Gakhar
- Protein Crystallography Facility and Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Kyle Wickersham
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Kevin Francis
- Department of Chemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Alexandra Vardi-Kilshtain
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University , Ramat-Gan 5290002, Israel
| | - Dan T Major
- 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
| |
Collapse
|
13
|
Dixit M, Das S, Mhashal A, Eitan R, Major D. Practical Aspects of Multiscale Classical and Quantum Simulations of Enzyme Reactions. Methods Enzymol 2016; 577:251-86. [DOI: 10.1016/bs.mie.2016.05.046] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
14
|
Wang S, MacKay L, Lamoureux G. Development of Semiempirical Models for Proton Transfer Reactions in Water. J Chem Theory Comput 2015; 10:2881-90. [PMID: 26588263 DOI: 10.1021/ct500164h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
This letter presents a method for the parametrization of semiempirical models for proton transfer reactions in water clusters. Two new models are developed: AM1-W, which is a reparameterization of the classic AM1 model, and AM1PG-W, which is a modified AM1-like model including a pairwise correction to the core repulsion function. Both models show good performance on hydrogen-bonding energies and on proton transfer energy profiles, which are of great importance for proton transfer reactions in large water clusters and in proteins. The parametrization method introduced is general and can be used to develop any other system-specific semiempirical models.
Collapse
Affiliation(s)
- Shihao Wang
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| | - Laurent MacKay
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| | - Guillaume Lamoureux
- Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM) and ‡Department of Physics, Concordia University , Montréal, Canada
| |
Collapse
|
15
|
Vardi-Kilshtain A, Nitoker N, Major DT. Nuclear quantum effects and kinetic isotope effects in enzyme reactions. Arch Biochem Biophys 2015; 582:18-27. [DOI: 10.1016/j.abb.2015.03.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2015] [Revised: 03/02/2015] [Accepted: 03/03/2015] [Indexed: 11/28/2022]
|
16
|
Engel H, Eitan R, Azuri A, Major DT. Nuclear quantum effects in chemical reactions via higher-order path-integral calculations. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.01.001] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
17
|
Choe H, Ha JM, Joo JC, Kim H, Yoon HJ, Kim S, Son SH, Gengan RM, Jeon ST, Chang R, Jung KD, Kim YH, Lee HH. Structural insights into the efficient CO2-reducing activity of an NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA. ACTA ACUST UNITED AC 2015; 71:313-23. [PMID: 25664741 DOI: 10.1107/s1399004714025474] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2014] [Accepted: 11/20/2014] [Indexed: 11/11/2022]
Abstract
CO2 fixation is thought to be one of the key factors in mitigating global warming. Of the various methods for removing CO2, the NAD-dependent formate dehydrogenase from Candida boidinii (CbFDH) has been widely used in various biological CO2-reduction systems; however, practical applications of CbFDH have often been impeded owing to its low CO2-reducing activity. It has recently been demonstrated that the NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA (TsFDH) has a higher CO2-reducing activity compared with CbFDH. The crystal structure of TsFDH revealed that the biological unit in the asymmetric unit has two conformations, i.e. open (NAD(+)-unbound) and closed (NAD(+)-bound) forms. Three major differences are observed in the crystal structures of TsFDH and CbFDH. Firstly, hole 2 in TsFDH is blocked by helix α20, whereas it is not blocked in CbFDH. Secondly, the sizes of holes 1 and 2 are larger in TsFDH than in CbFDH. Thirdly, Lys287 in TsFDH, which is crucial for the capture of formate and its subsequent delivery to the active site, is an alanine in CbFDH. A computational simulation suggested that the higher CO2-reducing activity of TsFDH is owing to its lower free-energy barrier to CO2 reduction than in CbFDH.
Collapse
Affiliation(s)
- Hyunjun Choe
- Department of Chemical Engineering, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Jung Min Ha
- Department of Bio and Nano Chemistry, Kookmin University, Seoul 136-702, Republic of Korea
| | - Jeong Chan Joo
- Department of Chemical Engineering, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Hyunook Kim
- Department of Chemistry, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Hye-Jin Yoon
- Department of Chemistry, College of Natural Sciences, Seoul National University, Seoul 151-742, Republic of Korea
| | - Seonghoon Kim
- Department of Chemistry, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Sang Hyeon Son
- Department of Bio and Nano Chemistry, Kookmin University, Seoul 136-702, Republic of Korea
| | - Robert M Gengan
- Department of Chemistry, Faculty of Applied Sciences, Durban University of Technology, Durban, South Africa
| | - Seung Taeg Jeon
- Department of Bio and Nano Chemistry, Kookmin University, Seoul 136-702, Republic of Korea
| | - Rakwoo Chang
- Department of Chemistry, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Kwang Deog Jung
- Energy Research Center, Korea Institute of Science and Technology, Seoul, Republic of Korea
| | - Yong Hwan Kim
- Department of Chemical Engineering, Kwangwoon University, Seoul 139-701, Republic of Korea
| | - Hyung Ho Lee
- Department of Chemistry, College of Natural Sciences, Seoul National University, Seoul 151-742, Republic of Korea
| |
Collapse
|
18
|
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
|
19
|
Carvalho ATP, Barrozo A, Doron D, Kilshtain AV, Major DT, Kamerlin SCL. Challenges in computational studies of enzyme structure, function and dynamics. J Mol Graph Model 2014; 54:62-79. [PMID: 25306098 DOI: 10.1016/j.jmgm.2014.09.003] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2014] [Revised: 09/13/2014] [Accepted: 09/16/2014] [Indexed: 01/23/2023]
Abstract
In this review we give an overview of the field of Computational enzymology. We start by describing the birth of the field, with emphasis on the work of the 2013 chemistry Nobel Laureates. We then present key features of the state-of-the-art in the field, showing what theory, accompanied by experiments, has taught us so far about enzymes. We also briefly describe computational methods, such as quantum mechanics-molecular mechanics approaches, reaction coordinate treatment, and free energy simulation approaches. We finalize by discussing open questions and challenges.
Collapse
Affiliation(s)
- Alexandra T P Carvalho
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Alexandre Barrozo
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden
| | - Dvir Doron
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Alexandra Vardi Kilshtain
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry Bar-Ilan University, Ramat-Gan 52900, Israel.
| | - Shina Caroline Lynn Kamerlin
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC Box 596, S-751 24 Uppsala, Sweden.
| |
Collapse
|
20
|
Zinovjev K, Tuñón I. Exploring chemical reactivity of complex systems with path-based coordinates: Role of the distance metric. J Comput Chem 2014; 35:1672-81. [DOI: 10.1002/jcc.23673] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2014] [Revised: 05/16/2014] [Accepted: 06/18/2014] [Indexed: 01/06/2023]
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física; Universitat de València; 46100 Burjassot Spain
| | - Iñaki Tuñón
- Departament de Química Física; Universitat de València; 46100 Burjassot Spain
| |
Collapse
|
21
|
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
|
22
|
Doron D, Weitman M, Vardi-Kilshtain A, Azuri A, Engel H, Major DT. Multiscale Quantum-Classical Simulations of Enzymes. Isr J Chem 2014. [DOI: 10.1002/ijch.201400026] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
23
|
Gerber RB, Shemesh D, Varner ME, Kalinowski J, Hirshberg B. Ab initio and semi-empirical Molecular Dynamics simulations of chemical reactions in isolated molecules and in clusters. Phys Chem Chem Phys 2014; 16:9760-75. [DOI: 10.1039/c3cp55239j] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Recent progress in “on-the-fly” trajectory simulations of molecular reactions, using different electronic structure methods is discussed, with analysis of the insights that such calculations can provide and of the strengths and limitations of the algorithms available.
Collapse
Affiliation(s)
- R. B. Gerber
- Institute of Chemistry and The Fritz Haber Research Center
- The Hebrew University of Jerusalem
- Jerusalem 91904, Israel
- Department of Chemistry
- University of California
| | - D. Shemesh
- Institute of Chemistry and The Fritz Haber Research Center
- The Hebrew University of Jerusalem
- Jerusalem 91904, Israel
| | - M. E. Varner
- Department of Chemistry
- University of California
- Irvine 92697, USA
| | - J. Kalinowski
- Department of Chemistry
- University of Helsinki
- , Finland
| | - B. Hirshberg
- Institute of Chemistry and The Fritz Haber Research Center
- The Hebrew University of Jerusalem
- Jerusalem 91904, Israel
| |
Collapse
|
24
|
Liang S, Roitberg AE. AM1 Specific Reaction Parameters for Reactions of Hydroxide Ion with Halomethanes in Complex Environments: Development and Testing. J Chem Theory Comput 2013; 9:4470-80. [DOI: 10.1021/ct400471m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shuai Liang
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| | - Adrian E. Roitberg
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|