5.1. Food model preparation
α-𝐷-Glucose (Crystals; Sigma Aldrich, St. Louis, USA) and WPI (Powder; carbohydrates or lipids as impurities < 10 ; Mullins Whey Inc., Mosinee, USA) were used to composite the food model. The glucose and WPI solutions (20, w/w) were prepared separately in deionized water and subsequently mixed to obtain solutions at different mass ratios (7:3, 1:1, 3:7, and 0:1; w/w). Further, the 5 mL of prepared solutions, loaded in preweighed 20 mL glass vials (semi-closed), were frozen in a still-air freezer (DW-HL240, Zhongkemeiling Co., Ltd., China) at − 20°C for 24 h, and then, were subsequently tempered at − 80°C for 3 h prior to lyophilization. The amorphous samples were obtained until the chamber pressure in a laboratory-used freeze dryer (10N/B, Scientz, Ningbo, China) was below 2 bar. It should note that lyophilizing glucose is extremely difficult due to its low Tg and high solubility nature21. In this study, therefore, the amorphous glucose was obtained experimentally via a modified quench-cooling approach reported by Simperler and Others23, in which approximately 1 g of glucose crystals was cooled to − 30°C and melted at 160°C, and then, quench-cooled again to − 30°C. Three units of each amorphous sample were stored in vacuum desiccators over desiccant (P2O5; Sigma-Aldrich, St. Louis, USA) to avoid water sorption and reach equilibrium at 30°C for further analysis.
5.2. Water sorption testing
The amorphous samples were weighed to monitor water sorption behavior as a function of time (24 h intervals until 120 h) over saturated solutions of LiCl, CH3COOK, MgCl2, K2CO3, Mg(NO3)2, NaNO2, and NaCl (Sigma Aldrich, St. Louis, Mo., USA) at respective aw of 0.11, 0.20, 0.31, 0.43, 0.53, 0.65, and 0.75 at storage temperature of 30°C, in vacuum desiccators. The Guggenheim-Anderson-de Boer (GAB) model (Eq. 1) was applied to fit the water sorption data of each sample, where the m and m0 referred to the weighted water content and monolayer water content; CGAB and KGAB were constants24.
$$\frac{\text{m}}{{\text{m}}_{\text{0}}}\text{=}\frac{{\text{C}}_{\text{GAB}}{\text{K}}_{\text{GAB}}{\text{a}}_{\text{w}}}{\text{(1}-{\text{K}}_{\text{GAB}}{\text{a}}_{\text{w}}\text{)(1}-{\text{K}}_{\text{GAB}}{\text{a}}_{\text{w}}\text{+}\text{Ck}{\text{a}}_{\text{w}}\text{)}}$$
1
5.3. Thermal Analysis
The thermal properties, including the onset-Tg value for each sample, were determined using a differential scanning calorimeter (DSC; Mettler-Toledo, Schwerzenbach, Switzerland). About 15 mg of prepared samples were transferred into a pre-weighed 50 mL aluminum pan and hermetically sealed before measurement. An empty punctured pan was used as a reference to minimize the systematic error caused by water vapor. Samples were scanned from − 20°C to over the Tg region at 5°C /min and then cooled at 10°C /min to the initial temperature. A second heating scan was run to well above the Tg at 5 ℃/min. The onset-Tg derived from second heating scans were recorded using STARe software (Version 8.10, Mettler-Toledo, Schwerzenbach, Switzerland). The Gordon-Taylor (GT) equation (Eq. 2) had proven to fit experimental onset-Tg data of glucose/WPI solid matrices, where w1 and w2 were the mass fractions of amorphous sample and water, Tg1 and Tg2 were their values, and kGT was a constant and its thermodynamic meaning discussed later.
$$\frac{{\text{W}}_{\text{2}}}{{\text{T}}_{\text{g}}-{\text{T}}_{\text{g2}}}\text{ = }{\text{k}}_{\text{GT}}{\text{w}}_{\text{1}}\left({\text{T}}_{\text{g}}-{\text{T}}_{\text{g1}}\right)$$
2
5.4. Dynamic-Mechanical Analysis
The mechanical properties of prepared samples were studied by using a dynamic-mechanical analyzer (DMA; Mettler-Toledo, Schwerzenbach, Switzerland). The loss modulus (E”) of materials as a function of temperature at different frequencies (0.5, 1, 3, 5, and 10 Hz) were determined in this study. Before starting an experiment, the instrument was balanced and set at zero to determine the zero-displacement position and return the force to the zero position. Approximately 100 mg samples of ground materials were spread on a titanium pocket-forming sheet. The length, width, and thickness (~ 2 mm) of the sample pocket between the clamps were measured. Samples were scanned from − 20 ℃ to over the Tg region with a cooling rate of 5 ℃/min and a heating rate of 2 ℃/min using the single cantilever bending mode to obtain E” values using DMA software (Version 1.43.00, Mettler-Toledo Schwerzenbach, Switzerland). During heating, the samples were analyzed for Tα values determined from the peak temperature of E” 25.
5.5. Molecular mobility determination
The temperature difference (Tα−Tg), at which relaxation times (τ) exceed time factors critical to the characteristics of the materials, was used to calculate S values, which can represent the extent of molecular mobility as noted above. The τ and the temperature of Tα above Tg were modeled and analyzed using the William-Landel-Ferry (WLF) equation (Eq. 3), where T, Tg, τ, τg, refers to the experiment temperature, onset-Tg, experimental a-relaxation time (oscillation frequency set in DMA measurement, τ = ½πf), and relaxation time in glass state (≈ 100 s). The WLF model constants C1 and C2 can be derived from a plot of 1/log(τ/τg) against 1/(Tα−Tg) using experimental τ with the assumption of τg = 100 s at the onset-Tg20. Moreover, the S value of the system is determined by Eq. (4), where C1, C2, and ds refers to the material-special WLF constants and the decrease in the number of logarithmic decades for τ (which ds = 4), respectively.
$$\text{Log}\left(\frac{\text{ τ}}{{\text{ τ}}_{\text{ }\text{g}}}\right)\text{= }\frac{-{\text{C}}_{\text{1}}\text{(}\text{T}-{\text{T}}_{\text{g}}\text{)}}{{\text{C}}_{\text{2}}\text{+}\text{(}\text{T}-{\text{T}}_{\text{g}}\text{)}}$$
3
$$\text{S= }\frac{{\text{d}}_{\text{s}}{\text{C}}_{\text{2}}}{-{\text{C}}_{\text{1}}-{\text{d}}_{\text{s}}}$$
4
Previous studies reported that the compositional dependent of S in non-crystalline sugar/protein solids could be represented by Eq. (5). Where w1 and w2 referred to the mass fractions of dry solids and water, ksp was a partition constant of molecular mobility, Sd1 and Sd2 represented the S value for anhydrous solids and amorphous water26.
$${\text{S}}_{\text{p}}\text{= }\frac{{\text{w}}_{\text{1}}{\text{S}}_{\text{d1}}\text{+}{\text{k}}_{\text{s}\text{p}}{\text{w}}_{\text{2}}{\text{S}}_{\text{d2}}}{{\text{w}}_{\text{1}}\text{+}{\text{k}}_{\text{s}}{\text{w}}_{\text{2}}}$$
5
5.6. Microbial response determination
Yeast Activation The D. hansenii (ACCC 20010; Xuanya Biotechnology Co. Ltd., China), isolated from the natural microflora, was chosen as a targeting microorganism because of its xerotolerant nature. D. hansenii was activated prior to inoculation on the bases of the method reported by Sharma and Others27. The lyophilized strains were dissolved and inoculated in a glassy tube containing 0.5 mL liquid YM agar (2.0% glucose, 0.5% yeast extract, 1.0% NaCl, 0.23% NaH2PO4, 0.5% (NH4)2SO4, and 1.8% agar; Sartorius Stedim Biotech, Globaltec Corp., Germany) at 30°C on a rotational shaker (200 rpm) for 24 h, and the successful activation achieved when the single colony was obtained.
Sample Inoculation The lyophilized sample were stored in vacuum desiccators over P2O5 as a desiccant to avoid water sorption and used UV light was for 24 h to eliminate environmental effects prior to inoculation. A tiny quantity of yeast-containing solution (~ 0.2 µl) was streaked on amorphous glucose/WPI solid matrices at mass rations of 1:0, 7:3, 1:1, 3:7, and 0:1. The inoculated glucose/WPI solid matrices were rehumidified over a saturated solution of NaCl, KCl, and K2SO4 (Sigma Aldrich, St. Louis, Mo., USA) at respective aw of 0.75, 0.83, and 0.92 aw at 30°C, respectively. Previous studies have verified that the low proportion of water incorporated with the inoculation did not raise the water content of the sample at aw > 0.758. Other equilibrated samples were not inoculated, but rather placed in a closed container over studied storage aw ranges and temperature for blank control. It is important to know that the whole inoculation was implemented in a clean bench, which can maintain a sterilized condition to avoid contamination from the surrounding ambient.
Growth Characterization Scanning electron microscopy (SEM; Phenom Pro, Phenom World. BV, Holland) was used to observe the morphology of microbes in glucose/WPI solid matrices at an acceleration voltage of 10 KV. The studied samples were coated using a gold-palladium alloy coater (Baltec Co., Man- chester, NH) and observed at 8000 × magnification. Built-in instrument software (SEM Center, JEOL, Japan) was used for image collection. The growth characters of D. hansenii in the glucose/WPI solid matrices were determined by an ATP fluorescence detector (Pi-102, Hygiena, USA) at 3 h intervals for 36 h and plotted the growth curve thereafter. The specific growth rate value of the D. hansenii in each system was determined by Eq. (6), and the cell doubling time of each system could be determined by Eq. (7). Where µ was the growth rate (h− 1), g was cell doubling time (h), N0 is the number of microbial cells at the beginning, Nt is the number of microbial cells at any time, and t is time (h).
$$\mu =\frac{(ln{N}_{t}-ln{N}_{0})}{t}$$
6
$$\text{g=}\frac{ln\left(2\right)}{\mu }$$
7
5.7. Statistical Analysis
The GAB isotherms and GT equation, S parameter, and microbial growth characteristics of triplicate measurements were analyzed and plotted in Microsoft Excel (2019, Microsoft, Inc., USA). The average values with a standard deviation of triplicate measurements were calculated. Additionally, the error bars and significance analysis were implemented in the confidence interval of 95% to represent the variability of data.