The Effect of Terpenoid Compounds on the Formation of Advanced Glycation Endproducts (AGEs) in Model Systems

: Background: Terpenoid compounds, despite their established antioxidant ability, are neglected as potential glycation regulators. Methods: In-vitro model systems of lysine (0.1 M) with glucose (0.1 M and 1 M) were incubated at 80 ◦ C and 100 ◦ C for 3 h in the presence of aniseed oil, thymol and linalool (2–8 µ M). Color development, absorbance at UV-Vis (280 nm, 360 nm, 420 nm), ﬂuorescence intensity ( λ exc = 370 nm, λ emm = 440 nm) and lysine depletion (HPLC-FL) were measured to monitor the progress of the Maillard reaction. Response Surface Methodology was used to analyze the impact of the ﬁve experimental conditions on the glycation indices. Results: All terpenoid compounds promoted color development and did not affect lysine depletion. The choice of terpenoid compound impacted glycation at 280 nm, 360 nm and 420 nm ( p < 0.02). The effect was stronger at lower temperatures ( p < 0.002) and 0.1 M glucose concentrations ( p < 0.001). Terpenoid concentration was important only at 360 nm and 420 nm ( p < 0.01). No impact was seen for ﬂuorescence intensity from the choice of terpenoid compounds and their dose ( p = 0.08 and p = 0.44 respectively). Conclusion: Terpenoid compounds show both anti- and proglycative activity based on the incubation conditions. Thymol showed the largest antiglycative capacity, followed by aniseed oil and linalool. Maximal antiglycative capacity was seen at 0.1 M glucose, 2 µ M terpenoid concentration, 80 ◦ C and 1 h incubation. a lysine


Introduction
Essential oils (EOs) constitute a large class of secondary plant metabolites with important defensive and attractive role between plants and their environment [1]. Liquid at ambient temperature, their volatile fraction is composed principally by mono-, sesqui-and di-terpenes alongside some volatile phenolic substances like phenylpropanoids [2]. Their ability to act as hydrogen donators or radical scavengers has been linked with inhibition of lipid peroxidation [3]. This ability to terminate the lipid peroxidation chain reaction is the main application of EOs as natural food antioxidants and it is linked to their terpenoid content [4].
Among terpenoids, those with phenolic ring like thymol and carvacrol possess the highest antioxidant activity followed by the monoterpene hydrocarbons with strongly activated methylene group like terpinolene, α-terpinene and γ-terpinene. Linalool, a tertiary allylic alcohol, shows a pro-oxidant effect while benzene derivatives (other than thymol and carvacrol) like estragol, a phenylopropanoid, are more efficient in preventing the formation of primary lipid oxidation products [5].
On the other hand, little is known about the impact of terpenoid compounds on other oxidative reaction in foods. Glycation is an important chemical reaction occurring mainly during food processing when a free amino group of protein interacts with the carbonyl group of glucose or other reducing sugars [6][7][8]. Among the complex mixture of compounds produced, several carbonyl derivatives such as Strecker aldehydes and dicarbonyls (glyoxal, methylglyoxal, 3-deoxyglucosone) are formed [9,10]. Those highly reactive carbonyl derivatives can react with the ε-amino group of lysine or the guanidino group of arginine to generate products known as Advanced Glycation End-products (AGEs). AGEs are a very heterogeneous group of compounds and can be classified into fluorescent (such as pentosidine) and nonfluorescent AGEs (such as Nε-(carboxymethyl)lysine CML) as well as cross-linking and non-cross-linking AGEs. Some AGEs exhibit both fluorescent and cross-linking properties, while others show none of the mentioned characteristics [7,11].
Mainly produced through the thermal processing of foods, mitigation of AGEs formation in foods is often dealt through modifications in heating time and temperature [12]. Alternatively, molecules that exhibit enhanced antioxidant activity, sugar autoxidation inhibition, reactive dicarbonyls trapping, and amino group binding, could be used as AGEs inhibitors [13,14]. Many synthetic inhibitors, have been used to control AGEs production in food systems and in humans (aminoguanidine, metformin, irbesartan, pyridoxamine, etc.) but their severe side effects highlight the need for effective and safe alternatives like natural plant extract [8]. Phenolic compounds with potent antioxidant and free radical scavenging activity and the ability to trap reactive carbonyl species are commonly studied [11,13]. Phenolic acids such as ferulic and isoferulic acid [15,16], cinnamic [17], chlorogenic [14], vanillic, syringic, sinapic, p-coumaric and caffeic acids [18,19] have showed satisfactory AGEs inhibition in vitro. In food systems, flavonoids such as rutin, quercetin, genistein and narigenin, stilbenes (resveratrol) [13,20] and tannins (punicalagin) [21] have been identified as potent AGEs inhibitors. In addition, several herbal infusions of medicinal plants have been successfully screened as inhibitors [22]. Hydroxytyrosol, the main phenolic compound in olive oil, was also recently identified as a potential AGEs inhibitor through the binding of methylglyoxal [23].
Despite their antioxidant and radical scavenging properties, EOs have not been fully studied as AGEs inhibitors. In a PatchDock study, limonene was shown to inhibit albumin glycation through hydrophobic binding of its monoterpene hydrocarbon [24]. In another study, the essential oil of caraway which consists mainly of S-carvone (a monoterpene aldehyde), prevented the formation of early Maillard intermediates. Interestingly, the antiglycative capacity of caraway oil was not linked to the carbonyl group in the S-carvone structure. The same study, also showed that anethole and thymol the major constituents in sweet cumin, were more likely to act as glycation promoters than inhibitors [25].
In the present study, linalool and thymol, two monoterpene compounds as well as aniseed oil, a commercial product purchased by Sigma Aldrich were selected to be studied, in terms of their role in lysine glycation. The main essential oil compound in aniseed oil is anethole (about 80-90%) together with small amounts of estragol, an isomer of anethole. The terpenoid compounds' choice was based on the fact that they are main essential oil compounds of herbs used in European cuisine. The objective of this study was to test the impact of these terpenoid compounds on lysine glycation models representative of thermal food processing. Their role as potential Maillard modulators was tested for different terpenoid concentrations, temperatures, heating times and glucose concentrations.

Reagents
D-Glucose, L-Lysine monohydrochloride and linalool were purchased from Acros Organics, thymol from Applichem and anise oil from Sigma Aldrich (CAS: 8007-70-3, nature identical, synonym: aniseed oil). All solutions (glucose, lysine and the test compounds) were prepared in phosphate buffer 0.1 M, pH 7.2 ± 0.1 (PBS, Bectron, Dickinson and Company, Bergen County, NJ, USA) and Tween 80 (Serva Electrophoresis GmbH, Heidelberg, Germany) which helps to incorporate aniseed oil and thymol in phosphate buffer. As linalool possesses a very good aqueous solubility at pH 7.0, the use of Tween 80 was not necessary [26].

Experimental Model of Protein Glycation
A glucose:lysine model with or without (control) the presence of terpenoid compounds were used as the protein glycation model. Two different glucose:lysine concentrations (10:1 and 1:1 molar ratios) were used. Each reaction mixture consisted of 3 mL of glucose either 1 M or 0.1 M, 3 mL of lysine 0.1 M and 6.25, 12.5, 18.75 and 25 µL of terpenoid compounds working solutions 2 mM, resulting to a final concentration of 2, 4, 6 and 8 µM. The models were added to glass vials with ribbed stopers and heated using a Thermo Shaker LLG, Unithermix 1 Pro at 80 • C and 100 • C at a fixed 600 rpm rotational speed for a total of 3 h. All working solutions were freshly prepared at 100 mL final volume and in the case of aniseed oil and thymol they were mixed under magnetic stirring for 1 min followed by homogenization using an UltraTurrax appliance for 1 min to achieve better dispersion of the lipophilic droplets into the aqueous phase. All experiments were performed in duplicate. Samples and controls were taken at, 1, 2 and 3 h of heating. Aliquots of 200 µL of the reaction mixtures and controls were transferred into wells of a 96-well plate and kept at −20 • C until analysis.

Color Development Assay
The color of the studied model systems was determined using a Lovibond Spectrocolorimeter LC100. CIE L*a*b* values were measured using daylight (D 65/10 • ) with brightness. Color development was quantified using controls at t = 1, 2 and 3 h as the reference value and via the measurement of lightness, redness and yellowness values following the equation: L = lightness, a = redness, b = yellowness and L • , α • , and b • are the color of the control at t = 1, 2 and 3 h of heating respectively [27].

Assessment of the Maillard Reaction Evolution by Ultraviolet-Visible (UV-Vis) Spectroscopy
UV-Vis spectroscopy (Epoch 2, BioTek Instruments, Inc., Winooski, VT, USA) in three distinct wavelengths was used to measure the concentration of Maillard products at each experimental point expressed in arbitrary units. Absorbance at 280 nm was used for the quantification of early-stage products such as Schiff bases and Amadori rearrangement products. Intermediate products such as unsaturated carbonyls, dicarbonyls, ROS and furans that come from the dehydration, oxidation and fragmentation of initial products were measured at 360 nm, and the final products of non-enzymatic browning such as cross-linked products, aromatic compounds and colored polymers were determined at 420 nm [28].
The ability of each terpenoid compound to impact the formation of early, intermediate and advanced glycation products was calculated by the following equation: where Abs 280_360_420_sample is the absorbance of the reaction mixtures with the presence of the test molecules and Abs control is the absorbance of the reaction mixtures without the presence of the test molecules.

Measurement of Fluorescent AGEs Formation
Following incubation at 80 • C and 100 • C for 1, 2 and 3 h, fluorescence intensity was measured using a fluorescent spectrometer (Victor X2 Fluorimeter 2030, Perkin-Elmer-Massachusetts, USA) at an excitation wavelength of λ exc = 370 nm and an emission wavelength of λ emm = 440 nm. The capacity (%) of each terpenoid compound to impact the formation of fluorescent AGEs was based on [22] and calculated by the following equation: where Fluo sample is the fluorescence intensity of each tested molecule with glycation models; Fluo control is the fluorescence intensity of the glycation models.

Preparation of Solutions
Standard stock solutions of lysine (10 µmol/mL) were prepared in 0.1 M HCl. For derivatization, it is essential that the pH value of the standard is in the range of pH 7. Therefore, the standard was first diluted 1:1 with 0.1 M NaOH. The stock solution was then diluted with water to prepare the standard working solutions ranging from 0.1 to 0.5 µmol/mL for calibration. Solutions were stored at 4 • C.
The derivatization reagent was prepared by mixing 1 mL of OPA solution (10 mg/mL, dissolved in methanol) in 10 mL of borate buffer (consisting of 0.8 M borate buffer in KCl and 0.8 M NaOH with a final pH of 9.9) and 16 µL MCE. The derivatization reagent is stable within 9 days [29].
For the preparation of the samples, 50 µL of the heated model system were transferred to a 10 mL volumetric flask and then filled with analytical grade water.

Derivatization Procedure
The derivatization took place manually by mixing 250 µL of the diluted sample or the standard solution with 150 µL of derivatization reagent. After 2 min, 50 µL of the mixture were added to 200 µL of neutralization solution (as described in Section 2.4.3) [30] and 20 µL of the final solution were injected into the HPLC system.

Statistical Analysis and Response Surface Methodology
Descriptive statistics, specifically box and whiskers plots, were used to estimate the data distribution and to explore potential differences between the terpenoid compounds and their impact on amino acid glycation.
A D-Optimal design was combined with Response Surface Methodology (RSM) as a post hoc analysis [32]. A reduced quadratic randomized model was retrofitted to the experimental conditions. Time, temperature, glucose concentration and terpenoid compounds were used as numerical values and forced to predefined values as described above (fixed values), while the choice of terpenoid compounds was set as a categorical value. A total of 33 runs were determined by a selection criterion chosen during the experimental design ( Table 1). Color Development, Absorbance at 280 nm, 360 nm and 420 nm and Fluorescence Intensity expressed as (%) difference from the control were used for response prediction. For this purpose, a randomly selected sample was used for response prediction. The model's fitness was confirmed by analysis of variance (ANOVA) and the determination coefficient (R 2 ) using p-values. Dependent variables were analyzed either via a Linear or 2FI model. Model choice was based on the following criteria: test of normality for residuals, residuals versus predicted and residuals versus runs plots, Predicted R 2 and Lack of Fit test. Box-Cox analysis were used to identify the most suitable data transformations. Namely, Color development underwent an inverse square root transform, and Absorbance at 420 nm underwent a power transform with a λ = 2.56 and k = 66.583. Color Development was analyzed with a linear model and Absorbances at 280 nm, 360 nm and 420 nm and Fluorescence Intensity were analyzed with a 2FI model. Significance level was set at α = 0.05. All analyses were carried out using the IBM SPSS Statistics 25 Software and Design-Expert 11.0.5.0 (Stat-Ease, Inc., Minneapolis, MN, USA).

Color Development
The addition of terpenoid compounds during the heat treatment of the glucose/lysine solutions resulted to positive color development. All terpenoid compounds tested were associated with greater color development at all time-points. Not only color development was larger at 100 • C as indicated by the greater median values in all systems, but all systems also showed greater variability at 100 • C as indicated by the large interquartile ranges ( Figure 1). Large variabilities were considered a positive finding in this study, as it indicated potential differences in the Maillard reaction process among the different study conditions.

Color Development
The addition of terpenoid compounds during the heat treatment of the glucose/lysine solutions resulted to positive color development. All terpenoid compounds tested were associated with greater color development at all time-points. Not only color development was larger at 100 °C as indicated by the greater median values in all systems, but all systems also showed greater variability at 100 °C as indicated by the large interquartile ranges ( Figure 1). Large variabilities were considered a positive finding in this study, as it indicated potential differences in the Maillard reaction process among the different study conditions. Based on the RSM analysis, temperature was the only factor significantly and positively associated with color development (p < 0.001). No differences were seen among Based on the RSM analysis, temperature was the only factor significantly and positively associated with color development (p < 0.001). No differences were seen among the different terpenoid compounds (Figure 1) or their concentration (p = 0.10 and p = 0.38 respectively).

Assessment of the Maillard Reaction Evolution by Ultraviolet-Visible (UV-Vis) Spectroscopy
The measurement of the absorbance in specific wavelengths associated with different stages of the Maillard reactions, however, highlighted differences among the terpenoid compounds behavior. As shown in Figure 2, the terpenoid compounds tested have different effects on the absorbance measurements at 280 nm, both among them and in different incubation conditions. In 80 • C systems, aniseed oil (median-77.7%) and thymol (median-80.5%) were associated with 100% lower concentrations of early glycation products in 0.1 M glucose in the first hour of incubation (Supplementary Table S1) which were reduced over time. Linalool (median-39.6%) showed a constant association with reduced levels of early glycation products throughout the 3 h of incubation. In 1 M glucose at 80 • C, aniseed oil (median-9.35%) and thymol (median-5.35%) showed reduced inhibitory activity, if any, while linalool (median 10.25%) was associated with promotion of the early glycation product formation as early as the first hour of incubation (Supplementary Table S1). For incubations at 100 • C the effects of terpenoid compounds on early glycation products formation indicated almost no effect at 1 M glucose (for aniseed median-1.0%, for thymol median-0.15% and for linalool median-0.6%) whereas at 0.1 M glucose, thymol (median-8.7%) and aniseed oil (median-1.8%) showed small inhibitory activity. In addition, the inhibitory activity of aniseed oil and thymol was seen after 3 h of incubation at all glucose concentration (Supplementary Table S1). At 100 • C, all systems showed a much smaller variability in absorbance values at 280 nm compared to the incubations at 80 • C. At 80 • C and low glucose concentration the variability was larger as indicated by the interquartile range which was more pronounced for aniseed oil and thymol than linalool. Interestingly at 80 • C and 0.1 M glucose the net effect of all terpenoid compounds was towards reduced absorbance (inhibitory effect) ( Figure 2). different effects on the absorbance measurements at 280 nm, both among them and in different incubation conditions. In 80 °C systems, aniseed oil (median-77.7%) and thymol (median-80.5%) were associated with 100% lower concentrations of early glycation products in 0.1 M glucose in the first hour of incubation (Supplementary Table S1) which were reduced over time. Linalool (median-39.6%) showed a constant association with reduced levels of early glycation products throughout the 3 h of incubation. In 1 M glucose at 80 °C, aniseed oil (median-9.35%) and thymol (median-5.35%) showed reduced inhibitory activity, if any, while linalool (median 10.25%) was associated with promotion of the early glycation product formation as early as the first hour of incubation (Supplementary Table S1). For incubations at 100 °C the effects of terpenoid compounds on early glycation products formation indicated almost no effect at 1 M glucose (for aniseed median-1.0%, for thymol median-0.15% and for linalool median-0.6%) whereas at 0.1 M glucose, thymol (median-8.7%) and aniseed oil (median-1.8%) showed small inhibitory activity. In addition, the inhibitory activity of aniseed oil and thymol was seen after 3 h of incubation at all glucose concentration (Supplementary Table  S1). At 100 °C, all systems showed a much smaller variability in absorbance values at 280 nm compared to the incubations at 80 °C. At 80 °C and low glucose concentration the variability was larger as indicated by the interquartile range which was more pronounced for aniseed oil and thymol than linalool. Interestingly at 80 °C and 0.1 M glucose the net effect of all terpenoid compounds was towards reduced absorbance (inhibitory effect) ( Figure 2). The RSM model confirmed that, the choice of terpenoid compound (p = 0.001), glucose concentration (p < 0.001), temperature (p < 0.001) and time (p = 0.02) were associated with different effects on absorbance at 280 nm. Significant effects were also seen for the interactions between glucose concentration and temperature (p < 0.001), The RSM model confirmed that, the choice of terpenoid compound (p = 0.001), glucose concentration (p < 0.001), temperature (p < 0.001) and time (p = 0.02) were associated with different effects on absorbance at 280 nm. Significant effects were also seen for the interactions between glucose concentration and temperature (p < 0.001), temperature and terpenoid compound (p = 0.004), glucose concentration and terpenoid compound concentration (p = 0.004) and terpenoid compound and its concentration (p = 0.046). At both temperatures, in low glucose levels (0.1 M) higher terpenoid compound concentrations were associated with reduced inhibitory capacity overtime, while the opposite effect was seen in high glucose levels (1 M) (Figure 3). These were specifically true for linalool and aniseed oil.
On the contrary overtime higher thymol concentrations were associated with higher inhibitory capacity in all temperatures and glucose levels (data not shown).
In terms of their effect on intermediate glycation products measured at 360 nm unlike what was seen at 280 nm, all systems showed a smaller variability in absorbance measurements at 360 nm at 80 • C compared to 100 • C, with the greatest variability (interquartile ranges) seen at 100 • C and 1 M glucose (Figure 4).
Aniseed oil and thymol showed 8-14% inhibitory capacity on average at 80 • C (Supplementary Table S2), with greater inhibition in lower glucose concentration ( Figure 4A). At 100 • C the same terpenoid compounds showed inhibitory capacity at 2 h and 3 h of incubation (Supplementary Table S2), but only in low glucose concentrations (for aniseed median-13.35% and for thymol median-9.3%) ( Figure 4B). Aniseed oil (median 18.75%) and thymol (median 24.55%) promoted the production of intermediate glycation products at 100 • C and high glucose concentration ( Figure 4B). Linalool (median-13.85%) successfully inhibited the production of intermediate glycation products at higher glucose concentrations and that only for the 2 h and 3 h timepoints (Supplementary Table S2). temperature and terpenoid compound (p = 0.004), glucose concentration and terpenoid compound concentration (p = 0.004) and terpenoid compound and its concentration (p = 0.046). At both temperatures, in low glucose levels (0.1 M) higher terpenoid compound concentrations were associated with reduced inhibitory capacity overtime, while the opposite effect was seen in high glucose levels (1 M) (Figure 3). These were specifically true for linalool and aniseed oil. On the contrary overtime higher thymol concentrations were associated with higher inhibitory capacity in all temperatures and glucose levels (data not shown).   According to the RSM analysis, all factors were associated with differences in absorbance at 360 nm (p < 0.02). Significant effects were also seen for the interactions between time and temperature (p = 0.03), glucose concentration and terpenoid compound concentration (p = 0.03), terpenoid compound and its concentration (p = 0.02), glucose concentration and temperature (p = 0.007), glucose concentration and terpenoid compound (p = 0.02) as well as terpenoid compound and temperature (p = 0.02). In terms of their effect on intermediate glycation products measured at 360 nm unlike what was seen at 280 nm, all systems showed a smaller variability in absorbance measurements at 360 nm at 80 °C compared to 100 °C, with the greatest variability (interquartile ranges) seen at 100 °C and 1 M glucose (Figure 4). Aniseed oil and thymol showed 8-14% inhibitory capacity on average at 80 °C (Supplementary Table S2), with greater inhibition in lower glucose concentration ( Figure 4Α). At 100 °C the same terpenoid compounds showed inhibitory capacity at 2 h and 3 h of incubation (Supplementary Table S2), but only in low glucose concentrations (for aniseed median-13.35% and for thymol median-9.3%) ( Figure 4Β). Aniseed oil (median 18.75%) and thymol (median 24.55%) promoted the production of intermediate glycation products at 100 °C and high glucose concentration ( Figure 4B). Linalool (median-13.85%) successfully inhibited the production of intermediate glycation products at higher glucose concentrations and that only for the 2 h and 3 h timepoints (Supplementary Table S2).
According to the RSM analysis, all factors were associated with differences in absorbance at 360 nm (p < 0.02). Significant effects were also seen for the interactions between time and temperature (p = 0.03), glucose concentration and terpenoid compound concentration (p = 0.03), terpenoid compound and its concentration (p = 0.02), glucose concentration and temperature (p = 0.007), glucose concentration and terpenoid compound (p = 0.02) as well as terpenoid compound and temperature (p = 0.02).
At both temperatures, there was inhibition of the formation of intermediate glycation products at 0.1 M glucose at all time points, however there was stronger inhibitory effects with increasing time at higher glucose levels (1 M) ( Figure 5). At 1 M glucose, increasing terpenoid concentration was associated with lower inhibitory effect ( Figure 5). The later was seen in all terpenoid compounds tested. Interestingly, at 0.1 M glucose aniseed oil and linalool showed a positive association between terpenoid concentration and inhibitory effect (data not shown). The relationship between thymol concentration and its inhibitory effect remained negative at 0.1 M glucose as well.  (Figure 5). At 1 M glucose, increasing terpenoid concentration was associated with lower inhibitory effect ( Figure 5). The later was seen in all terpenoid compounds tested. Interestingly, at 0.1 M glucose aniseed oil and linalool showed a positive association between terpenoid concentration and inhibitory effect (data not shown). The relationship between thymol concentration and its inhibitory effect remained negative at 0.1 M glucose as well.
At 420 nm, at 80 • C thymol showed no inhibition (at 1 M glucose median-0.5%, at 0.1 M glucose median-0.95%) compared to the other two terpenoid compounds in both glucose systems (for aniseed oil at 0.1 M median-4.25%, at 1 M glucose median-11.55% and for linalool −3.65 and −14.25, respectively). At 100 • C, linalool showed higher inhibitory effect in both glucose systems (at 0.1 M glucose median-17.5% and at 1 M glucose median-21.2%) ( Figure 6B) compared to the other two terpenoids compounds tested. In the higher glucose concentration thymol promoted the production of advanced glycation products in high temperature (median 26.8%) throughout the heating time. Similar effects were seen for aniseed oil at higher glucose concentration (median 11.1%) (Supplementary Table S3). Once again, larger variability in absorbance measurements at 420 nm was seen at higher temperatures ( Figure 6).
According to the RSM model, absorbance at 420 nm was influenced by glucose concentration, incubation temperature, the choice of terpenoid compound and its concentration (p < 0.003 for all). Interactions between terpenoid compound choice and its concentration, terpenoid compound and glucose concentration and terpenoid compound and temperature also significantly impacted absorbance at 420 nm (p < 0.003).
More specifically, linalool showed the greatest inhibitory capacity followed by aniseed oil, while thymol did not show any inhibitory capacity (data not shown). As far as the interactions between terpenoid concentration, glucose concentration and temperature are concerned, for all compounds tested increasing concentrations were associated with lower inhibitory capacity at 1 M glucose concentration in both temperatures (Figure 7). At 0.1 M glucose, linalool and aniseed oil showed a positive relationship between terpenoid concentration and inhibitory capacity (data not shown). However, thymol showed decreasing inhibitory capacity with higher concentration in both glucose levels (data not shown). Appl. Sci. 2022, 12, x FOR PEER REVIEW 10 of 19 At 420 nm, at 80 °C thymol showed no inhibition (at 1 M glucose median-0.5%, at 0.1 M glucose median-0.95%) compared to the other two terpenoid compounds in both glucose systems (for aniseed oil at 0.1 M median-4.25%, at 1 M glucose median-11.55% and for linalool −3.65 and −14.25, respectively). At 100 °C, linalool showed higher inhibitory effect in both glucose systems (at 0.1 M glucose median-17.5% and at 1 M glucose median-21.2%) ( Figure 6B) compared to the other two terpenoids compounds tested. In the higher glucose concentration thymol promoted the production of advanced glycation products in high temperature (median 26.8%) throughout the heating time. Similar effects were seen for aniseed oil at higher glucose concentration (median 11.1%) (Supplementary Table S3). Once again, larger variability in absorbance measurements at 420 nm was seen at higher temperatures ( Figure 6).

Measurement of Fluorescent AGEs Formation
At 80 • C terpenoid concentration was negatively associated with the inhibitory effects specifically. aniseed oil had a similar inhibitory effect at all concentrations, while for thymol and linalool showed reduced inhibitory capacity with increasing concentrations (at 80 • C).
More specifically aniseed oil showed the greatest inhibitory capacity of all three terpenoid compounds tested at 80 • C both at 0.1 M (median-25.45%) and 1 M glucose (median-14.38%), followed by thymol (for 0.1 M median-18.01% and for 1 M median-10.67%) ( Figure 8A). Linalool only inhibited the production of fluorescent glycation products at 80 • C and low glucose concentrations (median-25.45%) and that after 2 h of heating (Supplementary Table S4). As far as variability is concerned, in the case of fluorescence intensity there is no clear pattern between temperatures and glucose concentrations. The highest interquartile ranges were observed at 100 • C incubation but for linalool that effect was seen at 1 M glucose, while for aniseed oil and thymol at 0.1 M glucose ( Figure 8B). According to the RSM model, absorbance at 420 nm was influenced by glucose concentration, incubation temperature, the choice of terpenoid compound and its concentration (p < 0.003 for all). Interactions between terpenoid compound choice and its concentration, terpenoid compound and glucose concentration and terpenoid compound and temperature also significantly impacted absorbance at 420 nm (p < 0.003).
More specifically, linalool showed the greatest inhibitory capacity followed by aniseed oil, while thymol did not show any inhibitory capacity (data not shown). As far as the interactions between terpenoid concentration, glucose concentration and temperature are concerned, for all compounds tested increasing concentrations were associated with lower inhibitory capacity at 1 M glucose concentration in both temperatures (Figure 7). At 0.1 M glucose, linalool and aniseed oil showed a positive relationship between terpenoid concentration and inhibitory capacity (data not shown). However, thymol showed decreasing inhibitory capacity with higher concentration in both glucose levels (data not shown). highest interquartile ranges were observed at 100 °C incubation but for linalool that effect was seen at 1 M glucose, while for aniseed oil and thymol at 0.1 M glucose ( Figure 8B). The RSM showed that time, glucose concentration and temperature were significantly associated with differences in fluorescence intensity (p < 0.006 for all). The interactions between terpenoid compound and its concentration and terpenoid compound and temperature were also significant factors of fluorescence intensity (p = 0.004 and p = 0.02 respectively). A positive relationship between terpenoid concentration and inhibition of fluorescent product formation at both glucose concentration is seen at 100 °C where at 80 °C seems to be independent as increasing the heating time. Overall, inhibition of fluorescent product formation was only seen at 80 °C and not at 100 °C (Figure 9).  The RSM showed that time, glucose concentration and temperature were significantly associated with differences in fluorescence intensity (p < 0.006 for all). The interactions between terpenoid compound and its concentration and terpenoid compound and temperature were also significant factors of fluorescence intensity (p = 0.004 and p = 0.02 respectively). A positive relationship between terpenoid concentration and inhibition of fluorescent product formation at both glucose concentration is seen at 100 • C where at 80 • C seems to be independent as increasing the heating time. Overall, inhibition of fluorescent product formation was only seen at 80 • C and not at 100 • C (Figure 9).

Specific Markers of Maillard Reaction-HPLC Determination of Lysine
As terpenoid compound concentration had a significant impact on the Maillard reaction products only in specific wavelengths and for the extreme terpenoid concentrations, lysine depletion was studied only in the case when 8 µM of each was added into the model systems.
In terms of the overall intensity of the Maillard reactions in the system tested, the addition of aniseed oil, linalool and thymol showed an increase in the amount of lysine lost during the three hours of heating at both 80 • C and 100 • C. As shown in Figure 10, this effect was the same among the compounds tested and the controls (p > 0.05) where the depletion of lysine reaches a plateau after the first hour of heating. In the high glucose concentration system at 100 • C lysine depletion was the highest among the conditions tested reached more than 70% after the first hour of heating.

Specific Markers of Maillard Reaction-HPLC Determination of Lysine
As terpenoid compound concentration had a significant impact on the Maillard reaction products only in specific wavelengths and for the extreme terpenoid concentrations, lysine depletion was studied only in the case when 8 μM of each was added into the model systems.
In terms of the overall intensity of the Maillard reactions in the system tested, the addition of aniseed oil, linalool and thymol showed an increase in the amount of lysine lost during the three hours of heating at both 80 °C and 100 °C. As shown in Figure 10, this effect was the same among the compounds tested and the controls (p > 0.05) where the depletion of lysine reaches a plateau after the first hour of heating. In the high glucose concentration system at 100 °C lysine depletion was the highest among the conditions tested reached more than 70% after the first hour of heating.

Discussion
This is one of the few studies that investigate the role of terpenoid compounds on the evolution of Maillard reaction. The study highlighted a general promotion of glycation by

Discussion
This is one of the few studies that investigate the role of terpenoid compounds on the evolution of Maillard reaction. The study highlighted a general promotion of glycation by terpenoid compounds as measured by color development and lysine depletion. However, measurement of absorbance and fluorescence intensity showed both antiglycative and proglycative capacity, depending on terpenoid compound concentration but most importantly the temperature and glucose concentration of the system. In detail, all terpenoid showed some antiglycative capacity at low glucose levels (0.1 M) and low temperatures (80 • C). In terms of effect sizes, thymol and aniseed oil showed stronger inhibitory capacity reducing the concentration of early glycation products by approximately 100% compared to the control, at least for the first hour of heating, while linalool showed much smaller inhibitory capacity if any. The RSM model allowed to test for the most desirable experimental condition in order to minimize the extent of the Maillard reactions and the production of AGEs in all five parameters tested. According to this model, for aniseed oil a system with 0.1 M glucose, 4 µM aniseed oil heated at 80 • C for 1 h would be the optimal solution (Desirability = 0.845). For thymol and linalool, a system with 0.1 M glucose, 2 µM aniseed oil heated at 80 • C for 1 h would be the optimal solution (Desirability = 0.949 and 0.713, respectively). According to the RSM analysis, the overall antiglycative capacity of the terpenoid compounds tested could be ranked Thymol > Aniseed Oil > Linalool.
Such a behavior could be attributed to the high antioxidant activity of thymol and aniseed oil. Research has highlighted that not only phenols (like thymol and carvacrol) but also compounds with methylene groups show high antioxidant activity [5]. Anethole and estragole, the main compounds of the aniseed oil, possess methylene groups on their structure which may be responsible for the observed inhibition. A pro-oxidant behavior for linalool was also observed in the same research article. On the other hand, the presence of anethole and thymol in sweet cumin were positively linked with the formation of early radical intermediates on experiment where equimolar amounts of glucose and glycine with varying amounts of pure terpenoid solutions were refluxed for 1 h [25]. The dual role of terpenoid compounds to act as antioxidant and prooxidant depending on their concentration was underlined by many researchers [4].
At the same time, a dose response relationship was rarely observed and when observed the direction of the relationship was unusual. In the cases of a dose response, it was the extreme terpenoid concentrations, of 2 µM and 8 µM, which showed the antiglycative activity and the middle terpenoid concentrations that showed a proglycative activity. Although the increasingly stronger antiglycative capacity with increasing terpenoid concentrations could be physiologically relevant, the fact that the smallest concentration tested in this experiment (2 µM) showed the same antiglycative capacity as 8 µM terpenoid concentrations was an unexpected finding. Previous research on the influence of terpenoid compounds on AGEs inhibition revealed that higher concentrations are more effective for such a behavior. Particularly cymene, an aromatic monoterpene, inhibited by 56.6% the total AGEs fluorescence when used at a concentration of 100 µM [24]. In the same study, the inhibitory capacity of cymene was attributed to the stabilization effect of the compound on the secondary structure of BSA. Based on our findings we could only stipulate those different mechanisms of action are responsible for the antiglycative activity at extreme terpenoid concentrations. One possible mechanism could be that the molecules antioxidant activity is linked with its antiglycative capacity at high terpenoid concentrations, whereas for the low terpenoid concentrations, mechanisms like hydrophobic binding on glycation sites should be further tested in the future [33].
Inhibition of the formation of fluorescent AGEs was only seen at 80 • C. This temperature is representative of a gentle cooking procedure as that obtained in traditional earthenware pots, consisting integral part of the Mediterranean culinary tradition [34], and it is concluded that it facilitates the action of terpenoid compounds towards inhibition of AGEs formation. On the contrary, cooking at 100 • C, representative of boiling procedure, promoted the formation of AGEs in both glucose levels.
Thymol and aniseed oil seems to have a uniform behavior toward the formation of the Amadori products during the early stage of the glycation in the equimolar system perhaps through their antioxidant and antiradical properties. Such a behavior was not observed in the higher glucose concentration suggesting that the compounds under investigation were not involved in glucose autoxidation. Previous research has showed that molecules with more than one phenolic rings such as hydroxytyrosol, and flavonoids possess the ability to act during the intermediate stages of glycation by trapping glyoxal and methylglyoxal and inhibiting AGEs formation [6,13,23]. Based on the later, thymol and aniseed oil molecules with one and none phenolic rings, respectively, are not expected to exhibit dicarbonyl trapping activities and hence would have limited impact on glucose autoxidation driven reactions.
As with many in vitro trials, this study possesses a series of questions for future research. This initial screening showed that the terpenoid compounds have both antiglycative and proglycative capacity depending on the glycation conditions and that these compounds may be modulating glycation in more than one mechanism. Although these findings are interesting, the study of mechanisms of action as well as the translation of the current findings to actual foods would require the experiments to be repeated in model systems using intact protein molecules like albumin. Especially the possibility of terpenoid-amino acid interactions that allow for binding and the creation of a physical barrier around glycation sites could only be tested in systems involving intact protein molecules.
Furthermore, if dose-response is a factor of interest, then future research should focus either on larger terpenoid concentrations than the ones tested herein, or even more to study terpenoid concentrations potentially smaller than 2 µM and understand their impact on protein glycation. As essential oils are rarely added in foods, small dose of these terpenoid compounds have a larger relevance to the conditions of thermal processing of foods. As terpenoid rich herbs and spices are often added during the thermal processing of foods, the final food is enriched in the terpenoid compounds of the herb/spice but during cooking the transfer of lipophilic compounds from the herb/spice to the food system is likely to be small. The addition of more objective markers of glycation like the measurement of dicarbonyl compounds and cross-link agents, even the production of reactive oxygen species would be an important addition to better describe the chemical transformations taking place in these systems.
Despite those limitations, this study is still valuable as it is to our knowledge the first study to simultaneously study the impact of five experimental conditions (glucose, temperature, time, terpenoid compounds and terpenoid compound dose) on lysine glycation. It is this design that allowed the detection of the opposing behavior of terpenoid compounds on glycation markers based on the experimental conditions and to pose novel research questions that could further explore the potential of terpenoid compounds to be used as antiglycative agents during thermal processing.

Conclusions
Terpenoid compounds commonly found in essential oils have the capacity to regulate Maillard reactions observed as a promotion of color development and lysine depletion in model systems. The qualitative impact of those compounds on specific stages of the Maillard reaction, however, is more complex including both anti-and pro-glycative effects which are dependent on the incubation conditions. Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/app12020908/s1, Table S1: Raw values of absorbance at 280 nm at 1, 2, 3 h for model systems with various concentrations of terpenoid compounds, in two different glucose concentrations and in two different temperatures; Table S2: Raw values of absorbance at 360 nm at 1, 2, 3 h for model systems with various concentrations of terpenoid compounds, in two different glucose concentrations and in two different temperatures; Table S3: Raw values of absorbance at 420 nm at 1, 2, 3 h for model systems with various concentrations of terpenoid compounds, in two different glucose concentrations and in two different temperatures; Table S4: Raw values of fluorescence intensity at λexc = 370 nm, λemm = 440 nm at 1, 2, 3 h for model systems with various concentrations of terpenoid compounds, in two different glucose concentrations and in two different temperatures.