1
|
Shin W, Ran X, Yang ZJ. Accelerated Entropic Path Sampling with a Bidirectional Generative Adversarial Network. J Phys Chem B 2023; 127:4254-4260. [PMID: 37133810 DOI: 10.1021/acs.jpcb.3c01202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
The role of entropy in mediating the dynamic outcomes of chemical reactions remains largely unknown. To evaluate the change of entropy along post-transition state paths, we have previously developed entropic path sampling that computes configurational entropy from an ensemble of reaction trajectories. However, one major caveat of this approach lies in its high computational demand: about 2000 trajectories are needed to converge the computation of an entropic profile. Here, by leveraging a deep generative model, we developed an accelerated entropic path sampling approach that evaluates entropic profiles using merely a few hundred reaction dynamic trajectories. The new method, called bidirectional generative adversarial network-entropic path sampling, can enhance the estimation of probability density functions of molecular configurations by generating pseudo-molecular configurations that are statistically indistinguishable from the true data. The method was established using cyclopentadiene dimerization, in which we reproduced the reference entropic profiles (derived from 2480 trajectories) using merely 124 trajectories. The method was further benchmarked using three reactions with symmetric post-transition-state bifurcation, including endo-butadiene dimerization, 5-fluoro-1,3-cyclopentadiene dimerization, and 5-methyl-1,3-cyclopentadiene dimerization. The results indicate the existence of a "hidden entropic intermediate", which is a dynamic species that binds to a local entropic maximum where no free energy minimum is formed.
Collapse
Affiliation(s)
- Wook Shin
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Xinchun Ran
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Zhongyue J Yang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Vanderbilt Institute of Chemical Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, United States
- Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
| |
Collapse
|
2
|
Tao Z, Qiu T, Subotnik JE. Symmetric Post-Transition State Bifurcation Reactions with Berry Pseudomagnetic Fields. J Phys Chem Lett 2023; 14:770-778. [PMID: 36652556 DOI: 10.1021/acs.jpclett.2c02668] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We investigate how the Berry force (i.e., the pseudomagnetic force operating on nuclei as induced by electronic degeneracy and spin-orbit coupling (SOC)) might modify a post-transition state bifurcation (PTSB) reaction path and affect product selectivity for situations when multiple products share the same transition state. To estimate the magnitude of this effect, Langevin dynamics are performed on a model system with a valley-ridge inflection (VRI) point in the presence of a magnetic field (that mimics the Berry curvature). We also develop an analytic model for such selectivity that depends on key parameters such as the surface topology, the magnitude of the Berry force, and the nuclear friction. Within this dynamical model, static electronic structure calculations (at the level of generalized Hartree-Fock with spin-orbit coupling (GHF+SOC) theory) suggest that electronic spin induced Berry force effects may indeed lead to noticeable changes in methoxy radical isomerization.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
Carpenter BK. Prediction of Kinetic Product Ratios: Investigation of a Dynamically Controlled Case. J Phys Chem A 2023; 127:224-239. [PMID: 36594780 PMCID: PMC9841574 DOI: 10.1021/acs.jpca.2c08301] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Of the various factors influencing kinetically controlled product ratios, the role of nonstatistical dynamics is arguably the least well understood. In this paper, reactions were chosen in which dynamics played a dominant role in product selection, by design. Specifically, the reactions studied were the ring openings of cyclopropylidene to allene and tetramethylcyclopropylidene to tetramethylallene (2,4-dimethylpenta-2,3-diene). Both reactions have intrinsic reaction coordinates that bifurcate symmetrically, leading to products that are enantiomeric once the atoms are uniquely labeled. The question addressed in the study was whether the outcomes─that is, which product well on the potential energy surface was selected─could be predicted from their initial conditions for individual trajectories in quasiclassical dynamics simulations. Hybrid potentials were developed based on cooperative interaction between molecular mechanics and artificial neural networks, trained against data from electronic structure calculations. These potentials allowed simulations of both gas-phase and condensed-phase reactions. The outcome was that, for both reactions, prediction of initial selection of product wells could be made with >95% success from initial conditions of the trajectories in the gas phase. However, when trajectories were run for longer, looking for "final" products for each trajectory, the predictability dropped off dramatically. In the gas-phase simulations, this drop off was caused by trajectories hopping between product wells on the potential energy surface. That behavior could be suppressed in condensed phases, but then new uncertainty was introduced because the intermolecular interactions between solute and bath, necessary to permit intermolecular energy transfer and cooling of the hot initial products, often led to perturbations of the initial directions of trajectories on the potential energy surface. It would consequently appear that a general ability to predict outcomes for reactions in which nonstatistical dynamics dominate remains a challenge even in the age of sophisticated machine-learning capabilities.
Collapse
|
4
|
Loco D, Chataigner I, Piquemal J, Spezia R. Efficient and Accurate Description of Diels-Alder Reactions Using Density Functional Theory. Chemphyschem 2022; 23:e202200349. [PMID: 35696652 PMCID: PMC9796631 DOI: 10.1002/cphc.202200349] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 06/11/2022] [Indexed: 01/01/2023]
Abstract
Modeling chemical reactions using Quantum Chemistry is a widely used predictive strategy capable to complement experiments in order to understand the intrinsic mechanisms guiding the chemicals towards the most favorable reaction products. However, at this purpose, it is mandatory to use reliable and computationally tractable theoretical methods. In this work, we focus on six Diels-Alder reactions of increasing complexity and perform an extensive benchmark of middle- to low-cost computational approaches to predict the characteristic reactions energy barriers. We found that Density Functional Theory, using the ωB97XD, LC-ωPBE, CAM-B3LYP, M11 and MN12SX functionals, with empirical dispersion corrections coupled to an affordable 6-31G basis set, provides quality results for this class of reactions, at a small computational effort. Such efficient and reliable simulation protocol opens perspectives for hybrid QM/MM molecular dynamics simulations of Diels-Alder reactions including explicit solvation.
Collapse
Affiliation(s)
- Daniele Loco
- Sorbonne Université, Laboratoire de Chimie ThéoriqueUMR 7616 CNRS4 Place Jussieu75005ParisFrance
- Qubit PharmaceuticalsIncubateur Paris Biotech Santé24 rue du Faubourg Saint Jacques75014ParisFrance
| | - Isabelle Chataigner
- Sorbonne Université, Laboratoire de Chimie ThéoriqueUMR 7616 CNRS4 Place Jussieu75005ParisFrance
- Normandie Univ.UNIROUENCNRS, INSA Rouen, COBRA76000RouenFrance
| | - Jean‐Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie ThéoriqueUMR 7616 CNRS4 Place Jussieu75005ParisFrance
| | - Riccardo Spezia
- Sorbonne Université, Laboratoire de Chimie ThéoriqueUMR 7616 CNRS4 Place Jussieu75005ParisFrance
| |
Collapse
|
5
|
Meuwly M. Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
6
|
Tremblay MT, Yang ZJ. The effect of zero‐point energy in simulating organic reactions with post‐transition state bifurcation. J PHYS ORG CHEM 2022. [DOI: 10.1002/poc.4322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Zhongyue J. Yang
- Department of Chemistry Vanderbilt University Nashville Tennessee USA
- Center for Structural Biology Vanderbilt University Nashville Tennessee USA
- Vanderbilt Institute of Chemical Biology Vanderbilt University Nashville Tennessee USA
- Data Science Institute Vanderbilt University Nashville Tennessee USA
| |
Collapse
|
7
|
Zhang X, Lefebvre PL, Harvey JN. Effect of solvent motions on the dynamics of the Diels-Alder reaction. Phys Chem Chem Phys 2021; 24:1120-1130. [PMID: 34928279 DOI: 10.1039/d1cp05272a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
How solvent motions affect the dynamics of chemical reactions in which the solute undergoes a substantial shape change is a fundamental but elusive issue. This work utilizes reactive simulation and Grote-Hynes theory to explore the effect of solvent motions on the dynamics of the Diels-Alder reaction (in the reverse direction, this reaction involves very substantial solute expansion) in aprotic solvents. The results reveal that the solvent environment is not sufficiently constraining to influence transition state passage dynamics, with the calculated transmission coefficients being close to unity. Even when solvent motions are suppressed or artificially slowed down, the solvent only affects the reaction dynamics in the transition state region to a very small extent. The only notable effect of solvent occurs far from the transition state region and corresponds to caging of the reactants within the reactant well.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium.
| | - Pierre-Louis Lefebvre
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium. .,Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Jeremy N Harvey
- Theoretical and Computational Chemistry, Department of Chemistry, KU Leuven Celestijnenlaan 200F, 3001, Leuven, Belgium.
| |
Collapse
|
8
|
Wang X, Zhang C, Jiang Y, Wang W, Zhou Y, Chen Y, Zhang B, Tan RX, Ge HM, Yang ZJ, Liang Y. Influence of Water and Enzyme on the Post-Transition State Bifurcation of NgnD-Catalyzed Ambimodal [6+4]/[4+2] Cycloaddition. J Am Chem Soc 2021; 143:21003-21009. [PMID: 34851644 DOI: 10.1021/jacs.1c10760] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The enzyme NgnD catalyzes an ambimodal cycloaddition that bifurcates to [6+4]- and [4+2]-adducts. Both products have been isolated in experiments, but it remains unknown how enzyme and water influence the bifurcation selectivity at the femtosecond time scale. Here, we study the impact of water and enzyme on the post-transition state bifurcation of NgnD-catalyzed [6+4]/[4+2] cycloaddition by integrating quantum mechanics/molecular mechanics quasiclassical dynamics simulations and biochemical assays. The ratio of [6+4]/[4+2] products significantly differs in the gas phase, water, and enzyme. Biochemical assays were employed to validate computational predictions. The study informs how water and enzyme affect the bifurcation selectivity through perturbation of the reaction dynamics in the femtosecond time scale, revealing the fundamental roles of condensed media in dynamically controlling the chemical selectivity for biosynthetic reactions.
Collapse
Affiliation(s)
- Xin Wang
- State Key Laboratory of Coordination Chemistry, Jiangsu Key Laboratory of Advanced Organic Materials, Chemistry and Biomedicine Innovation Center, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China.,Henan-Macquarie University Joint Centre for Biomedical Innovation, School of Life Sciences, Henan University, Kaifeng, Henan 475004, China
| | - Chun Zhang
- School of Pharmaceutical Sciences, Jiangnan University, Wuxi 214122, China
| | - Yaoyukun Jiang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Wen Wang
- State Key Laboratory of Pharmaceutical Biotechnology, Institute of Functional Biomolecules, Chemistry and Biomedicine Innovation Center, School of Life Sciences, Nanjing University, Nanjing 210023, China
| | - Yuan Zhou
- Henan-Macquarie University Joint Centre for Biomedical Innovation, School of Life Sciences, Henan University, Kaifeng, Henan 475004, China
| | - Yu Chen
- State Key Laboratory of Coordination Chemistry, Jiangsu Key Laboratory of Advanced Organic Materials, Chemistry and Biomedicine Innovation Center, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Bo Zhang
- State Key Laboratory of Pharmaceutical Biotechnology, Institute of Functional Biomolecules, Chemistry and Biomedicine Innovation Center, School of Life Sciences, Nanjing University, Nanjing 210023, China
| | - Ren Xiang Tan
- State Key Laboratory of Pharmaceutical Biotechnology, Institute of Functional Biomolecules, Chemistry and Biomedicine Innovation Center, School of Life Sciences, Nanjing University, Nanjing 210023, China
| | - Hui Ming Ge
- State Key Laboratory of Pharmaceutical Biotechnology, Institute of Functional Biomolecules, Chemistry and Biomedicine Innovation Center, School of Life Sciences, Nanjing University, Nanjing 210023, China
| | - Zhongyue J Yang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States.,Vanderbilt Institute of Chemical Biology, Vanderbilt University, Nashville, Tennessee 37235, United States.,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States.,Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Yong Liang
- State Key Laboratory of Coordination Chemistry, Jiangsu Key Laboratory of Advanced Organic Materials, Chemistry and Biomedicine Innovation Center, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
9
|
Qin ZX, Tremblay M, Hong X, Yang ZJ. Entropic Path Sampling: Computational Protocol to Evaluate Entropic Profile along a Reaction Path. J Phys Chem Lett 2021; 12:10713-10719. [PMID: 34709848 DOI: 10.1021/acs.jpclett.1c03116] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Fleeting intermediates constitute dynamically stepwise mechanisms. They have been characterized in molecular dynamics trajectories, but whether these intermediates form a free energy minimum to become entropic intermediates remains elusively defined. We developed a computational protocol known as entropic path sampling to evaluate the entropic variation of reacting species along a reaction path based on an ensemble of trajectories. Using cyclopentadiene dimerization as a model reaction, we observed an entropy maximum along the reaction path which originates from an enhanced conformational flexibility as the reacting species enter into a flat energy region. As the reacting species further approach product formation, unfavorable entropic restriction fails to offset the potential energy drop, resulting in no free energy minimum along the post-TS pathway. Our results show that cyclopentadiene dimerization involves an entropy maximum that leads to dynamic intermediates with elongated lifetimes, but the reaction does not involve entropic intermediates.
Collapse
Affiliation(s)
- Zhi-Xin Qin
- Center of Chemistry for Frontier Technologies, Department of Chemistry, State Key Laboratory of Clean Energy Utilization, Zhejiang University, Hangzhou 310027, China
| | - Matthew Tremblay
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Xin Hong
- Center of Chemistry for Frontier Technologies, Department of Chemistry, State Key Laboratory of Clean Energy Utilization, Zhejiang University, Hangzhou 310027, China
| | - Zhongyue J Yang
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee 37235, United States
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Vanderbilt Institute of Chemical Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
- Data Science Institute, Vanderbilt University, Nashville, Tennessee 37235, United States
| |
Collapse
|
10
|
Zhang X, Vázquez SA, Harvey JN. Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials. J Chem Theory Comput 2021; 17:1277-1289. [PMID: 33550803 DOI: 10.1021/acs.jctc.0c01059] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d3 (CD3CN) and dichloromethane-d2 (CD2Cl2). However, insights into how CD2Cl2 facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD2Cl2. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVBCI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD2Cl2 interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVBCI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD2Cl2 as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD2Cl2. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD2Cl2 interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| | - Saulo A Vázquez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Jeremy N Harvey
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| |
Collapse
|
11
|
Scivetti I, Sen K, Elena AM, Todorov I. Reactive Molecular Dynamics at Constant Pressure via Nonreactive Force Fields: Extending the Empirical Valence Bond Method to the Isothermal-Isobaric Ensemble. J Phys Chem A 2020; 124:7585-7597. [PMID: 32820921 DOI: 10.1021/acs.jpca.0c05461] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Empirical Valence Bond (EVB) method offers a suitable framework to obtain reactive potentials through the coupling of nonreactive force fields. In this formalism, most of the implemented coupling terms are built using functional forms that depend on spatial coordinates, while parameters are fitted against reference data to model the change of chemistry between the participating nonreactive states. In this work, we demonstrate that the use of such coupling terms precludes the computation of the stress tensor for condensed phase systems and prevents the possibility to carry out EVB molecular dynamics in the isothermal-isobaric (NPT) ensemble. Alternatively, we make use of coupling terms that depend on the energy gaps, defined as the energy differences between the participating nonreactive force fields, and derive a general expression for the EVB stress tensor suitable for computation. Implementation of this new methodology is tested for a model of a single reactive malonaldehyde solvated in nonreactive water. Mass densities and probability distributions for the values of the energy gaps computed in the NPT ensemble reveal a negligible role of the reactive potential in the limit of low concentrated solutions, thus corroborating for the first time the validity of approximations based on the canonical NVT ensemble, customarily adopted for EVB simulations. The presented formalism also aims to contribute to future implementations and extensions of the EVB method to research the limit of highly concentrated solutions.
Collapse
Affiliation(s)
- Ivan Scivetti
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K.,Department of Chemistry, University of Liverpool, Liverpool L69 3BX, U.K
| | - Kakali Sen
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Alin M Elena
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| | - Ilian Todorov
- Daresbury Laboratory, Sc. Tech., Keckwick Lane, Daresbury, Warrington WA4 4AD, U.K
| |
Collapse
|
12
|
Roytman VA, Jin S, Nguyen VT, Nguyen VD, Haug GC, Larionov OV, Singleton DA. Bond Memory in Dynamically Determined Stereoselectivity. J Am Chem Soc 2019; 142:85-88. [PMID: 31852185 DOI: 10.1021/jacs.9b12227] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The carboborative ring contraction of cyclohexenes exhibits an abnormal selectivity pattern in which a formally concerted double migration gives rise to predominant but not exclusive inversion products. In dynamic trajectories, the inversion and retention products are formed from the same transition state, and the trajectories accurately account for the experimental product ratios. The unusual origin of the selectivity is the dynamically retained non-equivalence of newly formed versus pre-existing bonds after the first bond migration.
Collapse
Affiliation(s)
- Vladislav A Roytman
- Department of Chemistry , Texas A&M University , P.O. Box 30012, College Station , Texas 77842 , United States
| | - Shengfei Jin
- Department of Chemistry , The University of Texas at San Antonio , San Antonio , Texas 78249 , United States
| | - Vu T Nguyen
- Department of Chemistry , The University of Texas at San Antonio , San Antonio , Texas 78249 , United States
| | - Viet D Nguyen
- Department of Chemistry , The University of Texas at San Antonio , San Antonio , Texas 78249 , United States
| | - Graham C Haug
- Department of Chemistry , The University of Texas at San Antonio , San Antonio , Texas 78249 , United States
| | - Oleg V Larionov
- Department of Chemistry , The University of Texas at San Antonio , San Antonio , Texas 78249 , United States
| | - Daniel A Singleton
- Department of Chemistry , Texas A&M University , P.O. Box 30012, College Station , Texas 77842 , United States
| |
Collapse
|
13
|
Hare SR, Bratholm LA, Glowacki DR, Carpenter BK. Low dimensional representations along intrinsic reaction coordinates and molecular dynamics trajectories using interatomic distance matrices. Chem Sci 2019; 10:9954-9968. [PMID: 32055352 PMCID: PMC6991188 DOI: 10.1039/c9sc02742d] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Accepted: 08/23/2019] [Indexed: 01/22/2023] Open
Abstract
Most chemical transformations (reactions or conformational changes) that are of interest to researchers have many degrees of freedom, usually too many to visualize without reducing the dimensionality of the system to include only the most important atomic motions. In this article, we describe a method of using Principal Component Analysis (PCA) for analyzing a series of molecular geometries (e.g., a reaction pathway or molecular dynamics trajectory) and determining the reduced dimensional space that captures the most structural variance in the fewest dimensions. The software written to carry out this method is called PathReducer, which permits (1) visualizing the geometries in a reduced dimensional space, (2) determining the axes that make up the reduced dimensional space, and (3) projecting the series of geometries into the low-dimensional space for visualization. We investigated two options to represent molecular structures within PathReducer: aligned Cartesian coordinates and matrices of interatomic distances. We found that interatomic distance matrices better captured non-linear motions in a smaller number of dimensions. To demonstrate the utility of PathReducer, we have carried out a number of applications where we have projected molecular dynamics trajectories into a reduced dimensional space defined by an intrinsic reaction coordinate. The visualizations provided by this analysis show that dynamic paths can differ greatly from the minimum energy pathway on a potential energy surface. Viewing intrinsic reaction coordinates and trajectories in this way provides a quick way to gather qualitative information about the pathways trajectories take relative to a minimum energy path. Given that the outputs from PCA are linear combinations of the input molecular structure coordinates (i.e., Cartesian coordinates or interatomic distances), they can be easily transferred to other types of calculations that require the definition of a reduced dimensional space (e.g., biased molecular dynamics simulations).
Collapse
Affiliation(s)
- Stephanie R Hare
- University of Bristol School of Chemistry , Cantock's Close , Bristol , UK BS8 1TS
- University of Bristol School of Mathematics , University Walk , Bristol , UK BS8 1TW
| | - Lars A Bratholm
- University of Bristol School of Chemistry , Cantock's Close , Bristol , UK BS8 1TS
- University of Bristol School of Mathematics , University Walk , Bristol , UK BS8 1TW
| | - David R Glowacki
- University of Bristol School of Chemistry , Cantock's Close , Bristol , UK BS8 1TS
- University of Bristol School of Computer Science , Merchant Venturers Building, Woodland Road , Bristol , UK BS8 1UB
| | - Barry K Carpenter
- Cardiff University School of Chemistry , Main Building, Park Place , Cardiff , UK CF10 3AT .
| |
Collapse
|
14
|
Yang Z, Jamieson CS, Xue XS, Garcia-Borràs M, Benton T, Dong X, Liu F, Houk K. Mechanisms and Dynamics of Reactions Involving Entropic Intermediates. TRENDS IN CHEMISTRY 2019. [DOI: 10.1016/j.trechm.2019.01.009] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
15
|
Shannon RJ, Hornung B, Tew DP, Glowacki DR. Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field. J Phys Chem A 2019; 123:2991-2999. [PMID: 30793911 DOI: 10.1021/acs.jpca.8b12006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.
Collapse
Affiliation(s)
- R J Shannon
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Mechanical Engineering , Stanford University , 452 Escondido Mall , Stanford , California 94305 , United States
| | - B Hornung
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D P Tew
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D R Glowacki
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Computer Science , University of Bristol , Bristol BS8 1UB , U.K
| |
Collapse
|
16
|
Amabilino S, Bratholm LA, Bennie SJ, Vaucher AC, Reiher M, Glowacki DR. Training Neural Nets To Learn Reactive Potential Energy Surfaces Using Interactive Quantum Chemistry in Virtual Reality. J Phys Chem A 2019; 123:4486-4499. [DOI: 10.1021/acs.jpca.9b01006] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Silvia Amabilino
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Lars A. Bratholm
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Simon J. Bennie
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Alain C. Vaucher
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | | |
Collapse
|
17
|
Unke OT, Brickel S, Meuwly M. Sampling reactive regions in phase space by following the minimum dynamic path. J Chem Phys 2019; 150:074107. [DOI: 10.1063/1.5082885] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Oliver T. Unke
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Sebastian Brickel
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
18
|
Nikitin K, O'Gara R. Mechanisms and Beyond: Elucidation of Fluxional Dynamics by Exchange NMR Spectroscopy. Chemistry 2019; 25:4551-4589. [PMID: 30421834 DOI: 10.1002/chem.201804123] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2018] [Indexed: 12/31/2022]
Abstract
Detailed mechanistic information is crucial to our understanding of reaction pathways and selectivity. Dynamic exchange NMR techniques, in particular 2D exchange spectroscopy (EXSY) and its modifications, provide indispensable intricate information on the mechanisms of organic and inorganic reactions and other phenomena, for example, the dynamics of interfacial processes. In this Review, key results from exchange NMR studies of small molecules over the last few decades are systemised and discussed. After a brief introduction to the theory, the key types of dynamic processes are identified and fundamental examples given of intra- and intermolecular reactions, which, in turn, could involve, or not, bond-making and bond-breaking events. Following that logic, internal molecular rotation, intramolecular stereomutation and molecular recognition will first be considered because they do not typically involve bond breaking. Then, rearrangements, substitution-type reactions, cyclisations, additions and other processes affecting chemical bonds will be discussed. Finally, interfacial molecular dynamics and unexpected combinations of different types of fluxional processes will also be highlighted. How exchange NMR spectroscopy helps to identify conformational changes, coordination and molecular recognition processes as well as quantify reaction energy barriers and extract detailed mechanistic information by using reaction rate theory in conjunction with computational techniques will be shown.
Collapse
Affiliation(s)
- Kirill Nikitin
- School of Chemistry, University College Dublin, Belfield, Dublin, Ireland
| | - Ryan O'Gara
- School of Chemistry, University College Dublin, Belfield, Dublin, Ireland
| |
Collapse
|
19
|
Zhang X, Harvey JN. EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions. Phys Chem Chem Phys 2019; 21:14331-14340. [DOI: 10.1039/c8cp06686h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.
Collapse
Affiliation(s)
- Xiaoyong Zhang
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry
- KU Leuven
- B-3001 Leuven
- Belgium
| | - Jeremy N. Harvey
- Department of Chemistry and Division of Quantum Chemistry and Physical Chemistry
- KU Leuven
- B-3001 Leuven
- Belgium
| |
Collapse
|
20
|
The influence of the solvent’s mass on the location of the dividing surface for a model Hamiltonian. Chem Phys Lett 2019. [DOI: 10.1016/j.cpletx.2019.100030] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
21
|
Kulkarni Y, Kamerlin SCL. Computational physical organic chemistry using the empirical valence bond approach. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2019. [DOI: 10.1016/bs.apoc.2019.07.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
22
|
Hare SR, Li A, Tantillo DJ. Post-transition state bifurcations induce dynamical detours in Pummerer-like reactions. Chem Sci 2018; 9:8937-8945. [PMID: 30627409 PMCID: PMC6296359 DOI: 10.1039/c8sc02653j] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2018] [Accepted: 09/25/2018] [Indexed: 12/28/2022] Open
Abstract
A post-transition state bifurcation (PTSB) involved in a Pummerer-type rearrangement is characterized using density functional theory (DFT) calculations on potential energy stationary points and direct dynamics simulations. A sensitivity of the ratio of products produced via this PTSB to solvent dielectric constant is revealed and implications of such a dependence for selectivity control of organic reactions are discussed.
Collapse
Affiliation(s)
| | - Ang Li
- Shanghai Institute of Organic Chemistry , China
| | | |
Collapse
|
23
|
Chen S, Yu P, Houk KN. Ambimodal Dipolar/Diels–Alder Cycloaddition Transition States Involving Proton Transfers. J Am Chem Soc 2018; 140:18124-18131. [DOI: 10.1021/jacs.8b11080] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Shuming Chen
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Peiyuan Yu
- Energy Storage and Distributed Resources Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - K. N. Houk
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| |
Collapse
|
24
|
Garcia-Meseguer R, Carpenter BK. Re-Evaluating the Transition State for Reactions in Solution. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800841] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Barry K. Carpenter
- School of Chemistry; Cardiff University; CF10 3AT Cardiff United Kingdom
| |
Collapse
|
25
|
Meuwly M. Reactive molecular dynamics: From small molecules to proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1386] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Markus Meuwly
- Department of Chemistry University of Basel Basel Switzerland
- Department of Chemistry Brown University Providence Rhode Island
| |
Collapse
|
26
|
Athokpam B, Ramesh SG. Alkyl hydrogen atom abstraction reactions of the CN radical with ethanol. J Chem Phys 2018; 148:134503. [PMID: 29626852 DOI: 10.1063/1.5021634] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a study of the abstraction of alkyl hydrogen atoms from the β and α positions of ethanol by the CN radical in solution using the Empirical Valence Bond (EVB) method. We have built separate 2 × 2 EVB models for the Hβ and Hα reactions, where the atom transfer is parameterized using ab initio calculations. The intra- and intermolecular potentials of the reactant and product molecules were modelled with the General AMBER Force Field, with some modifications. We have carried out the dynamics in water and chloroform, which are solvents of contrasting polarity. We have computed the potential of mean force for both abstractions in each of the solvents. They are found to have a small and early barrier along the reaction coordinate with a large energy release. Analyzing the solvent structure around the reaction system, we have found two solvents to have little effect on either reaction. Simulating the dynamics from the transition state, we also study the fate of the energies in the HCN vibrational modes. The HCN molecule is born vibrationally hot in the CH stretch in both reactions and additionally in the HCN bends for the Hα abstraction reaction. In the early stage of the dynamics, we find that the CN stretch mode gains energy at the expense of the energy in CH stretch mode.
Collapse
Affiliation(s)
- Bijyalaxmi Athokpam
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| | - Sai G Ramesh
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
27
|
Essafi S, Harvey JN. Rates of Molecular Vibrational Energy Transfer in Organic Solutions. J Phys Chem A 2018; 122:3535-3540. [DOI: 10.1021/acs.jpca.7b12563] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stéphanie Essafi
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Jeremy N. Harvey
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| |
Collapse
|
28
|
Yang Z, Dong X, Yu Y, Yu P, Li Y, Jamieson C, Houk KN. Relationships between Product Ratios in Ambimodal Pericyclic Reactions and Bond Lengths in Transition Structures. J Am Chem Soc 2018; 140:3061-3067. [DOI: 10.1021/jacs.7b13562] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Zhongyue Yang
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Xiaofei Dong
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Yanmin Yu
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
- Beijing
Key Laboratory for Green Catalysis and Separation, Department of Chemistry
and Chemical Engineering, Beijing University of Technology, Beijing 100124, China
| | - Peiyuan Yu
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
- Energy
Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yingzi Li
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Cooper Jamieson
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - K. N. Houk
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| |
Collapse
|
29
|
Shannon R, Glowacki DR. A Simple “Boxed Molecular Kinetics” Approach To Accelerate Rare Events in the Stochastic Kinetic Master Equation. J Phys Chem A 2018; 122:1531-1541. [DOI: 10.1021/acs.jpca.7b12521] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Robin Shannon
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
| | - David R. Glowacki
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, BS8 1UB, United Kingdom
| |
Collapse
|
30
|
Bailey JO, Singleton DA. Failure and Redemption of Statistical and Nonstatistical Rate Theories in the Hydroboration of Alkenes. J Am Chem Soc 2017; 139:15710-15723. [PMID: 29029560 PMCID: PMC5831723 DOI: 10.1021/jacs.7b07175] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Our previous work found that canonical forms of transition state theory incorrectly predict the regioselectivity of the hydroboration of propene with BH3 in solution. In response, it has been suggested that alternative statistical and nonstatistical rate theories can adequately account for the selectivity. This paper uses a combination of experimental and theoretical studies to critically evaluate the ability of these rate theories, as well as dynamic trajectories and newly developed localized statistical models, to predict quantitative selectivities and qualitative trends in hydroborations on a broader scale. The hydroboration of a series of terminally substituted alkenes with BH3 was examined experimentally, and a classically unexpected trend is that the selectivity increases as the alkyl chain is lengthened far from the reactive centers. Conventional and variational transition state theories can predict neither the selectivities nor the trends. The canonical competitive nonstatistical model makes somewhat better predictions for some alkenes but fails to predict trends, and it performs poorly with an alkene chosen to test a specific prediction of the model. Added nonstatistical corrections to this model make the predictions worse. Parametrized Rice-Ramsperger-Kassel-Marcus (RRKM)-master equation calculations correctly predict the direction of the trend in selectivity versus alkene size but overpredict its magnitude, and the selectivity with large alkenes remains unpredictable with any parametrization. Trajectory studies in explicit solvent can predict selectivities without parametrization but are impractical for predicting small changes in selectivity. From a lifetime and energy analysis of the trajectories, "localized RRKM-ME" and "competitive localized noncanonical" rate models are suggested as steps toward a general model. These provide the best predictions of the experimental observations and insight into the selectivities.
Collapse
Affiliation(s)
- Johnathan O. Bailey
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| | - Daniel A. Singleton
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| |
Collapse
|
31
|
Hare SR, Tantillo DJ. Post-transition state bifurcations gain momentum – current state of the field. PURE APPL CHEM 2017. [DOI: 10.1515/pac-2017-0104] [Citation(s) in RCA: 105] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
AbstractThe existence of post-transition state bifurcations on potential energy surfaces for organic and biological reaction mechanisms has been known for decades, but recently, new reports of bifurcations have been occurring at a much higher rate. Beyond simply discovering bifurcations, computational chemists are developing techniques to understand what aspects of molecular structure and vibrations control the product selectivity in systems containing bifurcations. For example, the distribution of products seen in simulations has been found to be extremely sensitive to the local environment of the reacting system (i.e. the presence of a catalyst, enzyme, or explicit solvent molecules). The outlook for the future of this field is discussed, with an eye towards the application of the principles discussed here by experimental chemists to design a reaction setup to efficiently generate desired products.
Collapse
Affiliation(s)
- Stephanie R. Hare
- University of California, Davis, 1 Shields Avenue, Davis, CA 95616, USA
| | - Dean J. Tantillo
- University of California, Davis, 1 Shields Avenue, Davis, CA 95616, USA
| |
Collapse
|
32
|
Hare SR, Pemberton RP, Tantillo DJ. Navigating Past a Fork in the Road: Carbocation-π Interactions Can Manipulate Dynamic Behavior of Reactions Facing Post-Transition-State Bifurcations. J Am Chem Soc 2017; 139:7485-7493. [PMID: 28504880 DOI: 10.1021/jacs.7b01042] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Dynamics calculations are described for carbocation rearrangements involving product-forming pathways with post-transition-state bifurcations. We show that noncovalent interactions with associated benzene rings (a simple model of aromatic amino acid side chains) can switch inherent dynamical tendencies for competing modes of disrotation, establishing that meaningful changes in dynamically controlled product selectivity can be achieved with few weak noncovalent interactions.
Collapse
Affiliation(s)
- Stephanie R Hare
- Department of Chemistry, University of California, Davis , One Shields Avenue, Davis, California 95616, United States
| | - Ryan P Pemberton
- Department of Chemistry, University of California, Davis , One Shields Avenue, Davis, California 95616, United States
| | - Dean J Tantillo
- Department of Chemistry, University of California, Davis , One Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
33
|
Zhan S, Mårtensson D, Purg M, Kamerlin SCL, Ahlquist MSG. Capturing the Role of Explicit Solvent in the Dimerization of RuV
(bda) Water Oxidation Catalysts. Angew Chem Int Ed Engl 2017; 56:6962-6965. [DOI: 10.1002/anie.201701488] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Revised: 03/23/2017] [Indexed: 11/06/2022]
Affiliation(s)
- Shaoqi Zhan
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| | - Daniel Mårtensson
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| | - Miha Purg
- Department of Cell and Molecular Biology; Uppsala University; 75124 Uppsala Sweden
| | - Shina C. L. Kamerlin
- Department of Cell and Molecular Biology; Uppsala University; 75124 Uppsala Sweden
| | - Mårten S. G. Ahlquist
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| |
Collapse
|
34
|
Zhan S, Mårtensson D, Purg M, Kamerlin SCL, Ahlquist MSG. Capturing the Role of Explicit Solvent in the Dimerization of RuV
(bda) Water Oxidation Catalysts. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201701488] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Shaoqi Zhan
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| | - Daniel Mårtensson
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| | - Miha Purg
- Department of Cell and Molecular Biology; Uppsala University; 75124 Uppsala Sweden
| | - Shina C. L. Kamerlin
- Department of Cell and Molecular Biology; Uppsala University; 75124 Uppsala Sweden
| | - Mårten S. G. Ahlquist
- Division of Theoretical Chemistry & Biology; School of Biotechnology; KTH Royal Institute of Technology; 10691 Stockholm Sweden
| |
Collapse
|
35
|
Spezia R, Martínez-Nuñez E, Vazquez S, Hase WL. Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:20170035. [PMID: 28320909 PMCID: PMC5360905 DOI: 10.1098/rsta.2017.0035] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 01/24/2017] [Indexed: 06/06/2023]
Abstract
In this Introduction, we show the basic problems of non-statistical and non-equilibrium phenomena related to the papers collected in this themed issue. Over the past few years, significant advances in both computing power and development of theories have allowed the study of larger systems, increasing the time length of simulations and improving the quality of potential energy surfaces. In particular, the possibility of using quantum chemistry to calculate energies and forces 'on the fly' has paved the way to directly study chemical reactions. This has provided a valuable tool to explore molecular mechanisms at given temperatures and energies and to see whether these reactive trajectories follow statistical laws and/or minimum energy pathways. This themed issue collects different aspects of the problem and gives an overview of recent works and developments in different contexts, from the gas phase to the condensed phase to excited states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.
Collapse
Affiliation(s)
- Riccardo Spezia
- Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, CEA CNRS Université Paris Saclay, 91025 Evry, France
- LAMBE, Université d'Evry, 91025 Evry, France
| | - Emilio Martínez-Nuñez
- Departamento de Química Física and Centro Singular de Investigación en Química, Biológica y Materiales Moleculares (CIQUS), Campus Vida, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - Saulo Vazquez
- Departamento de Química Física and Centro Singular de Investigación en Química, Biológica y Materiales Moleculares (CIQUS), Campus Vida, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - William L Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX 79409, USA
| |
Collapse
|
36
|
Quapp W, Bofill JM, Ribas-Ariño J. Analysis of the Acting Forces in a Theory of Catalysis and Mechanochemistry. J Phys Chem A 2017; 121:2820-2838. [PMID: 28338327 DOI: 10.1021/acs.jpca.7b00022] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The theoretical description of a chemical process resulting from the application of mechanical or catalytical stress to a molecule is performed by the generation of an effective potential energy surface (PES). Changes for minima and saddle points by the stress are described by Newton trajectories (NTs) on the original PES. From the analysis of the acting forces we postulate the existence of pulling corridors built by families of NTs that connect the same stationary points. For different exit saddles of different height we discuss the corresponding pulling corridors; mainly by simple two-dimensional surface models. If there are different exit saddles then there can exist saddles of index two, at least, between. Then the case that a full pulling corridor crosses a saddle of index two is the normal case. It leads to an intrinsic hysteresis of such pullings for the forward or the backward reaction. Assuming such relations we can explain some results in the literature. A new finding is the existence of roundabout corridors that can switch between different saddle points by a reversion of the direction. The findings concern the mechanochemistry of molecular systems under a mechanical load as well as the electrostatic force and can be extended to catalytic and enzymatic accelerated reactions. The basic and ground ansatz includes both kinds of forces in a natural way without an extra modification.
Collapse
Affiliation(s)
- Wolfgang Quapp
- Mathematisches Institut, Universität Leipzig , PF 100920, D-04009 Leipzig, Germany
| | | | | |
Collapse
|
37
|
Diamantis P, Unke OT, Meuwly M. Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N. PLoS Comput Biol 2017; 13:e1005450. [PMID: 28358830 PMCID: PMC5391117 DOI: 10.1371/journal.pcbi.1005450] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2016] [Revised: 04/13/2017] [Accepted: 03/13/2017] [Indexed: 11/18/2022] Open
Abstract
In heme proteins, the efficient transport of ligands such as NO or O2 to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it. Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.
Collapse
Affiliation(s)
| | - Oliver T. Unke
- Department of Chemistry, University of Basel, Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Basel, Switzerland
- * E-mail:
| |
Collapse
|
38
|
Gillick-Healy MW, Jennings EV, Müller-Bunz H, Ortin Y, Nikitin K, Gilheany DG. Two Independent Orthogonal Stereomutations at a Single Asymmetric Center: A Narcissistic Couple. Chemistry 2016; 23:2332-2339. [DOI: 10.1002/chem.201604080] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2016] [Indexed: 01/31/2023]
Affiliation(s)
| | | | | | - Yannick Ortin
- School of Chemistry; University College Dublin; Ireland
| | | | | |
Collapse
|
39
|
Carpenter BK, Harvey JN, Orr-Ewing AJ. The Study of Reactive Intermediates in Condensed Phases. J Am Chem Soc 2016; 138:4695-705. [DOI: 10.1021/jacs.6b01761] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Barry K. Carpenter
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, U.K
| | - Jeremy N. Harvey
- Department
of Chemistry, KU Leuven, Celestijnen Laan 200F, B-3001 Heverlee, Belgium
| | - Andrew J. Orr-Ewing
- School
of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, U.K
| |
Collapse
|
40
|
Glowacki DR, Orr-Ewing AJ, Harvey JN. Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD3CN treated with a parallel multi-state EVB model. J Chem Phys 2015; 143:044120. [DOI: 10.1063/1.4926996] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David R. Glowacki
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
- PULSE Institute and Department of Chemistry, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | | | - Jeremy N. Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Heverlee, Belgium
| |
Collapse
|
41
|
Molecular Dynamics of the Diels–Alder Reactions of Tetrazines with Alkenes and N2 Extrusions from Adducts. J Am Chem Soc 2015; 137:4749-58. [DOI: 10.1021/jacs.5b00014] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
42
|
Abstract
Combining QMDFF with EVB allows to generate reactive force fields of useful quality for widely varying reactions with minimal effort.
Collapse
Affiliation(s)
- Bernd Hartke
- Institute for Physical Chemistry
- Christian-Albrechts-University
- D-24118 Kiel
- Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry
- Institut für Physikalische und Theoretische Chemie
- Universität Bonn
- 53115 Bonn
- Germany
| |
Collapse
|