In Silico Prediction of Medium Effects on Esterification Equilibrium Using the COSMO‐RS Method

This paper presents a new approach for predicting solvent effects on esterification reactions of industrial importance in the field of biocatalysis. The COSMO‐RS method has been used to calculate the activity coefficients of the chemical species involved in various reactions, carried out in different solvents. For comparison we also used the traditional UNIFAC method. Three lipase‐catalyzed esterifications were considered: (1) 1‐dodecanoic acid with menthol in n‐hexane, n‐heptane, cyclohexane, 2,2,4‐trimethylpentane, toluene, acetonitrile, and 2‐methyl‐2‐butanol; (2) 1‐dodecanoic acid and 1‐dodecanol in n‐hexane, n‐heptane, cyclohexane, 2,2,4‐trimethylpentane, and toluene; and (3) glycerol and n‐octanoic acid in acetonitrile, benzene, and toluene and in the neat reaction mixture (without any solvent). Predicted activities were used to calculate the thermodynamic equilibrium ratio. This should be independent of medium, and the variation in COSMO‐RS values is at most 9‐fold (corresponding to a ΔG° of about 5.5 kJ/mol, which would still be a very useful prediction) and often only 2‐fold (corresponding to less than 2 kJ/mol or 0.5 kcal/mol, therefore comparable with experimental error). UNIFAC is weaker, especially when important roles are played by conformational freedom, intramolecular interactions, strong polar effects, and charge distribution of molecules in the mixture. The relative percent deviations from the mean of equilibrium constants in different solvents range between 17 and 49 for COSMO‐RS versus 32 to 65 for UNIFAC. The COSMO‐RS method opens up new perspectives for the development of theoretical models for solvent selection with general applicability.

[1]  H. Hasse,et al.  Thermodynamics of phase and chemical equilibrium in a strongly nonideal esterification system , 2005 .

[2]  A. Klibanov,et al.  Insights into the Solvent Dependence of Chymotryptic Prochiral Selectivity , 1998 .

[3]  Andreas Klamt.  COSMO and COSMO-RS , 2002 .

[4]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[5]  A Zaks,et al.  Industrial biocatalysis. , 2001, Current opinion in chemical biology.

[6]  A. Klamt,et al.  Prediction of infinite dilution activity coefficients using COSMO-RS , 2003 .

[7]  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 .

[8]  A. Klibanov,et al.  Prediction of the Solvent Dependence of Enzymatic Prochiral Selectivity by Means of Structure-Based Thermodynamic Calculations , 1996 .

[9]  Alexander M. Klibanov,et al.  Predicting the solvent dependence of enzymatic substrate specificity using semiempirical thermodynamic calculations , 1993 .

[10]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

[11]  D. Bunk,et al.  Thermodynamics of the lipase-catalyzed esterification of glycerol and n-octanoic acid in organic solvents and in the neat reaction mixture , 2001 .

[12]  A. Klibanov.  Improving Enzymes by Using them in Organic Solvents , 2001 .

[13]  Y. Dudal,et al.  Influence of water activity on the synthesis of triolein catalyzed by immobilized Mucor miehei lipase , 1995, Biotechnology and bioengineering.

[14]  Y. B. Tewari.  Thermodynamics of the lipase-catalyzed transesterification of (−)-menthol and dodecyl dodecanoate in organic solvents , 2000 .

[15]  A. Klamt.  Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .

[16]  C. Ebert,et al.  Influence of organic solvents on enzyme chemoselectivity and their role in enzyme-substrate interaction , 1996 .

[17]  Stanley I. Sandler,et al.  A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model , 2002 .

[18]  S. Grimme,et al.  An improved method for density functional calculations of the frequency-dependent optical rotation , 2002 .

[19]  D. Vanderah,et al.  Thermodynamics of the Lipase-Catalyzed Esterification of 1-Dodecanoic Acid with (−)-Menthol in Organic Solvents , 1998 .

[20]  A. Klamt,et al.  Refinement and Parametrization of COSMO-RS , 1998 .

[21]  C. Akoh,et al.  Enzymatic synthesis of geraniol and citronellol esters by direct esterification in n-hexane , 1993, Biotechnology Letters.

[22]  A. Kiener,et al.  Industrial biocatalysis today and tomorrow , 2001, Nature.

[23]  F. Molinari,et al.  Production of geranyl acetate and other acetates by direct esterification catalyzed by mycelium of Rhizopus delemar in organic solvent , 2004, Biotechnology Letters.

[24]  P. Halling,et al.  Solvent selection for biocatalysis in mainly organic systems: Predictions of effects on equilibrium position , 1990, Biotechnology and bioengineering.

[25]  P. Halling,et al.  Solvation of CBZ–amino acid nitrophenyl esters in organic media and the kinetics of their transesterification by subtilisin , 1994, Biotechnology and bioengineering.

[26]  A. Janssen,et al.  Solvent effects on lipase‐catalyzed esterification of glycerol and fatty acids , 1993, Biotechnology and bioengineering.

[27]  Maurizio Fermeglia,et al.  Prediction of phase equilibria for binary mixtures by molecular modeling , 2001 .

[28]  A. Klamt,et al.  Fast solvent screening via quantum chemistry: COSMO‐RS approach , 2002 .

[29]  J. Jongejan,et al.  Do organic solvents affect the catalytic properties of lipase? Intrinsic kinetic parameters of lipases in ester hydrolysis and formation in various organic solvents , 1995, Biotechnology and bioengineering.

[30]  V. Nagarajan,et al.  Biocatalysis: applications and potentials for the chemical industry. , 2002, Trends in biotechnology.

[31]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[32]  P. Halling,et al.  Solvent effects on biocatalysis in organic systems: Equilibrium position and rates of lipase catalyzed esterification , 1991, Biotechnology and bioengineering.