# Impact of Flaxseed Gums on the Colloidal Changes and In Vitro Digestibility of Milk Proteins

Thierry Hellebois, Jennyfer Fortuin, Claire Gaiani, Christos Soukoulis

► To cite this version:

Thierry Hellebois, Jennyfer Fortuin, Claire Gaiani, Christos Soukoulis. Impact of Flaxseed Gums on the Colloidal Changes and In Vitro Digestibility of Milk Proteins. Foods, 2022, 11 (24), pp.4096. 10.3390/foods11244096. [hal-04117162](https://hal-univ-lorraine.fr/hal-04117162)

HAL Id: hal-04117162

https://hal.univ-lorraine.fr/hal-04117162

Submitted on 5 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Impact of Flaxseed Gums on the Colloidal Changes and In Vitro Digestibility of Milk Proteins

Thierry Hellebois <fnref n="1" />,<fnref n="2" /> ID, Jennyfer Fortuin <fnref n="1" />,<fnref n="3" />, Claire Gaiani <fnref n="2" /> ID and Christos Soukoulis <fnref n="1" />,*

**Abstract:** Flaxseed (*Linum usitatissimum* L.) mucilage is one of the most studied plant seed gums in terms of its techno-functional and health-promoting properties. Nonetheless, the interplay of flaxseed gum (FG) with other food biopolymers, such as milk proteins, under in vitro digestion conditions remains underexplored. The aim of the present work was to investigate the colloidal interplay between flaxseed gum (golden or brown) and milk proteins (sodium caseinate or whey protein isolate) under simulated in vitro digestion conditions and its relationship with the attained in vitro protein digestibility. The presence of flaxseed gum in the milk protein food models and in the oral food boluses obtained was associated with the occurrence of segregative microphase separation. Flaxseed gum exhibited a prominent role in controlling the acid-mediated protein aggregation phenomena, particularly in the sodium caseinate gastric chymes. The addition of FG in the food models was associated with a higher amount of intact total caseins and β-lactoglobulin at the end of the gastric processing step. Monitoring of the intestinal processing step revealed a very advanced cleavage of the whey proteins (>98%) and caseins (>90%). The degree of the milk protein hydrolysis achieved at the end of the intestinal processing was significantly higher in the systems containing flaxseed gum (i.e., 59–62%) than their gum-free protein counterparts (i.e., 46–47%). It was postulated that the electrostatic milk protein complexation capacity and, to a lesser extent, the thickening effect of flaxseed gum influenced the in vitro digestibility of the milk proteins.

**Keywords:** *Linum usitatissimum* L.; mucilage; in vitro digestion; whey protein; sodium caseinate; protein–polysaccharide interactions

## 1. Introduction

Flax (*Linum usitatissimum* L.) is one of the most important industrial crops cultivated for its fibres and seeds. Although flax fibres are prevalently exploited in the manufacturing and composite biomaterials domain [1], flaxseed is considered a rapidly emerging food biomass in the production of sustainable dietary relevant ingredients such as soluble dietary fibres (mucilage), proteins, and lipids [2]. Its mucilaginous matter, located in the outermost seed layer, accounts for 3–9% of the weight of the whole seed. From an osidic composition viewpoint, flaxseed gum (FG) comprised three major polysaccharidic fractions: a low molecular weight rhamnogalacturonan-I (RG-I) rich ($M_w$ ~300–700 kDa), a high molecular weight arabinoxylan (AX) rich ($M_w$ ~2000–4000 kDa), and a medium molecular weight ($M_w$ ~700–1500 kDa composite heteropolymeric (RG-I and AX) fraction [3]. Other than the phenotypic and genotypic characteristics of flaxseed, the aqueous extraction conditions, i.e., temperature, pH, ionic strength, seed-to-water ratio, the implementation of extraction-aiding physical processing (e.g., ultrasound or microwave), and the gum purification and dehydration methods determine the structure conformation and consequently, the inherent techno-functional properties (e.g., thickening, gelling, stabilising) of flaxseed gums [4].

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Hellebois, T.; Fortuin, J.; Gaiani, C.; Soukoulis, C. Impact of Flaxseed Gums on the Colloidal Changes and In Vitro Digestibility of Milk Proteins. *Foods* **2022**, *11*, 4096. https://doi.org/10.3390/foods11244096

Academic Editors: Peng Zhou and Song Miao

Received: 15 November 2022

Accepted: 13 December 2022

Published: 18 December 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Environmental Research and Innovation (ERIN) Department, Luxembourg Institute of Science and Technology (LIST), L-4362 Esch-sur-Alzette, Luxembourg</footnote>

<footnote n="2">² Laboratoire d'Ingénierie des Biomolécules (LIBio), Université de Lorraine, F-54000 Nancy, France</footnote>

<footnote n="3">³ Food Quality and Design Group, Wageningen University and Research, 6708 NL Wageningen, The Netherlands</footnote>

<footnote>* Correspondence: christos.soukoulis@list.lu; Tel.: +35-227-5888-5106</footnote>
Milk proteins i.e., caseins ($\alpha_{s-}$, $\beta$- and $\kappa$-) and whey proteins (primarily $\beta$-lactoglobulin, $\alpha$-lactalbumin, lactoferrin, immunoglobulins, and proteose peptones), possess a substantial techno-functional and dietary role in food product design [5]. From a nutritional viewpoint, milk proteins are widely recognised as natural sources of essential amino acids and bioactive peptides [6]. Due to their inherent structure conformational dissimilarity, milk proteins exert a diversified colloidal response to intragastric conditions (i.e., acid and pepsin-induced coagulation), which in turn can impact their resistance to the peptic cleavage that takes place during gastrointestinal transit. Due to their loose and more flexible structure, caseins undergo extensive intragastric coagulation compared to the well-defined globular structure conformation of whey proteins [6,7]. On the one hand, $\alpha_{s-}$, $\beta$- and $\kappa$-caseins are less resistant to peptic cleavage than $\alpha$-lactalbumin and $\beta$-lactoglobulin as a result of the selectivity of pepsin for the looser structure of the former. On the other hand, it is well documented that the ingestion of caseins favours the formation of large aggregates in the stomach, impeding gastric emptying, and thus, caseins exhibit a more sustained postprandial aminoacidemic response than whey proteins [8]. Other parameters associated with food processing, such as homogenisation, and enzymatic or thermal pretreatment, also have been reported to influence the digestibility of milk proteins [6,9–11].

Polysaccharides play a substantial role in the food industry, serving as natural sources of dietary fibre and as structuring, texturising, thickening, and stabilising agents [12]. It is well documented that mixing milk proteins with polysaccharides may result in associative (complex coacervation) or non-associative segregative phase separation phenomena [13]. The colloidal interplay between milk proteins and polysaccharides impacts the textural and structural aspects of the food matrix and modulates its colloidal configuration throughout oro-gastrointestinal transit [14]. In this context, the presence of polysaccharides in proteinaceous-rich food matrices may sustain their disintegration rate and suppress the susceptibility of milk proteins to peptic cleavage [15–20]. Borreani and co-workers [15] demonstrated that regardless of the protein source (i.e., skim milk powder, whey protein isolate, or sodium caseinate), anionic polysaccharides (sodium alginate) exerted a better ability to hinder the intragastric pepsinolytic activity than the non-ionic exemplars (konjac gum). Studying binary blends of sodium alginate with microparticulated whey protein, Koutina and co-workers [17] reported that the degree of the peptic cleavage of whey proteins was adversely associated with the molecular weight of sodium alginate, which was ascribed to the formation of protein-polysaccharide coacervates. Implementing a semi-dynamic in vitro digestion model, Markussen and co-workers [20] demonstrated that the occurrence of segregative phase separation and complex coacervation phenomena are the major drivers of the intragastric-induced peptic cleavage of milk protein concentrate (MPC). In a recent study, Hellebois and co-workers [21] demonstrated that alfalfa seed galactomannan modulated the acid-induced aggregation of milk proteins, leading to higher rates of free amino group release in the intestinal digesta.

Notwithstanding the broad techno-functionality of flaxseed gums, their interplay with proteins under in vitro digestion-simulating conditions remains unexplored. The present study aimed to understand the impact of two different types of flaxseed gums (brown and golden) on the in vitro digestibility of two milk protein (whey protein isolate and sodium caseinate) liquid food models. From a mechanistic standpoint, it is hypothesised that flaxseed gum, due to its anionic character, can modulate the peptic cleavage of milk proteins through its ability to control the acid-induced protein aggregation phenomena and hinder the proteases diffusivity sterically at the solid-liquid interface boundaries.

## 2. Materials and Methods

### 2.1. *Materials*

Whey protein isolate (WPI) PRODIET 90S, with a protein content of 87.7% wt., was kindly donated by Ingredia (Arras, France). Sodium caseinate (NaCN, protein purity of 89.4% wt.) was obtained from Sigma-Aldrich (Leuven, Belgium). The exact protein content of the milk proteins was determined using the Dumas method (Elementar Vario Cube,
Langensenbold, Germany). The gum extraction and isolation from golden and brown
flaxseed (GF and BF, respectively) was conducted at pH = 8 following the method explicitly
described by Hellebois and co-workers [3]. All the other chemicals used were purchased
from Sigma-Aldrich and were of analytical grade.

## 2.2. *Preparation of the Milk Protein Food Models*

The appropriate amount of protein powder was dissolved in ultrapure water (18.2 mΩ, Merck Millipore, Burlington, MA, United States) to achieve a protein concentration of 10% wt. and was stirred until complete dissolution. Following the complete dissolution of the protein, the solution was heat treated at 80 °C for 20 min in a shaking water bath (SW22, Julabo, Seelbach, Germany). The solution obtained was then homogenised three times using a high-pressure homogeniser (Panda plus 2000, GEA, Düsseldorf, Germany) at 500 bar before being centrifuged at 18,500 × g for 10 min to remove the non-solubilised particles. The homogeneity of the protein suspension was confirmed using laser light scattering particle size analysis (Mastersizer 3000, Malvern Instruments, Worcestershire, United Kingdom). To prevent bacterial growth, 0.02% wt. sodium azide (NaN₃) was added. BF or GF was then added to the solution to achieve 0.1, 0.5, and 1.0% wt. and the solution was solubilised overnight. Both WPI and NaCN solutions were prepared using this method.

## 2.3. *In Vitro Gastrointestinal Processing of the Milk Protein Food Models*

The in vitro digestion experiment was based on the INFOGEST protocol [22], as previously employed in protein-polysaccharide biopolymer systems. Briefly, ten millilitres of the initial matrices (WPI/BF or GF; NaCN/BF or GF) were diluted 1:1 with simulated saliva fluids containing 75 U mL⁻¹ of α-amylase and incubated for 3 min in a shaking water bath (SW22, Julabo, Seelbach, Germany) at 37 ± 0.1 °C, 100 min⁻¹. The oral bolus was then diluted 1:1 with simulated gastric fluids containing 2000 U mL⁻¹ of pepsin, adjusted with HCl 1 M at pH = 2.5 and incubated for 2 h at 37 °C ± 0.1 °C, f = 100 min⁻¹. The intestine phase was initiated by diluting 1:1 the gastric chymes with simulated intestine fluids containing the bile salts and 100 U mL⁻¹ of pancreatin and was then adjusted to pH = 7 with NaOH 1 M. Samples were collected in the initial food matrix at the end of the oral phase and at t = 0, 5, 10, 20, 30, 60, 90, and 120 min in the gastric and intestinal phases. The proteases activity was ceased by the addition of a protease inhibitor mix (Cytiva, Marlborough, MA, USA), and the samples were either kept in ice for direct analysis or quench frozen in liquid nitrogen and kept at -80 °C until use.

## 2.4. *Steady State Flow Rheological Measurements*

The initial food matrices and the oral, gastric, and intestinal phases were rheologically assessed. The flow behaviour of the systems was measured at a shear rate from 0.1 to 1000 s⁻¹ at 37 °C using a double gap (DG 26.7) geometry of 27.1 mm diameter mounted on an oscillatory rheometer (MCR 302, Anton Paar, Graz, Austria). The apparent viscosity of the food matrices and digesta was measured at 50 s⁻¹. The contribution of the sedimentable aggregates formed during the gastric phase on the rheological properties was assessed by measuring the soluble fraction after a centrifugation step at 4500 × g for 5 min.

## 2.5. *Confocal Laser Scanning Microscopy (CLSM)*

The microstructure of the mixed biopolymer systems was characterised using a CLSM microscope (LSM 880 with Airy scan, Zeiss, Jena, Germany) equipped with a ×10 objective. The detailed sample preparation and analysis method is described by Hellebois and co-workers [21].

## 2.6. *Particle Size Distribution Measurements*

The volume-weighted particle size distributions and the volume-weighted mean particle diameters (d₄,₃) of aliquots obtained from the initial food matrix and the oral,
gastric, and intestinal phases obtained were determined using static laser light scattering (Mastersizer 3000, Malvern Panalytical, Malvern, United Kingdom). The refractive indexes of whey protein, sodium caseinate, and water were set at 1.45, 1.34, and 1.33, respectively.

## 2.7. Protein Digestibility

The cleavage of the intact milk proteins into peptides during the in vitro gastrointestinal processing steps was monitored by SDS-PAGE analysis adopting the protocol described by Hellebois and co-workers [23]. Prior to the analysis, all systems were diluted to 5, 5, 10, and 20 µg of proteinaceous matter per well for the initial food matrix, and oral, gastric, and intestinal chymes, respectively. The samples, XT sample buffer (Bio-Rad, Hercules, CA, United States), and 0.6 XT reducing agent (Bio-Rad) were mixed at a 1:0.25:0.05 ratio, heat treated at 95 °C for 5 min, and then transferred into the SDS-PAGE gel wells.

## 2.8. Protein Hydrolysis Quantification

The amino acid release throughout the in vitro digestibility was assessed by the o-phthal dialdehyde (OPA) method carefully described by Hellebois and coworkers [21]. The degree of hydrolysis of the protein, after considering each dilution step, was determined as follows:

$$DH (\%) = \frac{NH_{2\text{digested}} - NH_{2\text{FM}}}{NH_{2\text{total}} - NH_{2\text{FM}}} \times 100 \quad (1)$$

where DH is the degree of hydrolysis, and NH₂digested, NH₂FM, and NH₂total denote the primary amino group content of the digested samples, food matrix, and hydrolysed FM aliquots, respectively.

## 2.9. Statistical Analyses

The normal distribution of the experimental data was confirmed by implementing the Shapiro–Wilk test and Q–Q plot representation normality tests. One-way ANOVA coupled with Tukey’s post-hoc means comparison test was performed to assess the significance of the tested parameters (milk protein type, FG type and content) on the in vitro digestibility data. All analysed were carried out using Origin 2019b software (OriginLab Inc, Northampton, MA, United States).

# 3. Results and Discussion

## 3.1. Rheological Behaviour and Colloidal Changes during the In Vitro Simulated Oral Processing

In keeping with our previous findings [21], the milk protein food models exhibited a Newtonian flow behaviour with the apparent viscosities of 23.8 and 10.8 mPa s for NaCN and WPI, respectively. Upon FG addition, a proportional increase to the FG concentration increase was observed in the apparent viscosities of the food matrices (Figure 1). No significant differences were found in the average macroviscosity values of the NaCN-based food models as a function of the FG phenotype. On average, the GF-stabilised WPI-based food models exhibited significantly (p < 0.05) higher viscosities than their BF counterparts (i.e., 914 and 644 mPa s, respectively). The rheological behaviour of colloidal dispersions in the presence of polysaccharides is known to be influenced by the molecular properties of the latter, such as the molecular weight, intrinsic viscosity, hydrodynamic radius, and surface charge density [24,25]. According to our previous observations [3], GF gum was characterised by a significantly higher molar mass and intrinsic viscosity (Mw = 1340 kDa, [η] = 6.6 dL g⁻¹) than BF gum (Mw = 1147 kDa, [η] = 5.1 dL g⁻¹). Except for the direct contribution of the polysaccharidic fraction to the viscosity of the bulk aqueous phase, the protein–polysaccharide molecular interactions e.g., electrostatic forces (electrostatic complexation or repulsion) and hydrogen bonding can also result in antagonistic or synergistic effects on the bulk rheology [26]. Although flaxseed gum is recognised as an anionic biopolymer, differences in the absolute surface charge density have been reported, i.e., |ζ| = 25–30 and 32–34 mV (pH = 7) for GF and BF gum, respectively [3].
Close to pH neutrality conditions (i.e., food models and oral boluses) implies that the observed reduction in the macroviscosity of the BF-stabilised milk protein solutions is due to the stronger repulsive forces between the BF and milk protein polymer chains. It should be noted that the macroviscosity dependence on the FG concentration followed a linear trend, with the highest slopes detected in the WPI-based solutions. Contrary to non-ionic polysaccharides, e.g., alfalfa galactomannan [21], where the macroviscosity of the milk protein solutions was radically changed at $c > c^*$, in the FG-stabilised protein food models, it was not possible to experimentally justify such a dependence on the critical concentration ($c^*$) i.e., $c^*_{GF} = 0.53$ and $c^*_{BF} = 0.76\%$ [3].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Impact of the golden (GF) and brown (BF) flaxseed gums (FG) on the apparent viscosities of the whey protein isolate (WPI) (A) and sodium caseinate (NaCN) (B) based food models and the oral, gastric, and intestine phases obtained. The term “gastric centrifuged” refers to the samples obtained following centrifugation of the gastric chymes at 4500 g for 5 min. a–d, A–D Different letters between the bars denote a significant difference among the samples differing in FG content and phenotype (uppercase) or in vitro digestion step (lowercase).</figcaption></figure>

While the solute dilution degree plays a central role in the viscosimetric response of the oral boluses, other parameters, such as the changes in counterion composition and ionic strength upon mixing with the artificial saliva, may be equally influential. As illustrated in Figure 1, admixing the initial protein food models with the artificial saliva was associated with an approximately one-order decrease in the macroviscosity values. The milk protein type and FG phenotype had a significant impact ($p < 0.05$) on the macroviscosity of the oral food boluses obtained. Hence, the WPI systems stabilised by BF experienced the highest decrease in macroviscosity upon oral bolus formation. On the other hand, the impact of the FG concentration on the macroviscosity of the model food boluses was unclear ($p > 0.05$).

To gain insight into the colloidal changes occurring during the oral processing, as influenced by the FG type and concentrations, the initial food models and the oral food
boluses obtained were analysed by means of CLSM (Figures 2 and 3) and static laser light scattering (SLS) (Figure 4). FG in the protein food models was associated with segregative microphase separation phenomena as a cause of the intermolecular repulsive forces between the FG and milk proteins [27]. From a mechanistic point of view, segregative phase separation in polysaccharide-protein aqueous systems is associated with either the Flory-Huggins or depletion flocculation theories [28]. In colloidal dispersions (i.e., biopolymer microspheres dispersed into the semi-dilute polysaccharide-rich aqueous phase), segregative phase separation is driven primarily by the depletion-flocculation mechanism, where the milk proteins behave as colloidal microspheres. Indeed, it has previously been demonstrated that both WPI and NaCN may exhibit a microsphere-like behaviour stemming from either their nano-aggregated (WPI) or rod-like self-assembled (NaCN) structure conformation [13,29].

As illustrated in Figure 2, in the WPI-based food models, a similar segregative phase separation pattern was observed between BF and GF. At $c_{FG} = 0.1\%$ wt., the WPI food models exhibited a biphasic w/w emulsion-like conformation where FG-rich water microdroplets were uniformly distributed into the continuous protein-rich aqueous phase. At $c_{FG} \ge 0.5\%$ wt., the FG-stabilised WPI food models had a heterogenous microstructure as evidenced by the presence of protein microaggregates dispersed into the continuous phase. Given that the aggregative phenomena took place at $c \sim c^*$ ($c^*_{GF} = 0.53$ and $c^*_{BF} = 0.76\%$ wt. [3]), it is assumed that the transition from the dilute to the semi-dilute regime promotes the noncovalent self-association of the whey proteins [27]. To a lesser extent, electrostatic complexation between FG and WPI proteins may also occur. By adding the artificial saliva, a decisive reduction in the prominence of the demixing phenomena in the systems containing $0.1\%$ wt. was observed (Figure 2 (B1,B2,E1,E2)), which can be ascribed to the shifting of the thermodynamic equilibrium close to the binodal curve [27]. At $c_{FG} \ge 0.5\%$ wt., the dilution of the initial food matrices with the oral fluids resulted in a sparsely heterogeneous microstructure (Figure 2 (C1,C2,F1,F2)).

In contrast to the WPI exemplars, the NaCN-based food models exerted a demixed water-in-water emulsion-like microstructure for the whole range of $c_{FG}$, with the BF-stabilised systems being more pronounced (Figure 3(A1-G1)). Interestingly, the emulsion-like segregated microstructure was detected in the entity of the NaCN-based oral boluses, but to a lesser extent (Figure 3(A2-G2)). No clear interrelationship between the FG concentration and the magnitude of the segregative phenomena was noted. According to the SLS measurements (Figure 4), the presence of soluble protein aggregates ($d_{4,3} = 258$ and 225 nm for WPI and NaCN, respectively) was confirmed. The impact of FG on the particle size of the protein aggregates was dependent on the protein type. As for the WPI food models, a proportional increase to the FG concentration increase was observed in the mean size of the soluble protein aggregates. Although the coacervation between FG and WPI at pH > $pH_c$, i.e., ~5.2–5.4 [30] is considered unlikely to occur as both biopolymers were negatively charged, the formation of soluble protein aggregates under weakly acidic conditions (pH > 5.5) has been reported, as a result of the intermolecular association between FG and WPI through oppositely charged patch junction zones [31].
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Impact of flaxseed gums on the in vitro oro-gastric induced colloidal changes of the WPI-based food models adopting the INFOGEST 2.0 in vitro digestion model. **1:** Food matrix, **2:** Oral phase; **3:** gastric (start); **4:** gastric (end); (**A:** WPI only; (**B**–**D:** WPI + BF; (**E**–**G:** WPI + GF. Scale bar = 100 µm.)</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** CLSM assessment of the impact of flaxseed gums on the in vitro oro-gastric induced col- loidal changes of the NaCN-based food models adopting the INFOGEST 2.0 in vitro digestion model. **1:** Food matrix, **2:** Oral phase; **3:** gastric (start); **4:** gastric (end); (**A:** NaCN only; (**B**–**D:** NaCN + BF; (**E**–**G:** NaCN + GF. Scale bar = 100 µm.)</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Changes in the protein aggregate particle size as they occurred throughout the simulated in vitro oro-gastrointestinal processing of the WPI-(A) and NaCN-(B) based food models. <sup>a-c</sup> Different letters between the bar phases denote a significant difference among the samples differing in the amount of gum for the same digestive.</figcaption></figure>

Nonetheless, the FG phenotype was not influential on the mean size of the soluble protein aggregates. The addition of FG into the NaCN food models did not modify their colloidal properties for the entire range of gum concentration. In mild acidic conditions, NaCN is characterised by a significantly higher negative surface charge density than WPI [21], and hence it is assumed that the strong electrostatic repulsion between the FG and NaCN molecules hampered the formation of soluble aggregates. On admixing with the artificial saliva, none of the protein food models tested experienced any significant change in terms of the protein aggregate mean size, which implies that the changes in the ionic strength and counterion composition of the food boluses did not modify the mechanism of the FG–milk protein interactions.

## 3.2. Rheological Behaviour and Colloidal Changes under In Vitro Gastro-Intestinal Conditions

To offset the technical limitations of the static in vitro digestion models associated with the dynamics of the physiology of digestion [32], the gastric fluids were added at a constant rate and under constant stirring to achieve the target gastric pH<sub>gastric</sub> = 2.5 within 5 min (t = 0 min). The enzymatic activity of pepsin was predetermined, and the amount of the HCl required for achieving the pH<sub>gastric</sub> = 2.5 was customised per food model to minimise the buffering effect of the milk proteins and FG. In addition, the peptic cleavage process was hindered using a protease inhibitor cocktail, whilst a real-time analysis of the samples, stored in an ice bath, was conducted on the same day.

All gastric phases exhibited a suspension-like behaviour owing to the formation of acid-mediated protein aggregates. The viscosimetric characterisation of the gastric phases (before and after centrifugation) revealed that the NaCN-based food models experienced the highest increase in apparent viscosity at 50 s<sup>-1</sup> compared to the WPI exemplars (Figure 1). Interestingly, centrifugation resulted in an abrupt (almost four-order magnitude) decline
in the macroviscosity values of the NaCN-based gastric chymes, which can be attributed to the depletion of their sedimentable protein aggregate matter (Figure 1B). On the other hand, no significant differences in the macroviscosity of the actual and centrifuged WPI gastric phases were observed, indicating the presence of finer and more uniform protein aggregates (Figure 1A).

As well illustrated in the acquired CLSM micrographs (Figure 3) and SLS analysis data (Figure 4), the NaCN-based gastric phases comprised large (d$_{4,3}$ = 91 µm) acid-induced protein aggregates. Although caseins may also undergo pepsin-induced clotting, in the case of NaCN, this is considered unlikely owing to its nonmicellar structure conformation [6]. The bimodal pattern of the particle size distribution curves of the NaCN-based gastric chymes is generally confirmatory of the data in the literature [15,20,21]. At the end of the gastric processing step, a reduction in the d$_{4,3}$ values (i.e., d$_{4,3}$ = from 91 to 71 µm) throughout digestion was observed as a result of the peptic cleavage and the in vitro simulating mechanical (antral) forces. The presence of increasing FG content was associated with a proportional reduction in the mean size of the NaCN aggregates, as demonstrated by the decay in the intensity of the above-micron particle population peak (Figure 5). Nonetheless, no significant differences in the d$_{4,3}$ values were observed at the start and endpoint of the in vitro gastric processing (i.e., 15.7 and 17.3 µm, respectively). Interestingly, the decay in the d$_{4,3}$ values followed an exponential trend as follows:

$$d = d_0 + A \cdot \exp(-k_c \cdot c_{FG}) \quad (2)$$

where d, d$_{0}$ denote the d$_{4,3}$ values (in µm) of the NaCN aggregates in the absence (d$_{0}$) and presence of FG (d), A is a constant, and k$_{c}$ is the NaCN aggregate size reduction rate (%$^{-1}$) as a function of the FG content.

<figure><img src="image_8.png" /><figcaption>Figure 5. Particle size distribution profiles of the initial and end NaCN-based gastric phases (1: t = 0 min and 2: t = 120 min) as influenced by the phenotype and concentration of flaxseed gum. (A1,A2) = golden flaxseed and (B1,B2) = brown flaxseed.</figcaption></figure>
According to our calculations, the $k_c$ values at $t = 0$ min were $11.02 \pm 0.82$ and $11.42 \pm 0.99\%^{-1}$ for GF and BF, respectively, which were significantly higher than the corresponding ones at $t = 120$ min (i.e., $4.54 \pm 0.51$ and $4.18 \pm 0.73\%^{-1}$ for GF and BF, respectively). This suggests that the FG impeded the acid-induced aggregation respon- siveness of NaCN to the gastric fluids resulting in smaller and more uniform particle size distributions. Indeed, this was confirmed by employing the CLSM (Figure 3) and SLS measurements (Figure 5). Although the mechanism was not evidenced, it can be assumed that the colloidal changes observed were governed by the ability of FG to adsorb onto the NaCN particles via electrostatic complexation upon acidification, and its direct effect on the continuous aqueous phase macroviscosity of the gastric chymes, which diminished the diffusion rate of pepsin into the solid (protein aggregates)-water interface [15,20,21].

In opposition to the aforementioned findings, the initial WPI-based gastric phases exhibited a relatively uniform microstructure ($d_{4,3} = 0.69$ µm), in agreement with our previous observations [21]. In the absence of FG, a bimodal particle size distribution pattern was confirmed in the WPI exemplars. A significant increase in WPI acid-formed aggregates ($d_{4,3} = 4.9$ µm) was detected at the end of the gastric processing step. Markussen and co- workers [20] have reported a similar behaviour in semi-dynamically digested milk protein concentrate (MPC)-based food models. Upon mixing with the gastric fluids, all gastric phases containing FG were characterised by highly polydisperse colloidal suspensions (Figure 6). Despite the positive correlation between the $d_{4,3}$ values and the $c_{FG}$, the mean size of the particles was substantially changed only at $c_{FG} = 1\%$ wt. (Figure 4). GF addition resulted in the formation of the largest protein aggregates, which may be associated with its higher negative surface charge density in highly acidic conditions (i.e., pH = 2–4) than BF [3], favouring its electrostatic bridging with whey proteins. Throughout gastric processing, a re-organisation of the morphological features of the WPI-FG gastric phases was observed, as illustrated in the CLSM micrographs (Figure 2). As seen in Figure 6, the submicron fraction of the particles (representing the WPI nanoaggregates colloidal dispersion), was decreased at the expense of the fraction larger than a micron, i.e., 1–50 µm, when $c_{FG} > 0.5\%$ wt. This implies that the formation of WPI–FG coacervates took place throughout gastric processing leading to the progressive depletion of the colloidal WPI soluble aggregates and, eventually, to the increase of the sedimentable biopolymer matter. Indeed, the centrifugation of the gastric endpoint phases resulted in a significant decrease in their apparent viscosities (Figure 1). Omitting the discrepancies due to ionic strength, the macroviscosity of the centrifuged gastric phases (i.e., 1.6–2.6 mPa s) was significantly lower than that of the dilute flaxseed gum solutions (i.e., 9.5–17.6 mPa s [3]). Hence, the colloidal changes in the WPI-FG gastric phases are driven primarily by an electrostatic complexation mechanism.

### 3.3. *In Vitro Digestibility of the Milk Proteins*

The peptic cleavage of milk proteins and the degree of proteolysis were measured employing SDS-PAGE (Figures 7 and 8) and OPA (Figure 9) assay, respectively. To gain insight into the kinetics of the peptic cleavage of the milk proteins, the SDS-PAGE images obtained (Figure 7) were analysed densitometrically using ImageJ software (Figure 8), as previously reported by Hellebois and co-workers [21]. Different aliquots of the oral, gastric, and intestinal phases (i.e., 5, 10, and 20 µg of proteinaceous matter per well) were used to improve the detection of the intact and fragmented proteins. Three significant bands identified at 25–35, 18.2, and 14.2 kDa allotted to total caseins ($a_s$ –, $\beta$ –, and $\kappa$–casein), $\beta$-lactoglobulin ($\beta$-Lg), and $\alpha$-lactalbumin ($\alpha$-La) were semi-quantified and the densitometric data were fitted to an exponential decay model [21,33,34], as follows:
<figure><img src="image_9.png" /><figcaption>Figure 6. Particle size distribution profiles of the initial and end WPI-based gastric phases (1: t = 0 min and 2: t = 120 min) as influenced by the phenotype and concentration of flaxseed gum. (A1,A2) = Golden flaxseed and (B1,B2) = Brown flaxseed.</figcaption></figure>

$$c = c_{120} + (c_0 - c_{120})\exp(-kt) \quad (3)$$

where $c_0$ and $c_{120}$ denote the normalised percentage of the residual intact protein at the beginning or end of the gastric or intestinal processing step, t is the gastric or intestinal processing time (in h), and k ($h^{-1}$) is the peptic cleavage rate of proteins to polypeptides. The time $\tau$ (in h) required to cleave the 50% of the initial proteins was calculated according to the half-time (for first-order kinetics) equation:

$$\tau = \frac{\ln 2}{k} \quad (4)$$

In the absence of FG (Figure 8D), β-lactoglobulin underwent a faster peptic cleavage during the gastric processing step than α-lactalbumin and total caseins ($k = 14.9 \pm 0.9$, $2.7 \pm 0.7$ and $1.9 \pm 0.8 \text{ min}^{-1}$, respectively). Besides the structure conformational aspects of the proteins [35], several processing parameters associated with the implementation and severity of heat treatment, pressure processing, enzymatic pre-treatment, etc., may affect the sensitivity of milk proteins to peptic cleavage [9,10,36,37]. According to Figure 8A–C, the presence of FG was associated with a higher amount of intact β-lactoglobulin (42–57%) and total caseins (22–27%) and a lower amount of intact α-lactalbumin (16–28%), at the end of the gastric processing step.
<figure><img src="image_10.png" /><figcaption>Figure 7. *Cont.*</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 7.** SDS-PAGE electropherogram of the gastric (1) and intestinal (2) phases of the NaCN-(A–G) and WPI-(H–N) based food models as influenced by the phenotype and concentration (0.1–1% wt.) of flaxseed gum. An SDS-PAGE analysis was conducted on gastric and intestine aliquots sampled at predetermined time points (0, 5, 10, 20, 30, 60, 90, and 120 min). Abbreviations used: β-Lg: β-lactoglobulin, α-La: α-lactalbumin.</figcaption></figure>
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 8.** SDS-PAGE densitometric analysis illustrating the kinetics of β-lactoglobulin (A), α-lactalbumin (B), and total caseins (C) fate in the gastric (t = 0–120 min) and intestinal (t = 120–240 min) phases as influenced by the phenotype and concentration of flaxseed gum. The calculated rates for β-Lg, α-La and casein are depicted in (D) (gastric rates) and (E) (intestinal rates).

<sup>a-c</sup> Different letters for the same protein type denote a significant difference among the samples. Abbreviations used: β-Lg: β-lactoglobulin, α-La: α-lactalbumin.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 9.** Degree of protein hydrolysis in the WPI-(A) and NaCN-(B) based gastric and intestinal phases as influenced by the phenotype and concentration (0.1–1% wt.) of flaxseed gum. The protein hydrolysis rates of WPI and NaCN during gastric (k<sub>GASTRIC</sub>) and intestinal (k<sub>INTESTINAL</sub>) steps are displayed in (C–F). <sup>a–d</sup>, <sup>A–C</sup> Different letters between the samples denote a significant difference.</figcaption></figure>

It is presumed that the differences in the pepsin-induced cleavage of whey proteins may stem from the differences in their secondary structure conformation and surface charge density. β-Lactoglobulin, due to its higher isoelectric point (pI<sub>β-Lg</sub>~5.2, pI<sub>α-La</sub>~4.6), can undergo electrostatic complexation with FG earlier than α-lactalbumin, leading to the steric hampering of the peptic cleavage in the former case. From a kinetic perspective, α-lactalbumin and total caseins underwent pepsin-induced peptic cleavage at a slower pace (τ = 15.2 and 21.4 min, respectively) than β-lactoglobulin (τ = 2.8 min, Figure 8D). In the presence of FG, the time required for cleaving the 50% of the initial intact proteins was shortened (τ = 2.1, 6.1, and 5.7 min for β-lactoglobulin, α-lactalbumin, and total caseins, respectively), although its phenotype did not play any significant role (Figure 8D). The ability of the polysaccharides to modify the kinetics of the pepsin cleavage has been reported in previous studies [15–17,20,21]. For anionic polysaccharides, the major mechanistic pathway associated with the reduced digestibility of milk proteins concerns their direct molecular interactions that sterically impede the access of proteases to the cleaving sites of the protein substrate [15,20]. In addition, some anionic polysaccharides i.e., sodium alginate, can selectively inhibit the enzymatic activity of proteases (i.e., pepsin) [38]. Although these results could lead us to speculate to some extent on the pepsin-inhibiting capacity of FG, further experiments are required in order to confirm this hypothesis.

By monitoring the temporal changes in the intact protein residue in the intestinal phases, an almost complete cleavage of the milk proteins was detected as less than 2–3%, and 10% of intact whey proteins and total caseins remained at the end of the intestinal processing. Based on the calculated kinetic parameters (Figure 8E), no significant differences in the time required for cleaving the 50% of the intragastric residual intact protein matter
were found i.e., $\tau = 10.3, 9.3$, and $13.7$ min for β-lactoglobulin, α-lactalbumin and total caseins, respectively. The presence of FG did not modify significantly the peptic cleaving rate of α-lactalbumin (i.e., $\tau = 11.3$ min). On the other hand, the pancreas mediated cleaving rates of β-lactoglobulin and total caseins were increased (i.e., $\tau = 4.7$ and $7$ min, respectively); however the differences observed were nonsignificant.

In order to better understand the impact of the FG phenotype and content on the extent of proteolysis throughout the in vitro simulated gastric and intestinal processing steps, the data obtained from the OPA assay were fitted to the following mathematical model, as previously reported by Le Feunteun and co-workers [33]:

$$ \textit{Gastric phases : } DH_{\text{gastric}} = DH_{\text{gastric}, \infty} \exp(1 - k_{\text{gastric}} t) \quad (5) $$

where the $DH_{\text{gastric},\infty}$ denotes the degree of hydrolysis (%) at the end of the gastric processing and $k_{\text{gastric}}$ is the rate of hydrolysis (in $h^{-1}$), and t the digestion time.

$$ \textit{Intestinal phases : } DH_{\text{intestinal}} = DH_{\text{intestinal},120} + (DH_{\text{intestinal},0} - DH_{\text{intestinal},120})\exp(-k_{\text{intestinal}}t) \quad (6) $$

where the $DH_0$, and $DH_{120}$ denote the degree of hydrolysis (%) achieved at the beginning and end of the intestine processing step, $k_{\text{intestinal}}$ is the rate of hydrolysis (in $h^{-1}$), and t the digestion time.

As displayed in Figure 9A,B, the extent of hydrolysis (% DH) achieved at the end of the gastric processing step was at most 9–11%, regardless of the milk protein type and FG presence. It should be noted that in all gastric chymes, the calculated hydrolysis rates were highly dependent on the changes in the amount of the free amino groups during the very early stages of the gastric processing i.e., within the first 30 min. Similar behaviour in the digestion of food biomacromolecules as influenced by the food matrix characteristics has been reported [39]. In line with previous studies [21,36,40], the rate of hydrolysis in the WPI only gastric phases was substantially higher than their NaCN only counterparts (Figure 9C,D). This is generally associated with the restricted colloidal responsiveness of WPI to pH and ionic strength changes that sterically favour the peptic cleaving activity of pepsin.

Incorporating FG into the initial food models induced a significant increase in the hydrolysis rates achieved during the gastric processing step, i.e., $k_{\text{gastric}} = 10.5 \pm 1.1$ and $21.1 \pm 3.5 h^{-1}$ (for WPI-based chymes in the absence and presence of FG, respectively) and $1.5 \pm 1.1$ and $19.4 \pm 6.4 h^{-1}$ for their NaCN counterparts. This can be ascribed to the ability of FG to control the acid-induced aggregation of the milk proteins and, therefore, the accessibility of pepsin to the protein substrate. The pepsinolysis-modulating role of FG was mainly flaxseed phenotype-driven (Figure 9C,D). In this context, the GF-stabilised food models underwent faster pepsinolysis ($k_{\text{gastric}} = 22 \pm 1.9$ and $23.9 \pm 5.7 h^{-1}$ for the WPI and NaCN-based food models, respectively) than their BF exemplars ($k_{\text{gastric}} = 20.3 \pm 3.9$ and $15 \pm 1.0 h^{-1}$ for the WPI- and NaCN-based food models, respectively). However, only in the case of NaCN the differences in the rates of hydrolysis were significant ($p < 0.05$). Note-worthily, the $k_{\text{gastric}}$ rates obtained hereby were significantly higher than those observed in the case of non-ionic polysaccharides (i.e., alfalfa galactomannan) [21]. Hence, it appears that the ability of FG to restrict the ripening of the acid-mediated protein aggregates via the formation of electrostatic complexes of controlled mean size is significantly more influential than its direct contribution to the solid-liquid interface mass transfer phenomena. In keeping with the latter, the modulating role of FG concentration in the induced pepsinolysis was unclear as it was only in the case of the GF-stabilised NaCN food models that there was a proportional increase to the gum content increase in hydrolysis rates.

On the completion of the intestinal digestion, a substantial increase in the DH of the proteins was observed i.e., ca. 47% for individual milk proteins and 59–62% when FG was incorporated in the milk protein-based food models (Figure 9A,B). Nonetheless, the pancreas-mediated hydrolysis rates were significantly ($p < 0.05$) lower in the intestinal digesta containing FG (Figure 9E,F). Although the hydrolysis rates were higher in the
intestinal phases containing BF ($k_{\text{intestinal}}$ = 1.37 and 0.86 $h^{-1}$ for the WPI- and NaCN-based food models, respectively) than their GF counterparts ($k_{\text{intestinal}}$ = 1.05 and 0.63 $h^{-1}$ for the WPI- and NaCN-based food models, respectively), the differences were nonsignificant. Our findings postulate that the equilibrium between the coacervates and precipitates, free proteins, and polyanions had a substantial impact on the extent of pancreas-induced hydrolysis [41,42]. Upon admixing of the gastric chymes with the simulating intestinal fluids, the concomitant increase in the pH (from 2.5 to 7.0), ionic strength, and cationic species ($Na^+$, $K^+$, and $Mg^{2+}$) resulted in the progressive dissolution of the milk protein-FG electrostatic complexes and the release of free proteins that can be cleaved more easily from the pancreas. At this stage, the solubilisation kinetics of the coacervates, together with the impact of the thickening capacity of the FG on the interphase mass transfer kinetics, became more influential on the rate of the polypeptide peptic cleavage.

**4. Conclusions**

In the present work, the ability of FG to modulate the colloidal transformation of milk proteins and, consequently, their cleavage throughout in vitro digestion was elucidated. Regardless of its phenotype, FG was associated with segregative phase separation in both tested milk protein-based food models. Upon oral processing, only the WPI-based food boluses were colloidally responsive to FG due to the formation of soluble protein-gum complexes. Under intragastric simulating conditions, both milk proteins underwent acid-induced aggregation with the colloidal changes (mean particle size and morphological aspects) in the NaCN-based systems being critically influenced by the presence of FG; however, the FG phenotype played a minor role. On the completion of the gastric processing, the presence of FG was associated with a higher amount of residual intact total caseins and $\beta$-lactoglobulin but it did not significantly impact the intact $\alpha$-lactalbumin residue. No remarkable differences in the amount of the residual intact proteins at the end of the intestinal digestion were observed. Although the FG did not remarkably influence the degree of hydrolysis during the gastric processing step, FG resulted in a significant increase in the degree of hydrolysis in the intestinal digesta obtained. From a mechanistic point of view, the equilibrium between the FG-milk protein coacervates and free polyelectrolytes, as influenced by the physicochemical changes occurring during the gastrointestinal transit, is the principal modulator of the in vitro digestibility of the milk proteins. To a lesser extent, the ability of FG to increase the microviscosity in the solid-liquid interphase can also affect the sensitivity of milk proteins to peptic cleavage, particularly in pH conditions where the formation of insoluble protein-polysaccharide complexes is not favoured.

**Author Contributions:** Conceptualisation, T.H. and C.S.; methodology, T.H. and C.S.; validation, T.H. and C.S.; formal analysis, T.H., J.F. and C.S.; investigation, T.H. and J.F.; writing—original draft preparation, T.H.; writing—review and editing, T.H., J.F., C.G. and C.S.; visualisation, T.H.; supervision, C.G. and C.S.; project administration, C.S.; funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Luxembourg National Research Fund, Grant number CORE/2018/SR/12675439. The APC was funded by the Luxembourg National Research Fund.

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

**Acknowledgments:** The authors would like to acknowledge Manon Hiolle (Ingredia, France) for generously providing the WPI.

**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.
## References

1. Bailey, C.; Gomina, M.; Breard, J.; Bourmaud, A.; Davies, P. Variability of Mechanical Properties of Flax Fibres for Composite Reinforcement. A Review. *Ind. Crops Prod.* **2020**, *145*, 111984. [CrossRef]

2. Liu, J.; Shim, Y.Y.; Tse, T.J.; Wang, Y.; Reaney, M.J.T. Flaxseed Gum a Versatile Natural Hydrocolloid for Food and Non-Food Applications. *Trends Food Sci. Technol.* **2018**, *75*, 146–157. [CrossRef]

3. Hellebois, T.; Fortuin, J.; Xu, X.; Shaplov, A.S.; Gaiani, C.; Soukoulis, C. Structure Conformation, Physicochemical and Rheological Properties of Flaxseed Gums Extracted under Alkaline and Acidic Conditions. *Int. J. Biol. Macromol.* **2021**, *192*, 1217–1230. [CrossRef] [PubMed]

4. Soukoulis, C.; Gaiani, C.; Hoffmann, L. Plant Seed Mucilage as Emerging Biopolymer in Food Industry Applications. *Curr. Opin. Food Sci.* **2018**, *22*, 28–42. [CrossRef]

5. Walstra, P.; Walstra, P.; Wouters, J.T.M.; Geurts, T.J. *Dairy Science and Technology*, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2005; ISBN 978-0-429-11614-8.

6. Dupont, D.; Tomé, D. Chapter 20—Milk Proteins: Digestion and Absorption in the Gastrointestinal Tract. In *Milk Proteins*, 3rd ed.; Boland, M., Singh, H., Eds.; Academic Press: Cambridge, MA, USA, 2020; pp. 701–714. ISBN 978-0-12-815251-5.

7. Goulding, D.A.; Fox, P.F.; O'Mahony, J.A. Chapter 2—Milk Proteins: An Overview. In *Milk Proteins*, 3rd ed.; Boland, M., Singh, H., Eds.; Academic Press: Cambridge, MA, USA, 2020; pp. 21–98. ISBN 978-0-12-815251-5.

8. Boirie, Y.; Dangin, M.; Gachon, P.; Vasson, M.-P.; Maubois, J.-L.; Beaufrère, B. Slow and Fast Dietary Proteins Differently Modulate Postprandial Protein Accretion. *Proc. Natl. Acad. Sci. USA* **1997**, *94*, 14930–14935. [CrossRef] [PubMed]

9. Barbé, F.; Ménard, O.; Le Gouar, Y.; Buffière, C.; Famelart, M.-H.; Laroche, B.; Le Feunteun, S.; Dupont, D.; Rémond, D. The Heat Treatment and the Gelation Are Strong Determinants of the Kinetics of Milk Proteins Digestion and of the Peripheral Availability of Amino Acids. *Food Chem.* **2013**, *136*, 1203–1212. [CrossRef] [PubMed]

10. Böttger, F.; Dupont, D.; Marcinkowska, D.; Bajka, B.; Mackie, A.; Macierzanka, A. Which Casein in Sodium Caseinate Is Most Resistant to in Vitro Digestion? Effect of Emulsification and Enzymatic Structuring. *Food Hydrocoll.* **2019**, *88*, 114–118. [CrossRef]

11. van Lieshout, G.A.A.; Lambers, T.T.; Bragt, M.C.E.; Hettinga, K.A. How Processing May Affect Milk Protein Digestion and Overall Physiological Outcomes: A Systematic Review. *Crit. Rev. Food Sci. Nutr.* **2020**, *60*, 2422–2445. [CrossRef]

12. Williams, P.A.; Phillips, G.O. 1—Introduction to Food Hydrocolloids. In *Handbook of Hydrocolloids*, 2nd ed.; Phillips, G.O., Williams, P.A., Eds.; Woodhead Publishing Series in Food Science, Technology and Nutrition; Woodhead Publishing: Sawston, UK, 2009; pp. 1–22. ISBN 978-1-84569-414-2.

13. Turgeon, S.L.; Beaulieu, M.; Schmitt, C.; Sanchez, C. Protein–Polysaccharide Interactions: Phase-Ordering Kinetics, Thermodynamic and Structural Aspects. *Curr. Opin. Colloid Interface Sci.* **2003**, *8*, 401–414. [CrossRef]

14. Ye, A. Gastric Colloidal Behaviour of Milk Protein as a Tool for Manipulating Nutrient Digestion in Dairy Products and Protein Emulsions. *Food Hydrocoll.* **2021**, *115*, 106599. [CrossRef]

15. Borreani, J.; Llorca, E.; Larrea, V.; Hernando, I. Adding Neutral or Anionic Hydrocolloids to Dairy Proteins under in Vitro Gastric Digestion Conditions. *Food Hydrocoll.* **2016**, *57*, 169–177. [CrossRef]

16. David, S.; Wojciechowska, A.; Portmann, R.; Shpigelman, A.; Lesmes, U. The Impact of Food-Grade Carrageenans and Consumer Age on the in Vitro Proteolysis of Whey Proteins. *Food Res. Int.* **2020**, *130*, 108964. [CrossRef] [PubMed]

17. Koutina, G.; Ray, C.A.; Lametsch, R.; Ipsen, R. The Effect of Protein-to-Alginate Ratio on in Vitro Gastric Digestion of Nanoparticulated Whey Protein. *Int. Dairy J.* **2018**, *77*, 10–18. [CrossRef]

18. Liu, L.; Kong, F. Influence of Nanocellulose on in Vitro Digestion of Whey Protein Isolate. *Carbohydr. Polym.* **2019**, *210*, 399–411. [CrossRef] [PubMed]

19. Ma, R.; Lin, Z.; Wu, Y.; Gao, Z.; Hu, B.; Xu, L.; Fang, Y.; Nishinari, K. Modulating the in Vitro Gastric Digestion of Heat-Induced Beta-Lactoglobulin Aggregates: Incorporation with Polysaccharide. *Food Chem.* **2021**, *354*, 129506. [CrossRef]

20. Markussen, J.Ø.; Madsen, F.; Young, J.F.; Corredig, M. A Semi Dynamic in Vitro Digestion Study of Milk Protein Concentrate Dispersions Structured with Different Polysaccharides. *Curr. Res. Food Sci.* **2021**, *4*, 250–261. [CrossRef]

21. Hellebois, T.; Gaiani, C.; Soukoulis, C. Impact of Alfalfa (Medicago Sativa L.) Galactomannan on the Microstructural and Physicochemical Changes of Milk Proteins under Static in-Vitro Digestion Conditions. *Food Chem. X* **2022**, *14*, 100330. [CrossRef]

22. Brodkorb, A.; Egger, L.; Alminger, M.; Alvito, P.; Assunção, R.; Ballance, S.; Bohn, T.; Bourlieu-Lacanal, C.; Boutrou, R.; Carrière, F.; et al. INFOGEST Static in Vitro Simulation of Gastrointestinal Food Digestion. *Nat. Protoc.* **2019**, *14*, 991–1014. [CrossRef]

23. Hellebois, T.; Gaiani, C.; Paris, C.; Planchon, S.; Renaut, J.; Soukoulis, C. Data on the In-Vitro Digestibility of Acid Gels Prepared from Brewers' Spent Grain Protein Isolates. *Data Brief* **2021**, *37*, 107160. [CrossRef]

24. Bai, L.; Liu, F.; Xu, X.; Huan, S.; Gu, J.; McClements, D.J. Impact of Polysaccharide Molecular Characteristics on Viscosity Enhancement and Depletion Flocculation. *J. Food Eng.* **2017**, *207*, 35–45. [CrossRef]

25. Ghosh, A.K.; Bandyopadhyay, P.; Ghosh, A.K.; Bandyopadhyay, P. *Polysaccharide-Protein Interactions and Their Relevance in Food Colloids*; IntechOpen: London, UK, 2012; ISBN 978-953-51-0819-1.

26. Musampa, R.M.; Alves, M.M.; Maia, J.M. Phase Separation, Rheology and Microstructure of Pea Protein–Kappa-Carrageenan Mixtures. *Food Hydrocoll.* **2007**, *21*, 92–99. [CrossRef]

27. Hellebois, T.; Gaiani, C.; Cambier, S.; Noo, A.; Soukoulis, C. Exploration of the Co-Structuring and Stabilising Role of Flaxseed Gum in Whey Protein Isolate Based Cryo-Hydrogels. *Carbohydr. Polym.* **2022**, *289*, 119424. [CrossRef] [PubMed]
28. Doublier, J.-L.; Garnier, C.; Renard, D.; Sanchez, C. Protein-Polysaccharide Interactions. *Curr. Opin. Colloid Interface Sci.* **2000**, *5*, 202–214. [CrossRef]

29. Farrer, D.; Lips, A. On the Self-Assembly of Sodium Caseinate. *Int. Dairy J.* **1999**, *9*, 281–286. [CrossRef]

30. Liu, J.; Shim, Y.Y.; Shen, J.; Wang, Y.; Reaney, M.J.T. Whey Protein Isolate and Flaxseed (*Linum Usitatissimum* L.) Gum Electrostatic Coacervates: Turbidity and Rheology. *Food Hydrocoll.* **2017**, *64*, 18–27. [CrossRef]

31. Soukoulis, C.; Cambier, S.; Serchi, T.; Tsevdou, M.; Gaiani, C.; Ferrer, P.; Taoukis, P.S.; Hoffmann, L. Rheological and Structural Characterisation of Whey Protein Acid Gels Co-Structured with Chia (*Salvia Hispanica* L.) or Flax Seed (*Linum Usitatissimum* L.) Mucilage. *Food Hydrocoll.* **2019**, *89*, 542–553. [CrossRef]

32. Xavier, A.A.; Mariutti, L.R. Static and Semi-Dynamic in Vitro Digestion Methods: State of the Art and Recent Achievements towards Standardization. *Curr. Opin. Food Sci.* **2021**, *41*, 260–273. [CrossRef]

33. Le Feunteun, S.; Verkempinck, S.; Floury, J.; Janssen, A.; Kondjoyan, A.; Marze, S.; Mirade, P.-S.; Pluschke, A.; Sicard, J.; van Aken, G.; et al. Mathe matical Modelling of Food Hydrolysis during in Vitro Digestion: From Single Nutrient to Complex Foods in Static and Dynamic Conditions. *Trends Food Sci. Technol.* **2021**, *116*, 870–883. [CrossRef]

34. Verkempinck, S.H.E.; Guevara-Zambrano, J.M.; Infantes-Garcia, M.R.; Naranjo, M.C.; Soliva-Fortuny, R.; Elez-Martínez, P.; Grauwet, T. Gastric and Small Intestinal Lipid Digestion Kinetics as Affected by the Gradual Addition of Lipases and Bile Salts. *Food Biosci.* **2022**, *46*, 101595. [CrossRef]

35. Huppertz, T.; Chia, L.W. Milk Protein Coagulation under Gastric Conditions: A Review. *Int. Dairy J.* **2021**, *113*, 104882. [CrossRef]

36. Halabi, A.; Croguennec, T.; Bouhallab, S.; Dupont, D.; Deglaire, A. Modification of Protein Structures by Altering the Whey Protein Profile and Heat Treatment Affects in Vitro Static Digestion of Model Infant Milk Formulas. *Food Funct.* **2020**, *11*, 6933–6945. [CrossRef] [PubMed]

37. Macierzanka, A.; Böttger, F.; Lansonneur, L.; Groizard, R.; Jean, A.-S.; Rigby, N.M.; Cross, K.; Wellner, N.; Mackie, A.R. The Effect of Gel Structure on the Kinetics of Simulated Gastrointestinal Digestion of Bovine β-Lactoglobulin. *Food Chem.* **2012**, *134*, 2156–2163. [CrossRef] [PubMed]

38. Chater, P.I.; Wilcox, M.D.; Brownlee, I.A.; Pearson, J.P. Alginate as a Protease Inhibitor in Vitro and in a Model Gut System; Selective Inhibition of Pepsin but Not Trypsin. *Carbohydr. Polym.* **2015**, *131*, 142–151. [CrossRef] [PubMed]

39. Pälchen, K.; Bredie, W.L.P.; Duijsens, D.; Isaac Alfie Castillo, A.; Hendrickx, M.; Van Loey, A.; Raben, A.; Grauwet, T. Effect of Processing and Microstructural Properties of Chickpea-Flours on in Vitro Digestion and Appetite Sensations. *Food Res. Int.* **2022**, *157*, 111245. [CrossRef] [PubMed]

40. Rinaldi, L.; Gauthier, S.F.; Britten, M.; Turgeon, S.L. In Vitro Gastrointestinal Digestion of Liquid and Semi-Liquid Dairy Matrixes. *LWT—Food Sci. Technol.* **2014**, *57*, 99–105. [CrossRef]

41. Chai, C.; Lee, J.; Huang, Q. The Effect of Ionic Strength on the Rheology of PH-Induced Bovine Serum Albumin/κ-Carrageenan Coacervates. *LWT—Food Sci. Technol.* **2014**, *59*, 356–360. [CrossRef]

42. Comert, F.; Malanowski, A.J.; Azarikia, F.; Dubin, P.L. Coacervation and Precipitation in Polysaccharide-Protein Systems. *Soft Matter* **2016**, *12*, 4154–4161. [CrossRef] [PubMed]