Journal of biological and health sciences http://biotecnia.unison.mx

Universidad de Sonora

ISSN: 1665-1456

Original Article

Synthesis of acylglycerols structured with n-3 PUFA in a packed bed reactor with recirculation: Op-timization by response surface

Síntesis de acilgliceroles estructurados con AGPI n-3 en un reactor de lecho empacado con recirculación: Optimización por superficie de respuesta

Juan Antonio Noriega Rodríguez1*

, Cristian Correa-Leyva1 , Esther Carrillo Pérez1 and Hugo S. García2

1 Universidad de Sonora. Sonora, México.

2 Instituto Tecnológico de Veracruz. México.

ABSTRACT

The enzymatic esterification of n-3 PUFA and glycerol (G), catalyzed by a Candida antarctica lipase, was studied in a recirculating packed-bed reactor for the synthesis of acyl- glycerols. An n-3 PUFA concentrate was prepared by chemi- cal hydrolysis of Menhaden’s oil followed by urea treatment. A rotatable central composite design was used to evaluate the effect of molar ratio (0.47 – 5.52 mol n-3 PUFA / mol G), temperature (28.14 – 71.86°C) and time (0.24 – 2.76 h), on the production structured acylglicerols. The analysis of variance shows that all principal factors have a significant effect (p < 0.05). It was determined through response surface methodo- logy that around 80 % of global esterification can be reached, by operating with a molar ratio of 0.5 mol PUFA/mol G, 70°C and 2.75 h. In general, models had a good adjust to experi- mental data (R2 > 0.92), and the optimal operating conditions for the formation of the different acylglycerols produced during the esterification (monocylglycerols, diacylglycerols and triacyglicerols) can be established.

Keywords: Esterification, n-3 polyunsaturated fatty acids,

immobilized lipase, HPTLC, regression models.

RESUMEN

Se estudió la esterificación enzimática de AGPI n-3 y glicerol

(G) catalizada por una lipasa de Candida antarctica en un reactor de lecho empacado con recirculación para la síntesis de acilgliceroles estructurados. Se preparó un concentrado de AGPI n-3 por hidrólisis química de aceite de Menhaden seguido de un tratamiento con urea. Se utilizó un diseño central compuesto rotable para evaluar el efecto de la rela- ción molar (0.47 – 5.52 mol AGPI n-3 / mol G), la temperatura (28.14 – 71.86 °C) y el tiempo (0.24 – 2.76 h) en la producción de acilgliceroles. El análisis de varianza mostró que todos los factores principales de los modelos presentaron un efecto significativo en las diferentes respuestas (p < 0.05). Mediante metodología de superficie de respuesta se determinó que se puede alcanzar una esterificación global cercana a 80 % al operar con una relación molar de 0.5 mol AGPI/mol G, 70

°C y 2.75 h de reacción. En general, los modelos tuvieron un buen ajuste a los datos experimentales (R2 > 0.92) y permiten establecer las condiciones de operación óptimas para la for- mación de los diferentes acilgliceroles producidos durante la esterificación (monocilgliceroles, diacilgliceroles y triacilgli- ceroles).

*Author for correspondence: Juan Antonio Noriega Rodríguez e-mail: juan.noriega@unison.mx

Received: August 12,2024

Accepted: November 21, 2024

Published: December 17, 2024

Palabras clave: Esterificación, ácidos grasos poliinsaturados n-3, lipasa inmobilizada, HPTLC, modelos de regresión.

INTRODUCTION

N-3 polyunsaturated fatty acids (PUFA), usually present in fish oil, play several important health-related functions and are considered essential because the human body cannot synthesize them (Khan et al., 2023, Hernández et al., 2016; Baeza et al., 2014). Different scientific studies indicate that the consumption of n-3 PUFA reduces the risk for coronary heart disease, neuromuscular diseases, allergies, diabetes, cancer and depression (Kołodziej et al., 2023; Shahidi and Ambigaipalan, 2018; Rimm et al., 2018; Damasco et al., 2018; Hallahan et al., 2016; Wang and Huang, 2015; Castellanos and Rodríguez, 2015). The great benefits of n-3 PUFA consump- tion have motivated the study of the formation of structured acylglycerols (AG) enriched with these fatty acids for nutra- ceutical and pharmaceutical purposes. Enzymatic processes have multiple advantages over conventional chemical processes, such as operation under mild temperature and pH conditions, reduced energy cost, high catalytic efficiency, and a wide range of fatty acid selectivity of lipases leading to specific and pure products. One method for the production of FAs enriched with n-3 PUFAs is through the direct enzy- matic esterification reaction of PUFAs and glycerol (Correa et al., 2017). During the stages of the esterification reaction, different products with high added value are produced such as monoacylglycerols (MAG), diacylglycerols (DAG) and tria- cylglycerols (TAG). MAGs are recognized for being excellent emulsifiers used for food, cosmetic and pharmaceutical applications. The consumption of oils consisting mainly of DAG contributes to weight loss, reduction of body fat, lower liver fat content, improvement of lipid metabolism by favo- ring oxidation and decreased storage of TAG in the adipose tissue (Nurul Nadiath et al., 2020).

1

The operation in packed bed reactors (PBR) presents mul- tiple advantages compared to batch reactors, for example, it facilitates the contact and separation of the enzyme with the substrates (Zhang et al., 2022), allows the reuse of the enzy- me, can be operated in continuous mode, and is usually more economical (Lee et al., 2016). Furthermore, immobilizing the enzyme provides it with greater catalytic and thermal stabi- lity, allowing a longer useful activity time. The immobilized

Volume XXVI

DOI: 10.18633/biotecnia.v26.2414

lipase from Candida antarctica has shown good stability at high temperatures and long operating times (Monteiro et al., 2021). The development of a continuous process in an PBR is the most appropriate to carry out reactions catalyzed by lipa- ses involving fats and oils for different industrial applications (Hong, 2014). However, when the substrates injected into the PBR are immiscible, such as n-3 PUFA and glycerol (G), the efficiency of the continuous process decreases considerably. When there is little miscibility in the substrates, an operation mode with recirculation can be chosen for the effective exe- cution of the process (Zhao et al., 2012).

Regression analysis using an experimental design, is the construction of predictive models that adequately describe the relationship between the dependent and independent variables of a study process. When the variables of a bioprocess are measurable, continuous and controllable during the experimental runs, a response surface methodology (RSM) study can be carried out considering a rotatable central composite design (RCCD) to determine the effect of the variables and optimize bioprocess responses (Filli et al., 2010). RSM is a collection of statistical and mathematical techniques useful for developing, improving and optimizing processes. The most extensive applications are particularly in situations where several input variables influence performance measures or quantity or quality characteristics of a product or a process (Montgomery, 2019).

In the present work, the effect of temperature, molar ratio and time on the direct enzymatic esterification reaction of n-3 PUFA and G, catalyzed by the immobilized lipase of Candida antarctica (NV-435) in a PBR with recirculation, was studied by means of a series of experiments based on a RCCD and analyzed using response surface methodology to optimize the production of structured acyglycerols.

MATERIAL AND METHODS

Reagents

Menhaden fish oil (Sigma-Aldrich) and food grade glycerol (J. T. Baker) were used for the reaction medium; tertbutylhydro- xyquinone (TBHQ, Química Dresen) as an antioxidant, and di- melsulfoxide (DMSO) as an emulsifier. The NV-435 lipase from Candida antarctica (Novo Nordisk Denmark) immobilized in macroporous acrylic resin was used as catalyst. For high- performance thin layer chromatography (HPTLC) analysis, 10x10 cm plates (Uniplate) were used with analytical grade silica gel, hexane and ethyl ether (Sigma-Aldrich), glacial acetic acid (Fermont) and reagent grade iodine (Fermont).

Isolation and concentration of PUFA

To isolate and concentrate PUFAs from Menhaden oil, the procedure described by Correa et al. (2017) was used. One hundred g of oil were weighed and TBHQ (0.02 % w/w) was added as an antioxidant. The mixture was saponified with 200 mL of 7 M KOH in 70 % ethanol by heating at reflux for one hour. A liquid extraction was carried out with distilled water (240 mL) and hexane (200 mL), the aqueous phase was

then acidified to pH 1.0 with HCl to subsequently extract the free fatty acids (FFA) with hexane (200 mL). The hexane was removed by evaporation in a rotovapor (40 rpm, 40°C, vacuum 25 mmHg), and the remaining moisture was remo- ved with anhydrous sodium sulfate. In a flask, each 15 g of FFA were mixed with 25 g of urea in 100 mL of 95 % ethanol and heated until the mixture turned a homogeneous yellow color. The solution was transferred to centrifuge tubes and quickly cooled in an ice water bath for subsequent cooling (4°C, 8 h). The crystals formed were separated by centrifuga- tion at 6000xg (15 min, 0°C). The supernatant was acidified to pH=4 and transferred to a separatory funnel to extract the PUFA with a 1:1 hexane:water mixture, stirring for 30 min. Finally, the organic phase was evaporated in a rotary evapo- rator (40 rpm, 40°C, vacuum 25 mmHg) to obtain the PUFA concentrate.

Enzymatic esterification in PBR with recirculation

The experimental runs to study the enzymatic esterification reaction were carried out based on a RCCD. Mixtures of 5 g with different molar proportions of n-3 PUFA and G (0.48, 1.5, 3, 4.5, 5.52 mol/mol), TBHQ (0.02 % w/ p) as antioxidant, and DMSO (1 % w/w) as emulsifier, were prepared in glass vials (20 mL) with magnetic stirrer. Each vial was placed in- side a double arm cooling jacket reactor (Wheaton Celstir®) placed on a magnetic stirring plate at 800 rpm. After stirring the mixture for 30 minutes, a peristaltic pump was used to recirculate the medium at a flow rate of 1.65 mL/min through an PBR (length = 5.104 cm and diameter = 0.717 cm) with 0.9115 g of immobilized lipase NV- 435 of Candida antarctica. A water recirculator bath was used to maintain an isother- mal process, experimental runs were carried out at different operating temperatures (28.14, 37, 50, 63, 71.86 °C) and 35 µL samples were taken at different reaction times (t = 0.24 , 0.75, 1.5, 2.25, 2.76 h). A representative diagram of the process is shown in Figure 1.

Analysis of acylglycerols (AG) and global esterification The analysis was carried out in a HPTLC equipment (CA- MAG). Aliquots (5 µL) of the sample were diluted in a chloroform:methanol solution (2:1 v/v) to obtain a final concentration of 20 µg AG/ µL sample. One microliter of each diluted sample was applied on a silica gel plate 1 cm from the edge and with a separation of 1 cm between each sample. The plate was saturated with the mobile phase of hexane:ethyl ether (38:62 v/v) with 1 % glacial acetic acid in a vertical double-channel glass chamber for thin layer chroma- tography by 40 minutes. The plate was placed into the mobile phase for chromatographic development until a path of 9 cm was reached. The plate was then removed from the chamber and the solvent was evaporated into the air. The plate was developed by atomizing a solution of iodine in ethanol (0.5

% w/v), until achieving a light yellow color on the surface of the plate and the solvent was cold evaporated for 30 min. The developed plate was scanned, and shadow area analysis of the spot bands was performed using ImageJ 1.45 software

cording to the significative variables (α = 0.05) resulted from the analysis of variance. Statistical analysis was performed in a JMPpro v17.0 ® platform.

RESULTS AND DISCUSSION

‌Table 1 records the experimental results of the enzymatic es- terification of n-3 PUFA and G from the different experimental runs established for the RCCD. The coefficient of determina- tion (R2) is a value that varies from 0 to 1 that expresses the prediction of a model to the experimental data, with a value close to 1 indicating a better prediction, however, in these

adj

models the value of R 2 is usually considered, which is a

corrected measure of goodness of fit to determine the preci- sion of the model that identifies the percentage of variance in the target field that is explained by the model inputs. Gu- tierrez and De la Vara (2012) suggest that an R2 greater than

0.75 is adequate to determine a good fit of the model to the

experimental data. In this study, R2 greater than 0.92 and R

2

adj

Figure 1. Operation diagram of the enzymatic esterification reaction in the packed bead reactor (PBR) with recirculation. 1: jacketed PBR, 2: jacketed double sidearm Celstir®, 3: magnetic stirrer, 4: peristaltic pump and 5: recir- culator water bath with temperature control.

Figura 1. Diagrama de operación de la reacción de esterificación enzimática en el reactor de lecho empacado (RLE) con recirculación. 1: RLE enchaqueta- do, 2 Celstir ® de doble brazo enchaquetado, 3: agitador magnético, 4: bom- ba peristáltica y 5: baño recirculador de agua con control de temperatura.

(NIH; Rasband, 2018). For the quantification of MAG, DAG and TAG, a standard curve was generated using lipid standards (Sigma-Aldrich) containing a mixture of cis-9-monoolein, cis-

values greater than 0.74 were obtained. Different bioprocess investigations have obtained R2 values similar to those repor- ted in this study (Pereira et al., 2018). The statistical analysis showed that the models had a significant effect with respect to the different responses (p < 0.05).

.

The F statistic was performed in the ANOVA, high F values were obtained, higher than the critical values, which indicate that the models are statistically significant (Nahe- miah et al., 2015). The results of the statistical analysis of the models are recorded in Table 2. To determine the fit model for each of the responses, the least significant parameters (p >

9-1,2-diolein, cis-9-1,3-diolein and cis-9-triolein.

The G concentration was calculated using a stoichio- metric mass balance considering the concentrations of

0.2) were eliminated to optimize the value of R

Global esterification (GE)

2

adj

MAG, DAG and TAG, determined by HPTLC analysis using the following equation:

𝐶𝐶𝐺𝐺 = 𝐶𝐶𝑜𝑜𝐺𝐺 − (𝐶𝐶𝑀𝑀𝐴𝐴𝐺𝐺 + 𝐶𝐶𝐷𝐷𝐴𝐴𝐺𝐺 + 𝐶𝐶𝑇𝑇𝐴𝐴𝐺𝐺) (1)

To determine the global esterification (GE), 30 µL ali- quots were taken and diluted in ethanol. The samples were titrated with a 0.05 N NaOH solution using phenolphthalein as an indicator.

The parameters with the greatest effect on the GE are the linear and quadratic terms of the molar ratio and reaction time (p < 0.005). In the analysis of variance, the effect of temperature did not have a representative effect on GE (p = 0.3954), but it was included in the equation to have a hierarchical model because the interaction term between temperature and time had a significant effect (p < 0.1). The model to describe the GE obtained through second-order polynomial regression is the following:

Experimental design and statistic analysis

A RCCD was selected to evaluate the effect of the molar ratio

GE = 3.019 M2 - 6.157 t2 - 2.647 M t + 0.223 T t - 15.369 M

+ 21.593 t - 0.180 T + 48.124

(3)

between PUFA and G (x1=M), temperature (x2=T) and time (x3=t), on GE (y1), MAG (y2); DAG (y3) and TAG (y4) production.

To predict the dependent variables (yi) a third order

polynomial regression model was considered:

The RSM graphs show that a GE close to 80 % can be achieved by operating with a molar ratio of 0.5 mol PUFA/ mol G and a temperature of 70°C for 3 h (Fig. 2a). It is im- portant to clarify that by decreasing the molar ratio (M), n-3

𝑖𝑖

𝑘𝑘

𝑦𝑦 = 𝛽𝛽

+ ∑ 𝛽𝛽 𝑥𝑥𝑖𝑖 + ∑ 𝛽𝛽 𝑥𝑥2 + ∑ 𝛽𝛽

0

𝑖𝑖=1

𝑖𝑖

𝑘𝑘

𝑖𝑖=1

𝑖𝑖𝑖𝑖 𝑖𝑖

𝑘𝑘

𝑖𝑖=1

𝑖𝑖𝑖𝑖𝑖𝑖 𝑖𝑖

𝑘𝑘−1 𝑘𝑘

𝑥𝑥3 + ∑ ∑ 𝛽𝛽 𝑥𝑥𝑖𝑖𝑥𝑥𝑗𝑗 + 𝜀𝜀

𝑖𝑖=1 𝑖𝑖<𝑗𝑗

𝑖𝑖𝑗𝑗

(2)

PUFAs have more binding sites in the G structure due to the surplus of this substrate considering that the stoichiometric ratio is 3:1 mol PUFA/mol G, therefore, it is normal to obtain higher degrees of esterification. Linder et al. (2005) studied

Where β0, βi, βii, βiii, and βij represent the interaction co-

efficients in the regression, the linear effect, quadratic effect,

‌cubic effect and the combined effects, respectively.

A simplified model was considered for each response ac-

the enzymatic esterification reaction of PUFA and G in a batch reactor, and obtained a higher GE at a molar ratio of

0.83 mol PUFA/mol G at 46°C. It is important to highlight the high reaction rate that occurs in the PBR, since a high GE was

Table 1. Experimental and predicted results of the enzymatic esterification of n-3 PUFAs and G in the PBR with recirculation.

Tabla 1. Resultados experimentales y predichos de la esterificación enzimática de los AGPI n-3 y G en el RLE con recirculación.

M= molar ratio (mol PUFA / mol G), T= Temperature (°C), t= time (h), GE = Global esterification (%), MAG % = Monoacylglycerols ( % mol), DAG % = Diacylglycerols (% mol), TAG % = Triacylglycerols (% mol).

* The percentage molar fraction was calculated based on the total fatty acids.

Table 2. Statistical terms obtained in the statistical analysis of the GE, MAG

%, DAG % and TAG % models.

Tabla 2. Términos estadísticos obtenidos en el análisis estadístico de los mo- delos de GE, % MAG, % DAG y % TAG.

GE = Global esterification (%), MAG % = Monoacylglycerols (% mol), DAG %

n-3 PUFA and G at different temperatures in a batch reactor, and it was observed that an increase in temperature raised the GE reached during the reaction relatively little.

Synthesis of monoacylglycerols

MAGs are the first product formed in the enzymatic esterifi- cation reaction, and as the reaction proceeds, the formation of DAG consumes the MAGs present in the reaction medium. The parameters with the greatest effect on the percentage molar fraction of MAG (% MAG) are the linear term of the molar ratio, and the quadratic terms of the molar ratio and temperature (p < 0.005). The statistical model that describes the MAG % is:

= Diacylglycerols (% mol), TAG % = Triacylglycerols (% mol).

* The percentage molar fraction was calculated based on the total fatty acids.

MAG % = 1.876 M2 + 0.018 T2 + 0.093 M T - 0.299 T t - 19.099 M

-1.585 T + 14.251 t + 107.631

(4)

obtained at relatively short operating times, unlike different studies in intermittent reactors where long reaction times are required to achieve a significant GE (Noriega et al., 2013; Wang et al., 2011). The RSM of temperature versus molar ratio provides very interesting information (Fig. 2b). Firstly, the model suggests that at higher molar ratios (M > 4) the GE can increase a little, however, it is not advisable to operate under these conditions because PUFA are the limiting reagent in es- terification, and these are the most expensive substrate of the process due to its laborious obtaining process (Torres et al., 2014). Secondly, it can be seen that temperature is not as sig- nificant a variable in the GE as time and molar ratio, obtaining similar GE at different operating temperatures in relation to the proportions of the substrates in the medium. Correa et al. (in press) studied the kinetics of enzymatic esterification of

Figure 3 shows the RSM graphs of the MAG %. The model suggests that the MAG % can be optimized, reaching 80 %, by operating with a molar ratio of 0.5 mol PUFA/mol G and 30°C for 2.75 h. Byun et al. (2007) studied the forma- tion of MAG in the enzymatic esterification of fish oil, and determined that, at low operating temperatures, greater MAG formation occurs. Similarly, Zhao et al. (2011) sought to optimize the formation of MAG in an enzymatic esterification reaction without solvents, and obtained conversions of 49.6

% at a molar ratio of 0.16 mol AGL/mol G and a temperature of 50°C. Figure 3b shows the response surface of the molar ratio against temperature (t = 0.5 h), observing that opera- ting with higher molar ratios can slightly increase the MAG %, this result is similar to that obtained by Noriega et al. (2013), in which the enzymatic esterification reaction of n-3 PUFA and G was studied in an intermittent reactor.

Figure 2. Response surface for global esterification (GE) in the PBR with recirculation. a) versus molar ratio and time at T=70°C, b) versus temperature and molar ratio at t = 1.5 h.

Figura 2. Superficie de respuesta para la esterificación global (EG) en el RLE con recirculación. a) contra la relación molar y el tiempo a T=70°C, b) contra la temperatura y la relación molar a t = 1.5 h.

Figure 3. Response surface for monoacylglyerols (MAG) in the PBR with recirculation. a) versus molar ratio and time at T = 28°C, b) versus temperature and molar ratio at t = 0.5 h.

Figura 3. Superficie de respuesta para los monoacilgliceroles (MAG) en el RLE con recirculación. a) contra la relación molar y el tiempo a T = 28°C, b) contra la temperatura y la relación molar a t = 0.5 h.

Synthesis of diacylglycerols

DAGs are an intermediate product in the enzymatic esteri- fication reaction, however, these are lipids with high added

DAG % = -1.265 M2 - 0.016 T2 - 5.297 t2 + 0.799 M t + 0.126 T t

+ 7.045 M + 1.497 T + 7.166 t - 18.672

(5)

value due to their multiple benefits. The parameters with the greatest effect in the model are the quadratic effects of time, molar ratio and temperature, presenting values p < 0.001. The model that determines the percentage mole fraction of DAG (% DAG) is the following:

The RSM analysis performed on temperature versus time (Fig. 4a) shows that at a temperature of 50.33°C, molar ratio of 3.26 mol PUFA/mol G and reaction time of 1.52 h, the production of DAG can be maximized in the enzymatic es- terification reaction (DAG %=36). This behavior is due to the fact that, in the reaction mechanism, at one of the first stages

of esterification, the MAG react with PUFA to form DAG, but in the next stage these are consumed for the formation of TAG. Different investigations have published kinetic profiles of DAG formation similar to those obtained in this investigation (Correa et al., 2017; Bornadel et al., 2013). Duan et al. (2010) sought to optimize the formation of 1,3-DAG from an enzy- matic esterification reaction in a batch reactor. They obtained optimal formation conditions at 60°C and a molar ratio of 2.5 mol FFA/ mol G, reaching a concentration of 40 %. Lo et al. (2007) studied the enzymatic esterification of palmitic and oleic acid catalyzed by the lipozyme lipase RM IM, obtaining optimal conditions for the formation of DAG at 66.29°C, ope- rating with a molar ratio of 2.14 mol FFA/mol G. Although the results obtained differ slightly from those reported by other researchers, the enzymatic esterification reaction presents many factors that can influence the performance and beha- vior of the process, such as the composition of the substrate (Rivero et al., 2020; Kahveci and Xu, 2011), the lipase used and the water present in the reaction medium.

Synthesis of triacylglycerols

Due to the low goodness of fit of a second-order polynomial model (R2=0.74) to describe the percentage mole fraction of TAG, it was decided to consider a third-order model that in- cluded the combined effect of all factors and their respective cubic terms. The terms with the greatest effect in the model are the quadratic and cubic terms of the molar ratio (p < 0.01). The equation that describes the % TAG is:

% TAG = -0.844 M3 - 3.695 t3 + 0.059 M T t + 7.550 M2 -0.009 T2

‌TAGs are the products generated at the last stage of the enzymatic esterification reaction; this is reflected in their low concentration obtained in the experimental results due to the reversibilities present during the reaction. The model ob- tained through RSM indicates that higher molar fractions of TAG are obtained by increasing the reaction time and with n molar ratio of approximately 4.5 mol PUFA/mol G, reaching a molar fraction close to 26 % (Fig. 5a). The low TAG conversions achieved in the experimental runs may also be due to the reversible effects caused by the water produced during the esterification reaction, for this reason, many investigations seek to eliminate or isolate the water present in the medium to optimize the process and obtain better performances in the formation of TAG (Rodrigues and Fernandez, 2010). No- riega et al. (2013)studied the yield of TAG formation in the enzymatic esterification of PUFA catalyzed by the NV-435 lipase. After a 24 h reaction they only reached yields close to 12 %, so they added molecular meshes in order to isolate the water in the reaction medium and after another 24 h they achieved yields close to 80 %; this implies that the water caused strong reversible effects in the formation of TAG. The MSR graph of temperature versus molar ratio (Fig. 5b), at t=1.5 h, reflects that higher TAG conversions can be obtained at temperatures close to 50°C.

It is recommended that the model obtained be con- sidered only within the operating range established in the experimental design because, although a good fit (R2=0.92) is observed, certain inconsistencies occur when extrapolating

+ 17.97 t2 - 0.153 M T - 4.139 M t - 8.384 M + 1.042 T

-22.073 t + 5.029

(6)

the study range. Studying a process through an experimental design is very effective to contemplate the effect of factors

Figure 4. Response surface for diacylglycerols (DAG) in the PBR with recirculation. a) versus temperature and time at molar ratio =3.26, b) versus temperature and molar ratio at t=1.5 h.

Figura 4. Superficie de respuesta para diacigliceroles (DAG) en el RLE con recirculación. a) contra la temperatura y el tiempo a una relación molar =3.26, b) contra la temperatura y la relación molar a un t=1.5 h.

Figure 5. Response surface for triacylglycerols (TAG) in the PBR with recirculation. a) versus molar ratio and time at T = 50°C, b) versus temperature and molar ratio at t = 1.5 h.

Figura 5. Superficie de respuesta de triacilgliceroles (TAG) en el RLE con recirculación. a) contra la relación molar y el tiempo a T = 50°C, b) contra la tempera- tura y la relación molar a t = 1.5 h.

on the responses analyzed, but the models are usually limited to values established at the treatment levels (Montgomery, 2019).

CONCLUSIONS

‌The operation of the PBR with recirculation presented mul- tiple benefits compared to an intermittent reactor, such as high degrees of esterification at relatively short operating ti- mes, the possibility of reusing the immobilized lipase which, throughout the experiments, did not present apparent acti- vity loss and in addition, resulted in an easy products sepa- ration. The experimental design considered to evaluate the effect of the ratio of substrates, temperature and operation time on the synthesis of structured acyglycerols allows the optimization of reaction conditions depending on the desi- red product. The analysis of variance showed that the main factors of the models presented a significant effect on the different responses (p < 0.05), and the established models had a good fit to the experimental data (R2 > 0.92). With RSM, it was determined that a GE close to 80 % can be achieved when operating with a molar ratio of 0.5 mol PUFA/mol G, 70°C and 2.75 h of reaction. To maximize the production of MAG, it is recommended to operate with a molar ratio of 0.5 mol PUFA/mol G and 30°C, for 2.5 h of reaction, obtaining 80 % mol MAG. The maximum formation of DAG, close to 36 mol %, is reached with a ratio of 3.26 mol PUFA/mol G, 50.3

°C and 1.5 h of reaction. As for TAG, a molar fraction of 26 % was achieved when operating with a ratio of 4.5 mol PUFA/ mol G, 50 °C and 2 h of reaction. The water produced during the reaction caused considerable reversible effects during the reaction, which may favor the production of MAG, but

impair the formation of TAG. This study can be considered as a reference for the design and scaling of similar bioprocesses.

AG Acylglycerol ANOVA Analysis of variance

DAG % Percent molar fraction of diacylglycerols DAG Diacylglycerol

FFA Free fatty acids

GE Global esterification

G Glycerol

MAG % Percent molar fraction of monocylglycerols MAG Monoacylglycerol

PBR Packed bed reactor

PUFA Polyunsaturated fatty acid

RCCD Rotatable central composite design RSM Response surface methodology

TAG % Percent molar fraction of triacylglycerols

   TAG                                    Triacylglycerol                      

Abbreviations

Acknowledgments

Author Correa-Leyva gratefully acknowledges financial sup- port for his PhD studies (scholarship 590674) from CONAH- CYT (Mexico).

CONFLICTS OF INTEREST

The authors declare no conflicts of interest.

REFERENCES

Baeza, R., López, L.X. and García, H.S. 2014. Biocatalytic modification of food lipids: reactions and applications. Revista Mexicana de Ingeniería Química 13 (1), 29-47. https://www.redalyc.org/articulo.oa?id=62031166003

Bornadel, A., Orellana, A.C., Adlercreutz, P. and Hatti K.R. 2013. Kinetic modeling of lipase-catalyzed esterification reaction between aleic acid and trimethylolpropane: a simplified model for multi-substrate multi-product ping–pong mechanisms. Biotechnology Progress 29, 1422-1429. https://doi.org/10.1002/btpr.1806

Castejón Natalia & Francisco Javier Señoráns. (2021). Integrated green and enzymatic process to produce omega-3 acylglycerols from Echium plantagineum using immobilized lipases. J Am Oil Chem Soc 98: 341–352. DOI 10.1002/ aocs.12464.

Castellanos, T.L. and Rodríguez, D.M. 2015. El efecto de omega 3 en la salud humana y consideraciones en la ingesta. Revista Chilena de Nutrición 42 (1), 90-95. http://dx.doi.org/10.4067/ S0717-75182015000100012

Correa, C., Tejeda, A., Martin, A.R., García, H.S. and Rodríguez, J.A. 2017. Cinética de esterificación enzimática de ácidos grasos poliinsaturados n-3 a glicerol: Mecanismo ping-pong multi- sustrato multi-producto. Revista Mexicana de Ingeniería Química 16 (3), 805-812. https://www.redalyc.org/articulo. oa?id=62053304010

Damasco, E., Ventura, J.L., Guevara, M. and Zentella, A. 2018. Omega-3 polyunsaturated fatty acids supplemented diet and its preventive effect on tumor growth in nude mice. Nutrición Clínica y Dietética Hospitalaria 38 (2), 16-21. DOI: 10.12873/382damasco

Duan, Z.Q., Du, W. and Liu, D.H. 2010. Novozym 435-catalyzed 1, 3-diacylglycerol preparation via esterification in t-butanol system. Process Biochemistry 45 (12), 1923-1927. https:// doi.org/10.1016/j.procbio.2010.03.007

Filli, K.B., Nkama, I., Abubakar, U.M. and Jideani, V.A. 2010. Influence of extrusion variables on some functional properties of extruded millet-soybean for the manufacture of ‘fura’: A Nigerian traditional food. African Journal of Food Science 4 (6), 342- 352. Recuperado de: bit.ly/3dgjrV6

Gutierrez Pulido, H., and De la Vara, S.R. 2012. Análisis y Diseño de Experimentos. 3er Edición. Ed. Mc Graw Hill. México.

Hallahan, B., Ryan, T., Hibbeln, J.R., Murray, I.T., Glynn, S., Ramsden, C.E., SanGiovanni, J.P. and Davis, J.M. 2016. Efficacy of omega-3 highly unsaturated fatty acids in the treatment of depression. The British Journal of Psychiatry 209 (3), 192-

201. https://doi.org/10.1192/bjp.bp.114.160242[ Hernández Rodas, M.C., Morales, P.J., Valenzuela, B.R., Morales,

I.G. and Valenzuela, B.A. 2016. Beneficios de los ácidos grasos poliinsaturados de cadena larga n-3 en la enfermedad por hígado graso no alcohólico. Revista Chilena de Nutrición 43 (2), 196-205. http://dx.doi.org/10.4067/S0717-

75182016000200013

Hong, S.I., Ma, N., No, D.S., Choi, N., Baik, J.Y., Kim, C.T., Kim, Y.,

Chang, E. and Kim, I.H. 2014. Enrichment of DHA from tuna oil in a packed bed reactor via lipase-catalyzed esterification. Journal of the American Oil Chemists’ Society 91 (11), 1877- 1884. DOI 10.1007/s11746-014-2536-9

Kahveci, D. and Xu, X. 2011. Repeated hydrolysis process is effective for enrichment of omega 3 polyunsaturated fatty acids in salmon oil by Candida rugosa lipase. Food Chemistry 129 (4), 1552-1558. https://doi.org/10.1016/j. foodchem.2011.05.142

Khan, I., Hussain, M., Jiang, B., Zheng, L., Pan, Y., Hu, J., Khan, A., Ashraf, A. and Zou, X. 2023. Omega-3 long-chain polyunsaturated fatty acids: Metabolism and health

implications. Progress in Lipid Research, 101255. https://doi. org/10.1016/j.plipres.2023.101255

Kołodziej, Ł., Czarny, P.L., Ziółkowska, S., Białek, K., Szemraj, J., Gałecki, P., Su, K.P. and Śliwiński, T. 2023. How fish consumption prevents the development of major depressive disorder? A comprehensive review of the interplay between n-3 PUFAs, LTP and BDNF. Progress in Lipid Research, 101254. https://doi.org/10.1016/j.plipres.2023.101254

Lee, E.J., Lee, M.W., Kim, H.J., Oh, S.W., Kim, Y. and Kim, I.H. 2016.

Preparation of high purity docosahexaenoic acid from microalgae oil in a packed bed reactor via two-step lipase- catalysed esterification. Journal of Functional Foods 21, 330-337. https://doi.org/10.1016/j.jff.2015.12.020

Montgomery, D. C. 2019. Response surface methods and designs. In Design and Analysis of Experiments, 10th ed. John Wiley & Sons.

Nahemiah, D., Nkama, I. and Badau, M.H. 2015. Application of response surface methodology (RSM) and central composite design (CCD) to optimize minerals composition of rice- cowpea composite blends during extrusion cooking. Int J Food Sci Nutr Eng 5, 40-52. DOI: 10.5923/j.food.20150501.06 Noriega Rodríguez, J.A., Carrillo Pérez, E., Gámez Meza, N., Medina Juárez, L.A., Baeza Jiménez, R. y García, H. S.

2013. Optimization of the lipase catalyzed production of structured acylglycerols with polyunsaturated fatty acids isolated from sardine oil. Journal of Food Research 2 (6), 97-

105. http://dx.doi.org/10.5539/jfr.v2n6p97

Nurul Nadiah Abd Razak, Yolande Pérès, Lai Ti Gew, Patrick Cognet, and Mohamed Kheireddine Aroua. 2020. Effect of reaction medium mixture on the lipase catalyzed synthesis of diacylglycerol. Industrial & Engineering Chemistry Research 59 (21), 9869-9881. DOI: 10.1021/acs.iecr.0c00298

Pereira, G.N., Holz, J.P., Giovannini, P.P., Oliveira, J.V., de Oliveira,

D. and Lerin, L.A. 2018. Enzymatic esterification for the synthesis of butyl stearate and ethyl stearate. Biocatalysis and Agricultural Biotechnology 16, 373-377. https://doi. org/10.1016/j.bcab.2018.09.008

Piñeiro, G., Lago, N. y Culebras, J.M. 2013. Papel de los ácidos grasos omega-3 en la prevención de enfermedades cardiovasculares. Revisión. Nutrición Hospitalaria 28, 1-5. http://dx.doi.org/10.3305/nh.2013.28.1.6312

Rasband, W.S. 1997-2018. ImageJ, U.S. National Institutes of Health, Bethesda, Maryland, USA, imagej.nih.gov/ij/.

Rimm, E.B., Appel, L.J., Chiuve, S.E., Djoussé, L., Engler, M.B., Kris- Etherton, P.M., Mozaffarian, D., Siscovick, D. and Lichtenstein,

A.H. 2018. Seafood long-chain n-3 polyunsaturated fatty acids and cardiovascular disease: a science advisory from the American Heart Association. Circulation, 138 (1), e35-e47. https://doi.org/10.1161/CIR.0000000000000574

Rivero, F., Padial, M., Guadix, E.M. and Morales, R. 2020. Novozyme 435 and lipozyme RM IM preferably esterify polyunsaturated fatty acids at the sn‐2 position. European Journal of Lipid Science and Technology 122 (10), 2000115. https://doi.org/10.1002/ejlt.202000115

Monteiro, R.R.C., Virgen-Ortiz, J.J., Berenguer-Murcia, A., da Rocha, T.N., dos Santos, J.C.S., Alcántara, A.R. and Fernandez- Lafuente, R. 2021. Biotechnological relevance of the lipase A from Candida antarctica. Catalysis Today, 362,141-154.

Rodrigues, R.C. and Fernandez, R. 2010. Lipase from Rhizomucor miehei as a biocatalyst in fats and oils modification. Journal of Molecular Catalysis B: Enzymatic 66 (1-2), 15-32. https:// doi.org/10.1016/j.molcatb.2010.03.008

Shahidi F. and Ambigaipalan P. 2018. Omega-3 Polyunsaturated fatty acids and their health benefits. Annual Review of Food Science and Technology 9, 345-81. https://doi.org/10.1146/ annurev-food-111317-095850

Torres, P., Pedrali, A., Bavaro, T., Ambrosini, S., Marrubini, G., Pappalardo, V.M., Massolini, G., Terreni, M. and Ubiali, D. 2014. Preparation of PUFA concentrates as acylglycerols via enzymatic hydrolysis of hempseed oil (Cannabis sativa L.) in a homogeneous low‐water medium. European Journal of Lipid Science and Technology 116 (11), 1496-1504. https:// doi.org/10.1002/ejlt.201400201

Wang, L., Wang, Y., Hu, C., Cao, Q., Yang, X. and Zhao, M. 2011. Preparation of diacylglycerol-enriched oil from free fatty acids using lecitase ultra-catalyzed esterification. Journal of the American Oil Chemists’ Society 88 (10), 1557-1565. DOI 10.1007/s11746-011-1821-0

Wang, Y. and Huang, F. 2015. N-3 polyunsaturated fatty acids and inflammation in obesity: local effect and systemic benefit. BioMed research international, 1-16. https://doi. org/10.1155/2015/581469

Xie, D., Chen, Y., Yu, J., Yang, Z., Wang, X. & Wang, X. (2023). Progress in enrichment of n-3 polyunsaturated fatty acid: A review. Critical Reviews in Food Science and Nutrition, 63 (32), 11310-11326. https://doi.org/10.1080/10408398.2022.

2086852.

Xu, X. 2000. Production of specific‐structured triacylglycerols by lipase‐catalyzed reactions: a review. European journal of lipid science and technology 102 (4), 287-303. https://

doi.org/10.1002/(SICI)1438-9312(200004)102:4<287::AID- EJLT287 > 3.0.CO;2-Q

Yuan, Q., Ramprasath, V., Harding, S., Rideout, T. and Chan, Y. 2010. Diacylglycerol oil reduces body fat but does not alter energy or lipid metabolism in overweight, hypertriglyceridemic women. The Journal of Nutrition 140, 1122-1126. https:// doi.org/10.3945/jn.110.121665

Zhao, T., Kim, B.H., Hong, S.I., Yoon, S.W., Kim, C.T., Kim, Y. and Kim, I. H. 2012. Lipase‐catalyzed production of pinolenic acid concentrate from pine nut oil using a recirculating packed bed reactor. Journal of Food Science 77 (2), 267-271. https://doi.org/10.1111/j.1750-3841.2011.02562.x

‌Zhen Zhang, Fang Liu, Xiang Ma, Huihua Huang, and Yong Wang. 2018. Two-stage enzymatic preparation of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) enriched fish oil triacylglycerols. Journal of Agricultural and Food Chemistry 66 (1), 218-227. DOI: 10.1021/acs.jafc.7b04101

Zhang, Z., Lee, W.J., Sun, X., and Wang Y. 2022. Enzymatic interesterification of palm olein in a continuous packed bed reactor: Effect of process parameters on the properties of fats and immobilized Thermomyces lanuginosus lipase. LWT- Food Science and Technology 162 :113459. https://doi. org/10.1016/j.lwt.2022.113459

Zhu, Y., Feng, Y., Wang, J., Yuan, Z., Miao, Y., Miao, T., Gao, B. and Zhang, L. 2024. Selective esterification design of lipases for TAG synthesis based on the unique structure of curved DHA. ACS Food Science & Technology, 4 (7), 1722-1730. https:// doi.org/10.1021/acsfoodscitech.4c00174