Parameterization of the ReaxFF reactive force field for a proline‐catalyzed aldol reaction

A parameterization of the ReaxFF reactive FF is performed using a Monte Carlo Simulated Annealing procedure for the modeling of a proline‐catalyzed aldol reaction. Emphasis is put on the accurate reproduction of the relative stabilities of several key intermediates of the reaction, as well as, on the description of the reaction path bridging these intermediates based on quantum mechanical calculations. Our training sets include new criteria based on geometry optimizations and short Molecular Dynamics simulations to ensure that the trained ReaxFF potentials adequately predict the structures of all key intermediates. The transferability of the sets of parameters obtained is assessed for various steps of the considered aldol reaction, as well as for different substrates, catalysts, and reagents. This works indeed highlights the challenge of reaching transferable parameters for several reaction steps. © 2016 Wiley Periodicals, Inc.

[1]  F. Giacalone,et al.  New Simple Hydrophobic Proline Derivatives as Highly Active and Stereoselective Catalysts for the Direct Asymmetric Aldol Reaction in Aqueous Medium , 2008 .

[2]  Sudhir B. Kylasa,et al.  The ReaxFF reactive force-field: development, applications and future directions , 2016 .

[3]  William A Goddard,et al.  Development of a ReaxFF reactive force field for glycine and application to solvent effect and tautomerization. , 2011, The journal of physical chemistry. B.

[4]  F. Giacalone,et al.  Supported proline and proline-derivatives as recyclable organocatalysts. , 2008, Chemical Society reviews.

[5]  Paul Ha-Yeon Cheong,et al.  Computational prediction of small-molecule catalysts , 2008, Nature.

[6]  Dmitry Bedrov,et al.  Reactions of singly-reduced ethylene carbonate in lithium battery electrolytes: a molecular dynamics simulation study using the ReaxFF. , 2012, The journal of physical chemistry. A.

[7]  S. Stuart,et al.  A reactive potential for hydrocarbons with intermolecular interactions , 2000 .

[8]  Paul Ha-Yeon Cheong,et al.  Theory of asymmetric organocatalysis of Aldol and related reactions: rationalizations and predictions. , 2004, Accounts of chemical research.

[9]  Markus Hütter,et al.  Parameterization of a reactive force field using a Monte Carlo algorithm , 2013, J. Comput. Chem..

[10]  M. Paddon-Row,et al.  Theory and modeling of stereoselective organic reactions. , 1986, Science.

[11]  Richard A. Lerner,et al.  Proline-Catalyzed Direct Asymmetric Aldol Reactions , 2000 .

[12]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[13]  A. V. Duin,et al.  ReaxFF: A Reactive Force Field for Hydrocarbons , 2001 .

[14]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[15]  K N Houk,et al.  Quantum mechanical predictions of the stereoselectivities of proline-catalyzed asymmetric intermolecular aldol reactions. , 2003, Journal of the American Chemical Society.

[16]  Shashank Chaturvedi,et al.  Determination of best-fit potential parameters for a reactive force field using a genetic algorithm , 2012, Journal of Molecular Modeling.

[17]  A. V. van Duin,et al.  Carbon cluster formation during thermal decomposition of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-triamino-2,4,6-trinitrobenzene high explosives from ReaxFF reactive molecular dynamics simulations. , 2009, The journal of physical chemistry. A.

[18]  Adi M. Treasurywala,et al.  Are Oxazolidinones Really Unproductive, Parasitic Species in Proline Catalysis? – Thoughts and Experiments Pointing to an Alternative View , 2007 .

[19]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[20]  Michael Page,et al.  Following steepest descent reaction paths. The use of higher energy derivatives with ab initio electronic structure methods , 1990 .

[21]  J. Tersoff,et al.  New empirical approach for the structure and energy of covalent systems. , 1988, Physical review. B, Condensed matter.

[22]  Pierre O. Hubin,et al.  Accessing the free energy profile of a ring closure in a proline-catalyzed reaction using a reactive force field , 2015, Theoretical Chemistry Accounts.

[23]  Betsy M Rice,et al.  Parameterizing complex reactive force fields using multiple objective evolutionary strategies (MOES). Part 1: ReaxFF models for cyclotrimethylene trinitramine (RDX) and 1,1-diamino-2,2-dinitroethene (FOX-7). , 2015, Journal of chemical theory and computation.

[24]  Susan B. Sinnott,et al.  Charge optimized many-body potential for the Si/SiO2 system , 2007 .

[25]  R. B. Sunoj,et al.  Stereocontrol in proline-catalyzed asymmetric amination: a comparative assessment of the role of enamine carboxylic acid and enamine carboxylate. , 2011, Chemical communications.

[26]  Vincenzo Barone,et al.  Exploring the conformational and reactive dynamics of biomolecules in solution using an extended version of the glycine reactive force field. , 2013, Physical chemistry chemical physics : PCCP.

[27]  A. V. van Duin,et al.  Simulations on the thermal decomposition of a poly(dimethylsiloxane) polymer using the ReaxFF reactive force field. , 2005, Journal of the American Chemical Society.

[28]  Arieh Warshel,et al.  An empirical valence bond approach for comparing reactions in solutions and in enzymes , 1980 .

[29]  A. V. van Duin,et al.  ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. , 2008, The journal of physical chemistry. A.

[30]  D. MacMillan,et al.  The advent and development of organocatalysis , 2008, Nature.

[31]  S. Papson.  “Model” , 1981 .

[32]  Julian Müller,et al.  Efficient global optimization of reactive force‐field parameters , 2015, J. Comput. Chem..

[33]  Dandan Jiang,et al.  Development, applications and challenges of ReaxFF reactive force field in molecular simulations , 2016, Frontiers of Chemical Science and Engineering.

[34]  Denis Jacquemin,et al.  Ab initio quantum chemical and ReaxFF-based study of the intramolecular iminium–enamine conversion in a proline-catalyzed reaction , 2012, Theoretical Chemistry Accounts.

[35]  R. B. Sunoj,et al.  Enamine versus oxazolidinone: what controls stereoselectivity in proline-catalyzed asymmetric aldol reactions? , 2010, Angewandte Chemie.

[36]  K N Houk,et al.  Kinetic and stereochemical evidence for the involvement of only one proline molecule in the transition states of proline-catalyzed intra- and intermolecular aldol reactions. , 2003, Journal of the American Chemical Society.

[37]  Henry S. Rzepa,et al.  The Houk–List transition states for organocatalytic mechanisms revisited , 2014 .

[38]  K N Houk,et al.  Transition states of amine-catalyzed aldol reactions involving enamine intermediates: theoretical studies of mechanism, reactivity, and stereoselectivity. , 2001, Journal of the American Chemical Society.

[39]  Adri C. T. van Duin,et al.  Delft molecular mechanics: a new approach to hydrocarbon force fields. Inclusion of a geometry-dependent charge calculation , 1994 .

[40]  Bernd Hartke,et al.  Global optimization of parameters in the reactive force field ReaxFF for SiOH , 2013, J. Comput. Chem..

[41]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[42]  A. V. van Duin,et al.  ReaxFF molecular dynamics simulation of thermal stability of a Cu3(BTC)2 metal-organic framework. , 2012, Physical chemistry chemical physics : PCCP.

[43]  W. Goddard,et al.  General Multiobjective Force Field Optimization Framework, with Application to Reactive Force Fields for Silicon Carbide. , 2014, Journal of chemical theory and computation.

[44]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .