Abstract
The ease of separation by gas chromatography between petroselinic and oleic acids depends on the alcohol moieties of their esters. The esters of higher molecular weight alcohols tend to be better separated on a 90%-biscyanopropyl-10%-cyanopropylphenyl polysiloxane capillary column (30 m × 0.25 mm i.d.). By analysis of free energies contribution from different parts of the molecules, it is tentatively concluded that the interaction between the double bond and the column stationary phase is interfered by the bulky alkyl group, and it is the major driving force for the separation of the two fatty acids.
1. Introduction
Petroselinic acid (18:1Δ6cis) is a major fatty acid in the seed oils of most of the members of the Apiaceae (Umbelliflerae), Araliaceae, Garryaceae species, Acer (Aceraceae), and Geraniaceae [1–4]. Petroselinic acid is an interesting oleochemical for food, cosmetics, and pharmaceutical industries [5]. Coriander seed oil, which contained about 70% petroselinic acid, showed a significant hypolipidomic effect on animals [6]. High petroselinic acid diets found to lower the n-6 long-chain fatty acids of the phospholipids [7]. Hence, petroselinic has become more and more important in industry and human health. In 1997, Santinelli and Damiani [8] pointed out that petroselinic acid was always accompanied by some oleic acid (18:1Δ9cis), and there was no single chromatographic technique suitable for their separation and quantitation. However, methods for its separation from oleic acid have gradually been improved. Lin et al. [9] and Liu and Hammon [10] reported the separation of these two positional isomers by reversed phase liquid chromatography. Kim et al. [11] were able to separate and quantitate petroselinic and oleic acids in coriander tissues as their methyl esters by gas chromatography (GC) on a 5%-phenyl-95%-methylpolysiloxane capillary column. In addition, Thies [12] reported the use of butyl esters for the separation of these two positional isomers at 220°C with a total run time of 13 min. Thies saw good resolution of the butyl esters under these conditions, but they were not baseline separated. On the other hand, with a highly polar column (100%-cyanopropyl polysiloxane) of 50 meters, the two isomers could almost be baseline separated as their isopropyl esters [13, 14]. These two isomers in coriander seed oil could also be separated and quantitated on a 30 m of 100%-cyanopropyl polysiloxane with two steps temperature program [15], but how well they were separated was not reported. Near-baseline separations of methyl petroselinate and oleate on 100%-cyanopropyl polysiloxane columns from two different manufacturers at 180°C were reported [16]. A recent report [17] showed that methyl petroselinate and oleate could be separated on ionic liquid stationary phase [1,9-di(3-vinyl-imidazolium) nonane bis-(trifluoromethyl) sulfonyl imidate] columns, with the lengths of 12 (the microbore capillary), 30, and 100 meters.
A detailed study on the effects of alcohol chain lengths and branched chain on the separations of petroselinic and oleic acid esters was reported by Isbell [18]. Analyses were conducted on a 30 m of 100%-cyanopropyl polysiloxane column. The authors found that 2-ethyl-1-hexyl esters were promising derivatives for quantitative analysis of the oils containing the mixture of petroselinic and oleic acids. With the 2-ethyl-1-hexyl esters, the interaction of the polar carboxyl group with the stationary phase was sterically masked. Therefore, the main interaction would be that between the small polarity differences of the olefins and the polar stationary phase of the column, which caused the separation.
However, methyl esters of these two acids could also be separated on a 100 m of 100%-cyanopropyl polysiloxane capillary column [19].
In this study, oils containing petroselinic and oleic acids are transesterified with alcohols of different chain lengths. The gas chromatographic separations of these two positional isomers are discussed in term of their standard free energy of transfer from solution to gas (). Thus, a general knowledge of GC related to is briefly reviewed.
In gas-liquid chromatography, the solute (at infinite dilution) partitioning between the stationary and the mobile phases is assumed to be at equilibrium or very close to equilibrium state. Thus, the change in directly relates to its equilibrium constant (). The , on the other hand, can be divided into and [20] as described in the following: where is the number of carbon atoms. If , = , that is is the free energy contributed from the hypothetical molecule of zero carbon atom or may be simply understood as the functional group. Differentiating (1) with respect to , one obtains = . Thus, is the change in per carbon atom. In addition, equation (1) can be expanded as described below [21].
From basic thermodynamics, where and are the changes in standard enthalpy and entropy, respectively. is the absolute temperature, is the retention factor, β is column phase ratio, and R is the universal gas constant. Equation (3) is obtained by substitution of (2) into (1), or where and are the changes in standard enthalpy and entropy of transfer from solution to gas of an hypothetical molecule of , respectively. and are the increments in standard enthalpy and entropy of transfer from solution to gas per carbon atom; and are retention times of the solute and unretained gas, respectively
At constant , equation(4) is reduced to (6) or Vant’ Hoff’s equation. where It was reported that all the four numeric values of (4) (, , , and ) remained unchanged when the column was broken or cut into different lengths [22]. They were still valid for predicting of fatty acid methyl esters (FAMEs) with good accuracy. On the other hand, for columns of the same stationary phase but differing in the inside diameters, only the numeric value of was changed while the other three remained unchanged [23]. The β value is approximately equal to column inside diameter/(4 × film thickness), while the other coefficients are thermodynamic parameters and depend on type of solute and solvent (stationary phase) pair.
For linear aliphatic solutes (RH and RX), having the same number of carbon atoms but differing in the X functional groups, the difference in the retention indices of RzX and RzH is called homomorphic factor [24]. When the solutes are fatty acid alkyl esters (FAAEs) of different alcohols, the general structure would be R1COORy, where Rz = R1COO- and X = Ry. In this study, the homomorphic factors of fatty esters of different alcohols are evaluated in term of thermodynamic parameters. Furthermore, the free energy of transfer from solution gas of the double bond () is postulated, and the relative change in is responsible for the separations of 18:1∆6 and 18:1∆9 esters of different alcohols.
2. Experimental
2.1. Materials
Coriander seeds were obtained from a plant grower (Bangkok, Thailand). Olive oil was from local supermarket in Bangkok (Thailand). Fatty acid methyl esters (C16–C20) and oleic acid methyl ester were obtained from Sigma-Aldrich Chem. Co. (St. Loius, MO). C1–C3 alcohols and solvents were reagent grade obtained from Lab Scan Co. (Bangkok, Thailand). Higher alcohols (C4–C8) were from Sigma-Aldrich Chem. Co. (St. Loius, MO).
2.2. Extraction of Coriander Seed Oil
One gram of coriander seed kernel was ground in a mortar. The paste was transferred to a 4 mL vial, and 1 mL of toluene was added. The mixture was intermittently vortexed for 1 min. The toluene layer was decanted into a new vial and transesterified as described in the next section.
2.3. Transesterification
Transesterification of FAMEs with other alcohols (C2–C8) was carried out via a microreactor as described by Kaewkool et al. [25]. A disposable syringe (3 mL) was plugged with a small piece of cotton wool. NaOH (0.5 g) was rapidly ground and packed into the disposable syringe.
Five hundred microliters of a standard mixture of FAMEs (C16–C20) in toluene, about 2 mg/mL each, were mixed with 0.5 mL of the desired alcohol (C2–C8). Transesterification was started by passing the mixture through the microreactor gravitationally. Elution rate was controlled manually by the plunger, such that the mixture was eluted out in 30–45 s. Another 1 mL of the toluene-alcohol (1 : 1 v/v) mixture was added to wash the microreactor. The wash time was about 20 s. The combined eluent was acidified with 0.1 mL glacial acetic acid and washed with 1 mL of water. The organic layer was separated and dried over anhydrous Na2SO4. Transesterification and soap formation were checked by a liquid chromatography-based method according to Kittiratanapiboon and Krisnangkura [26].
For coriander seed oil an equal amount of olive oil, was mixed in order to increase the oleate peak in the chromatogam. The oil was similarly transesterifed as described above. For secondary octanols, alkaline catalyst was not suitable due to its slow reaction rate and soap formation. Acid catalysis was performed according to Kalayasiri et al. [27].
2.4. Gas Chromatography (GC)
Gas chromatographic analysis was carried out on a Shimadzu gas chromatograph model 2010 (Shimadzu Inc., Tokyo, Japan) and equipped with an FID, a split-splitless injector, a data processor (CBM 102), and a 90%-biscyanopropyl-10%-cyanopropylphenyl polysiloxane (Rtx-2330) capillary column (30 m × 0.25 mm i.d., 0.25 μm film thickness; from Restek International, Bellefonte, PA). The chromatographic conditions were as follows: Nitrogen carrier gas flow, 1 mL/min; nitrogen makeup gas flow, 30 mL/min; detector/injector, 230°C; split ratio, about 50 : 1. Oven temperatures were set isothermally at 5°C intervals between 180°C and 200°C, for determination of the column coefficients (, , , and ) and 180°C for the separation of petroselinic and oleic esters.
2.5. Determination of the Four Coefficients of (4) for Fatty Acid Butyl Esters (FABEs)
FABEs, which are approximately half way between methyl and octyl esters, were chosen as the references. The four numeric coefficients for FABEs were determined as described by Krisnangkura et al. [21].
At constant absolute temperature (), equation (4) is reduced to (8), the well-known Martin’s equation where Thus, plotting against would result in a straight line with a slope of and an intercept of . Plotting against would also result in a straight line with a slope of and an intercept of .
The mixture of FABEs (C16–C20) was chromatographed at 5°C intervals between 180°C and 200°C. The four column coefficients were solved by Microsoft Office Excel 2007.
2.6. Determination of the Four Coefficients of (4) for Fatty Acid Esters of Other Alcohols
The difference in retention indices of RX and RH (R = alkyl, X is functional group) is called homomorphic factor. According to (1), is the free energy contributed from the characteristic functional group, X. is the free energy contributed from the hydrocarbon chain of carbon atoms. Thus, the difference in between RX and RH would simply arise from the of the two molecules, which have the same alkyl, R, group. Similarly, two fatty acid esters having the same but differing in the alcohol moieties, the difference in of the two fatty acid esters would be the difference in . When (1) is expanded to (4), the numeric values of and would remain unchanged. Only the and , derived from , are changed. Therefore, the numeric values and of FABEs are assigned for other fatty acid esters. The values of and were obtained by solving two simultaneous equations, at different temperatures. The numeric values of and for each fatty acid ester are averaged from 25 data (5 different temperature pairs and 5 fatty acids).
3. Results and Discussion
3.1. The Homomorphic Factor of Fatty Acid Esters
As pointed out in Section 2.6, the homomorphic factor is the difference in retention indices of RX and RH. According to (1), two FAAEs (R1COORy) having the same R1 but differing in the alcohol moieties (Ry), would be characteristic of the homomorphic factor. In this section, attempts are made to prove that FAAEs of different alcohols would differ in , while the values are unchanged.
3.1.1. The Numeric Values for FAAEs
The four numeric coefficients for FAAEs were determined as described in Section 2.6 and solved by Microsoft Office Excel 2007. The numeric values and the standard deviations (sd) of 95% confidence for , , , and are listed in Table 1. Although, only 4 significant numbers are shown in Table 1, it should be pointed out that the higher decimal number would give better agreement between the predicted and experimental values. Hence, unrounded values of the 4 coefficients were used in the calculations.
The average values of and for FAAEs and their sd are summarized in Table 1. The values of are more negative as the carbon numbers of the alcohols are increased. The plot of and carbon numbers of alcohols yield a straight line with the slope of −0.580 and intercept of −6.265. The , sd of intercept, and slope are 0.941, 0.298, and 0.059, respectively. On the other hand, the values of increase positively as the carbon numbers of the alcohols are increased. The plot of and carbon numbers of alcohols give a straight line with the slope and intercept, of 382.1 and 416.7, respectively. The , sd of intercept and slope are 0.961, 157.7, and 31.22, respectively.
is the difference between and . The values of for each fatty acid ester at any temperature can be estimated from the coefficients and (in Table 1, with the nominal β value of 250). The value of , for each fatty acid ester, at constant temperature, increases negatively as the carbon numbers of the alcohols are increased (Figure 1). However, the transformation of a solute in solution to gas involves two physical phenomena. The interaction between the solute and the column stationary phase is measured by free energy of solution, (), and the vaporization of the solute is measured by free energy of vaporization, , [28, 29]. Therefore, the increment in with Ry may involve or or both of them, and the separation between 18:Δ6cis and 18:Δ9cis cannot be ascribed to at this moment.
3.1.2. Evaluation of Coefficients in (4)
The four numeric values (, , , and ) for FAAEs of different alcohols are summarized in Table 1. The validations of these values were indirectly verified by comparing the calculated retention times () with the experimental values (). Both and values are summarized in Table 2. All the are very close to the values. The greatest difference is 1.56% for methyl arachidate. There are only 3 values which differed greater than ±1.00% from the values. Forty one data points (21.0%) have the differences between 0.50 and 1.00. The rest (155 data points or 77.5%) have the differences less than 0.50%. The good agreement between the and values suggests that all the coefficients are valid for prediction retention times of FAAEs. The slopes, intercepts, and standard deviations of the plots between and for each FAAE are summarized in Table 3. The , values of all the FAAE are greater than 0.997, and the slopes (the ratios of ) are very close to 1.00, except that of FAME. The anomaly of the C1 series was observed on other phenomenon, such as gas holdup time [30, 31]. Hence, they are additional information, which supports that all the column coefficients are valid for predicting the retention times of FAAEs. Therefore, it may conclude that when the alcohols are changed, the numeric values of of and remain unchanged, accordingly the value does not change. It would also imply that the interaction between the saturated hydrocarbon chain and the stationary phase is not affected by the alcohol moieties. Thus, the homomorphic factor is verified for fatty acid esters of different alcohols.
3.2. GC Separation of Octadecenoic Acid Alkyl Esters of Different Alcohols and Free Energy Contribution from the Double Bond
Figure 2 is the aligned chromatograms of FAAEs (with different alcohols) of a mixture of coriander seed and olive oils, separated on a 90%-biscyanopropyl-10%-cyanopropylphenyl polysiloxane (Rtx-2330) capillary column at 180°C. The methyl and ethyl esters of petroselinic and oleic are coeluted. Partial separations of these two acid esters are observed when the alkyl groups have 3–6 carbon atoms. FAAEs of higher alcohols (FAHpE, FAOE, and FAEHE) are baseline resolved. Results agree well with the report of Isbell et al. [18]. In Figure 2, the octyl esters show approximately the same resolution as those of the 2-ethyl-1-hexyl esters in separations of 18:1∆6 and 18:1∆9. However, the 2-ethyl-1-hexyl esters are eluted faster than the octyl esters. The retention times of the 2-ethyl-1-hexyl esters are comparable to those of the heptyl esters. The bulky 2-ethyl-1-hexyl group was ascribed for lowering the interaction between the polar carboxyl and the polar stationary phase. Therefore, the interaction between the small polarity differences of the olefin and the polar stationary phase became apparent and caused the separation [18]. Although the discussion of Isbell et al. [18] was theoretically sound, it lacked the supported evidence. On the contrary, results in Figure 1 show that values increase as the alcohol groups are larger, which contradicted to the explanation of Isbell et al. [18]. Therefore, a new free energy term is proposed, and it is expected that it would be able to find the cause of the separation between petroselinic and oleic esters.
Martin [20] were the first to divide a molecule into different parts, and a free energy term was assigned to each of them as shown in (1). Subsequently, different free energy terms were assigned to describe interactions between column stationary phase and different functional groups of the molecule [32, 33]. Hence, the free energy concept of Kollie and Poole [32] and Golovnya [33] was extended, in this work, to investigate the interaction between the double bond and the column stationary phase. The new free energy is called free energy of transfer from solution to gas of the double bond (). Equation (10) is obtained by inserting the into(1) where the subscripts 1 and 2 stand for saturated and unsaturated fatty acids, respectively.
Equation (10) can also be derived as follows: Equation (12) is the difference between (11). Equation (10) is obtained by substituting (1) into (12).
is the difference between of 18:1Δ9 and of 18:1Δ6 as shown in . The will be used as a tool to investigate the interaction of the double bond with the stationary phase of the column. The of petroselinic and oleic (and vaccenic 18:1Δ11cis) acid esters at 180°C are summarized in Table 4. The C1–C3 esters are not well separated from each other, and they are not listed in the table. The values of these three esters at 180°C are approximately −151, −133, and 128 cal·mol−1, respectively. Partial separations are observed for FABEs. Separations of FAAE are gradually improved as the carbon numbers of the alcohols are increased.
According to Isbell et al. [18], the bulky alkyl group interfered the interaction of the carboxyl group to the column stationary phase. The interference did not extend to the hydrocarbon chain. Hence, the relative changes in should be the same for both saturated and unsaturated FAAEs. However, the values of petroselinic, oleic, and vaccenic esters (Table 4) become less negative as the alkyl sizes are increased, suggesting that the alkyl groups interfere the interaction between the hydrocarbon chain and the stationary phase, and the degree of interferences increases as the carbon numbers of the alcohols are increased.
The values of the butyl esters of 18:1∆6, ∆9, and ∆11 are −102.4, −120.1, and −145.4 calmol −1, respectively. The values become less negative as the alkyl groups are larger. 2-Ethyl-1-hexyl esters have the lowest compared with other alkyl esters. Their values for 18:1Δ6,Δ9, and Δ11 are −36.7, −61.4, and −90.5 cal·mol −1, respectively. Therefore, the interaction of the double bond with the column stationary phase increases when the double bond is far away from the carboxyl group, that is, interference from the bulky alkyl group decreases along the distant.
The difference in of 18:1Δ9 and 18:1Δ6 () of butyl, pentyl, hexyl, heptyl, octyl and 2-ethylhexyl esters is −17.7, −21.8, −23.5, −23.3, −24.5, and −24.7 cal·mol−1, respectively. The increases in as the alkyl groups increased suggest that the relative interactions of the double bonds of the two acid esters with the stationary phase are increased. Thus, it may conclude that is the driving force for the separation of the two isomers. However, the values level after C6, implying that the resolution of these two positional isomers may not be further increased by simply increasing the size of the alcohol.
The of 2-ethyl-1-hexyl esters is −24.7 cal·mol −1. It is estimated that a column of about 65,000 required plate numbers () would be able to baseline separate the two positional isomers, while the butyl esters, whose is −17.7 cal·mol −1, would require up to 180,000 plate numbers. The of 18:1∆9 and ∆11 have similar trend as the 18:1∆6 and 18:1∆9 pair, and discussion is not necessary.
Petroselinic acid has a double bond closer to the carboxyl group than oleic acid. The interference from the bulky alkyl group should be stronger than that to the oleic acid. Similar observation was reported by Kuningas et al. [34] that the separation of geometric isomers (alkenes) depended on the differences in free energy.
The the adjacent pairs of octadecenoic esters would reflect on the ease of separation for each pair. Results in Table 4 show that the separation of 18:1∆9 and 18:1∆11 is much simpler than the separation of 18:1∆6 and 18:1∆9.
3.3. Separation of Octadecenoic Acid Octyl Esters
As discussed in Section 3.2, the lowering in was ascribed to the reduction in the interaction between the double bond in the FAAEs and the column stationary phase. 1-Octyl and 2-ethy-1-hexyl esters are two of the best derivatives for separation of petroselinic and oleic acids. However, esters of secondary octyl alcohols have not been reported (to our knowledge). Thus, they were investigated for their interferences to the interaction of the double bond with the stationary phase. The retention times, retention factors, , and resolution (Rs) between secondary octyl petroselinate and oleate are summarized in Table 5. The values of secondary octyl esters of petroselinate and oleate range from −24.5 to −25.8 cal·mol−1. The resolutions of these two acid secondary octyl esters (1.15-1.26) are close to that of the 1-octyl ester (1.26). However, more symmetric or compact alcohol had the shorter retention time. The order of elution for the esters is 4-octyl, 3-octyl, 2-octyl, and 1-octyl, respectively.
4. Conclusion
When the alcohol moieties of FAAE are changed, the does not change. It implies that the interaction between the saturated hydrocarbon chain and the stationary phase is not affected by the alcohol moieties. On the other hand, petroselinic and oleic have one cis-double bond in the hydrocarbon chain; its interaction with the column stationary phase would differ from that of the saturated hydrocarbon. Therefore, is introduced and it would provide an insight into the interaction of the double with the column stationary phase. The bulky alkyl group attached to the carboxyl group affects the interaction between the double bond and the stationary phase. The steric effect decreases along the distance. In addition, for FAAEs having the same carbon numbers, the more compact is alcohol moiety, the shorter is the retention time, while the resolutions of the adjacent positional isomers are approximately the same. For petroselinic and oleic acids, the driving force for their separation is the difference in of the two geometrical isomers.
5. Symbols and Abbreviation
Fatty Acids
18:1Δ6cis: | Petroselinic acid |
18:1Δ9cis: | Oleic acid |
18:1Δ11cis: | Vaccenic acid. |
Thermodynamics
: | Standard free energy of transfer from solution to gas |
: | Standard enthalpy of transfer from solution to gas |
: | Standard entropy of transfer from solution to gas |
: | Free energy of transfer from solution to gas of a hypothetical molecule of zero carbon atom or simply called free energy of the functional gr |
: | Enthalpy of transfer from solution to gas of a hypothetical molecule of zero carbon atom |
: | Entropy of transfer from solution to gas of a hypothetical molecule of zero carbon atom |
: | Free energy of transfer from solution to gas of a double bond |
: | Change in free energy of transfer from solution to gas per carbon atom |
: | Change in enthalpy of transfer from solution to gas per carbon atom |
: | Change in entropy of transfer from solution to gas per carbon atom |
: | The difference in free energy of transfer from solution to gas of two unsaturated FAAE (Δ9 and Δ6). |
General
β: | Column phase ratio |
FAAE: | Fatty acid alkyl ester |
FABE: | Fatty acid butyl ester |
FAEE: | Fatty acid ethyl ester |
FAEHE: | Fatty acid 2-ethylhexyl ester |
FAHpE: | Fatty acid heptyl ester |
FAHxE: | Fatty acid hexyl ester |
FAME: | Fatty acid methyl ester |
FAPnE: | Fatty acid pentyl ester |
FAPrE: | Fatty acid propyl ester |
FAOE: | Fatty acid octyl ester |
: | Retention factor |
: | Universal gas constant |
Rs: | Resolution |
sd: | Standard deviation |
: | Retention time of a solute |
: | Hold up time |
: | Absolute temperature |
: | Carbon number. |
Acknowledgment
This work was supported by Thailand Research Fund.