Solvent Effects on Equilibrium Position and Initial Rate of Lipase-catalyzed Esterification Reactions in Organic Solvents: Experimental Results and Prediction Capabilities

The solvent effect on the equilibrium position and the initial rate of esterification of 1-hexanol with acetic acid catalyzed by a lipase has been experimentally investigated. A variety of non-polar and polar solvents have been considered and the results obtained indicate that the solvent effect on the equilibrium conversion is very important compared to that for transesterification reactions. A theoretically sound methodology using the group-contribution UNIFAC model for the prediction of solvent effects on the equilibrium position of enzymatic reactions is presented and it is applied to the reaction of 1-hexanol with acetic acid as well as to a similar reaction from the literature. The results obtained are better than those from empirical methods proposed in the literature such as correlations with the octanol-water partition coefficient of the solvent, as well as the solubility of water in the solvent. Moreover, the proposed methodology can be used for the determination of the equilibrium constant of the reaction. For the prediction of the solvent effect on the initial rate of enzymatic reactions it is found that it is more accurately determined using the product of the activities of the reactants, which can be predicted by the UNIFAC model, than the octanol-water partition coefficient of the solvent or the solubility of water in the solvent.     

Solvent Effects on Equilibrium Position and Initial Rate of Lipase-catalyzed Esterification Reactions in Organic Solvents: Experimental Results and Prediction Capabilities

The solvent effect on the equilibrium position and the initial rate of esterification of 1-hexanol with acetic acid catalyzed by a lipase has been experimentally investigated. A variety of non-polar and polar solvents have been considered and the results obtained indicate that the solvent effect on the equilibrium conversion is very important compared to that for transesterification reactions. A theoretically sound methodology using the group-contribution UNIFAC model for the prediction of solvent effects on the equilibrium position of enzymatic reactions is presented and it is applied to the reaction of 1-hexanol with acetic acid as well as to a similar reaction from the literature. The results obtained are better than those from empirical methods proposed in the literature such as correlations with the octanol-water partition coefficient of the solvent, as well as the solubility of water in the solvent. Moreover, the proposed methodology can be used for the determination of the equilibrium constant of the reaction. For the prediction of the solvent effect on the initial rate of enzymatic reactions it is found that it is more accurately determined using the product of the activities of the reactants, which can be predicted by the UNIFAC model, than the octanol-water partition coefficient of the solvent or the solubility of water in the solvent.     

Combining solvent engineering and thermodynamic modeling to enhance selectivity during monoglyceride synthesis by lipase-catalyzed esterification

Monoglyceride synthesis by Rhyzomucor miehei lipase was investigated via direct esterification between glycerol (adsorbed onto silica gel) and oleic acid in organic solvents. The main difficulty is to avoid the unwanted production of di- and tri-glycerides. It was demonstrated that an increase in solvent polarity, using mixtures of n-hexane and 2-methyl-2-butanol (2M2B), improves drastically the selectivity toward monoglyceride formation. In pure n-hexane, the monoglyceride represents only 6 molar % of the total products at the thermodynamic equilibrium (34 and 60% for di- and tri-glyceride respectively). Use of an equivolume mixture of n-hexane/2M2B enables a product mixture to be obtained containing 94% of monoglyceride at equilibrium (2.4 and 0% for di- and tri-glyceride respectively). This positive effect is counterbalanced by a decrease both in initial velocities and in substrate conversion at thermodynamic equilibrium.A modeling, able to predict the three thermodynamic equilibria governing the 3 consecutive reactions, based on activity coefficient calculations using the UNIFAC model, is proposed. It takes into account both the partition of water between solvent and immobilized catalyst, and the partition of glycerol between solvent and silica gel. A good correlation with experimental data obtained in n-hexane/2M2B mixtures was observed.     

The effect of organic solvents on the equilibrium position of enzymatic acylglycerol synthesis

The effect of organic solvents on the equilibrium position of lipase-catalyzed esterification of glycerol and decanoic acid has been investigated. The reaction is carried out in an aqueous-organic two-phase system. In polar solvents, high mole fractions of monoacylglycerol and low mole fractions of triacylglycerol and measured, while in nonpolar solvents, the measured differences in the mole fractions of monodi-, and triacylglycerols are less. There is a good correlation between the ester mole fractions at equilibrium and the log P of the solvent (partition coefficient in n-octanolwater), however, only if the group of tertiary alcohols is excluded. In the plot of the easter mole fractions as a function of the logarithm of hte solubility of water in the organic solvent, the tertiary alcohols can be included; however, in this case other deviations appear.For the prediction of the effect of organic solvents on the ester mole fractions at reaction equilibrium in nondilute reaction systems with a water activity below 1, the program TREP (Two-phase Reaction Equilibrium Prediction) is developed, which is based on the UNIFAC group contribution method. With this model the equilibrium data are essentially predicted from basic thermodynamic data. The required equilibrium constants are estimated from experiments without an organic solvent in the reaction medium. The mole fractions calculated by TREP show the same trends as the experimentally measured mole fractions; however, some variation is observed in the absolute values. These deviations may be due to inaccuracies in the UNIFAC group contribution method. TREP is found to be a correct method to predict within some limits the ester mole fractions at equilibrium for all mixtures of solvents, substrates, and products. The production of monoester can be enhanced in reaction system with a sufficient high concentration of a polar solvent. In experiments with a triglymeto-decanoic acid ratio of 5, almost no di-and triesters can be detected at equilibrium. (c) 1993 John Wiley & Sons, Inc.     

Enzymatic Reactions in Non-Conventional Media: Prediction of Solvent Water Content for Optimum Water Activity

Most enzymes provide their optimum performance at a given water activity (aw), which is generally solvent independent. For a given organic liquid solvent at a specific temperature or for a supercritical solvent at a specific temperature and pressure this corresponds to a water concentration in which water has the desired activity. We present here a methodology for predicting this water concentration thus reducing substantially the amount of experimental work needed to find the optimum solvent with respect to equilibrium conversion.

If the enzyme optimum water activity is known, the methodology predicts the required water content in the solvent to achieve this aw value. If, in addition, the enzyme water activity curve is available, this methodology provides the total water that must be added to the system (enzyme plus solvent) so that a specific water activity can be obtained.

The same methodology can also be applied to predict the effect of the total water content of the system (initial or initial plus produced) on the water activity values. It is shown that: (a) for esterification reactions taking place in hydrophobic organic solvents, the produced water can lead to a substantial change in water activity, but not for less hydrophobic solvents; (b) introduction of dry CO₂ into a system, pre-equilibrated to a certain water activity at atmospheric pressure, can lead to a substantial decrease in the water activity especially at temperatures just above the critical one of the solvent and pressures larger than that.     

Solvent selection for biocatalysis in mainly organic systems: predictions of effects on equilibrium position

Predictions may be made for the influence of solvent choice on the equilibrium position of biocatalyzed reactions, based on data for the liquid-liquid distribution of the reactants. The most reliable predictions are probably for dilute systems, based on partition coefficients or correlations derived from them. The effective equilibrium constant for esterification reactions is predicted to alter by more than four orders of magnitude on changing between different water-immiscible solvents. The equilibrium constant correlates well with the solubility of water in the solvent, and is most favorable for synthesis in the least polar solvents (aliphatic hydrocarbons). Similar effects seem to apply for other reactions, including oxidation of alcohols and hydrolysis of chlorides. Predictions can be made for nondilute systems using the UNIFAC system of group contributions, but the reliability of these is more questionable.     

Enzymatic Preparation of Mono- and Di-Stearin by Glycerolysis of Ethyl Stearate and Direct Esterification of Glycerol in the Presence of a Lipase from Candida Antarctica (Novozym 435)

Mixtures of 1(3)-monostearin and distearin were prepared by direct esterification of glycerol with stearic acid or transesterification using ethyl stearate as acyl donor in the presence of Candida antarctica lipase (Novozym 435) using a variety of solvents of differing polarity. In all cases, the transesterification resulted in higher product yields. In n-heptane as reaction medium the addition of water (3%) was essential for high product yields, with mono- and distearin being produced in almost equal amounts. Using more polar solvents as reaction media, such as acetonitrile or acetone, again the highest yields were obtained in the transesterification mode; employing these solvents the reactions were much more selective towards the formation of monostearin.     

The control of lipase-catalysed transesterification and esterification reaction rates. Effects of substrate polarity, water activity and water molecules on enzyme activity

The reaction rate of two lipase-catalysed reactions, esterification and transesterification, were studied in a liquid/solid two-phase system in order to investigate the effect of water partition between the enzyme preparation and the liquid phase composed of only the reactants, i.e. without the conventional solvents. Lipase from Candida cylindracea was used for these studies. The enzyme was inactive in dehydrated systems. In the case of monoester synthesis, the reaction rate increased with increasing water activity. The reaction rates of the non-specific C. cylindracea lipase-catalysed reactions were very sensitive to the nature of the substrates in this unusual system. For instance, the transesterification reaction rate of ethyl propionate was 48 times higher with nonanol than heptanol in the case of dehydrated substrates, but only 2.2 times higher in the case of water-saturated substrates. The results presented here demonstrate the absolute necessity to consider the polarity of every substrate, because of its ability to modify the water partition between the solid phase (enzyme preparation) and the liquid phase (substrate and product), which results in drastic changes in enzyme activity. Contrary to esterification, which is known to be activated by the water produced, the rate of transesterification remained constant at the beginning of the reaction. However, when transesterification and esterification were carried out in the same liquid phase, the transesterification reaction rate was controlled by the water produced by the concomitant esterification. Activation effects of the water molecules produced during the enzymatic reaction were of exactly the same order of magnitude for both reactions.     

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 DeltaG degrees 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.     

Enantioselective esterification of ibuprofen by a novel thermophilic Biocatalyst: APE1547

The enantioselective esterification of ibuprofen catalyzed by a novel thermophilic esterase (APE1547) from the archaeon Aeropyrum pernix K1 was successfully conducted in organic solvents. The effects of acyl acceptor, substrate ratio, organic solvent, temperature, and water activity were investigated. Under optimum conditions, the highest enantioselectivity (E = 38.1) was obtained with a higher enzyme activity (216.5 μmol/g/h). Celites were added into the reaction mixture to remove the water produced in the esterification. The reaction achieved its equilibrium in approximately 96 h with a conversion of 57 and 99% (ee) of the un-reacted (S)-ibuprofen obtained.     

Alginate as immobilization matrix and stabilizing agent in a two-phase liquid system: application in lipase-catalysed reactions.

Alginate was evaluated as an immobilization matrix for enzyme-catalyzed reactions in organic solvents. In contrast to most hydrogels, calcium alginate was found to be stable in a range of organic solvents and to retain the enzyme inside the gel matrix. In hydrophobic solvents, the alginate gel (greater than 95% water) thus provided a stable, two-phase liquid system. The lipase from Candida cylindracea, after immobilization in alginate beads, catalysed esterification and transesterification in n-hexane under both batch and continuous-flow conditions. The operational stability of the lipase was markedly enhanced by alginate entrapment. In the esterification of butanoic acid with n-butanol, better results were obtained in the typical hydrophilic calcium alginate beads than in less hydrophilic matrices. The effects of substrate concentration, matrix area, and polarity of the substrate alcohols and of the organic solvent on the esterification activity were examined. The transesterification of octyl 2-bromopropanoate with ethanol was less efficient than that of ethyl 2-bromopropanoate with octanol. By using the hydrophilic alginate gel as an immobilization matrix in combination with a mobile hydrophobic phase, a two-phase liquid system was achieved with definite advantages for a continuous, enzyme-catalysed process.     

Comparison of direct esterification and transesterification of fructose by Candida antartica lipase

Fructose oleates synthesis was performed in a batch reactor by trans- or direct esterification. An immobilized lipase from Candida antartica was used. When a solvent was used, 65% and 46% of conversion of fructose were obtained by transesterification and direct esterification, respectively. These two reactions were also compared in a solvent-free melt. Both in molten media and with cosolvent, two isomeric forms of fructose oleates were produced.     

A comparison of lipase-catalysed ester and lactone synthesis in low-water systems: analysis of optimum water activity.

We investigated the effects of the lyophilisation medium (enzyme plus buffer salt and additives) and of water activity (a(w)) on the catalytic properties of lipase from Chromobacterium viscosum (lipase CV) in organic solvents; catalysis of ester and lactone synthesis were compared and, despite the similarities of the reactive groups involved in these reactions, some interesting differences were observed. Including 2-[N-morpholino]ethanesulfonic acid (MES) buffer in the lyophilisation medium of lipase CV increased its catalytic activity in transesterification and lactonisation, although the buffer salt requirement for maximal activity differed between the two reactions. Sorbitol, glucose, lactose, 18-crown-6 (crown ether 18-C-6), beta-cyclodextrin and bovine serum albumin were employed as alternative additives in the transesterification reaction, but were not as effective as MES buffer. Salt hydrates were used to investigate the effect of a(w) on esterification and lactonisation reactions catalysed by lipase CV. The maximum rate of hexadecanolide synthesis in toluene occurred at a(w) = 0.48. The optimum a(w) for the transesterification reaction in heptane/alcohol mixtures depended on the alcohol substrate employed (1-heptanol, 2-heptanol, or 3-methyl-3-hexanol) but not on the acyl donor (p-NP acetate or caprylate). The optimum a(w) values for both reactions were unchanged when a common solvent system (toluene/1-heptanol) was employed, indicating that the dependence of enzyme activity on a(w) is an intrinsic property of the enzyme-catalysed reaction and not a function of the solvent or other additives.     

Computer-aided real-time estimation of reaction conversion for lipase-catalyzed esterification in solvent-free systems

Real-time conversion estimation through macroscopic balancing was investigated for enzymatic esterification reactions in a solvent-free system. In principle, the conversion of ester synthesis can be determined from the amount of water produced by the reaction because water is formed as a by-product in the same molar ratio as the product. In this study, we show that the water production rate, and thereby the reaction conversion, can be estimated on-line from measurements of the relative humidity of the inlet and outlet air and the material balances of water in the system. In order to test the performance of the real-time conversion estimation method, the lipase-catalyzed esterification reaction of n-capric acid and n-decyl alcohol in solvent-free media was conducted while controlling the water activity at various values. When the reaction conversions estimated on-line were compared with those analyzed off-line by gas chromatography, good agreement was obtained: the average mean absolute error was +/- 2.4% of the reaction conversion despite the simplicity of the method. The on-line estimation method presented here requires no expensive or complicated analytical instruments and no sampling of reaction medium. It can be used for monitoring nonaqueous enzymatic reactions where water is produced or consumed during reaction.     

Water activity control: a way to improve the efficiency of continuous lipase esterification

During continuous lipase-catalyzed oleic acid esterification by ethanol in n-hexane, the oleic acid conversion, initially at 95%, decreases to 20% after 2 h. This decrease is caused by the accumulation of the water produced in the course of the reaction in the packed-bed reactor (PBR). In order to improve the PBR efficiency, it is necessary to evacuate the water produced. In this study, different approaches have been tested to control the water content in the PBR during continuous esterification. The first approach consisted in improving the water solubility by increasing the reaction medium polarity. The addition of polar additives to n-hexane, the use of more polar solvents, and the use of solvent-free reaction medium were tested as a means to favor the water evacuation from the PBR. First of all, the use ofn-hexane supplemented with acetone (3 M) or 2-methyl-2-propanol (1 M) enabled the conversion to be maintained at higher values than those obtained in pure n-hexane. The replacement of n-hexane by a more polar solvent, like the 5-methyl-2-hexanone, resulted in the same effect. In all cases, conversions at steady-state were always less than 95%, as obtained in pure n-hexane. This is explained by a decrease in the enzyme activity due to the increase in the medium polarity. Nevertheless, an increase in enzyme quantity allowed 90% conversion to be maintained during 1 week using 3 M acetone amended n-hexane. Good results (a steady-state conversion of about 80%) were obtained when esterification was carried out in a solvent-free reaction medium containing 2 M 2-methyl-2-propanol as a polar additive. The second approach consisted in the evaporation of the accumulated water by use of an intermittent airflow. Although this process did not enable constant esterification rate to be maintained, it did enable the initial conversion (95%) to be restored intermittently.     

Solvent effects on the catalytic activity of subtilisin suspended in compressed gases

We studied a model transesterification reaction catalyzed by subtilisin Carlsberg suspended in carbon dioxide, propane, and mixtures of these solvents under pressure. To account for solvent effects due to differences in water partitioning between the enzyme and the bulk solvents, we measured water sorption isotherms for the enzyme in each solvent. We measured catalytic activity as a function of enzyme hydration and obtained bell-shaped curves with maxima at the same enzyme hydration (12%) in all the solvents. However, the activity maxima were different in all media, being much higher in propane than in either CO(2) or the mixtures with 50 and 10% CO(2). Considerations based on the solvation ability of the solvents did not offer an explanation for the differences in catalytic activity observed. Our results suggest that CO(2) has a direct adverse effect on the catalytic activity of subtilisin. (c) 1996 John Wiley & Sons, Inc.     

Effect of water activity on the lipase catalyzed esterification of geraniol in ionic liquid [bmim]PF6

Enzymatic reactions in non-aqueous media have been shown to be effective in carrying out chemical transformation where the reactants are insoluble in water or water is a byproduct limiting conversion. Ionic liquids, liquid organic salts with infinitesimal vapor pressure, are potentially useful alternatives to organic solvents. It is known that the thermodynamic water activity is an important variable affecting the activity of enzymes in non-aqueous solvents. This study investigated the influence of water activity on the esterification of geraniol with acetic acid in ionic liquid [bmim]PF6 catalyzed by immobilized Candida antarctica lipase B. The conversion of geraniol in [bmim]PF6 was significant although the reaction rate was slower than in organic solvents. The profile of initial reaction rate-water activity was determined experimentally, and differed from the data reported for other non-aqueous solvents. A maximum in the initial reaction rate was found at aw = 0.6. The pseudo reaction equilibrium constant, Kx, was measured experimentally for the reaction. The average value of Kx in [bmim]PF6 was 12, 20-fold lower than the value reported for the same system in hexane.     

Solvent effects on biocatalysis in organic systems: equilibrium position and rates of lipase catalyzed esterification

Porcine pancreatic lipase immobilized on celite particles has been employed as a catalyst for the esterification of dodecanol and decanoic acid in a predominantly organic system. Solvent influence on the equilibrium position and on the catalyst activity has been studied using 20 solvents, including aliphatic and aromatic hydrocarbons, ethers, ketones, nitro- and halogenated hydrocarbons, and esters. The equilibrium constant for esterification correlates well with the solubility of water in the organic solvent, which in turn shows a good relationship with a function of Guttman's donor number and the electron pair acceptance index number of the solvent. This may be rationalized in terms of the requirements for solvation of water and of the reactants. The catalyst activity, measured as the initial rate of the esterification reaction, is best correlated as a function of both n-octanol-water partition coefficient (log P) and either the electron pair acceptance index or the polarizability.     

Enzyme reaction kinetics in organic solvents: A theoretical kinetic model and comparison with experimental observations

A theoretical kinetic model has been developed in order to describe the enzyme reaction in organic solvents. In this model the hydration of the enzyme molecule was examined and the equilibrium kinetic constants expressed in terms of thermodynamic activity. Analysis of a proposed kinetic model shows that the enzyme reaction rate in organic solvents is determined by two factors: substrate solvation and enzyme hydration, which are determined by the activity coefficient of the substrate and the water activity of the reaction media, respectively. The activity coefficient of the substrate and the water activity have been calculated using the UNIFAC equation to analyze the effects of organic solvents on the rate of enzyme reaction, and the results were compared with experimental data. Predictions of the proposed model were found to be in good agreement with previous experimental observations.     

Enzymatic transesterification of olive oil and its precursors

The effect of different solvents and three different acyl acceptors on the transesterification of triolein (as a model compound) was investigated. The yield of biodiesel (methyl or ethyl ester) was monitored as a function of time. The yield of the product was also determined in a solvent-free system for two different modes of stirring. The results indicate that the highest yield is obtained in a solvent-free system with mechanical stirring. Methyl acetate is also effective as a solvent and acyl acceptor. Biodiesel was also produced by transesterification of triglycerides (triolein) present in olive oil with methanol and Novozym® 435. The effect of the molar ratio of methanol to triolein, mode of methanol addition, enzyme activity and reaction temperature on overall conversion and yield was determined. The final conversion and yield of biodiesel after a reaction time of 24 h were unaffected by changes in these parameters over the range studied. Preliminary findings indicate that the results obtained from small scale reactors and fresh oil can be extended to larger reactors and used oil.