Optimization of Pea Protein Isolate-Stabilized Oil-in-Water Ultra-Nanoemulsions by Response Surface Methodology and the Effect of Electrolytes on Optimized Nanoemulsions

: Nanoemulsions are optically transparent and offer good stability, bioavailability, and control over the targeted delivery and release of lipophilic active components. In this study, pea protein isolate (PPI)-stabilized O/W nanoemulsions were evaluated using response surface methodology to obtain optimized ultra-nanoemulsions of Sauter mean diameter ( D 3,2 ) < 100 nm using a high-pressure homogenizer (HPH). Furthermore, the effect of food matrix electrolytes, i.e., the pH and ionic strength, on the emulsion (prepared at optimized conditions) was investigated. The results revealed that the droplet size distribution of emulsions was mainly inﬂuenced by the PPI concentration and the interaction of oil concentration and HPH pressure. Moreover, a non-signiﬁcant increase in droplet size was observed when the nanoemulsions (having an initial D 3,2 < 100 nm) were stored at 4 ◦ C for 7 days. Based on the current experimental design, nanoemulsions with a droplet size < 100 nm can effectively be prepared with a high PPI concentration (6.35%), with less oil (1.95%), and at high HPH pressure (46.82 MPa). Such emulsions were capable of maintaining a droplet size below 100 nm even at ionic conditions of up to 400 mM NaCl and at acidic pH.


Introduction
Consumers' preferences have drawn substantial scientific and industrial interest for the preparation of more natural and healthier foods [1][2][3]. The presence of both an aqueous phase and a lipid phase in foods is considered healthier, as the aqueous and lipid phases are responsible for the uptake of hydrophilic and lipophilic micronutrients in the gastrointestinal tract, respectively. As the majority of foods are naturally rich in aqueous phases, the dispersion of a certain amount of lipids in foods (e.g., soups, sauces, smoothies, purees, sausages, beverages, etc.) can significantly enhance the nutritional and sensory quality of foods by providing a suitable solvent system for lipophilic active components for enhanced stability, controlled release, and bioaccessibility [2][3][4][5][6][7][8]. The multi-phase system in which a liquid is dispersed in the continuous phase of (at least partly) immiscible liquid is referred to as an emulsion [9,10]. Systems with lipids/oils dispersed in the continuous aqueous phase are therefore referred to as oil-in-water (O/W) emulsions.
Emulsions with nanoscale droplet sizes are often referred to in the literature as mini emulsions, nanoemulsions, ultrafine emulsions, sub-micron emulsions, etc., and typically, those in the range of 20-100 nm radii are preferably referred to as nanoemulsions or ultrafine emulsions [3,11,12]. These droplets (r < 100 nm) are smaller than the wavelength of visible light (λ-390-750 nm), tend to be transparent or translucent to the naked eye, and possess stability against sedimentation or creaming [11][12][13]. Hence, these emulsions offer a promising application in the incorporation of lipophilic active components such

Materials
Commercial food grade pea protein isolates (PPI) (80% protein) from MYVEGAN TM UAE were used in this study. Sunflower oil (refractive index = 1.473) was purchased from a local grocery store. Deionized water (refractive index = 1.33) was used throughout the study.

Experimental Design and Evaluation
A central composite rotatable design (CCRD) with six center points was used to determine the combinations of three independent variables that were each present at five levels (−α, −1, 0, 1, and +α), as shown in Table 1. The center points and respective factorial points for PPI (5 ± 2%) and oil (7 ± 3%) were selected considering the solid-not-fat and fat percentage in milk of different origin as a reference, while the HPH pressure range (30 ± 10 MPa) was based on the results of previous studies [26,40]. Twenty experimental runs were generated, and three responses of droplet size distribution (i.e., Sauter mean diameter, median droplet diameter, and the span) were measured. The obtained responses were evaluated using response surface methodology. A quadratic polynomial model was fitted to each response using the following equation: where y is the response, β 0 is a constant, β i is the linear coefficient, β ii is the quadratic coefficient, β ij is the interaction coefficient, k is the number of independent variables, β i indicates the independent variables, and X i and X j indicate the difference between the independent variables and their respective center points. The optimization was carried out using the optimization settings presented in Table 2. Using the optimized factor conditions, the selected response was predicted and validated. ) was dispersed in deionized water, magnetically stirred (650 rpm) at room temperature for 60 min, and left undisturbed at 4 • C overnight. PPI solution was centrifuged (10,000× g; 30 min at 4 • C) to obtain the supernatant; the supernatant was centrifuged again under the previous conditions to obtain the final supernatant, which was then used as the aqueous phase for coarse emulsion preparation.
Sunflower oil (1.95-12.05%), which was added to the aqueous phase and coarse emulsion, was prepared via homogenization at 10,000 rpm for 2 min using ULTRA-TURRAX (T25 D, IKA ® -Werke GmbH & Co. KG, Staufen, Germany) with an 18 G probe. The nanoemulsions were produced by passing the coarse emulsions through a high-pressure homogenizer (FT90, Armfield Ltd., UK) 3 times with a pressure of 13. 18-46.82 MPa. In addition to the samples required for CCRD, the nanoemulsions were also prepared under optimum process conditions with different ionic strengths (0-400 mM NaCl in the aqueous phase) and with different pH values (3, 5, and 7 pH of the aqueous phase, adjusted using 1 N solutions of HCl and NaOH). The aqueous phases with different ionic strengths and pH values were also kept overnight at 4 • C prior to emulsification.

Emulsion Characterization
All nanoemulsions were characterized for droplet size distribution using a laser diffraction particle size analyzer (Malvern Mastersizer 3000, Malvern Panalytical Ltd., Malvern, UK). The refractive indices of 1.46 and 1.33 were entered into the machine's SOP for sunflower oil and continuous phase, respectively. An emulsion sample was slowly added (using a 3 mL bubble pipette) into the wet sample dispersion unit (Hydro EV) containing deionized water until a 10% obscuration level was reached, and then, the measurement cycle was run to record the droplet size distribution. From each measurement (an average of three readings), the Sauter mean diameter (D 3,2 ), median diameter (D 50 ), and span value (δ) were taken for further analysis. The Sauter mean diameter (also known as surface-volume mean diameter) can be defined as where n i and D i are the number and diameters of emulsion droplets in a particular size fraction [41], whereas the median diameter corresponds to a size when the cumulative droplet size distribution reaches 50% ( Figure 1). Furthermore, the span value was used to describe the width of droplet size distribution (i.e., the polydispersity of the emulsion), which can be defined as where D 10 , D 50 , and D 90 are the droplet sizes corresponding to a cumulative droplet size distribution of 10%, 50%, and 90%, respectively ( Figure 1). where D10, D50, and D90 are the droplet sizes corresponding to a cumulative droplet size distribution of 10%, 50%, and 90%, respectively ( Figure 1). The emulsions were then stored at refrigeration temperature (~4 °C) in capped glass bottles (500 mL). In order to determine the stability of emulsions, the samples were analyzed again for droplet size distribution after a week using the same protocol as described above.

Statistical Analysis
JMP PRO 15 software (SAS Institute Inc., North Carolina, USA) was used to develop the experimental design and to evaluate independent variables via response surface methodology and optimization. The Tukey-HSD test was used for the mean comparison wherever applicable.

Response Surface Evaluation
Response surface evaluation for each of the emulsion characteristics viz. the Sauter mean diameter (D3,2), median diameter (D50), and the span value (δ) was carried out, and the response surface coefficients, coefficient of determination (R 2 ), and significance of regression (p-value) were determined ( Table 3). The evaluation revealed that the main effects of PPI (%) and oil (%), the quadratic effect of PPI (%), and interaction between the oil (%) and HPH pressure (MPa) had significant effects (p < 0.05) on D3,2, whilst only the main effects of PPI (%) and HPH pressure (MPa) had significant effects (p < 0.05) on D50. No significant effect (p > 0.05) of the evaluated factors was observed for the span; however, the quadratic effect of PPI (%) was significant at a 10% level of significance. The interaction plots for each of the emulsion characteristics are presented in Figure 2. The emulsions were then stored at refrigeration temperature (~4 • C) in capped glass bottles (500 mL). In order to determine the stability of emulsions, the samples were analyzed again for droplet size distribution after a week using the same protocol as described above.

Statistical Analysis
JMP PRO 15 software (SAS Institute Inc., Cary, NC, USA) was used to develop the experimental design and to evaluate independent variables via response surface methodology and optimization. The Tukey-HSD test was used for the mean comparison wherever applicable.

Response Surface Evaluation
Response surface evaluation for each of the emulsion characteristics viz. the Sauter mean diameter (D 3,2 ), median diameter (D 50 ), and the span value (δ) was carried out, and the response surface coefficients, coefficient of determination (R 2 ), and significance of regression (p-value) were determined ( Table 3). The evaluation revealed that the main effects of PPI (%) and oil (%), the quadratic effect of PPI (%), and interaction between the oil (%) and HPH pressure (MPa) had significant effects (p < 0.05) on D 3,2 , whilst only the main effects of PPI (%) and HPH pressure (MPa) had significant effects (p < 0.05) on D 50 . No significant effect (p > 0.05) of the evaluated factors was observed for the span; however, the quadratic effect of PPI (%) was significant at a 10% level of significance. The interaction plots for each of the emulsion characteristics are presented in Figure 2. Table 3. Response surface coefficients, coefficient of determination (R 2 ), and significance of models (p-value) for the Sauter mean diameter (D 3,2 ), median droplet diameter (D 50 ), and span (δ) of pea protein isolate (PPI)-stabilized nanoemulsions via high-pressure homogenization (HPH).

Effect of PPI Concentration
In general, droplets' sizes reduce with higher PPI concentrations, and a similar behavior was observed in the current study. In fact, a limited particle concentration resulted in the insufficient interfacial coverage of droplets, which led to bridge flocculation and ultimately to coalescence [31,[42][43][44]. This was also evident from a higher span value (due to a wider droplet size distribution, or in some cases due to a bi-modal distribution) for those emulsion samples containing lower PPI concentrations. At extreme poor coverage of the interface, emulsification failure occurred [42], although this was not observed in any of the emulsification runs in the present study.
With the increase in PPI concentration, a gradual decrease in the droplets' size was observed. The decrease in droplet size with the increase in particle concentration has been reported previously [31,[45][46][47] and is valid until the maximum critical particle concentration is reached, beyond which, the increase in particle concentration contributes to the viscosity and density and subsequent emulsion stability [2,42,43]. For low oil content and high HPH pressure, the maximum critical particle concentration was ~6%, i.e., a further increase in PPI concentration did not show a further decrease in droplet size. However, for high oil content and low HPH pressure, the increase in PPI (beyond 6%) resulted in a slight increase in the droplet size with a very high span. This could be associated with more than one factor, most likely due to an increased coalescence probability at higher oil content immediately after droplet disruption [48] and/or due to an incomplete homogenization as a result of (i) an overall increase in emulsion viscosity and (ii) a decrease in the flow velocity of the emulsion in the dispersing zone [49,50]. This was also indicated by a higher span value of the respective emulsion samples. The large-sized droplets at high PPI concentrations could also be partly associated with the depletion flocculation [31].

Effect of PPI Concentration
In general, droplets' sizes reduce with higher PPI concentrations, and a similar behavior was observed in the current study. In fact, a limited particle concentration resulted in the insufficient interfacial coverage of droplets, which led to bridge flocculation and ultimately to coalescence [31,[42][43][44]. This was also evident from a higher span value (due to a wider droplet size distribution, or in some cases due to a bi-modal distribution) for those emulsion samples containing lower PPI concentrations. At extreme poor coverage of the interface, emulsification failure occurred [42], although this was not observed in any of the emulsification runs in the present study.
With the increase in PPI concentration, a gradual decrease in the droplets' size was observed. The decrease in droplet size with the increase in particle concentration has been reported previously [31,[45][46][47] and is valid until the maximum critical particle concentration is reached, beyond which, the increase in particle concentration contributes to the viscosity and density and subsequent emulsion stability [2,42,43]. For low oil content and high HPH pressure, the maximum critical particle concentration was~6%, i.e., a further increase in PPI concentration did not show a further decrease in droplet size. However, for high oil content and low HPH pressure, the increase in PPI (beyond 6%) resulted in a slight increase in the droplet size with a very high span. This could be associated with more than one factor, most likely due to an increased coalescence probability at higher oil content immediately after droplet disruption [48] and/or due to an incomplete homogenization as a result of (i) an overall increase in emulsion viscosity and (ii) a decrease in the flow velocity of the emulsion in the dispersing zone [49,50]. This was also indicated by a higher span value of the respective emulsion samples. The large-sized droplets at high PPI concentrations could also be partly associated with the depletion flocculation [31]. Hence, the desired droplets of <100 nm could only be achieved at a high PPI concentration (~6%), with less oil content using high HPH pressure.

Effect of Oil Content
The increase in dispersed phase in general results in an increase in the droplets' size [51,52], which was also evident in our results. Such an increase in droplets' size can occur due to an increase in the interfacial area, limiting the stabilizer with an increase in the oil content, as also explained in the above section. In addition, the increase in oil content is followed by a decrease in the aqueous continuous phase, which results in the decrease in the total stabilizer concentration of the system, especially when stabilizers are dispersed in the aqueous phase. The reduction in the effective stabilizer concentration to oil ratio results in large droplets with wide distribution and low stability [53][54][55]. In this study, as the PPI was present in the aqueous phase, the increase in oil content led to the decrease in the effective PPI concentration in the total emulsion system, limiting the particles, which could explain the increase in droplets' size with the increase in oil content. However, in some cases, the increase in dispersed phase can reduce the droplets' size, for instance, when the stabilizers were present in the dispersed phase [56] or when the stabilizer is in sufficiently high concentration in the continuous phase [51,57]. However, at low HPH pressure, the increase in oil content resulted in a slight reduction in the droplets' size. This could possibly be associated with droplet disruption due to an increased probability of liquid-liquid interactions (indirect breaking) at low flow velocities [58]. However, this can also be regarded as an incomplete homogenization which is also evident from a high D 50 and a very high span value as compared to those at high HPH pressure.

Effect of HPH Pressure
In general, the increase in HPH pressure increases the degree of emulsification and results in smaller-sized droplets [9,59]. In this study, the increase in HPH pressure also resulted in the decrease in the droplets' size, but only when the PPI was present in a sufficiently high concentration or when the oil content was low. High HPH pressure corresponds to a high energy density in the dispersing zone, which results in a smaller droplet size. However, the newly created interfacial area should be immediately stabilized; otherwise, the disruption results will be superimposed by coalescence [50]. This explains the effectiveness of high HPH pressure to produce nanoemulsions at high PPI concentrations or low oil fractions, with the former being more important. At a low PPI concentration, the increase in HPH pressure had a negligible effect on D 3,2 , while D 50 followed the general trend. This could be associated with the inadequate stabilization of droplets, resulting in a bimodal distribution, as evidenced by a high span value. Furthermore, provided the sufficient availability of stabilizer, a higher oil fraction (up to 7%) can also be well homogenized into nanoemulsion (while considering droplet size and span value together) even at a moderate HPH pressure (~35 MPa).
The reduction in droplets' size with increasing HPH pressure could partly be associated with protein modification (although it was not the aim of present investigation). The state of a protein significantly influences its emulsifying properties [60]. Several studies have revealed that the increase in HPH pressure can physically modify proteins to improve their emulsifying and stabilizing properties and thus can stabilize smaller-sized droplets at the same particle concentration [26,[61][62][63].

Prediction and Validation
The coefficient of determination for the median droplet diameter (D 50 ) and the span value (δ) were below 80% and were more likely to deviate away from the predicted values. Increasing the number of runs or further customization of the model could potentially improve the coefficient of determination. However, a good coefficient of determination and a highly significant response surface model was obtained for D 3,2 . Therefore, optimization was carried out to minimize D 3,2 below 100 nm. The optimization setting and optimized factor values for the nanoemulsions are presented in Table 2. The actual value of 80.53 ± 2.45 nm was within the upper range of the prediction value of 53.87 ± 34.23. Hence, the response surface model can efficiently be used to optimize the PPI-stabilized nanoemulsions and to predict D 3,2 within the experimental range of the study.

Emulsion Stability
Out of 20 emulsification runs shown in Table 1, 10 emulsion samples (having D 3,2 ≤ 100 nm) were selected for a short-term stability study to further explain the suitability of emulsification conditions. The selected samples were stored at 4 • C, and then, their droplet size distributions were again analyzed after 7 days. The change in droplet size between day 0 and day 7 is graphically shown in Figure 3a, where each data point represents the average D 3,2 of the selected emulsification runs. It can be visualized that for the majority of emulsification conditions, a non-significant change (p > 0.05) in average D 3,2 was observed between day 0 and day 7. However, the data points for day 7 seem to be more scattered than those for day 0. This is further explained by plotting D 3,2 (for Day 7) as a function of the PPI-to-Oil ratio (P/O) and HPH pressure used during the emulsification of respective samples (Figure 3b). A clear trend in the reduction in D 3,2 can be seen with the increase in HPH pressure; and furthermore, a high P/O resulted in the better stability of the nanoemulsions, which is again in accordance with the discussion presented in Section 3.1. In addition to HPH pressure and P/O ratio, another aspect of nanoemulsions' stability is their extremely small droplet size, which is evident from some data points (D 3,2 ≤ 100 nm), but the samples were still stable at a relatively lower P/O ratio. In fact, the nanosized emulsion droplets usually have better stability against coalescence [2,64]; specifically, they have a low collision probability at a low dispersed phase fraction [65,66]. Therefore, the droplet sizes for these samples did not grow significantly even after 7 days. emulsification conditions. The selected samples were stored at 4 °C, and then, their d let size distributions were again analyzed after 7 days. The change in droplet size betw day 0 and day 7 is graphically shown in Figure 3a, where each data point represent average D3,2 of the selected emulsification runs. It can be visualized that for the maj of emulsification conditions, a non-significant change (p > 0.05) in average D3,2 wa served between day 0 and day 7. However, the data points for day 7 seem to be m scattered than those for day 0. This is further explained by plotting D3,2 (for Day 7) function of the PPI-to-Oil ratio (P/O) and HPH pressure used during the emulsificati respective samples (Figure 3b). A clear trend in the reduction in D3,2 can be seen wit increase in HPH pressure; and furthermore, a high P/O resulted in the better stabili the nanoemulsions, which is again in accordance with the discussion presented in Se 3.1. In addition to HPH pressure and P/O ratio, another aspect of nanoemulsions' stab is their extremely small droplet size, which is evident from some data points (D3,2 ≤ nm), but the samples were still stable at a relatively lower P/O ratio. In fact, the nanos emulsion droplets usually have better stability against coalescence [2,64]; specifically, have a low collision probability at a low dispersed phase fraction [65,66]. Therefore droplet sizes for these samples did not grow significantly even after 7 days.

Effect of Salt on Optimized Conditions
For the optimized set of factor values, the size distribution characteristics of e sions at different concentrations of salt in the range of most food applications (0-400 were evaluated and are presented in Figure 4a-c. The droplets' sizes for all emuls were <100 nm. When a small amount (25 mM) of salt was present in the aqueous p the size of emulsion droplets was reduced, but a further increase in the salt concentra (200-400 mM) resulted in the formation of droplets with slightly bigger diameters. A tionally, the span value decreased initially with the increase in NaCl from 0 to 50 beyond which, a significant increase in span was observed. Very low salt concentra

Effect of Salt on Optimized Conditions
For the optimized set of factor values, the size distribution characteristics of emulsions at different concentrations of salt in the range of most food applications (0-400 mM) were evaluated and are presented in Figure 4a-c. The droplets' sizes for all emulsions were <100 nm. When a small amount (25 mM) of salt was present in the aqueous phase, the size of emulsion droplets was reduced, but a further increase in the salt concentration (200-400 mM) resulted in the formation of droplets with slightly bigger diameters. Additionally, the span value decreased initially with the increase in NaCl from 0 to 50 mM, beyond which, a significant increase in span was observed. Very low salt concentrations are in general favorable to emulsions, while a further increase in the salt concentration can reduce the electrostatic repulsion between the particles, resulting in some particle aggregation, which can increase the effective particle size at the interface and thus the size of emulsion droplets [2,67]. Despite the reduction in electrostatic repulsions, no phase separation or flocs were seen, indicating that the main force involved in the stabilization of emulsions using PPI was steric repulsion rather than electrostatic repulsion [46,68]. The slight increase in droplets' diameter at 200-400 mM NaCl could also be due to the coalescence of smaller droplets under the effect of strong electrostatic screening [46,69]. It is, however, important to note that the optimum salt concentrations for an emulsion preparation and their effect could depend on the type of salt in itself and the salt susceptibility of the proteins at the given pH [2,70,71]. ration or flocs were seen, indicating that the main force involved in the stabilization of emulsions using PPI was steric repulsion rather than electrostatic repulsion [46,68]. The slight increase in droplets' diameter at 200-400 mM NaCl could also be due to the coalescence of smaller droplets under the effect of strong electrostatic screening [46,69]. It is, however, important to note that the optimum salt concentrations for an emulsion preparation and their effect could depend on the type of salt in itself and the salt susceptibility of the proteins at the given pH [2,70,71].

Effect of pH on Optimized Emulsion
For the optimized set of factor values, the emulsion droplet size characteristics of emulsions (D3,2, D50, and δ) at different pH values are presented in Figure 4d-f. The emulsions prepared at pH 3 and 7 resulted in very fine droplets, whilst near the isolectric point (pI), i.e., pH 5, relatively larger droplets were formed. Near the pI, the protein particles lose their charge and electrostatic repulsion, and thus tend to form aggregates [2]. The aggregated particles will eventually form larger-sized droplets, as the size of droplets is significantly influenced by the size of particles. The smallest size of droplet that a particle can stabilize is its own size [2], although some reports are available in which the droplets smaller than the initial particle size were stabilized; however, it was suggested that the sizes of particles were reduced during the emulsification, and the droplets' sizes were larger than the sizes of effective particles at the interface [72,73].
The emulsions prepared at pH 3 and 7 were far smaller than those at pH 5. At pH 7, PPI have a high negative charge, low hydrophobicity, and high solubility, which allow

Effect of pH on Optimized Emulsion
For the optimized set of factor values, the emulsion droplet size characteristics of emulsions (D 3,2 , D 50 , and δ) at different pH values are presented in Figure 4d-f. The emulsions prepared at pH 3 and 7 resulted in very fine droplets, whilst near the isolectric point (pI), i.e., pH 5, relatively larger droplets were formed. Near the pI, the protein particles lose their charge and electrostatic repulsion, and thus tend to form aggregates [2]. The aggregated particles will eventually form larger-sized droplets, as the size of droplets is significantly influenced by the size of particles. The smallest size of droplet that a particle can stabilize is its own size [2], although some reports are available in which the droplets smaller than the initial particle size were stabilized; however, it was suggested that the sizes of particles were reduced during the emulsification, and the droplets' sizes were larger than the sizes of effective particles at the interface [72,73].
The emulsions prepared at pH 3 and 7 were far smaller than those at pH 5. At pH 7, PPI have a high negative charge, low hydrophobicity, and high solubility, which allow them to associate better with the oil-water interface to lower interfacial tension [74]. In contrast, at pH 3, PPI possess high positive charge, high hydrophobicity, and high solubility, which leads to the formation of stronger interfacial viscoelastic films [30,74]. While some studies have revealed that the emulsions formed at pH 3 formed smaller droplets [30,75,76] and much better interfacial rheology [74], the emulsions prepared under both of these conditions (pH 3 and 7) were statistically similar. The better results of PPI at pH 3 were suggested to be due to the ability of PPI to possess a smaller hydrodynamic diameter (134-165 nm) [30]. In this study, emulsions with droplets' size <100 nm were prepared, which implies that the effective size of PPI particles is much smaller than the suggested size, possibly due to the impact of high HPH pressure on PPI size and emulsion stabilization properties [26,[61][62][63].

Conclusions
Response surface methodology can be effectively used to optimize PPI-stabilized O/W ultra-nanoemulsions of Sauter mean diameter (D 3,2 ) < 100 nm. The droplets' size was mainly influenced by the PPI concentration and the interaction of oil content and HPH pressure. The increase in PPI concentration and HPH pressure reduced the emulsions' droplets' size, while the increase in oil fraction had a slight tendency to increase the emulsions' droplets' size. The span values of the emulsion droplets were high for a very low PPI concentration and for the combination of high oil fraction and low HPH pressure. The prepared emulsions were stable, as no significant difference in the emulsions size distribution was observed when these emulsions were evaluated on day 7. The desired emulsions (<100 nm) were optimized at 6.45% PPI, with 1.95% oil and at 46.82 HPH pressure. The emulsions prepared under such conditions were capable of maintaining the droplets' size <100 nm even when the ionic conditions were altered with up to 400 mM NaCl and the pH was adjusted to 3 or 7. However, micro-emulsions were formed at pH 5. The encapsulation and bioaccessibility studies of lipophilic micronutrients could further elaborate upon the potentials of PPI for use in the industrial preparation of such nanoemulsions.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.