A free‐energy perturbation method based on Monte Carlo simulations using quantum mechanical calculations (QM/MC/FEP method): Application to highly solvent‐dependent reactions

This study describes the framework of the quantum mechanical (QM)/Monte Carlo (MC)/free‐energy perturbation (FEP) method, a FEP method based on MC simulations using quantum chemical calculations. Because a series of structures generated by interpolating internal coordinates between transition state and reactant did not produce smooth free‐energy profiles, we used structures from the intrinsic reaction coordinate calculations. This method was first applied to the Diels–Alder reaction between methyl vinyl ketone and cyclopentadiene and produced ΔG  sol‡ values of 20.1 and 21.4 kcal mol−1 in aqueous and methanol solutions, respectively. They are very consistent with the experimentally observed values. The other two applications were the free‐energy surfaces for the Cope elimination of N,N‐dimethyl‐3‐phenylbutan‐2‐amine oxide in aqueous, dimethyl sulfoxide, and tetrahydrofuran solutions, and the Kemp decarboxylation of 6‐hydroxybenzo‐isoxazole‐3‐carboxylic acid in aqueous, dimethyl sulfoxide, and CH3CN solutions. The calculated activation free energies differed by less than 1.8 kcal mol−1 from the experimental values for these reactions. Although we used droplet models for the QM/MC/FEP simulations, the calculated results for three reactions are very close to the experimental data. It was confirmed that most of the interactions between the solute and solvents can be described using small numbers of solvent molecules. This is because a few solvent molecules can produce large portions of the solute–solvent interaction energies at the reaction centers. When we confirmed the dependency on the droplet sizes of solvents, the QM/MC/FEP for a large droplet with 106 water molecules produced a ΔG  sol‡ value similar to the experimental values, as well as that for a small droplet with 34 molecules. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011

[1]  H. Bernhard Schlegel,et al.  Reaction Path Following in Mass-Weighted Internal Coordinates , 1990 .

[2]  Kimito Funatsu,et al.  A Novel Approach to Retrosynthetic Analysis Using Knowledge Bases Derived from Reaction Databases , 1999, J. Chem. Inf. Comput. Sci..

[3]  Kenzi Hori,et al.  Parallelized QM/MC simulations applied to the Diels-Alder reaction in aqueous, methanol and propane solutions , 2008 .

[4]  Peter A. Kollman,et al.  Solvation Free Energies of Amides and Amines: Disagreement between Free Energy Calculations and Experiment , 1995 .

[5]  D. Kemp,et al.  Physical organic chemistry of benzisoxazoles. IV. Origins and catalytic nature of the solvent rate acceleration for the decarboxylation of 3-carboxybenzisoxazoles , 1975 .

[6]  Ronald Breslow,et al.  Hydrophobic Effects on Simple Organic Reactions in Water , 1991 .

[7]  堀 憲次,et al.  Ethyl Methacrylate 及び Benzyl Methacylateへの適用 , 2007 .

[8]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[9]  Jan B. F. N. Engberts,et al.  Diels-Alder reactions in aqueous solutions. Enforced hydrophobic interactions between diene and dienophile , 1991 .

[10]  Kenzi Hori,et al.  Theoretical study of a reaction path via a hydrogen-bonded intermediate for the alkaline hydrolysis of esters in the gas phase , 1992 .

[11]  Manuel A. Aguilar,et al.  A new method to locate saddle points for reactions in solution by using the free‐energy gradient method and the mean field approximation , 2004, J. Comput. Chem..

[12]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[13]  William L Jorgensen,et al.  Influence of inter- and intramolecular hydrogen bonding on kemp decarboxylations from QM/MM simulations. , 2005, Journal of the American Chemical Society.

[14]  William L. Jorgensen,et al.  Solvent effects on a Diels-Alder reaction from computer simulations , 1991 .

[15]  F.J.Olivares del Valle,et al.  Geometry optimization of molecules in solution: Joint use of the mean field approximation and the free-energy gradient method , 2003 .

[16]  Jan B. F. N. Engberts.  The Role of the Solvent in Chemical Reactions , 2004 .

[17]  Ronald Breslow,et al.  Hydrophobic acceleration of Diels-Alder reactions , 1980 .

[18]  William L. Jorgensen,et al.  Investigation of solvent effects on pericyclic reactions by computer simulations , 1994 .

[19]  Kenzi Hori,et al.  Theoretical study on oxygen exchange accompanying alkaline hydrolysis of esters and amides. The role of water for the exchange reaction , 1999 .

[20]  Masataka Nagaoka,et al.  Structure optimization via free energy gradient method: Application to glycine zwitterion in aqueous solution , 2000 .

[21]  William L Jorgensen,et al.  Cope elimination: elucidation of solvent effects from QM/MM simulations. , 2006, Journal of the American Chemical Society.

[22]  徹 山口,et al.  量子化学計算を用いたモンテカルロシミュレーション（QM/MCシミュレーション）．酢酸メチルのアルカリ加水分解反応への適用 , 2008 .

[23]  Michel Dupuis,et al.  Two-dimensional free-energy surface on the exchange reaction of alkyl chloride/chloride using the QM/MM-MC method. , 2008, Physical chemistry chemical physics : PCCP.

[24]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[25]  Donald J. Cram,et al.  Room Temperature Wolff-Kishner Reduction and Cope Elimination Reactions , 1962 .

[26]  Kenji Hori,et al.  Hybridization-promoted and cytidine-selective activation for cross-linking with the use of 2-amino-6-vinylpurine derivatives. , 2005, The Journal of organic chemistry.

[27]  Donald J. Cram,et al.  Studies in Stereochemistry. XXXV. Mechanism of Ei Reaction of Amine Oxides , 1963 .

[28]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[29]  J. E. Mccarty,et al.  Studies in Sterochemistry. XXIV. The Preparation and Determination of Configuration of the Isomers of 2-Amino-3-phenylbutane, and the Steric Course of the Amine Oxide Pyrolysis Reaction in this System , 1954 .

[30]  Ronald Breslow,et al.  Diels-Alder reactions in nonaqueous polar solvents. Kinetic effects of chaotropic and antichaotropic agents and of .beta.-cyclodextrin , 1988 .

[31]  Hidetoshi Yamamoto,et al.  An Attempt Judging Possibility of Synthesis Routes from the TOSP Program Using DFT Calculations.: Application to Synthesis Routes for Substituted Fran-2,3-dione@@@置換フラン-2,3-ジオンへの適用 , 2005 .

[32]  William L. Jorgensen,et al.  PDDG/PM3 and PDDG/MNDO: Improved semiempirical methods , 2002, J. Comput. Chem..

[33]  K. Fukui.  The path of chemical reactions - the IRC approach , 1981 .

[34]  H. Bernhard Schlegel,et al.  An improved algorithm for reaction path following , 1989 .

[35]  Hendrik Zipse,et al.  A Quantum Mechanical and Statistical Mechanical Exploration of the Thermal Decarboxylation of Kemp's Other Acid (Benzisoxazole-3-carboxylic Acid). The Influence of Solvation on the Transition State Geometries and Kinetic Isotope Effects of a Reaction with an Awesome Solvent Effect , 1995 .

[36]  Masataka Nagaoka,et al.  Transition‐state optimization on free energy surface: Toward solution chemical reaction ergodography , 1998 .

[37]  Weitao Yang,et al.  QM/MM Minimum Free Energy Path: Methodology and Application to Triosephosphate Isomerase. , 2007, Journal of chemical theory and computation.

[38]  Hiroyuki Teramae,et al.  Theoretical study on the nature of the bonding and tautomerism in complexes with the tetrahydroborate anion , 1988 .

[39]  Jerry M. Parks,et al.  Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: sequential sampling and optimization on the potential of mean force surface. , 2008, The Journal of chemical physics.

[40]  Jacopo Tomasi,et al.  Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes , 1982 .

[41]  William L. Jorgensen,et al.  A Quantum Mechanical and Molecular Mechanical Method Based on CM1A Charges: Applications to Solvent Effects on Organic Equilibria and Reactions , 1998 .

[42]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[43]  Julian Tirado-Rives,et al.  Molecular modeling of organic and biomolecular systems using BOSS and MCPRO , 2005, J. Comput. Chem..

[44]  J. Stewart.  Optimization of parameters for semiempirical methods II. Applications , 1989 .

[45]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[46]  D. Kemp,et al.  Physical organic chemistry of benzisoxazoles. III. Mechanism and the effects of solvents on rates of decarboxylation of benzisoxazole-3-carboxylic acids , 1975 .