1
|
Chen BW. Equilibrium and kinetic isotope effects in heterogeneous catalysis: A density functional theory perspective. CATAL COMMUN 2023. [DOI: 10.1016/j.catcom.2023.106654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023] Open
|
2
|
Giese TJ, Zeng J, Ekesan Ş, York DM. Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions. J Chem Theory Comput 2022; 18:4304-4317. [PMID: 35709391 PMCID: PMC9283286 DOI: 10.1021/acs.jctc.2c00151] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a fast, accurate, and robust approach for determination of free energy profiles and kinetic isotope effects for RNA 2'-O-transphosphorylation reactions with inclusion of nuclear quantum effects. We apply a deep potential range correction (DPRc) for combined quantum mechanical/molecular mechanical (QM/MM) simulations of reactions in the condensed phase. The method uses the second-order density-functional tight-binding method (DFTB2) as a fast, approximate base QM model. The DPRc model modifies the DFTB2 QM interactions and applies short-range corrections to the QM/MM interactions to reproduce ab initio DFT (PBE0/6-31G*) QM/MM energies and forces. The DPRc thus enables both QM and QM/MM interactions to be tuned to high accuracy, and the QM/MM corrections are designed to smoothly vanish at a specified cutoff boundary (6 Å in the present work). The computational speed-up afforded by the QM/MM+DPRc model enables free energy profiles to be calculated that include rigorous long-range QM/MM interactions under periodic boundary conditions and nuclear quantum effects through a path integral approach using a new interface between the AMBER and i-PI software. The approach is demonstrated through the calculation of free energy profiles of a native RNA cleavage model reaction and reactions involving thio-substitutions, which are important experimental probes of the mechanism. The DFTB2+DPRc QM/MM free energy surfaces agree very closely with the PBE0/6-31G* QM/MM results, and it is vastly superior to the DFTB2 QM/MM surfaces with and without weighted thermodynamic perturbation corrections. 18O and 34S primary kinetic isotope effects are compared, and the influence of nuclear quantum effects on the free energy profiles is examined.
Collapse
Affiliation(s)
- Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Jinzhe Zeng
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
3
|
Passos MO, Lins IA, Venâncio MF, Alves TV. Differences in the torsional anharmonicity between reactant and transition state: the case of 3-butenal + H abstraction reactions. Phys Chem Chem Phys 2021; 23:25414-25423. [PMID: 34751697 DOI: 10.1039/d1cp03981d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients for the hydrogen-abstraction reactions of 3-butenal by a hydrogen atom were obtained applying multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). Torsional anharmonicity due to the hindered rotors was taken into account by calculating the rovibrational partition function using the extended two-dimensional torsional (E2DT) method. For comparison, rovibrational partition functions were also estimated using the multistructural method with torsional anharmonicity based on a coupled torsional potential (MS-T(C)). By contrast, with (E)-2-butenal reactions, the abstraction reactions of 3-butenal proceed via five reaction channels (R1)-(R5). In a conformational search, 45 distinguishable structures of transition states were found, including enantiomers, which were separated into six conformational reaction channels (CRCs). The individual reactive paths were constructed, the recrossing and semiclassical transmission coefficients estimated, and the multipath rate constants were obtained. High torsional barriers between the wells of CRC2/CRC6 indicate a harmonic behavior. Consequently, a difference between the torsional anharmonicity of 3-butenal and the transition states is responsible for the increase in the thermal rate constants for channel (R2). Analysis of the contributions of each conformer of the transition state shows an important contribution of the high-energy rotamers in the total flux of (R1)-(R5). After fitting the rate constants in a four-parameter equation, the activation energy estimation showed a strong temperature dependence.
Collapse
Affiliation(s)
- Maiara Oliveira Passos
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| | - Igor Araujo Lins
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| | - Mateus Fernandes Venâncio
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| | - Tiago Vinicius Alves
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| |
Collapse
|
4
|
Passos MO, Lins IA, Alves TV. Rate coefficients and product branching ratios for (E)-2-butenal + H reactions. Phys Chem Chem Phys 2020; 22:14246-14254. [PMID: 32555895 DOI: 10.1039/d0cp02142c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate constants for the hydrogen abstraction reactions of (E)-2-butenal by hydrogen atoms were calculated, for the first time, using the multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). After a torsional potential energy surface exploration, ten conformations of the transition states (including the mirror images) were found and separated into four conformational reaction channels (CRCs). Individual energy paths of each CRC were built, recrossing and quantum tunneling effects estimated, and the thermal rate constants obtained. Due to the hindered rotors, the torsional anharmonicity was incorporated in the rate coefficient through the calculations of the rovibrational partition functions using the extended two-dimensional torsional method (E2DT). For comparison, the one-well (1W-CVT/SCT) and harmonic multipath (MP-CVT/SCT) thermal rate constants were also estimated. In addition, kinetic Monte Carlo (KMC) simulations were performed to predict the product branching ratios. For all kinetic approaches, the formation of products of (R1) is predominant. Compared to the harmonic multipath estimation, the percentage of reaction (R4) increases by approximately 9% when the torsional anharmonicity is taken into account. For the reactions (R2) and (R3), the product branching ratio is slightly decreased when compared with the harmonic simulation.
Collapse
Affiliation(s)
- Maiara Oliveira Passos
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| | - Igor Araujo Lins
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| | - Tiago Vinicius Alves
- Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.
| |
Collapse
|
5
|
High-temperature rate constants for CH3OCOH + OH reactions: the effects of multiple structures and paths. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-02599-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
6
|
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
|
7
|
Ferro-Costas D, Martínez-Núñez E, Rodríguez-Otero J, Cabaleiro-Lago E, Estévez CM, Fernández B, Fernández-Ramos A, Vázquez SA. Influence of Multiple Conformations and Paths on Rate Constants and Product Branching Ratios. Thermal Decomposition of 1-Propanol Radicals. J Phys Chem A 2018; 122:4790-4800. [DOI: 10.1021/acs.jpca.8b02949] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
8
|
Bao JL, Xing L, Truhlar DG. Dual-Level Method for Estimating Multistructural Partition Functions with Torsional Anharmonicity. J Chem Theory Comput 2017; 13:2511-2522. [DOI: 10.1021/acs.jctc.7b00232] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Lili Xing
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|