Emulsification of Rosemary and Oregano Aqueous Extracts and Their In Vitro Bioavailability

Due to their richness in phenolic compounds, Mediterranean plants such as rosemary and oregano are increasingly recommended for consumption for their numerous health benefits. The pH shift and the presence of digestive enzymes significantly reduce the bioavailability of these biochemicals as they pass through the gastrointestinal tract. To prevent this degradation of phenolic compounds, methods such as emulsification of plant aqueous extracts are used. The aim of this study was to investigate the effects of emulsification conditions on the chemical properties (total polyphenolic content and antioxidant activity) of emulsified rosemary and oregano extracts. Response surface methodology was applied to optimize sunflower oil concentration, rotational speed, and emulsifier concentration (commercial pea protein). The emulsions prepared under optimal conditions were then used in bioavailability studies (in vitro digestion). The antioxidant activity of the emulsified rosemary/oregano extracts, measured by the DPPH method, remained largely stable when simulating in vitro digestion. Analysis of antioxidant activity after in vitro simulation of the gastrointestinal system revealed a higher degree of maintenance (up to 76%) for emulsified plant extracts compared to aqueous plant extracts. This article contributes to our understanding of how plant extracts are prepared to preserve their biological activity and their application in the food industry.


Introduction
Recently, the market has become more interested in aromatic herbs from the Mediterranean region, such as rosemary and oregano, due to their high content of bioactive compounds and consequent medicinal properties [1]. These plants have been used since ancient times as spices in the traditional Mediterranean diet and as medicines in folk medicine. Nowadays, they are considered important factors for food safety and quality because they inhibit various factors that affect food quality, such as oxidation and microbial spoilage [2,3].
Oregano (Origanum vulgare L.) and rosemary (Rosmarinus officinalis L.) have been the subject of numerous studies due to their significant antioxidant potential and their role in the prevention of diseases related to oxidative stress [4][5][6][7]. Various factors such as environmental conditions, food processing, storage, and digestive system conditions reduce polyphenol concentrations and result in a reduction or complete loss of their bioactivity [8]. In addition, many polyphenolic compounds from natural sources have limited water solubility and unpleasant taste [9]. Therefore, the application of phenolic compounds requires a protective formulation, such as encapsulation [9]. Emulsions are considered the most popular encapsulation system due to their high encapsulation efficiency, preservation of chemical stability of the encapsulated molecules, and controlled release of molecules [10]. For example, Tian et al. [11] presented the protective effect of emulsification with a xanthanlocustane gum mixture on tea polyphenols, while Ye et al. [12] prepared thermodynamically stable water-in-oil emulsions loaded with tea polyphenols using zein as a stabilizer. Bamba coverage of oil droplets, resulting in a decrease in the emulsion droplet size. For the emulsified rosemary extracts (Figure 2), the highest average Feret diameter (57.96 µ m) was determined for sample No. 10, while the lowest average Feret diameter was determined for sample No. 9 (25.73 µ m). For emulsified rosemary extracts with emulsifier concentrations of 0.1% and 0.5%, the lowest average Feret's diameter values were measured. From the obtained results, it could be concluded that emulsified oregano extracts had slightly higher average Feret diameter values compared to emulsified rosemary extracts.  For the emulsified oregano extracts (Figure 2), the highest average Feret diameter (97.97 µm) was determined for sample No. 9, while the lowest average Feret diameter was determined for sample No. 7 (28.56 µm). For emulsified oregano extracts, the lowest values of average Feret's diameter were obtained for emulsions containing 1% commercial pea protein, while the highest average Feret's diameter values were observed for emulsions containing 0.1% emulsifier concentration. Karadag et al. [30] and Sohn et al. [31] also reported a decrease in an emulsion droplet size with increasing emulsifier concentration. According to Sohn et al. [31], higher emulsifier concentration provides greater surface coverage of oil droplets, resulting in a decrease in the emulsion droplet size. For the emulsified rosemary extracts (Figure 2), the highest average Feret diameter (57.96 µm) was determined for sample No. 10, while the lowest average Feret diameter was determined for sample No. 9 (25.73 µm). For emulsified rosemary extracts with emulsifier concentrations of 0.1% and 0.5%, the lowest average Feret's diameter values were measured. From the obtained results, it could be concluded that emulsified oregano extracts had slightly higher average Feret diameter values compared to emulsified rosemary extracts.

Physical and Chemical Properties of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
Due to the growing interest in using medicinal plant extract formulations, it is necessary to determine the optimal process conditions for their preparation. In this work, the stability of prepared oil-in-water emulsions of rosemary/oregano extracts with the addition of commercial pea protein as an emulsifier was studied, and the effects of emulsifier concentration, rotational speed, and amount of oil in the emulsion were discussed. Emulsification of aqueous plant extracts was carried out based on the Box-Behnken experimental design. The physical and chemical properties of the obtained emulsified extracts are listed in Tables 1 and 2. The actual pea protein concentrations for rosemary and oregano extracts were determined using the Bradford method before the emulsification process began. For 0.1%, 0.5% and 1% pea protein suspensions in rosemary extract, the values were 74.91 mg L −1 , 484.81 mg L −1 and Average Feret diameters of oil-in-water emulsions containing rosemary extract and oil-inwater emulsions containing oregano extract. Different letters above the columns represent significant differences at p < 0.05.

Physical and Chemical Properties of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
Due to the growing interest in using medicinal plant extract formulations, it is necessary to determine the optimal process conditions for their preparation. In this work, the stability of prepared oil-in-water emulsions of rosemary/oregano extracts with the addition of commercial pea protein as an emulsifier was studied, and the effects of emulsifier concentration, rotational speed, and amount of oil in the emulsion were discussed. Emulsification of aqueous plant extracts was carried out based on the Box-Behnken experimental design. The physical and chemical properties of the obtained emulsified extracts are listed in Tables 1 and 2. The actual pea protein concentrations for rosemary and oregano extracts were determined using the Bradford method before the emulsification process began. For 0.1%, 0.5% and 1% pea protein suspensions in rosemary extract, the values were 74.91 mg L −1 , 484.81 mg L −1 and 767.49 mg L −1 , respectively, while for 0.1%, 0.5% and 1% pea protein suspensions in oregano extract, the values were 279.86 mg L −1 , 682.69 mg L −1, and 1283.39 mg L −1 , respectively.
Zeta potential or electrokinetic potential is a parameter that describes the electrochemical equilibrium on a solid surface and its liquid medium. It is of practical importance in controlling colloidal stability and flocculation processes. In the literature, emulsions with measured zeta potential values lower than −30 mV or higher than +30 mV are considered stable [32]. The higher the zeta potential, the stronger the electrostatic repulsion and the more stable the system [33], indicating the applicability of the zeta potential in predicting and controlling emulsion stability [34]. The measured zeta potential and conductivity values for the emulsified extracts are listed in Table 1. According to the results obtained for the emulsion of oil in aqueous rosemary extract, sample No. 8 (oil concentration 25%; emulsifier concentration 0.5%; rotational speed 35,000 rpm) had the lowest zeta potential value (−67.22 ± 2.24 mV), indicating the highest stability of the prepared oil-inwater emulsion containing rosemary extract. Sample No. 12 (oil concentration 15%; emulsifier concentration 1%; rotational speed 35,000 rpm) has the highest zeta potential value (−38.52 ± 3.28 mV), indicating the lowest stability of the prepared oil-in-water emulsion Average Feret diameters of oil-in-water emulsions containing rosemary extract and oil-inwater emulsions containing oregano extract. Different letters above the columns represent significant differences at p < 0.05. Table 1. Physical properties of oil-in-water emulsions containing rosemary/oregano extract with the addition of commercial pea protein as emulsifier. A,B The same superscript capital letters within a row denote no significant differences (p > 0.05) between the values obtained for the different emulsions under the same experimental conditions according to Tukey's ANOVA. a-j The same superscript lowercase letters within a column denote no significant differences (p > 0.05) between the values obtained for the emulsions under the different experimental conditions according to Tukey's ANOVA.

Sample
Zeta Potential (mV) Conductivity (mS cm −1 )  Zeta potential or electrokinetic potential is a parameter that describes the electrochemical equilibrium on a solid surface and its liquid medium. It is of practical importance in controlling colloidal stability and flocculation processes. In the literature, emulsions with measured zeta potential values lower than −30 mV or higher than +30 mV are considered stable [32]. The higher the zeta potential, the stronger the electrostatic repulsion and the more stable the system [33], indicating the applicability of the zeta potential in predicting and controlling emulsion stability [34]. The measured zeta potential and conductivity values for the emulsified extracts are listed in Table 1. According to the results obtained for the emulsion of oil in aqueous rosemary extract, sample No. 8 (oil concentration 25%; emulsifier concentration 0.5%; rotational speed 35,000 rpm) had the lowest zeta potential value (−67.22 ± 2.24 mV), indicating the highest stability of the prepared oil-in-water emulsion containing rosemary extract. Sample No. 12 (oil concentration 15%; emulsifier concentration 1%; rotational speed 35,000 rpm) has the highest zeta potential value (−38.52 ± 3.28 mV), indicating the lowest stability of the prepared oil-in-water emulsion with rosemary extract. The zeta potential of the emulsions is affected by both the concentration of the vegetable oil and the concentration of the surfactant. Based on the present results, it can be concluded that an increase in the oil concentration increases the zeta potential, indicating that the fatty acids present in the oil have an important influence on the zeta potential. The results obtained are in agreement with those of Rezvani et al. [35], where increasing the ratio of oil phase to water phase significantly increased the zeta potential of an orange oil-in-water emulsion. For oil-in-water emulsions containing oregano extract, the lowest zeta potential value was also measured for sample No. 8 (−57.18 ± 4.47 mV), indicating the highest stability of the emulsified oregano extract, while sample No. 7 (oil concentration 10%; emulsifier concentration 0.5%; rotational speed 35,000 rpm) showed the lowest stability (−37.71 ± 2.31 mV). The zeta potential range measured for all 17 experiments indicates the stability of all prepared emulsified rosemary and oregano extracts (values below −30 mV). These results are in agreement with the results obtained by Hinderink et al. [36], who measured the zeta potential of emulsions containing 1% purified pea protein. The results showed negative zeta potential values or less than −30 mV even after two weeks, indicating the stability of the prepared emulsions. Oil-in-water emulsions with rosemary/oregano extract prepared with a natural emulsifier (commercial pea protein) in this work showed higher stability than oil-in-water emulsions with mint extract prepared with commercial emulsifiers (polyethylene glycol 1500, 6000, and 20,000). The zeta potential values of the emulsions prepared with polyethylene glycol of different molecular masses did not exceed values less than −32.128 mV [33]. Conductivity values (Table 1), measured for emulsified rosemary extracts, ranged from 0.0043 ± 0.00 mS cm −1 (sample No. 11: oil concentration 15%; emulsifier concentration 0.1%; rotational speed 35,000 rpm) to 0.0453 ± 0.05 mS cm −1 (sample No. 8) while for emulsified oregano extracts, ranged from 0.0051 ± 0.00 mS cm −1 (sample No. 1: oil concentration 10%; emulsifier concentration 0.1%; rotational speed 25,000 rpm) to 0.0499 ± 0.05 mS cm −1 (sample No. 6: oil concentration 25%; emulsifier concentration 0.5%; rotational speed 15,000 rpm). The highest conductivity value was obtained for emulsified oregano extract (sample No. 6: 0.0499 ± 0.05 mS cm −1 ) and emulsified rosemary extract (sample No. 8: 0.0453 ± 0.05 mS cm −1 ).

Rosemary
Total polyphenolic content in the prepared emulsions was determined by the Folin-Ciocalteu method, while antioxidant activity was determined by the DPPH and FRAP methods ( Table 2). The highest TPC (53.93 ± 0.76 mg GAE g dw −1 ) was measured for the emulsified rosemary extract marked as sample No. 10 (oil concentration 15%; emulsifier concentration 1%; rotational speed 15,000 rpm), while the lowest TPC (1.35 ± 1.14 mg GAE g dw −1 ) was observed for sample No. 2 (oil concentration 25%; emulsifier concentration 0.1%; rotational speed 25,000 rpm). For the emulsified oregano extracts, the highest TPC (113.79 ± 1.90 mg GAE g dw −1 ) was measured for sample No. 3 (oil concentration 10%; emulsifier concentration 1%; rotational speed 25,000 rpm), while the lowest TPC (6.68 ± 0.00 mg GAE g dw −1 ) was measured for sample No. 11. Vallverdù-Queralt et al. [37] studied the polyphenolic profile of the most commonly used herbs as spices and reported the highest TPC values for rosemary, thyme, and oregano. The results of this study showed that the emulsified oregano extracts had higher TPC values than the emulsified rosemary extracts. According to the USDA National Nutrient Database [38], fresh rosemary leaves (nutritional values per 100 g) contain 21.8 mg of vitamin C and 3.31 g of protein, while dried oregano leaves (values per 100 g) contain 2.3 mg of vitamin C and 9 g of protein.
The differences in TPC values of the emulsified extracts could be due to the compounds such as ascorbic acid and proteins reacting with Folin-Ciocalteu reagent and interfering with the determination of TPC by the Folin -Ciocalteu assay [39,40]. Phenolic compounds in spicy vegetables are responsible for their antioxidant properties [41]. According to the results presented in Table 2, the highest antioxidant activity value determined by the DPPH method (0.35 ± 0.04 mmol Trolox g dw −1 ) was obtained for the emulsified rosemary extract No. 12, while the lowest DPPH value (0.06 ± 0.01 mmol Trolox g dw −1 ) was obtained for sample No. 11. Emulsified oregano extracts showed slightly higher DPPH values, compared to emulsified rosemary extracts. The highest DPPH value (0.83 ± 0.06 mmol Trolox g dw −1 ) was observed for emulsified oregano extract No. 8, while the lowest value (0.30 ± 0.03 mmol Trolox g dw −1 ) was found for sample No. 6. The highest FRAP value (0.49 ± 0.03 mmol FeSO4·7H2O g dw −1 ) was obtained for the emulsified rosemary extract, marked as sample No. 3, while the lowest FRAP value was recorded for sample No. 9 (oil concentration 15%; emulsifier concentration 0.1%; rotational speed 15,000 rpm). For emulsified oregano extracts, the highest FRAP (1.01 ± 0.08 mmol FeSO4·7H2O g dw −1 ) was measured for sample No. 10, while the lowest FRAP value (0.06 ± 0.01 mmol FeSO4·7H2O g dw −1 ) was obtained for sample No. 1. As for the antioxidant activity measured by DPPH method, the emulsified oregano extracts showed higher antioxidant activity measured by the FRAP method than the emulsified rosemary extracts. The antioxidant activity of the emulsified plant extracts depended on both the polyphenolic content of the aqueous extracts used [42] and the phenolic content of the oil phase used [43,44]. As noted by Zeb [45] the sunflower oil used for emulsification in this work is rich in linoleic acid and oleic acid as well as various antioxidants, including phenolic molecules. The positive effect of the oil phase on the antioxidant activity of the prepared elusions was also described by Thanonkaew et al. [46]. They analyzed the effect of the concentration of cold press rice bran oil (CPRBO) in emulsions and showed that increasing the oil concentration in the emulsion increased the antioxidant activity.

Response Surface Modeling of the Emulsification Process
Second-degree polynomials with interacting members were used to describe the experimental data. The estimated model parameters and results from analysis of variance  Tables 3 and 4, respectively, while the response surfaces of the significant variables affecting TPC, DPPH, and FRAP for both emulsified plant extracts are shown in Figures 3 and 4, respectively. ANOVA was applied to investigate the effects of process variables on the chemical properties of oil-in-water emulsions containing rosemary/oregano extract. The regression parameters are presented together with the coefficients of determination in Tables 3 and 4. For the emulsified rosemary extracts (Table 3), the emulsifier concentration and the rotational speed had a significant effect (p < 0.05) on the TPC. The obtained positive value for the linear parameter indicates that an increase in the commercial pea protein concentration causes an increase in the TPC value for emulsified rosemary extracts. A negative value for the linear parameter indicates that an increase in rotational speed results in a decrease in the TPC value of the emulsions. Oil concentration had no statistically significant effect on the TPC value (p > 0.05). From the calculated coefficient of determination (R 2 = 0.94; R 2 adj. = 0.90), it can be concluded that the influence of the process parameters on the TPC value can be excellently described by the developed predictive models. Table 3. Response surface models regression coefficients and determination coefficients for chemical properties of oil-in-water emulsions containing rosemary extract with the addition of commercial pea protein as an emulsifier (bold values: p < 0.05).

Analyzed Property Input Variable Intercept ± Standard Error
Regression Coefficients (L-Linear,

Optimization of the Emulsification Process with Respect to Chemical Properties of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract with the Addition of Emulsifier
Response Surface Methodology (RSM) was applied to define optimal conditions for the preparation of oil-in-water emulsions containing rosemary/oregano extract with the addition of commercial pea protein as an emulsifier. The influence of three independent variables (process parameters) such as oil concentration in emulsified plant extracts, emulsifier concentration, and rotational speed was analyzed at three levels: (oil concentration = 10%, 15%, 25%; emulsifier concentration = 0.1%, 0.5%, 1%; rotational speed = 15,000 rpm, 25,000 rpm, 35,000 rpm). The duration of emulsification was 4 min. Optimization was per- Regarding the antioxidant activity determined by the DPPH method, the influence of the process parameters on the DPPH value of the emulsified rosemary extracts cannot be described by the prediction models developed (R 2 = 0.35; R 2 adj . = 0.00). The studied process parameters had no significant effect (p > 0.05) on the DPPH values determined in the experiment. Regarding the antioxidant activity determined by the FRAP method, an increase in oil concentration and rotational speed leads to a decrease in the FRAP value of the emulsified rosemary extracts, while an increase in emulsifier concentration leads to an increase in the FRAP values. This can be explained by the antioxidant potential of the emulsifier used. This means that a higher concentration of the emulsifier leads to a higher antioxidant activity of the prepared emulsion. Based on the positive value obtained for the quadratic parameter, the rotational speed had a statistically significant effect (p < 0.05) on the FRAP, compared to the oil and emulsifier concentration (p < 0.05). The calculated coefficient of determination (R 2 = 0.76; R 2 adj. = 0.62) indicates a good relationship between the FRAP and the process variables.
For the oil-in-water emulsions containing oregano extract (Table 4), the studied process parameters had a significant effect (p < 0.05) on the obtained TPC values. An increase in the emulsifier concentration leads to a significant increase in the TPC value of the emulsified oregano extracts (the largest positive value of the linear regression coefficient was predicted). Based on the calculated coefficient of determination (R 2 = 0.99; R 2 adj. = 0.99), it can be concluded that the influence of the process parameters on the TPC value of the emulsions can be excellently described by the developed predictive models. The concentration of commercial pea protein and the rotational speed had a significant influence on the DPPH value of the emulsified oregano extracts (p < 0.05) compared to the influence of oil concentration (p > 0.05). From the positive values of the linear regression coefficients, it can be seen that an increase in emulsifier concentration and rotational speed leads to an increase in the DPPH value of the emulsions. According to the coefficient of determination (R 2 = 0.76; R 2 adj. = 0.60), there is a strong relationship between the experimental DPPH values and the parameters of the emulsification process. Regarding the antioxidant activity determined by the FRAP method, emulsifier and oil concentration significantly influenced the FRAP value of emulsified oregano extracts (p < 0.05) compared to the influence of rotational speed (p > 0.05). Increasing the concentration of the commercial pea protein led to an increase in the FRAP value of the emulsified extracts. The determination coefficient calculated for the FRAP value (R 2 = 0.91; R 2 adj. = 0.84) indicates a strong relationship between FRAP and the process parameters, which means that the predictive models excellently describe the experimental FRAP values of the emulsified oregano extracts.

Optimization of the Emulsification Process with Respect to Chemical Properties of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract with the Addition of Emulsifier
Response Surface Methodology (RSM) was applied to define optimal conditions for the preparation of oil-in-water emulsions containing rosemary/oregano extract with the addition of commercial pea protein as an emulsifier. The influence of three independent variables (process parameters) such as oil concentration in emulsified plant extracts, emulsifier concentration, and rotational speed was analyzed at three levels: (oil concentration = 10%, 15%, 25%; emulsifier concentration = 0.1%, 0.5%, 1%; rotational speed = 15,000 rpm, 25,000 rpm, 35,000 rpm). The duration of emulsification was 4 min. Optimization was performed simultaneously for the chemical properties of the emulsions based on the desirability profiles obtained from the RSM-predicted values. The desirability scale in the range of 0 (undesirable) to 1 (very desirable) was used [47].
Based on the obtained results from Figure 5, the optimal process conditions for the emulsification of plant extracts, concerning chemical properties, for the emulsified rosemary extracts, were: oil concentration = 16.17%, emulsifier concentration = 0.52% and rotational T with values of TPC = 18.59 mg GAE g dw −1 , DPPH = 0.2 mmol Trolox g dw −1 and FRAP = 0.37 mmol FeSO4·7H2O g dw −1 .
The optimal experimental value for DPPH was the same as the RSM-predicted value, while values for TPC and FRAP slightly differed from the RSM-predicted TPC and FRAP values. The model's predicted values were: TPC = 27.64 mg GAE g dw −1 , DPPH = 0.2 mmol Trolox g dw −1 , and FRAP = 0.31 mmol FeSO4·7H2O g dw −1 .
relationship between TPC and DPPH values and emulsifier concentration is proportional. The emulsifier concentration had a minor impact on FRAP results. While having no impact on DPPH, increasing rotational speed from 15,000 rpm to 35,000 rpm enhanced TPC. From the minimum to maximum rotational speed, the FRAP value was significantly reduced ( Figure 5).  Considering the emulsified oregano extracts (Figure 6), an increase in the oil concentration and rotational speed did not influence TPC, while the increase in the emulsifier concentration from 0.1% to 1% significantly increased TPC value. The increase in the oil concentration influenced the DPPH value, and the increase in the emulsifier concentration and rotational speed resulted in a significant DPPH increase. An increase in the oil concentration significantly decreased FRAP value, while an increase in the concentration of commercial pea protein significantly increased FRAP value. Rotational speed did not influence the FRAP value.
The investigations by Mehmood et al. [48] and Chong et al. [49] provided more evidence for the successful application of RSM in the assessment of parameters affecting emulsification. Mehmood and colleagues investigated the effects of olive oil concentrations, surfactant levels, and ultrasonic homogenization time on the physicochemical characteristics of generated olive-oil-based nanoemulsions for the encapsulation of α-tocopherol using a food-grade non-ionic surfactant. A significant influence of independent variables was observed on the response variables. They concluded that the particle size of nanoemulsions mainly depended on homogenization time, while surfactant concentration had a significant effect on antioxidant activity. Chong et al. [49] prepared The optimal process conditions for the emulsified oregano extracts were as follows ( Figure 6): oil concentration = 16.25%, emulsifier concentration = 0.52%, and rotational speed = 25,000 rpm with values of TPC = 57.55 mg GAE g dw −1 , DPPH = 0.31 mmol Trolox g dw −1 , and FRAP = 0.58 mmol FeSO4·7H2O g dw −1 . The optimal experimental values for TPC and FRAP were close to the RSM predicted values, while the optimal experimental value for DPPH significantly differed from the predicted value. Model predicted values were TPC = 60.23 mg GAE g dw −1 , DPPH = 0.57 mmol Trolox g dw −1 , and FRAP = 0.54 mmol FeSO4·7H2O g dw −1 . Based on the obtained results for the optimization of emulsification, increasing the oil concentration up to 25% did not affect TPC and DPPH values of the emulsified rosemary extracts but significantly decreased FRAP values. TPC and DPPH values were considerably affected by increasing the commercial pea protein content from 0.1% to 1%; the relationship between TPC and DPPH values and emulsifier concentration is proportional. The emulsifier concentration had a minor impact on FRAP results. While having no impact on DPPH, increasing rotational speed from 15,000 rpm to 35,000 rpm enhanced TPC. From the minimum to maximum rotational speed, the FRAP value was significantly reduced ( Figure 5).
Considering the emulsified oregano extracts (Figure 6), an increase in the oil concentration and rotational speed did not influence TPC, while the increase in the emulsifier concentration from 0.1% to 1% significantly increased TPC value. The increase in the oil concentration influenced the DPPH value, and the increase in the emulsifier concentration and rotational speed resulted in a significant DPPH increase. An increase in the oil concentration significantly decreased FRAP value, while an increase in the concentration of commercial pea protein significantly increased FRAP value. Rotational speed did not influence the FRAP value.
nanoemulsions consisting of palm oil, rich in tocopherol, with the addition of commercial emulsifiers (Tween 80, Span 80). The effects of emulsifier concentration, solvent amount, and homogenization pressure were investigated toward droplet size and polydispersity index. The actual values of the nanoemulsion prepared at optimal conditions were in good agreement with the predicted values obtained from RSM.

The In Vitro Digestion of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
The oil-in-water emulsions containing rosemary/oregano extract, prepared at the optimal conditions determined by the RSM were further used for in vitro simulation of the gastrointestinal system. After in vitro simulation, physical (TDS, conductivity) and chemical (TPC, DPPH, FRAP) properties of the plant extracts and emulsified plant extracts were measured. Before starting in vitro digestion, initial concentrations of the samples were obtained by measuring their physical and chemical properties (Table 5).
According to the definition, electrical conductivity is a measure of how easily an electric current can pass through water or some other solution or suspension. It is influenced by parameters such as the presence of weak electrolytes, pH, and temperature [26]. Conductivity could be used for determining the antioxidant activity of aqueous plant extracts [50,51]. According to the results presented in Table 5, initial conductivity and TDS values were higher for plant extracts compared to emulsions. The highest initial TDS and conductivity were measured for oregano extract (0.31 ± 0.01 g L −1 and 0.63 ± 0.00 mS cm −1 ) and rosemary extract (0.09 ± 0.00 g L −1 and 0.18 ± 0.00 mS cm −1 ) compared to emulsified oregano extract (0.04 ± 0.00 g L −1 and 0.01 ± 0.04 mS cm −1 ) and emulsified rosemary extract (0.02 ± 0.00 g L −1 and 0.04 ± 0.00 mS cm −1 ). The investigations by Mehmood et al. [48] and Chong et al. [49] provided more evidence for the successful application of RSM in the assessment of parameters affecting emulsification. Mehmood and colleagues investigated the effects of olive oil concentrations, surfactant levels, and ultrasonic homogenization time on the physicochemical characteristics of generated olive-oil-based nanoemulsions for the encapsulation of α-tocopherol using a food-grade non-ionic surfactant. A significant influence of independent variables was observed on the response variables. They concluded that the particle size of nanoemulsions mainly depended on homogenization time, while surfactant concentration had a significant effect on antioxidant activity. Chong et al. [49] prepared nanoemulsions consisting of palm oil, rich in tocopherol, with the addition of commercial emulsifiers (Tween 80, Span 80). The effects of emulsifier concentration, solvent amount, and homogenization pressure were investigated toward droplet size and polydispersity index. The actual values of the nanoemulsion prepared at optimal conditions were in good agreement with the predicted values obtained from RSM.

The In Vitro Digestion of Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
The oil-in-water emulsions containing rosemary/oregano extract, prepared at the optimal conditions determined by the RSM were further used for in vitro simulation of the gastrointestinal system. After in vitro simulation, physical (TDS, conductivity) and chemical (TPC, DPPH, FRAP) properties of the plant extracts and emulsified plant extracts were measured. Before starting in vitro digestion, initial concentrations of the samples were obtained by measuring their physical and chemical properties (Table 5). Table 5. Physical properties of rosemary/oregano extracts and oil-in-water emulsions containing rosemary/oregano extract, with the addition of commercial pea protein as emulsifier, before and during in vitro digestion. A,B The same superscript capital letters within a row denote no significant differences (p > 0.05) between the values obtained for the different plant aqueous extracts and emulsions at the same experimental conditions according to Tukey's ANOVA. a-d The same superscript lowercase letters within a column denote no significant differences (p > 0.05) between values during in vitro digestion according to Tukey's ANOVA.

Sample
Total Dissolved Solids (g L −1 ) Conductivity (mS cm −1 ) According to the definition, electrical conductivity is a measure of how easily an electric current can pass through water or some other solution or suspension. It is influenced by parameters such as the presence of weak electrolytes, pH, and temperature [26]. Conductivity could be used for determining the antioxidant activity of aqueous plant extracts [50,51]. According to the results presented in Table 5, initial conductivity and TDS values were higher for plant extracts compared to emulsions. The highest initial TDS and conductivity were measured for oregano extract (0.31 ± 0.01 g L −1 and 0.63 ± 0.00 mS cm −1 ) and rosemary extract (0.09 ± 0.00 g L −1 and 0.18 ± 0.00 mS cm −1 ) compared to emulsified oregano extract (0.04 ± 0.00 g L −1 and 0.01 ± 0.04 mS cm −1 ) and emulsified rosemary extract (0.02 ± 0.00 g L −1 and 0.04 ± 0.00 mS cm −1 ).

Rosemary
After in vitro digestion, TDS and conductivity of plant extracts and emulsions increased due to the addition of digestive enzymes, HCl, and bile salts throughout the whole simulation process of a gastrointestinal system [19] (Table 5). According to Figure 7, aqueous plant extracts had higher TPC values, concerning their emulsions, before in vitro simulation of gastrointestinal digestion. The initial TPC value for oregano extract was 100.27 ± 7.71 mg GAE g dw −1 , and for rosemary extract, the TPC value was 35.00 ± 0.19 mg GAE g dw −1 . The initial TPC value for emulsified oregano extract was 40.95 ± 5.19 mg GAE g dw −1 , while for the emulsified rosemary extract, the TPC value was 2.70 ± 0.00 mg GAE g dw −1 . The total polyphenolic content of aqueous plant extracts and emulsified plant extracts decreased during in vitro digestion. This TPC decrease was the most significant for the aqueous oregano extract and for the emulsified oregano extract (Figure 7b) compared to rosemary (Figure 7a). Regarding oregano extract and emulsified oregano extract, TPC decreased by approximately 95%, while for the emulsified rosemary extract, about 32% of TPC decreased.   TPC values for plant extracts and emulsified plant extracts, after the third step of in vitro digestion (duodenum), increased compared to their values measured in the second step (stomach) (Figure 7). Tarko and Duda-Chodak [52] reported a higher polyphenol value and antioxidant activity after finished digestion due to the action of digestive enzymes and pH changes in the stomach and duodenum, leading to polyphenols decomposition to components with higher antioxidant potential.
Results obtained for the antioxidant activity determined by the DPPH method showed higher initial values measured for plant extracts compared to emulsified plant extracts (Figure 8), which was expected concerning their initial TPC values (Figure 7). The highest DPPH decrease after the first digestion step was observed for oregano extract in comparison with its initial value (0.82 ± 0.06 mmol Trolox g dw −1 ), while DPPH value measured for the emulsified oregano extract was completely preserved (0.17 ± 0.01 mmol Trolox g dw −1 ). After the second digestion step, DPPH values started to increase, for all analyzed samples (Figure 8), as a result of pH influence and the difference in the ability of plants to reduce DPPH due to different portions of bioactive components [53]. The highest DPPH in the stomach was measured for oregano extract (0.25 ± 0.17 mmol Trolox g dw −1 ). After the third step, DPPH values for all samples decreased, with the exception of emulsified oregano extract (Figure 8b).    Generally speaking, DPPH values measured for oregano extract and emulsified oregano extract were higher compared to rosemary extract and emulsified rosemary extract. It is evident that plant emulsions prepared under optimal conditions showed better preservation of the antioxidant activity determined by the DPPH method compared to their aqueous extracts. Results obtained for the antioxidant activity determined by the FRAP method were not in agreement with the results obtained by the DPPH method during in vitro digestion. The reason could be the different profiles of polyphenolic components present in plants that can affect the ability to reduce certain radicals [47]. The higher initial FRAP values were measured for plant extracts (0.61 ± 0.02 mmol FeSO4·7H2O g dw −1 for rosemary and 1.38 ± 0.02 mmol FeSO4·7H2O g dw −1 for oregano, respectively). The initial FRAP values measured for emulsified plant extracts were 0.21 ± 0.01 mmol FeSO4·7H2O g dw −1 for rosemary and 0.56 ± 0.01 mmol FeSO4·7H2O g dw −1 for oregano. After the first digestion step (in the mouth), FRAP values significantly decreased for the rosemary extract and emulsified rosemary extract (Figure 9a). The same trend could also be observed for oregano samples (Figure 9b). After the second digestion step (in the stomach), FRAP values, for both extracts and emulsions, did not significantly change compared to the first step (Figure 9). At the end of in vitro simulation of gastrointestinal digestion (duodenum), the highest FRAP value was measured for rosemary extract (0.1 ± 0.00 mmol FeSO4·7H2O g dw −1 ), while the lowest/same was measured for emulsified rosemary and oregano extracts (0.03 ± 0.00 mmol FeSO4·7H2O g dw −1 ). The obtained results show that the antioxidant activity measured by the FRAP method remained mostly conserved in emulsified rosemary extract.

Bioavailability of Polyphenolic Compounds from Rosemary/Oregano Extracts and Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
After in vitro simulation of the gastrointestinal system, the bioavailability of total phenolics from plant extracts and emulsified plant extracts was determined, and results are presented in Table 6. Bioavailability was calculated as the ratio of TPC (DPPH, FRAP) after the end of in vitro digestion (duodenum) and the initial concentration (before starting in vitro digestion).
Regarding TPC, the highest bioavailability was calculated for the emulsified rosemary extract (67.78%), while other samples show a significant decrease in TPC during in vitro digestion (Figure 7). The results indicate a significant influence of pH changes on phenolic compounds during in vitro digestion [26,54]. Gutiérrez-Grijalva et al. [54] pointed out that the phenolic compounds' stability under various pH conditions, such as those present during gastrointestinal digestion, is influenced by the composition of the molecules and the distribution of the -OH radicals on their rings. They also stated that the loss of certain phenolic compounds following the GI process has been linked, in addition to pH changes, to potential interactions between polyphenols and other elements of the digestive fluids, such as enzymes and electrolytes, which requires additional research.
Regarding the antioxidant activity determined by the DPPH method, higher bioavailability values were observed for emulsified plant extracts compared to plant extracts

Bioavailability of Polyphenolic Compounds from Rosemary/Oregano Extracts and Oil-in-Water Emulsions Containing Rosemary/Oregano Extract
After in vitro simulation of the gastrointestinal system, the bioavailability of total phenolics from plant extracts and emulsified plant extracts was determined, and results are presented in Table 6. Bioavailability was calculated as the ratio of TPC (DPPH, FRAP) after the end of in vitro digestion (duodenum) and the initial concentration (before starting in vitro digestion).
Regarding TPC, the highest bioavailability was calculated for the emulsified rosemary extract (67.78%), while other samples show a significant decrease in TPC during in vitro digestion (Figure 7). The results indicate a significant influence of pH changes on phenolic compounds during in vitro digestion [26,54]. Gutiérrez-Grijalva et al. [54] pointed out that the phenolic compounds' stability under various pH conditions, such as those present during gastrointestinal digestion, is influenced by the composition of the molecules and the distribution of the -OH radicals on their rings. They also stated that the loss of certain phenolic compounds following the GI process has been linked, in addition to pH changes, to potential interactions between polyphenols and other elements of the digestive fluids, such as enzymes and electrolytes, which requires additional research. Table 6. Bioavailability of total polyphenols from rosemary/oregano extracts and oil-in-water emulsions containing rosemary/oregano extract after in vitro digestion, with the addition of commercial pea protein as emulsifier. Regarding the antioxidant activity determined by the DPPH method, higher bioavailability values were observed for emulsified plant extracts compared to plant extracts ( Table 6). This could be an indication that the desired effect of preserving antioxidant activity was achieved by the emulsification of plant extracts. The highest bioavailability was observed for emulsified oregano extract (76.47%), followed by emulsified rosemary extract (69.23%) and rosemary extract (42.11%). Results have shown that the antioxidant activity of emulsified plant extracts measured by the DPPH method remained mostly stable when simulating the in vitro digestion process (Figure 8). The bioavailability of the plant extracts and emulsified plant extracts for the FRAP method has shown significantly lower values compared to the DPPH values. Results have shown low preservation of the antioxidant activity determined by the FRAP method (Table 6, Figure 9) compared to the DPPH method.

Plant Materials, Sunflower Oil and Emulsifier
Dried rosemary (Rosmarinus officinalis L.) and organic pea protein powder were purchased from Nutrigold (Zagreb, Croatia), and dried oregano (Origanum vulgare L.) was purchased from SonnentoR (Sprögnitz, Austria). Edible sunflower oil (VitaDor, Germany) was obtained from a local supermarket. Plant materials were collected during the season of 2020 (oregano) and 2021 (rosemary), dried naturally, and stored at ambient conditions before use. The country of origin of rosemary is India, and the country of origin of oregano is Austria.

Preparation of Aqueous Plant Extracts
Dry rosemary/oregano (12 g) was mixed with 600 mL of deionized water, previously thermostated at 80 • C. Solid-liquid extraction was performed using Ika HBR4 digital oilbath (IKA-Werk GmbH & Co.KG, Staufen, Germany) at T = 80 • C and 250 rpm for 30 min. After extraction, samples were filtered, using a vacuum filtration system, through a 100% cellulose paper filter (LLG Labware, Meckenheim, Germany) with a 5-13 µm diameter pore size [50]. Aqueous plant extracts were further used for the preparation of oil-in-water emulsions containing plant extract with the addition of an emulsifier (organic pea protein powder).

Preparation of Pea Protein Suspensions in Aqueous Plant Extracts
Pea protein suspensions in aqueous plant extracts were prepared according to Sridharan et al. [55] with a slight modification of their method. In order to prepare 1% pea protein concentration, 4.5 g of pea protein powder from controlled organic cultivation was suspended in the 450 mL of aqueous rosemary/oregano extracts. Suspension was homogenized under magnetic stirring (500 rpm) at room temperature for several minutes. The pH value was adjusted to 7 by addition of NaOH (0.5 M). After pH adjustment, pea protein suspension was filtered using a filter paper with a cut off of 4-12 µm (Rundfilter, MN 640 m dia 11 cm, Macherey-Nagel, Düren, Germany) using vacuum filtration. Preparation of 0.1% and 0.5% pea protein concentrations in aqueous rosemary/oregano extracts was performed by diluting 1% pea protein concentration with distilled water.

Determination of the Concentration of Commercial Pea Protein According to the Bradford Method
The actual protein concentrations in 0.1%, 0.5%, and 1% pea protein suspensions were determined by the Bradford colorimetric method, using bovine serum albumin (BSA) as standard (1 mg mL −1 ) [56]. A 0.5 mL amount from each sample was mixed with 0.5 mL of Bradford reagent, and after 30 min incubation, the absorbance was measured at 595 nm using a UV-VIS spectrophotometer (Biochrom Libra S12, Cambridge, UK). Data were expressed as mg L −1 .

Oil-in-Water Emulsions Containing Plant Extract Preparation
Oil-in-water emulsions containing rosemary/oregano extract were prepared according to the conditions defined using Box-Behnken experimental design [57] (Table 7). A certain volume of aqueous rosemary/oregano extracts containing an appropriate emulsifier concentration (organic pea protein) was placed in a 15 mL falcon test tube with a certain volume of sunflower oil, in order to prepare 7 mL of emulsion. The mixture was homogenized using a homogenizer (OMNI TH220-PCRH homogenizator, Omni International, USA) for 4 min. The zeta potential and electrical conductivity of the prepared emulsions were measured immediately, while the rest of the samples were stored at 4 • C until further analysis for chemical properties (total polyphenolic content, antioxidant activity).
After the optimal conditions for the emulsification of aqueous plant extracts were determined, the same optimal conditions were used for preparation of oil-in-water emulsions containing rosemary/oregano extract, which were further used for in vitro simulation of the gastrointestinal system (GI).

Extraction Yield
Extraction yield was expressed as the dry matter content of the plant material. Dry matter was determined using a standard AOAC method 930.15 [58]. A certain amount of plant material samples (1-5 ± 0.0001 g) was dried in weighing dishes to constant weight at 105 • C, for 3 h. All measurements were performed in duplicate.

Determination of Physical Characteristics of Oil-in-Water Emulsions Containing Plant Extract Determination of Zeta Potential, Electrical Conductivity and Total Dissolved Solids
The surface charges of the emulsion droplets and electrical conductivity were measured using a Zetasizer Ultra instrument (Malvern Panalytical Limited ZSU3305, Malvern, UK). The emulsion samples were diluted two hundred times in distilled water. The diluted samples were then filled in a Malvern omega cuvette. The samples were placed in the instrument and equilibrated for 10 min before measurements were performed. All measurements were performed in triplicate.
Total dissolved solids (TDS) were determined using the SevenCompact conductometer (SevenCompact, Mettler Toledo, Switzerland). The probe was immersed into the sample (aqueous plant extracts or oil-in-water emulsions containing plant extract), and the values were read. Samples were homogenized before each measurement. All measurements were performed in duplicate.

Determination of the Average Feret's Droplet Diameter of Oil-in-Water Emulsions
Containing Rosemary/Oregano Extract According to Grgić et al. [59], Feret diameter is defined as the perpendicular distance between two tangents located on opposite sides of a particle. Samples of oil-in-water emulsions containing rosemary/oregano extract were photographed using a microscope equipped with a camera (BTC type LCD-35, Bresser, Germany) at 10x magnification. The average Feret diameter of the droplets was measured using the software tool ImageJ (v.1.8.0. National Institutes of Health, Bethesda, MD, USA). A total of 30 droplets were analyzed on each image. The average Feret's diameter was estimated as the average value of 30 measurements of the Feret's droplet diameter [59].
The reaction mixture was shaken and incubated for 5 min at 37 • C. Before the second step (stomach), it is necessary to adjust the pH value of the reaction mixture to 2 by adding HCl (37%). After pH adjustment, the second step starts with the addition of 60 mg of porcin-pepsin solution in 4 mL HCl (0.01 mol L −1 ). This mixture was shaken and incubated at 37 • C for 2 h. Before duodenal digestion, the pH value of the reaction mixture was adjusted to 6.5 by adding NaHCO 3 . The third step starts with the addition of pancreatin (0.08 g pancreatin dissolved in 10 mL phosphate buffer) and bile salts (0.5 g bile salts dissolved in 10 mL phosphate buffer). The reaction mixture was shaken and incubated at 37 • C for 2 h. Samples (1 mL) were taken after each digestion phase and analyzed for electrical conductivity, TDS, TPC, and the antioxidant activity (DPPH, FRAP) [26].

Determination of Chemical Characteristics of Aqueous Plant Extracts and Oil-in-Water Emulsions Containing Plant Extract Determination of Total Polyphenolic Content (TPC)
Total polyphenolic content (TPC) of aqueous rosemary/oregano extracts and oil-inwater emulsions containing rosemary/oregano extract was determined spectrophotometrically, according to Singleton and Rossi [61], using Folin-Ciocalteu reagent. Shortly, 7.9 mL of distilled water was mixed with 500 µL of Folin-Ciocalteu reagent (Folin-Ciocalteu reagent:water at a 1:2 ratio) and 100 µL sample. The reaction was started with addition of 1.5 mL of 20% Na 2 CO 3 solution. After 2 h of incubation in a dark place, the absorbance of the reaction mixture was measured at λ = 765 nm using spectrophotometer (Biochrom Libra S11, Cambridge, UK). Measurements were performed in duplicates and the results were derived from a calibration curve for gallic acid (GA) (0-500 mg L −1 ) and expressed as mg GA equivalents (GAE)/g dry weight (dw) plant material [50].

Determination of Antioxidant Activity by DPPH Method
Antioxidant activity measured by DPPH (1,1-diphenyl-2-picrylhydrazyl) scavenging method was performed as described by Brand-Williams et al. [62]. The reaction mixture consisted of 100 µL sample and 3.9 mL of DPPH radical (0.094 M) dissolved in methanol. After 30 min of incubation, the absorbance of the reaction mixture was measured at λ = 515 nm. Measurements were performed in duplicate, and the results were derived from a calibration curve for Trolox (0-1 mmol L −1 ) and expressed as mmol Trolox equivalents/g dry weight (dw) plant material [50].

Determination of Antioxidant Activity by FRAP Method
The FRAP (Ferric ion reducing antioxidant power) assay was carried out according to Benzie and Strain [63]. The reaction mixture consisted of 50 µL sample and 950 µL of FRAP reagent. The mixture was incubated for 4 min, and the absorbance was measured at λ = 593 nm. Measurements were performed in duplicate, and the results were derived from a calibration curve for FeSO 4 · 7H 2 O (0-1 mmol L −1 ) and expressed as mmol FeSO 4 · 7H 2 O equivalents/g dry weight (dw) plant material [50].

Response Surface Methodology (RSM)
Response Surface Methodology (RSM) was used to investigate the influence of three independent variables (oil concentration in emulsified plant extracts, emulsifier concentration, rotational speed) on the chemical characteristics (TPC, DPPH, FRAP) of emulsified extracts using software Statistica 14.0. (TIBCO ® Statistica, Palo Alto, CA, USA). The experiment was performed according to Box-Behnken experimental design [59]. Each independent variable was analyzed at three levels: oil concentration (10%, 15%, 25%), emulsifier concentration (0.1%, 0.5%, 1%), and rotational speed (15,000 rpm, 25,000 rpm, 35,000 rpm). The resulting experimental design comprised 17 experiments, as shown in Table 7. For the description of the experimental data obtained for chemical properties (TPC, DPPH, FRAP) of oil-in-water emulsions containing rosemary/oregano extract, prediction models (second-order polynomial models) were developed. In order to describe how well models fit a set of experimental data, the coefficient of determination (R 2 ) was calculated.

Statistical Analysis
The means of the results of physical and chemical properties of oil-in-water emulsions containing rosemary/oregano extract with the addition of commercial pea protein as emulsifier and physical and chemical properties of rosemary/oregano extracts and oil-inwater emulsions containing rosemary/oregano extract, with the addition of commercial pea protein as emulsifier, before and during in vitro digestion, were evaluated using analysis of variance (ANOVA) while Tukey's test was used to compare significant differences (p < 0.05) between the samples. The ANOVA was performed using Statistica 14.0. (TIBCO ® Statistica, Palo Alto, CA, USA). Differences between different plant aqueous extracts/emulsions at the same experimental conditions were analyzed, as well as the differences between different experimental conditions for the selected plant aqueous extracts/emulsions. Prior to statistical analysis, the normality of the data was evaluated using Kolmogorov-Smirnov test, implemented in Statistica software. p-values of the analyzed data sets were not significant, and therefore the assumption of normality can be accepted. Homogeneity of variance was assessed using Levene's test implemented in Statistica software. p-values of the analyzed data sets were higher than 0.05, indicating the homogeneity of the variances.

Conclusions
Phenolic compounds derived from plant sources have the potential to act as natural antioxidants and are relevant to food preservation and the prevention of human diseases. It is important to maintain the biological activity of these phytochemicals, and emulsification represents a good way of keeping their diverse biological effectiveness. To highlight the advantages of emulsion systems often used in the food industry, it is necessary to improve the process of emulsion preparation. The optimization of emulsification represents a good way to select corresponding process conditions in order to create emulsions with good nutritional, functional, and sensory properties, as well as biological activities. Additionally, the utilization of food-grade emulsions containing plant-based proteins as emulsifiers is also emerging since plants are recognized as rich sources of protein, intended for different purposes, with low production costs and efficient extraction protocols. Constant improvement in their functional properties (emulsification, foamability, solubility, gelation, hydration, miscibility, aroma trapping, and viscosity) is leading to their further applications in the food and pharmaceutical industries.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.