Fabricating Pea Protein Micro-Gel-Stabilized Pickering Emulsion as Saturated Fat Replacement in Ice Cream

Unsaturated fat replacement should be used to reduce the use of saturated fat and trans fatty acids in the diet. In this study, pea protein micro-gels (PPMs) with different structures were prepared by microparticulation at pH 4.0–7.0 and named as PPM (pH 4.0), PPM (pH 4.5), PPM (pH 5.0), PPM (pH 5.5), PPM (pH 6.0), PPM (pH 6.5), and PPM (pH 7.0). Pea protein was used as a control to evaluate the structure and interfacial properties of PPMs by particle size distribution, Fourier transform infrared spectroscopy (FTIR), free sulfhydryl group content, and emulsifying property. PPM (pH 7.0) was suitable for application in O/W emulsion stabilization because of its proper particle size, more flexible structure, high emulsifying activity index (EAI) and emulsifying stability index (ESI). The Pickering emulsion stabilized by PPM (pH 7.0) had a uniform oil droplet distribution and similar rheological properties to cream, so it can be used as a saturated fat replacement in the manufacture of ice cream. Saturated fat was partially replaced at different levels of 0%, 20%, 40%, 60%, 80%, and 100%, which were respectively named as PR0, PR20, PR40, PR60, PR80, and PR100. The rheological properties, physicochemical indexes, and sensory properties of low-saturated fat ice cream show that PPM (pH 7.0)-stabilized emulsion can be used to substitute 60% cream to manufacture low-saturated fat ice cream that has high structural stability and similar melting properties, overrun, and sensory properties to PR0. The article shows that it is feasible to prepare low-saturated fat ice cream with PPM (pH 7.0)-stabilized Pickering emulsion, which can not only maintain the fatty acid profile of the corn oil used, but also possess a solid-like structure. Its application is of positive significance for the development of nutritious and healthy foods and the reduction of chronic disease incidence.


Introduction
Chronic diseases, such as cardiovascular disease, have become the most serious diseases threatening humans in the world today, seriously endangering people's physical health.[1].Many studies have shown that the intake of saturated fats and trans fatty acids increases the risk of chronic disease, whereas the intake of unsaturated fats decreases it [2,3].Therefore, in recent years, the nutritional guidelines by various government agencies have recommended that unsaturated fat replacements could be used to reduce the use of saturated fat and trans fatty acids in the diet, so as to achieve a proportion of solid fat in total daily calorie intake of less than 10% [4].
Ice cream is a common oil-in-water system food.It is very popular among consumers because of its delicate, smooth taste and rich, mellow flavor.At present, ice cream processing generally continues the previous traditional ice cream processing methods, with water, milk, cream, sucrose, etc. as the main raw materials.It is characterized by high sugar, calories, and saturated fatty acid.Although consumers are already aware that excessive intake of saturated fat is not conducive to human health, ice cream requires saturated fat to provide structure and a creamy mouthfeel [5].Therefore, researchers have begun to investigate different techniques for reducing saturated fat to ensure that healthier, nutritious foods are obtained without degrading or destroying the structural and organoleptic properties of the food [6].
Ice cream is a favorite casual dessert, and the market for low-saturated fat ice cream is growing rapidly as consumers focus on healthy eating.However, simply reducing the saturated fat content of ice cream can affect many properties such as viscosity, hardness, melting, and flavor.The difficulty in low-saturated fat ice cream is how to reduce the saturated fat content of ice cream without affecting its sensory properties [7].Currently, the main approach to reduce saturated fat in food is to use saturated fat replacements (saturated fat mimics and saturated fat replacers like protein-based fat replacers, carbohydrate-based fat replacers, lipid-based fat replacers, and complex fat replacers) to substitute some or all of the saturated fat to ensure sensory properties while reducing the calories [8].However, these saturated fat replacements have not been commercialized due to a number of problems.The ester bonds of saturated fat replacements cannot be hydrolyzed by lipases and therefore cannot be digested and absorbed by the body [9].Thus, high intake of saturated fat replacements can lead to many intestinal problems.The application of protein matrix and carbohydrate matrix fat mimics would lead to loss of fat-soluble flavor substances, shortening the shelf life [10,11].Therefore, the development of a novel saturated fat replacement is of positive significance for the development of nutritious and healthy foods.
In recent years, consumers have gradually realized that saturated fatty acid intake increases the incidence of chronic diseases, so the demand for low-saturated fat food increases year-by-year.Protein-stabilized O/W emulsions, including pea protein-, soybean protein-, wheat gliadin-and corn gliadin-stabilized emulsions, have attracted more and more attention as new saturated fat replacements [12].On the one hand, protein-stabilized emulsions enable the solidification of vegetable oils.On the other hand, vegetable oil contained in plant protein-stabilized emulsions contains more unsaturated fat esters than fat, which is more nutritious [13].However, protein-stabilized conventional emulsions are sensitive to environmental factors, such as acid and heat, and tend to destabilize during processing and storage [14].Therefore, solving this destabilization problem by modifying the molecular structure of proteins is of increasing interest to the food industry.
The protein micro-gel particles produced by microparticulation require heating of the protein to form a gel structure while/later applying high shear to disrupt the gel network into particles [15].Protein micro-gel particles tend to have great interfacial properties because of the increasing hydrophobicity of the protein surface during microparticulation [16].Pickering emulsion stabilized by protein micro-gel particles has advantages such as good stability and biocompatibility [17].Thus, we hypothesized that a protein micro-gel particle-stabilized Pickering emulsion could be used to replace saturated fat in ice cream.
Pea protein has antioxidant, anti-hypertension, anti-inflammatory, cholesterol reduction, and regulation of intestinal bacterial activity, and it has hypo-sensitization [18].However, the poor solubility, emulsification, and emulsification stability of pea protein limit its application in food processing [19].Thus, in this study, pea protein was microparticulated to improve its processing performance, such as emulsification.Different structural and physicochemical properties of pea protein micro-gel particles were obtained by adjusting the preparation pH during microparticulation.Then, an O/W emulsion stabilized by pea protein micro-gel particles was designed to substitute saturated fat and was used as saturated fat replacement in low-saturated fat ice cream.And the saturated fat in ice cream was partially replaced at different rates of 20%, 40%, 60%, 80%, and 100%.The physicochemical properties and rheological properties of full-saturated fat and low-saturated fat ice cream were investigated to find the most suitable replacement ratio.

Preparation of Pea Protein Micro-Gel Particles
PPI (18 g/100 mL) dispersion was magnetically stirred for 2 h, and the initial pH of PPI dispersion was 7.22.Then, the pH of PPI dispersion was adjusted from 4.0 to 7.0 using 1 M HCl solution.Then, the PPI dispersions at pH 4.0 to 7.0 were respectively heated at 90 • C for 30 min.The PPI gels were prepared by cooling the heated PPI dispersions to room temperature and then refrigerating at 4 • C for 12 h.Then, the PPI gels were homogenized in a T18 Ultra-Turrax (IKA, Staufen, Germany) at 8000 rpm for 3 min to obtain PPI micro-gel particles.They were named pea protein micro-gel (PPM) (pH 4.0), PPM (pH 5.5), PPM (pH 5.0), PPM (pH 6.0), PPM (pH 6.5), and PPM (pH 7.0), respectively.

Particle Size Distribution of PPMs
Particle size distribution of the protein micro-gel was measured by a Battersize2000 particle size analyzer (Better, Dandong, China).The samples were added to a stirred measuring cell containing 800 mL of deionized water at room temperature, with deionized water as the dispersing medium.And the light shading percentage was between 10% and 15% in order to avoid the multiple scattering effects of proteins and to ensure the accuracy of the measurement.The relative refractive indices of the micro-gel particles and the dispersing medium (deionized water) were 1.52 and 1.33, respectively.Each sample with different treatments was measured three times and the average value was taken.

Free Sulfhydryl Determination of PPI and PPMs
Free sulfhydryl (SH F ) was measured according to Yang et al. [20], with minor modifications.Each sample was diluted to 10 mg/mL and then mixed with 3 mL of Tris-glycine buffer (0.086 M Tris-glycine, 0.004 M EDTA, pH 8.0) and 0.03 mL of Ellman's reagent (4.0 mg/mL DTNB) in the same buffer.After fully dissolving, the mixtures were centrifuged at 8000 rpm for 10 min.A Tris-glycine buffer and Ellman reagent were added into the supernatant and mixed evenly.A412 was determined by UV spectrophotometer (UV-2500, Shimadzu, Shanghai, China) after 25 min.In the experiment, the solution with the Ellman reagent without protein was used as the control.SH F content was calculated according to Equation (1), written as where D is dilution ratio and C is the protein concentration (mg/mL).

Fourier Transform Infrared Spectroscopy of PPI and PPMs
The wavenumber range of the iS10 FTIR spectrometer (Thermo Nicolet, Madison, WI, USA) of Nicole Power was 4000-400 cm −1 .The resolution of the spectrometer was 4 cm −1 .The signal to noise ratio was 50,000:1, and the scanning count was 32 times [21].

Emulsifying Properties of PPI and PPMs
Emulsifying activity index (EAI) and emulsion stability index (ESI) were measured according to Figueroa-González et al. [22], with minor modifications.A total of 2 mL corn oil and 8 mL protein solution were mixed.Then, the mixture was homogenized by using a T18 Ultra-Turrax (IKA, Staufen, Germany) at 10,000 rpm for 3 min.Immediately, 50 µL of the emulsion was diluted with 2.45 mL 0.1% SDS solution.The A0 and A10 of the diluted emulsion after 0 min and 10 min of resting at 500 nm were recorded.Then, the EAI and ESI values were calculated using Equations ( 2) and (3), written as where φ is the volume fraction of corn oil (v/v) (φ = 0.20), DF denotes dilution ratio (200), and C is the concentration of emulsifier (g/mL).

Emulsion Preparation
A 60 wt% oil phase and 40 wt% PPI or PPM (pH 7.0) were mixed to prepare the Pickering emulsion.The mixture was sheared at 10,000 rpm for 5 min by using the T18 Ultra-Turrax (IKA, Staufen, Germany).

Microstructure Observation
The microstructure of the emulsion droplets was observed using LSM800 CLSM (Zeiss, Bartenwueburg, Germany).A total of 20 µL of 0.1% (w/v) Nile red solution was used as a fluorescent dye of oil in the 1 mL emulsion, with excitation and emission wavelengths at 633 nm and 488 nm, respectively [23].

Rheological Properties Dynamic Oscillatory Measurements
The emulsion was characterized by a MARS 60 rheometer (Haake, Karlsruhe, Germany) at 25 • C [24].The storage modulus G ′ (Pa), loss modulus G ′′ (Pa), and tan δ of the emulsions, with frequency ranging from 0.01 to 10 Hz, were measured using a clamp with a conical (25 mm) diameter and a distance of 1 mm between parallel plates.G ′ , G ′′ , and the angular frequency ω were fitted using Equations ( 4) and (5), written as where K ′ (Pa•s n ′ ) and K ′′ (Pa•s n ′′ ) are constants and n ′ and n ′′ may be referred to as the frequency exponents.

Steady-State Shear Properties
The steady-state shear properties of the PPM (pH 7.0)-stabilized emulsion were measured by the MARS 60 rheometer (Haake, Karlsruhe, Germany) at 25 • C. Shear rate was set to increase from 1 s −1 to 300 s −1 in 400 s [25].Shear rate (γ, s −1 ) and apparent viscosity (η, Pa•s) were fitted by a power-law model (Equation ( 6)).Then, the consistency coefficient K (Pa•s n ) and flow behavior index n were obtained [10].

In Vitro Digestion In Vitro Digestion Model
The in vitro digestion model, including three digestion stages of saliva, stomach, and intestine, was created [26,27].The formulations of the stimulated saliva fluid (SSF), stimulated gastric fluid (SGF), and stimulated intestinal fluid (SIF) are shown in Table S1.
In general, for the oral phase, 5 g of the emulsions were mixed with 4 mL SSF electrolyte stock solution (contained 3.75 mg/mL mucin).Then, 0.975 mL of deionized water and 25 µL CaCl 2 solution (0.3 mol/L) were added to the mixture.After the pH was adjusted to 7.0, the mixtures were then stirred at 200 rpm for 2 min (37 • C).
For the gastric phase, 7.5 mL SGF, 1.6 mL porcine pepsin solution (25,000 U/mL), 5 µL of 0.3 mol/L CaCl 2 solution, and 148 µL of 6 mol/L HCl solution were added into 10 mL of the above mixture while stirring at 300 rpm for 2 h at 37 • C.
For the intestinal phase, after the pH was adjusted to 7.0 using a 1 mol/L NaOH solution, the gastric chyme for each sample was mixed with a 18.5 mL SIF and 40 µL of 0.3 mol/L CaCl 2 solution.The mixtures were stirred at 200 rpm for 2 h at 37 • C.

Particle Size Distribution and Microstructure of Emulsions during In Vitro Digestion
Particle size distribution and microscopic images of PPI-stabilized emulsion and PPMstabilized Pickering emulsion during in vitro digestion were measured according to the methods in Sections 2.3.1 and 2.5.2, respectively.

Free Fatty Acids (FFA) Release of Emulsions
During simulated intestinal digestion, lipids were continuously hydrolyzed to FFA by pancreatic lipase.The amount of FFAs released from the system was calculated based on the volume of the NaOH solution (0.1 mol/L) consumed in 2 h according to Equation ( 7), written as where M Lipid is the molar mass of corn oil (872 g/mol) and m Lipid is the weight of corn oil in the emulsions (g).

Rheological Properties
Rheological properties of ice cream slurries were measured according to the method of Section 2.5.3.

Melting Measurement
After extrusion of ice cream, about 30 g of ice cream (the mass was recorded as M 0 ) was weighed and placed on a metal net with a sieve size of 3 mm × 3 mm at 25 The dropping mass of ice cream was weighed every 10 min until it reached 50 min.Melting and melting rate were calculated according to Equations ( 8) and ( 9), written as Melting rate (g/min) = Mass of dropping Melting time (9) where M 50 is the mass of ice cream after 50 min of dropping.

Overrun Measurement
The overrun of ice cream was measured by recording per volume weight of ice cream slurry (m) and per volume weight of ice cream (m f ) according to Equation (10) [29].

.5. Sensory Evaluation
The sensory evaluation methods were based on the methods of Roberta et al. [30], with some modifications.According to the Chinese standard [31] (criterion for sensory evaluation of frozen drinks and ice cream), 20 reviewers were trained for 20 days to form the review panel.Samples were divided as separate portions in 30 mL plastic cups and were assessed in duplicate.The samples were coded, and their order randomized.Ice cream samples were prepared and stored at −18 • C. Before evaluation, the samples were left at room temperature for 5 min.The samples were evaluated by the reviewers in turn for the following five main aspects: color, morphology, organization, taste and odor, and impurity.Distilled water was required to eliminate the effects between samples after different flavors.A 100-point intensity scale (1-10 points for color, 6-30 points for morphology, 5-30 points for organization, 5-20 points for taste and odor, and 3-10 points for impurity) was used for each term, as detailed in Table 2. Sensory acceptance was obtained by adding up the scores of the color, morphology, organization, taste and odor, and impurity.

Statistical Analysis
All the experiments were repeated at least three times.All data were assed with one-way ANOVA followed with Tukey's multiple-range test.Differences were considered statistically significant at p < 0.05.The results were resented as mean ± standard deviations.

Particle Size Distribution of PPMs
The mean particle sizes of the PPMs prepared at different pH ranging from 4.0 to 7.0 are shown in Figure 1A.It shows that the pH has a significant effect on the particle size distribution of the prepared PPMs.During microparticulation, pea protein was thermally denatured to form a gel structure.After high-speed homogenization, the gel was crushed into a three-dimensional mesh structure with a diameter of 13.65-56.78µm.The mean particle size of PPM (pH 6.0) was the maximum and was 56.78 µm.A decline of mean particle size was observed when the microparticulation pH was far from pH 6.0.Mean particle sizes of PPM (pH 4.0), PPM (pH 4.5), and PPM (pH 5.0) were 20.69 µm, 18.79 µm, and 20.80 µm, respectively.There was no significant difference among the mean particle sizes of PPM (pH 4.0), PPM (pH 4.5), and PPM (pH 5.0).The sample with the smallest particle size was PPM (pH 7.0), with a mean particle size of 13.65 µm.

Statistical Analysis
All the experiments were repeated at least three times.All data were assed with oneway ANOVA followed with Tukey's multiple-range test.Differences were considered statistically significant at p < 0.05.The results were resented as mean ± standard deviations.

Particle Size Distribution of PPMs
The mean particle sizes of the PPMs prepared at different pH ranging from 4.0 to 7.0 are shown in Figure 1A.It shows that the pH has a significant effect on the particle size distribution of the prepared PPMs.During microparticulation, pea protein was thermally denatured to form a gel structure.After high-speed homogenization, the gel was crushed into a three-dimensional mesh structure with a diameter of 13.65-56.78µm.The mean particle size of PPM (pH 6.0) was the maximum and was 56.78 µm.A decline of mean particle size was observed when the microparticulation pH was far from pH 6.0.Mean particle sizes of PPM (pH 4.0), PPM (pH 4.5), and PPM (pH 5.0) were 20.69 µm, 18.79 µm, and 20.80 µm, respectively.There was no significant difference among the mean particle sizes of PPM (pH 4.0), PPM (pH 4.5), and PPM (pH 5.0).The sample with the smallest particle size was PPM (pH 7.0), with a mean particle size of 13.65 µm.Guo et al. [32], Cui et al. [33], and Schmitt et al. [34] confirmed that the surface charge is close to zero at the isoelectric point (pI).Thus, protein molecules aggregate and precipitate under hydrophobic interactions.However, when the pH is far away from the pI, the surface charge of protein molecules increases, which generates intermolecular repulsion and inhibits aggregation, leading to a decrease in protein particle size.

Content of Free Sulfhydryl Groups
Free sulfhydryl groups (SH F ) reflect the changes of the protein tertiary structure and played key role in the protein rigid structure formation.Free sulfhydryl groups can form disulfide bonds after oxidation, and disulfide bonds can be further reduced to free sulfhydryl groups through hydrogenation under the conditions of mutual conversion during heat treatment [35] and shear treatment.
The free sulfhydryl group content of different PPMs is shown in Figure 1B.SH F content of PPI is 5.45 µmol/g.The SH F content of PPM (pH 4.0), PPM (pH 4.5), PPM (pH 5.0), PPM (pH 5.5), PPM (pH 6.0), and PPM (pH 6.5) were 0.07, 0.09, 0.11, 0.19, 0.64, and 4.75 µmol/g, which were significantly lower than PPI; it was due to the highly concentrated, accessible, and reactive free sulfhydryl groups that were oxidized to form the disulfide bond [36].
The SH F content of PPMs gradually increased with the increasing microparticulation pH.The SH F content of PPM (pH 7.0) was the highest (11.10 µmol/g) and significantly higher than PPI.The increasing content of the free sulfhydryl groups may be due to the unfolding of the protein molecular structure and the break of disulfide bonds caused by high-temperature heating.Mession et al. [37] confirmed that heating causes the unfolding of protein molecules and exposing of embedded free sulfhydryl groups, disulfide bonds, and hydrophobic groups.The change of the functional group distribution can enhance the intermolecular interaction and then promote the formation of the gel network [38].
Guo et al. [10] found that the SH F content of protein molecules increases the promotion of the interaction between protein molecules and forms a tight and uniform network of gels.In summary, pH and heat can change the content of sulfhydryl bonds and affect the functional properties of proteins.

Fourier Transform Infrared Spectroscopy Analysis of PPI and PPMs
Figure 1C displays the spectra of PPI and PPMs in the 4000 to 400 cm −1 range.The OMNIC software was used to fit the peaks in the range of 1700 to 1600 cm −1 to calculate the content of each secondary structure of proteins [39].The peaks located at 1650-1658 cm −1 reflect the α-helix, 1610-1640 cm −1 reflect the β-sheet, 1660-1700 cm −1 reflect the β-turn, and the peaks at 1640-1650 cm −1 reflect the random coil [40]; the results are shown in Table 3.As can be seen in Table 3, the PPI contains 11.66% α-helix, 33.04% β-sheet, 46.94% β-turn, and 8.36% random coil.The β-sheet and β-turn are the dominant secondary structures of PPI and PPMs, followed by random coil and α-helix.Compared with PPI, the proportion of α-helix significantly decreased, which is consistent with the research of Beck et al. [41].The β-sheet content in PPM (pH 4.0) and PPM (pH 4.5) significantly increased to 41.64% and 39.73%, respectively; the β-turn content had significantly decreased to 39.53% and 41.76%, respectively; random coil content showed no significant difference with PPI.The decreasing α-helix contents and increasing random coil content and β-sheet content were observed in PPM (pH 5.0) and PPM (pH 6.0).The secondary structures of PPM (pH 5.5) and PPM (pH 6.5) shifted from β-sheet and α-helix to β-turn and random coil, which is consistent with the results of Sun et al. [16].It can be concluded that the reduction in α-helical structures might have been due to the fact that microparticulation destroyed the rigid structure of protein molecules, causing unfolding and rearrangement.The rearranged protein molecules tend to have a higher degree of intermolecular flexibility [39].There is no significant difference in the content of each secondary structure between PPM (pH 7.0) and PPI.Therefore, it can be concluded that the different pH during macroparticulation influenced the degree of protein unfolding, the distribution of amino acid side chains, and changes in protein surface charges and [42].Thus, flexible structures in the protein microgel particle formed after the disruption of the natural structures during microparticulation could lead to differences in their physiochemical properties with PPI, such as interfacial properties [16].The values that do not bear the same letter in the same column are significantly different (p < 0.05).

Emulsifying Properties
Figure 2A shows that the EAI of PPI increased gradually with the increase of pH from 5.5 to 7.0.Different pH preparation of PPMs also had a significant effect on EAI (Figure 2A).The EAI at pH near the pI (pH 4.0-6.0)was lower than the EAI at pH 6.0-7.0.The EAI values of the protein micro-gel particles prepared at pH 6.0-7.0 were higher than that of pI (pH 4.0-6.0),which was mainly due to the higher degree of protein molecules unfolding.The unfolded protein molecules could rapidly adsorb on the oil-water interface to form a denser viscoelastic interfacial layer, which stabilized the dispersed oil droplets.At neutral pH far from the isoelectric point, pea protein forms a stronger and denser viscoelastic network when adsorbed at the oil-water interface than that at pI (pH 4.0-6.0)[43].Near the isoelectric point, protein molecules have no surface charge, so low intermolecular repulsion leads to molecular aggregation, which caused the lower EAI of proteins [44].Another reason could be the proper three-phase contact angle of the micro-gel particle, which is attributed to balanced distribution of hydrophilic and hydrophobic groups of PPM prepared at pH 6.0-7.0 [33].The pH values of the dispersions have significant impacts on the ESI of PPI (Figure 2B).Among them, PPI has the better ESI, at pH 4.0 and pH 7.0, than others, with values of 89.10% and 92.40%, respectively.Except for PPM (pH 4.0), different preparation pH had no significant effect on the ESI values of other PPM samples.It can be seen that protein micro-gel particles once adsorbed on the oil-water interface, form a stable interfacial layer to stabilize the oil droplets [45].
The EAI and ESI of PPMs were both significantly higher than that of the PPI at corresponding pH (Figure 2), indicating that micro-gel particles have higher potential applications in multiphase systems than PPI.Natural pea proteins are globular proteins with well-defined molecular conformations, where hydrophilic groups are distributed on the surface of the protein molecule, while the hydrophobic groups are embedded within the molecule.During microparticulation, protein molecules partially unfolded, then recombined and cross-linked to form the gel network structure [36].Subsequently, high-speed shear treatment disrupts the gel maintained by non-covalent bonds to obtain micro-gel particles.Simultaneously, the surface hydrophilic and hydrophobic groups of protein molecules undergo rearrangement [45].In summary, compared to the rigid structure of PPI, more flexible structures were formed during microparticulation, which can be better adsorbed on the oil-water interface to form an interfacial layer that stabilizes the oil droplets.The higher EAI (56.56 m 2 /g) and ESI (95.10%) of PPM (pH 7.0) than other samples (Figure 2) indicated that PPM (pH 7.0) is more suitable for application in O/W emulsion stabilization.Therefore, PPM (pH 7.0) was chosen for subsequent emulsion preparation and research.

Micromorphology Observation
After staining with Nile red, the oil droplets showed red after excitation at 636 nm.It can be seen from Figure 3A that oil droplets exist in the emulsion in a round shape with uniform size.There is an interfacial protein layer, which can not only effectively reduce the surface tension, but also help to form uniform and stable emulsion [46].This also explains well the better ESI of PPM (pH 7.0).And it can be found that the droplets size of the emulsion was less than 10 µm, which indicated that the PPM (pH 7.0)-stabilized emulsion can simulate the smooth and delicate texture of saturated fat and can be used as a saturated fat replacement to replace thin cream to prepare low-saturated fat ice cream.

Rheological Properties
Dynamic rheological measurements can reflect the formation and destruction of composite gel matrix [24].The viscoelasticity of the PPM (pH 7.0)-stabilized Pickering emulsion can be characterized by storage modulus (G ′ ) and loss modulus (G ′′ ). Figure 3B shows that G ′ and G ′′ of the PPM (pH 7.0)-stabilized Pickering emulsion slightly increased with the increasing frequency.The slopes of G ′ (K ′ = 4753.64Pa•s n ′ ) were higher than that of G ′′ (K ′′ = 768.06Pa•s n ′′ ), which indicated that G ′ was more frequency dependent than G ′′ [10].G ′ of the PPM (pH 7.0)-stabilized Pickering emulsion was always larger than G ′′ and tan δ was less than one throughout the tested frequency range, indicating that elastic behavior dominated the rheological properties [47].
Figure 3C shows that the viscosity of the PPM (pH 7.0)-stabilized Pickering emulsion decreased with the increase in shear rate, indicating that the emulsion is a shear-thinning fluid.In PPM (pH 7.0)-stabilized Pickering emulsion, the droplets are close enough together to interact with each other which may lead to the formation of a three-dimensional network of aggregated droplets.As the shear rate is increased, the hydrodynamic forces cause aggregates to become deformed and eventually disrupted which results in a reduction in the viscosity.The flow behavior index n is 0.81, which is less than 1.It indicated that the PPM (pH 7.0)-stabilized Pickering emulsion is pseudoplastic fluid.According to Long et al. [48], cream is a pseudoplastic fluid and exhibits shear-thinning behavior.PPM (pH 7.0)-stabilized Pickering emulsion showed the same rheology properties as cream, which ensured its potential application as a saturated fat replacement for cream.

In Vitro Digestion
Figure 4A,B shows the particle size distribution and microstructure of PPI-stabilized emulsion and PPM (pH 7.0)-stabilized Pickering emulsion during simulated in vitro digestion.Figure 4A 1 ,B 1 shows that before digestion, the oil droplets in both emulsions exhibited uniformly distribution in the form of round balls.The particle sizes d 32 in the two emulsions are not significantly different.However, there is a very small distribution peak at 10-100 µm in the particle size distribution curve of the PPI-stabilized emulsion, indicating the presence of a few large oil droplets, which leads to a significantly higher d 43 (8.52After simulated gastric digestion, the particle size distribution and microscopic morphology of both emulsions changed dramatically (Figure 4A 3 ,B 3 ).The oil droplet size increased significantly (p < 0.05), which was attributed to the coalescence of the oil droplets.The gastric acid (pH 2.0) and pepsin in the SGF caused the protein at the oil-water interface in the emulsions to be denatured, losing the role of forming interfacial film to protect the oil droplets.It further led to flocculation and even coalescence of the oil droplets [49].Compared to PPM (pH 7.0)-stabilized Pickering emulsions, PPI-stabilized emulsions showed more flocculation of oil droplets after simulated gastric digestion (Figure 4A 3 ,B 3 ).This suggests that the micro-gel particles formed by pea proteins have an enhanced ability to resist both the acidic environment and pepsin digestion [50].
After simulated intestinal digestion, the number of oil drops in the two emulsions decreased significantly under the dilution of digestive fluid and the digestion of lipasecolipase (Figure 4A 4 ,B 4 ); and the oil droplet size significantly increased (p < 0.05).However, the oil droplet size in PPM (pH 7.0)-stabilized Pickering emulsion was significantly smaller than that in PPI-stabilized emulsion (p < 0.05).And the number of oil droplets was significantly greater than that in PPI stabilized emulsion (p < 0.05).It indicated that PPM (pH 7.0) has the ability to protect the integrity of the emulsion structure in the simulated gastrointestinal fluid and resist the digestion of simulated gastrointestinal fluid.This can also be confirmed from Figure 4C; it shows that the FFA release rate of both emulsions increases rapidly within 45 min of simulated intestinal digestion.After 45 min, the FFA release rate slowed down and eventually approached 0. Obviously, the FFA release rate of PPM (pH 7.0)-stabilized Pickering emulsion was significantly less than that of the PPIstabilized emulsion (p < 0.05).This indicates that the PPM (pH 7.0) on the oil-water interface has the ability to block the adsorption of bile salts and lipase, thus the delay of the lipid digestion of emulsions.The results of Sarkar et al. [51] also highlighted that the protein gel is able to slow the diffusion of both lipase-colipase and proteases, thus inhibiting lipid digestion of O/W emulsion.Therefore, PPM (pH 7.0)-stabilized emulsion can be used as compound fat replacement to reduce the lipid content of food and decline the lipid bioavailability.

Steady-State Shear Properties of Ice Cream Slurry
As can be seen from Figure 5, the viscosity of all ice cream slurries show a decreasing trend with the increasing shear rate and show a nonlinear relationship, indicating that all ice cream slurries exhibited shear-thinning behavior.This is due to the fact that at rest or low flow rates, the macromolecular chain-like structures are entangled with each other, and the viscosity of all ice cream slurries is high.With the increasing shear rate, due to the effect of the shear stress, the disordered macromolecular chain rolled and rotated and contracted into clusters, which reduced the intermolecular cross-linking, and thus the viscosity gradually decreased [52].According to the power law model with high regression coefficients (R 2 = 0.98-0.99), the fitting parameters were used to obtain the flow behavior index (n) and consistency coefficient (K) of all ice cream slurries.As can be seen in Figure 5, there were no significant differences in the n values of the ice cream slurries as the substitution rate increased from 0% to 100%.The index n is greater than zero and less than one, indicating that all ice cream slurries are pseudoplastic fluids.K values significantly increased with the increasing substitution rate, which indicated that the ice cream slurries showed more viscous properties due to the increase in apparent viscosity.Wang et al. [53] also found that the K values of ice cream slurries increased with the increase of soybean oil body substitution rate.They claimed that it was due to the fact that the increased content of Pickering emulsion introduces PPM biomolecules, whose hydrogen bonding or hydrophobic interactions between the nonpolar side segments of the carbon skeleton produce new ordered and improved structures, which lead to differences in rheological behavior.

Dynamic Oscillatory Measurements of Ice Cream Slurry
The dynamic oscillatory of ice cream slurries can reflect the physicochemical properties and flow behavior during chewing in the mouth.Storage modulus G ′ and loss modulus G ′′ were analyzed, shown in Figure 6A,B.G ′ and G ′′ of all ice cream slurries, except PR40 and PR60, increased with increasing frequency until 1.3 Hz, which indicated the destruction of the structures.G ′ and G ′′ of PR40 and PR60 remained stable until 2.8 Hz, indicating that the structure was still intact.This also indicates the enhanced structural stability of low-fat ice cream slurries at 40% and 60% substitution rates.The tan δ (G ′′ /G ′ ) of ice cream slurries is a key indicator of the phase melting and structural changes [54], which are shown in Figure 6C.The tan δ of all the samples in the measured frequency range are less than one, indicating that the G ′ values are higher than the G ′′ values, which suggests that all the ice cream slurries dominate the liquid characteristics, which are also known as viscoelastic fluids; viscoelastic fluid characteristics are the desired rheological properties of ice cream.It can be seen that the application of PPM (pH 7.0)-stabilized Pickering emulsion replacement of saturated fat in ice cream did not change its solid/liquid behavior.This also indicates PPM (pH 7.0)-stabilized Pickering emulsion has the potential to be applied in low-fat foods as compound fat replacements.

Physicochemical Properties and Sensory Evaluation of Ice Cream
Table 4 shows the effects of the substitution rates of the PPM (pH 7.0)-stabilized Pickering emulsion on the physicochemical properties and sensory evaluation of the ice cream.There was no significant difference in the melting of each ice cream as the substitution rate increased to 80%.This indicates that the substitution of saturated fat by a moderate amount of PPM (pH 7.0)-stabilized Pickering emulsion does not affect the ice phase organization.However, the melting of PR100 (1.345%) was significantly greater than the other samples, due to the fact that the ice phase arrangement of PR100 was different from the other samples.In the PR100 with a lower overrun and stability, there was not enough network space structure to prevent melting and complete collapse.According to Kaleda et al. [55], the difference in melting between the two types of ice cream can only be attributed to differences in ice phase arrangement.The values do not bear the same letter in the same column are significantly different (p < 0.05).
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Figure 1 .Figure 1 .
Figure 1.Particle size distribution (A), free sulfhydryl content (B) and FTIR spectrum (C) of pea protein micro−gels (PPMs) and pea protein isolate (PPI).The values that do not bear the same letter on the column are significantly different (p < 0.05).

Figure 2 .
Figure 2. Emulsifying activity index (A) and emulsion stability index (B) of PPI at different pH and PPMs.The values that do not bear the same letter on the column are significantly different (p < 0.05).

Figure 2 .
Figure 2. Emulsifying activity index (A) and emulsion stability index (B) of PPI at different pH and PPMs.The values that do not bear the same letter on the column are significantly different (p < 0.05).

Figure 4 .
µm) than that of the PPM (pH 7.0)-stabilized Pickering emulsion (4.45 µm).Comparing Figure4A2 ,B 2 and Figure4A 1,B 1 , it can be found that after a short period of oral digestion, the particle size distribution and microscopic morphology of the two emulsions did not change significantly (p > 0.05).This can also be confirmed by the non-significantly different values of d 32 and d 43 .dddCont.

Figure 5 .
Figure 5. Flow curves of ice cream slurries measured at 25 °C.a−f : The values that do not bear the same letter on the column are significantly different (p < 0.05).

Figure 5 .
Figure 5. Flow curves of ice cream slurries measured at 25 • C. a−f : The values that do not bear the same letter on the column are significantly different (p < 0.05).

Table 1 .
Composition of ice cream at different saturated fat replacer ratios.

Table 2 .
Sensory evaluation scoring index and criteria of ice cream.

Table 3 .
Content of secondary structures of PPI and PPMs.