1. Introduction
Bioactive compounds used to enrich foods and beverages, such as vitamins, proteins or antioxidants, are highly sensitive to light, temperature and oxygen [
1], undergoing degradation reactions (e.g., oxidation or pigment destruction) during processing that decrease their bioavailability [
2]. To preserve them from degradation, those compounds can be encapsulated in suitable substrates according to the chosen functionalities [
3] or required delivery rates; the release of a bioactive compound within the human body could be fast (mouth release) or prolonged over time (digestive tract release). The choice of the encapsulation technique is then key to preserving the biocompound and creating a suitable carrier microstructure—e.g., highly porous matrices can enhance mass transfer, leading to faster release rates. This makes freeze-drying a convenient technique for encapsulation of active biocompounds [
4,
5], as it helps with keeping the original porous structures of products, and its low temperature conditions also contribute to minimising degradation reactions [
6,
7].
One of the most versatile substrates employed in bioprocessing applications (i.e., food, pharma and healthcare technologies) is gellan gum gel. This is a non-toxic, biocompatible and biodegradable polymer [
8] that has been extensively used as (i) a texturiser and gelling agent [
9] in food applications; (ii) to formulate oral, nasal and ophthalmic formulations [
10,
11]; and (iii) as a scaffold for tissue regeneration [
12,
13]. A recent study [
14] focused on the development of dried-gel structures from hydrocolloids has revealed the potential of gellan gum gels to be used as “controllable” carrier, showing that is possible to modulate the freeze dried-gel properties (i.e., target different microstructures and therefore different drying and rehydration kinetics) by modifying the pH of the initial gel solution.
To explore this promising path, this work focuses on the characterisation of freeze-dried gellan gum gels at different pHs as vitamin delivery systems. Freeze-drying kinetics, as well as release mechanisms and rates have been investigated using both empirical and modelling approaches. Gellan gum gels were prepared at different pHs (i.e., 5.2, 4 and 2.5) and then loaded with riboflavin (vitamin B2) before freeze-drying. Experimental drying curves were fitted to five common food drying models [
15] (i.e., Newton, Page, Henderson and Pabis, logarithmic and Wang and Singh), and the effects of different pHs on freeze-drying kinetics were assessed. In addition, Information Theory criteria (Akaike and Bayesian information criteria) were used to discriminate the models [
16] attending to their accuracy and complexity (i.e., number of parameters involved). The effect of pH on vitamin release (at ambient temperature) has been also studied using the classical Korsmeyer–Peppas model [
17,
18], and corresponding delivery mechanisms revealed. The findings of this work can help in the design of targeted freeze-dried gellan gum microstructures for the controlled delivery of active biocompounds in functional foods applications.
2. Materials and Methods
2.1. Materials
Low acyl (LA) gellan gum powder was provided by CPKelco (CPKelco, UK). Citric acid (purity 99%) and riboflavin (purity 98%) were supplied by Sigma-Aldrich (Sigma-Aldrich, Dorset, UK). All materials were used as received.
2.2. Preparation of Riboflavin-Loaded Gellan Gum Gels
Gellan gum powder was dissolved in distilled water at a concentration of 2% (
w/w) and stirred at 85 °C for 2 h to ensure complete mixing [
19]. This resulted in a gel solution at pH 5.2 (natural pH) (Seven compact Benchtop pH meter, Mettler Toledo UK). The gellan solutions at pH 4 and pH 2.5 were obtained by adding increasing quantities of citric acid at a 0.3 mol/L concentration [
14]. All the solutions were placed in cylindrical moulds (diameter = 22 mm) for gelation [
19]. Once the gels were formed, the moulds were cut into regular pieces (height = 15 mm,) and the resulting samples were soaked in a 2.7 × 10
−4 mol/L riboflavin solution for 18 h. Finally, the loaded gels were washed with distilled water and blotted with paper to remove the vitamin settled on the surface.
2.3. Freeze Drying
The riboflavin-loaded gels were frozen at −20 °C for 24h and then freeze-dried at increasing processing times, from 2 h up to 18 h, in a bench top freeze dryer (SCANVAC Coolsafe™, model 110-4, Lynge, Denmark) with condenser temperature of −110 °C and chamber pressure of 10 Pa [
14]. The experiments were performed in triplicates, and each batch of freeze-dried samples was weighted to measure the water content.
2.4. Water Activity Analysis
Water activity (a
w) of wet and dried gels was measured using an AquaLab
® dew point water activity meter (model 4 TE, Decagon Devices Inc., Pullman, WA, USA). The temperature-controlled sample chamber was set to 25 °C [
14,
15]. All analyses were carried out in triplicate.
2.5. Vitamin Release Experiments
Vitamin delivery analyses were performed using a UV–Vis spectrophotometer (Orion Aquamate, Thermo Scientific, UK) at 444 nm. The loaded gels were placed in stirred (250 rpm) distilled water (200 mL) at room temperature. To measure vitamin release, aliquots of 3 mL were withdrawn from the release medium, analysed with the spectrophotometer and poured back into the release medium. The vitamin content in the release medium was then expressed as normalised vitamin release (NVR), a dimensionless quantity defined as:
All analyses were carried out in triplicate.
2.6. Drying Kinetics
The kinetics of moisture loss during freeze-drying of the loaded gellan gum gels were described by five empirical models commonly employed to characterise drying kinetics in foods [
15,
16]: Newton, Page, Henderson and Pabis, logarithmic and Wang and Singh.
Table 1 lists them all alongside their expressions.
To fit the models to the experimental freeze-drying curves, the (dimensionless) moisture ratios (MRs) of the samples were calculated first from the measured water content data as follows [
15]:
where
X(t) is the moisture content on a dry basis for the different processing times (h),
X0 is the initial moisture content (
w/w d.b.) and
Xeq is the equilibrium moisture content (
w/w d.b). The equilibrium moisture content for the dehydrated gels was calculated using the GAB model with measured water activities and parameters for gellan gums presented in [
9]:
For all the models in
Table 2, the unknow parameters (parameters
aj and
ki, with
j = 1,2 and
i = 1, …, (6) were estimated using regression analysis. The error (e) between the experimental (
θ) and predicted (i.e., fitted) MR values (
) [
15],
was minimised for all the
i measurements that formed the experimental data set of size
N using a nonlinear least squares method (implemented in Matlab with tolerance 10
−10).
The goodness-of-fit of each fitted model was then assessed using three statistical measures that take into account the complexity (i.e., number of parameters,
p) of each model [
20]. These were:
- -
The adjusted
R2 [
15]:
where
R2 is the regression coefficient of determination.
- -
The corrected Akaike information criterion (
AICC) [
15,
21]:
where
AIC is the Akaike information criterion [
21,
22].
- -
The Bayesian information criterion (
BIC) [
20]:
where
L is the maximum log-likelihood of the estimated model.
The goodness of the fit (or the likelihood) can be increased by adding more parameters to the model. However, this will increase complexity and might result in overfitting (i.e., more parameters than can be estimated with the available data), all which is penalized with higher AICC and BIC values [
15,
20]. Therefore, the model with best performance will be the one with higher
and lower
AICC and
BIC values [
20].
2.7. Kinetics and Mechanisms of Vitamin Release
The Korsmeyer–Peppas model [
17] has been used to describe release kinetics and identify delivery mechanisms. It is a semi-empirical power law that relates the fractional release of vitamin/drug to the release time [
17,
18]:
where the
and
are the cumulative amounts of drug released at time
t (measured in hours
h) and infinite time, respectively; the constant
krel (in
units) relates to the structure and geometry of the delivery matrix (in this case, the freeze-dried gels); and the dimensionless exponent
nrel is the release mechanism indicator. For cylindrical substrates
nrel 0.45 defines Fickian mechanisms, while anomalous/non-Fickian delivery is described by n
rel > 0.45 [
17,
18]. Experimental release curves were fitted to Equation (8) using a nonlinear least squares method [
15], and parameters
krel and
nrel were estimated within 95% confidence intervals.
3. Results and Discussion
3.1. Effect of pH on Moisture Losses during Freeze-Drying
Figure 1 shows a comparison of the drying curves, in terms of the moisture ratio (MR) and freeze-drying processing times (h), for gellan gum gels at pH 2.5, 4 and 5.2 (natural) loaded with riboflavin. Gels at pH 2.5 exhibit the fastest drying rates, with most of the moisture content (MR
) removed during the first 6 h of the freeze-drying process. On the other hand, samples at pH 5.2 and 4 followed a very similar drying trend up to the first 4 h of processing. From this time onwards, the samples at pH 4 present a significantly slower drying rate; i.e., MR
at t = 8 h compared to MR
for pH 5.2 at the same time. The three samples were completely dried (i.e., free water totally removed) by the end of the freeze-drying experiments at t= 18 h, independent of their pHs.
The fastest drying rates observed for the gels at lower pHs can be explained by the effects of acidifying the gel solution. As reported by Cassanelli et al., [
14], the acidification step both enhances ice crystal nucleation and weakens the gel structure at pH values as low as pH 2.5. The combination of these two effects might favour a more interconnected pore structure—i.e., more nuclei will lead to more crystals that will find lower resistance in the weak gel structure to form a network. This could lead to faster drying rates and also affect the strength of the rehydrated structure.
3.2. Freeze-Drying Kinetics: Parameter Estimation and Model Discrimination
Estimated parameters for all the drying models considered in this work (i.e., Newton, Page, Henderson and Pabis, logarithmic and Wang and Singh) are listed in
Table 2, together with the RMSE (root mean square error) for each fitting and the results corresponding to the goodness-of-fit of each model. According to these results, the models that provide more accurate descriptions for the drying kinetics are the Page (Table Equation 1) and the Wang and Singh (Table Equation 5) models, both presenting
(in average) for all pH values.
For samples at modified pHs (i.e., pH 2.5 and 4), the Wang and Singh model presents the lowest RMSE (3.52 × 10
−4 for pH 2.5) and the highest
, while for the freeze-dried gels at pH 5.2, the Page model is the model that presents the best fitting (RMSE = 0.040 and
). This is in agreement with fittings reported in Cassanelli et al. [
19], which showed the Page model as the best option to describe the freeze-drying kinetics of non-loaded gellan gels at natural pH.
The goodness-of-fit for all models is illustrated in
Figure 2, where experimental values are plotted against predicted moisture ratios for each drying model at all pH studied. This graph also shows the accuracy of the Page and Wang models, for which most of the predicted points lie on the correlation line (see
Figure 2b,d).
When comparing models with similar accuracies, the AICC criterion constitutes the best measure to discriminate models, with more negative AICC values indicating better model performances. According to this, if the Newton and Page models were compared at pH 4—the pH at which both models show very similar RMSE, and BIC—the lower AICC (−21.39) of the Newton model would make it the preferred one. This criterion is also an indicator of the complexity (e.g., number of parameters) of the assessed models—the Newton model involves a single parameter (k1), compared to the two needed in the Page model (k2, n). On the other hand, the logarithmic model (Table Equation 4), with the highest number of parameters considered (p = 3), presents the less negative AICC values at each pH.
Effect of pH on the Drying Kinetic Parameters
The effect of pH on drying kinetic parameters has been determined by analysing the values of the constants in Newton (Table Equation 1) and Henderson and Pabis (Table Equation 3) models. These two models are derived from Newton’s cooling law and Fick’s Second law [
16], respectively, so their constants enclose physical meaning—as opposite to the Page and Wang and Singh models that are purely empirical [
16].
Parameters
k1 (Newton) and
k3 (Henderson) in
Table 2, both time constants (
h−1), characterise the drying rates of the system, while
a1 (Henderson) is a dimensionless parameter related to the shape and structural properties of the samples [
16].
For gellan gum gels at different pHs, both Newton and Henderson rate parameters (i.e.,
k1 and
k3, respectively) show very similar trends. The higher values (
k1 = 0.225 h
−1 and
k3 = 0.232 h
−1) correspond to samples at pH 2.5, indicating a faster dehydration process. On the other hand, rate constants for samples at pH 4 are the lowest (
k1 = 0.157 h
−1 and
k3 = 0.158 h
−1), which relates to the slower drying rate of these samples. This is in agreement with differences on moisture ratios (MR) at different pHs discussed in
Section 3.1.
The values of constant
a1 are again similar for samples at pH 2.5 and 5.2 (
a1 = 1.036 and
a1 = 1.043), which suggests no significant structural differences at those pH values. However, the value of
a1 for freeze-dried gels at pH 4 (
a1 = 1.006) suggests changes in microstructure that might be behind the different drying rates observed at this particular pH. These findings are in agreement with the mechanical properties (i.e., higher gel strength and Young’s modulus) reported in Cassanelli et al. [
14] for gellan gels at pH 4 before and after freeze-drying—“stronger” gels might make ice nucleation and growth difficult, and therefore affect the freeze-dried microstructures of the gels.
3.3. Riboflavin Release from Freeze-Dried Gellan Gels at Different pHs
Figure 3 presents experimental riboflavin release curves from freeze-dried gels prepared at different pHs, plotted as normalised vitamin released (NVR) across time. Data in this graph show significant differences in release times: freeze-dried gels at pH 4 completed the vitamin release in approximately 9.5 h; gels at natural pH (5.2) were fully unloaded after 6h, and total vitamin delivery took 3h for gels at pH 2.5. Samples at pH 2.5 presented a weak structure—in accordance with strength at fresh and freezing stages—that lead to breakage during the release experiments. This increased the surface area of the gels in contact with the release medium, which explains the shorter delivery times.
The observed differences in the riboflavin release times to the medium can be related to the different microstructures and mechanical properties of the gels. Both Norton et al. [
23] and Cassanelli et al. [
14] reported that freeze-dried gellan gum gels at pH 4 exhibited an aggregated and rigid structure. This can impede mass transfer within the gel, increasing the time needed to release the vitamin completely from the substrate and leading to longer delivery processes. A much lower level of aggregation and very weak structures were reported for freeze-dried gels at pH 2.5 [
14,
23], which is also in agreement with our experimental observations. According to the same authors, unloaded freeze-dried gels at natural pH (pH 5.2) exhibit intermediate levels of aggregation [
14,
24], explaining the also intermediate release times shown in
Figure 3.
3.4. Delivery Mechanisms at Different pHs
To estimate the value of the dimensionless parameter
nrel in Equation (8), and therefore determine the release mechanism governing riboflavin delivery, the portion of the release curves (
Figure 3) corresponding to the first 60% of the total released vitamin—i.e., release curve portions such that
—were fitted to the Korsmeyer–Peppas model [
17,
18]. Samples at pH 2.5 were not considered in this analysis, as they broke into several pieces during the release tests, leading to delivery conditions out of the scope of this work.
Table 3 lists the parameters
krel and
nrel (95% CI) estimated at pH 5.2 and pH 4 (see
Table 3 for parameter units). These results are discussed next.
3.4.1. Release from Freeze-Dried Gellan Gels at pH 5.2
Riboflavin delivery from gels at natural pH (pH 5.2) is characterised by a shape constant
krel = 0.509 diffusional coefficient
nrel = 0.131 (see
Table 3 for the corresponding confidence intervals). According to the classification given in [
21,
22], this indicates that the governing release mechanism is purely Fickian, as
nrel = 0.131 < 0.45, which is the limiting value of the diffusional coefficient for Fickian transport mechanisms. Therefore, we can define an apparent diffusion coefficient
Dapp (m
2/s) for samples at pH 5.2 using a short-time approximation of Fick’s Second law [
18]:
As the aspect ratio of the cylindrical gels is approximately in the order of 1, the predictive capabilities of the short-time solution include up to the 85% of the total vitamin release [
22]. Thus, values such that
in the release curve at pH 5.2 were fitted to Equation (9). This gave an estimated
Dapp = 1.325 × 10
−9 m
2/s with 95% CI defined by [1.086 × 10
−9, 1.564 × 10
−9] m
2/s.
Figure 4 presents the comparison between the experimental and predicted release curve using
Dapp, showing a good agreement between the modelled Fickian delivery mechanism and the experimental data—and thus confirming Fickian transport for riboflavin released from freeze-dried gellan gels at pH 5.2
3.4.2. Release from Freeze-Dried Gellan Gels at pH 4
As shown in
Table 3, the estimated shape constant and diffusional exponent at pH 4 were
krel = 0.287 and
nrel = 0.472, respectively. These estimates (i) confirm the structural difference of samples at pH 4 discussed in previous subsections, as the value of
krel at pH 4 is almost half of the corresponding to pH 5.2, and (ii) indicate an anomalous delivery mechanism from freeze-dried gellan matrices at pH 4, since we estimated
nrel > 0.45.
Anomalous mass transport can be defined as a mix between Fickian and non-Fickian mechanisms, for which the general form of the Korsmeyer–Peppas model presented in Equation (8) can be split into two contributions [
24]:
The first one (
) represents the Fickian part, while the second term (
) accounts for the relaxational contribution [
24], which is related to stresses and state transitions of polymeric matrices. The dimensionless diffusional exponent
mrel in Equation (10) can be determined from aspect ratio (i.e., height/diameter of the sample) correlations. For the ratio characterising the cylindrical samples used here (
1.5),
mrel = 0.43 [
24].
Estimated values for both
0.264 (0.240, 0.289) and
= 0.022 (0.007, 0.038) were then obtained by fitting Equation (10) to the experimental release curves at pH 4 for NMC <0.60 with 95% confidence intervals (values in the parenthesis). Units for
and
are (
. These estimates give an idea of the relevance of each contribution. For samples at pH 4,
, showing that the release of riboflavin from the freeze-dried gels at pH 4 is mostly Fickian. This is confirmed by the vitamin release percentages due to each contribution calculated as Peppas and Sahlin [
24]:
They are presented in
Figure 5a. Overall, the Fickian contribution is predominant along the release process, i.e., overall
, with values closer to 90% at the initial times of the delivery test, while relaxation effects are more important towards the end of the experiment, as the delivery of riboflavin is closer to completion.
Given the relevance of Fickian transport during the release process at pH 4, and with the estimated diffusional exponent
nrel so close to the Fickian limiting value of 0.45—confidence intervals for this parameter are actually cross-boarding this limit, i.e., (0.441, 0.504) as shown in
Table 3—a hypothetical pure Fickian riboflavin delivery at pH 4 has been also assessed. Following the procedure explained in
Section 3.4.1, an apparent diffusion coefficient D
app = 5.626 × 10
−10 m
2/s was estimated within a 95% confidence interval (5.409 × 10
−10, 5.842 × 10
−10) m
2/s. This estimate together with the short time approximately described in Equation (9) was used to obtain a predicted release curve, which is presented in
Figure 5b alongside the experimentally obtained curve. As the comparison reveals, the hypothetical pure Fickian mechanism describes the behaviour observed during the release tests well, and it could be used to predict riboflavin delivery—neglecting relaxation effects—from freeze-dried gellan gels at pH 4 with high accuracy.
4. Conclusions
This work demonstrates the potential to control biocompound release from freeze-dried gellan gum gels by modifying the pH of the substrate during gel formation, and prior to the encapsulation stage. As an exemplar of a relevant biocompound, riboflavin (i.e., vitamin B2) was used.
Freeze-drying kinetics, as well as release mechanisms, were experimentally investigated and modelled. Five different drying kinetics models were discriminated by accuracy and goodness-of-fit using statistical measures (i.e., RMSE, , AICC and BIC). For samples at natural pH (pH 5.2), the Page model provided the most accurate description of freeze-drying kinetics, while the Wang and Singh model predicted more accurately, the kinetics at acidified pH (i.e., 4 and 2.5).
Results revealed consistent differences in the behaviour of substrates at pH 4. Such differences reflect slower drying and release kinetics, as well as a different delivery mechanism—samples at natural pH (pH 5.2) exhibited Fickian transport, while acidified samples (pH 4) were characterised by an anomalous release mechanism, but with a predominantly Fickian contribution (80–90%).
Overall, this work shows the potential of modified pH freeze-dried gellan gum gel matrices for controlled riboflavin release, demonstrating that:
- (i)
These hydrogels could be used in different enriched food and/or beverage products.
- (ii)
Model-based approaches like the one presented here represent useful tools for the design of novel food formulations.
Author Contributions
Conceptualization, V.P., E.L.-Q.; Experimental work and writing, V.P.; Modelling work and writing, E.L.-Q. Supervision and resources, I.T.N. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by EPSRC grant no. EP/S023070/1.
Acknowledgments
Authors acknowledge financial support from the EPSRC Centre for Doctoral Training in Formulation Engineering (grant number EP/S023070/1) as well as support from the Centre’s Director, Prof. Fryer.
Conflicts of Interest
The authors declare no conflict of interest
References
- Sheraz, M.A.; Kazi, S.H.; Ahmed, S.; Anwar, Z.; Ahmad, I. Photo, thermal and chemical degradation of riboflavin. Beilstein J. Org. Chem. 2014, 10, 1999–2012. [Google Scholar] [CrossRef] [PubMed]
- Moreno, J.J. Innovative Processing Technologies for Foods with Bioactive Compounds; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
- De Marco, I.; Reverchon, E. Starch aerogel loaded with poorly water-soluble vitamins through supercritical CO2 adsorption. Chem. Eng. Res. Des. 2017, 119, 221–230. [Google Scholar] [CrossRef]
- Fang, Z.; Bhandari, B. Spray drying, freeze drying and related processes for food ingredient and nutraceutical encapsulation. In Encapsulation Technologies and Delivery Systems for Food Ingredients and Nutraceuticals; Woodhead Publishing: Cambridge, UK, 2012. [Google Scholar]
- Oikonomopoulou, V.P.; Krokida, M.K. Freeze drying and microwave freeze drying as encapsulation methods. Therm. Nontherm. Encapsul. Methods 2017, 39–65. [Google Scholar] [CrossRef]
- Karam, M.C.; Petit, J.; Zimmer, D.; Baudelaire Djantou, E.; Scher, J. Effects of drying and grinding in production of fruit and vegetable powders: A review. J. Food Eng. 2016, 188, 32–49. [Google Scholar] [CrossRef]
- Reyes, A.; Mahn, A.; Cares, V. Analysis of dried onions in a hybrid solar dryer, freeze dryer and tunnel dryer. Chem. Eng. Trans. 2015, 43, 139–144. [Google Scholar]
- Mahajan, H.S.; Gattani, S.G. Gellan Gum Based Microparticles of Metoclopromide Hydrochloride for Intranasal Delivery: Development and Evaluation. Chem. Pharm. Bull. 2009, 57, 388–392. [Google Scholar] [CrossRef] [Green Version]
- Abramovič, H.; Klofutar, C. Water adsorption isotherms of some gellan gum samples. J. Food Eng. 2006, 77, 514–520. [Google Scholar] [CrossRef]
- Agnihotri, S.A.; Jawalkar, S.S.; Aminabhavi, T.M. Controlled release of cephalexin through gellan gum beads: Effect of formulation parameters on entrapment efficiency, size, and drug release. Eur. J. Pharm. Biopharm. 2006, 63, 249–261. [Google Scholar] [CrossRef]
- Osmałek, T.; Froelich, A.; Tasarek, S. Application of gellan gum in pharmacy and medicine. Int. J. Pharm. 2014, 466, 328–340. [Google Scholar] [CrossRef]
- Grenier, J.; Duval, H.; Barou, F.; Lv, P.; David, B.; Letourneur, D. Mechanisms of pore formation in hydrogel scaffolds textured by freeze-drying. Acta Biomater. 2019, 94, 195–203. [Google Scholar] [CrossRef]
- Kim, D.; Thangavelu, M.; Cheolui, S.; Kim, H.S.; Choi, M.J.; Song, J.E.; Khang, G. Effect of different concentration of demineralized bone powder with gellan gum porous scaffold for the application of bone tissue regeneration. Int. J. Biol. Macromol. 2019, 134, 749–758. [Google Scholar] [CrossRef] [PubMed]
- Cassanelli, M.; Prosapio, V.; Norton, I.; Mills, T. Acidified/basified gellan gum gels: The role of the structure in drying/rehydration mechanisms. Food Hydrocoll. 2018, 82, 346–354. [Google Scholar] [CrossRef]
- Onwude, D.I.; Hashim, N.; Janius, R.B.; Nawi, N.M.; Abdan, K. Modeling the Thin-Layer Drying of Fruits and Vegetables: A Review. Compr. Rev. Food Sci. Food Saf. 2016, 15, 599–618. [Google Scholar] [CrossRef]
- Lopez-Quiroga, E.; Prosapio, V.; Fryer, P.J.; Norton, I.T.; Bakalis, S. Model discrimination for drying and rehydration kinetics of freeze-dried tomatoes. J. Food Process Eng. 2019, e13192. [Google Scholar] [CrossRef] [Green Version]
- Peppas, N.A.; Narasimhan, B. Mathematical models in drug delivery: How modeling has shaped the way we design new drug delivery systems. J. Control. Release 2014, 190, 75–81. [Google Scholar] [CrossRef]
- Ritger, P.L.; Peppas, N.A. A simple equation for description of solute release I. Fickian and non-fickian release from non-swellable devices in the form of slabs, spheres, cylinders or discs. J. Control. Release 1987, 5, 23–36. [Google Scholar] [CrossRef]
- Cassanelli, M.; Norton, I.; Mills, T. Role of gellan gum microstructure in freeze drying and rehydration mechanisms. Food Hydrocoll. 2018, 75, 51–61. [Google Scholar] [CrossRef]
- Spiess, A.-N.; Neumeyer, N. An evaluation of R 2 as an inadequate measure for nonlinear models in pharmacological and biochemical research: A Monte Carlo approach. BMC Pharmacol. 2010, 10, 6. [Google Scholar] [CrossRef] [Green Version]
- Moxon, T.E.; Nimmegeers, P.; Telen, D.; Fryer, P.J.; Van Impe, J.; Bakalis, S. Effect of chyme viscosity and nutrient feedback mechanism on gastric emptying. Chem. Eng. Sci. 2017, 171, 318–330. [Google Scholar] [CrossRef]
- Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
- Norton, A.B.; Cox, P.W.; Spyropoulos, F. Acid gelation of low acyl gellan gum relevant to self-structuring in the human stomach. Food Hydrocoll. 2011, 25, 1105–1111. [Google Scholar] [CrossRef]
- Peppas, N.A.; Sahlin, J.J. A simple equation for the description of solute release. III. Coupling of diffusion and relaxation. Int. J. Pharm. 1989, 57, 169–172. [Google Scholar] [CrossRef]
© 2020 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 (http://creativecommons.org/licenses/by/4.0/).