1. Introduction
Foams are dispersed systems of gas in a liquid matrix. Foams are found in nature and are also widely used in various technologies; therefore, they are of continuous interest to science since the works of Plateau in the 19th century [1,2]. Nevertheless, a general theory to describe foam behavior is still lacking. This is because the formation of foam and its further existence involves a multitude of phenomena at different length scales, among which interfacial phenomena play crucial roles [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. The creation of foam is a highly hydrodynamic process [15] and requires the presence of surface-active agents, which can adsorb at the foam interfaces reducing their free energy, and thus, reducing the global free energy of such an interface-dominated system. Once the foam is being created, its evolution towards decaying properties is determined by different destabilizing mechanisms (syneresis, coarsening, coalescence [4,5,6,7]), but the foam can reach a quasi-static state under hydrostatic equilibrium (dry foam).
Studies on foam stability are rather subjective, due to lack of a universal standard for generating foam—stabilities of foams made of a given solution can be specific to the conditions under which the foams are produced, and one should always keep this in mind when comparing results for similar systems [10]. Several characteristic parameters (bubble size and volume fraction distributions, the height of a generated foam column, characteristic time parameters of foam decay, etc.) have been used to describe foamability and foam stability [2,3,4,5,6,16,17]. One of the most popular parameters for describing the decay of a foam column is the foam half-lifetime t½, i.e., the time when the foam reaches half of its initial height. However, this parameter cannot recognize the individual impacts from different destabilization mechanisms. In this respect, Lunkenheimer et al. [16] introduced two characteristic parameters: time of deviation tdev and the time of transition ttr, which indicate the onset and the maximum rate of bubble coalescence in the course of foam decay, respectively; thus, allowing for discrimination (at least to some extent) between the interrelated actions of syneresis and coalescence. Moreover, it was found for surfactant stabilized foams that there is a linear relation between t½ and ttr, but it was also suggested that ttr allows a better description of the foam behavior [16].
Even if the evolution of the macroscopic properties of a given foam can be well described, revealing the reasons for the observed behavior, and thus, designing strategies for foam control require further knowledge on the microscopic level. In this respect, multi-scale approaches have been used, where macroscopic foaming properties are compared with the microscopic (nanoscopic) properties of lamellae (foam films) and interfaces [18,19,20,21,22,23,24,25,26,27] as foam’s building blocks. Regarding interfaces, on the focus are the traditional adsorption dynamics studies, as well as rheological [9,10,11,12,13,14,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46], electrical [27,29,35,47,48,49,50,51], and structural [30,50,51] interfacial characteristics—it is emphasized the crucial role of the rate of adsorption, surface viscoelasticity (surface elasticity is the resistance force to mechanical disturbances), interface charging state and interfacial molecular organization, respectively, to foamability, foam stability and foam rheology. Pronounced viscoelastic characteristics of interfaces have been recognized to enhance foaming properties [10,11,14,28,31]. Natural and technical foams are usually polydisperse, which results in a spectrum of capillary pressures in the Plateau borders and the nodes that are counterbalanced by the action of the disjoining pressure in the foam films. Direct correlations between the foam and the corresponding foam films can be, however, difficult. While experiments with single foam films are usually performed under well-defined conditions, the interconnected lamellae in a macroscopic foam can experience the simultaneous action of several processes, e.g., lamella bursting can cause so-called ‘collective effects’ inducing pressure shocks, interfacial shear and dilational perturbations, surface tension (excess) gradients, local enrichment of surface-active material, etc. [4,10,18,52,53,54].
Proteins are broadly used as stabilizers of dairy foams. The complexity and high sensitivity of globular proteins to environmental conditions make it difficult to describe their impact on the mechanisms of foam (de)stabilization [8,9,11,12,14,26]. Moreover, protein aggregation (reversible or irreversible) can have complex (enhancing/inhibiting) influence on interfacial and foaming properties as compared to those of the respective native protein [8,19,39,40,41,42,43,44,49,55]. The behavior of proteins at interfaces and in corresponding foam films and foams is governed by the interplay of the protein concentration and the solvent conditions (salt type and concentration, pH), and the unique effect of pH is of special interest.
The whey globular protein β-lactoglobulin (BLG) is frequently used as a model protein in foam and interfacial studies. Adsorption dynamics and surface rheology of native BLG layers at static water/air interfaces [8,25,35,36,37,38,56,57,58,59,60], as well as foam film thickness [25,61], have been previously studied as affected by pH. It was proven by vibrational sum-frequency generation (SFG) spectroscopy that the pH-dependent net charge of BLG in solution determines the magnitude of the electric field at the solution’s interface and thereof, a charge reversal can occur when the isoelectric point of BLG moieties (pI ≈ 5.1) is crossed [35]. Under the latter conditions maxima in the pH-dependencies of the surface pressure and surface excess [35,58,60], as well as of the surface shear elastic and viscous moduli [56] are observed, while at the same time foam film thickness drops to a minimum, and foam Newton black films (NBF) are obtained [61]. A qualitative inference of the adsorption process at the interface of a rising bubble in BLG solution was also earlier obtained; the results, however, showed that the pH-dependent bubble motion contravened the above-discussed trends and the bubble rising velocity increase with increasing pH in the studied range 3–7 [62].
Nevertheless, these trends are not necessarily valid for the pH-dependent behavior of protein foams. For example, different stabilities have been observed for foams obtained from protein solutions that have similar surface pressures and dilational viscoelasticity moduli and were found to produce foam films with comparable thicknesses [22,25,61]. For BLG, Lech et al. [25] reported poor foamability and foam stability at pH 3, while those were much better for pH 7; at the same time, at pH 5 foam half-lifetime increased in a jump-like manner with the increase of the BLG concentration—again, foam stability could not be explained by the interfacial characteristics and the film thickness investigated in the same study. Previous studies showed that the stability against the rupture of single foam films could be well correlated with foam stability [18,27,52,53], Therefore, in the present study, we aim at elucidating the reasons of the peculiar pH-dependent behavior of BLG foams by investigating the stability of BLG foam films. We present experimental results on foam films and foams obtained from solutions with two different BLG concentrations, but at a constant buffer concentration of 10 mM (which induces the formation of black foam films [61]). Foam stability is described by the characteristic times t½, tdev and ttr. Single black foam films were investigated in a porous plate cell to determine their stability in terms of the critical disjoining pressure of film rupture, Πcr, [18,52,53]. In this fourth part of the series, we discuss the obtained results in the context of our three previous papers [58,59,60] from this series, as well as of the above-mentioned literature results.
2. Materials and Methods 2.1. Chemicals and Solutions
Native β-lactoglobulin (molecular weight ≈18.3 kDa) has been isolated and purified from whey protein isolate (WPI) as starting material (Davisco, Le Seur, Minnesota, USA, protein content of 94% based on dry matter, comprising 18% α-lactalbumin, 44% β-lg genetic variant A, 30% β-lactoglobulin genetic variant B, and 8% minor proteins) by a very mild selective thermal aggregation (50 ± 1 °C/held for 120 min) under milieu conditions (addition of citric acid and trisodium citrate to adjust pH 3.4 and to capture the calcium ion stabilizing α-Lactalbumin). Under these conditions, the other whey proteins (mainly α-lactalbumin, BSA and immunoglobulins) were destabilized, while BLG was not affected. By subsequent microfiltration/ultrafiltration in diafiltration mode and spray drying, a powder of native BLG at high yield (>95%) was obtained. The procedure was described in detail in [63]. Analysis of native whey proteins was carried out by reversed-phase RP-HPLC, also described in detail in [63]. The powder sample used in the present study contained total protein ≈ 98.9% (of which BLG content > 99%, ratio of BLG-A/BLG-B ≈ 1.22), salts ≈ 0.7%, and traces of lactose (<0.05%). Aqueous solutions were prepared in Na2HPO4/Citric acid/Milli-Q buffers by keeping a constant buffer concentration of Cbuff = 10 mM and varying the BLG concentration (CBLG) and pH. All experiments were performed at room temperature of ≈ 22°C.
2.2. Foams
A dynamic foam analyzer DFA100 (Krüss, Hamburg, Germany) was used for foam stability measurements. It consists of a glass column, which is equipped with a light sensor array and a LED light panel. The whole glass column is scanned in real-time, and the intensity of the light detected by the sensors records the height of the foam HF and the height of the solution HL at the bottom of the glass column [2]. In a typical experiment, the instrument was programmed to sparge air through the porous frit (pore size 10–16 µm) at a flow rate of 0.15 L/min until the foam reached a height of HF = 160 mm. Subsequently, the gas flow was stopped automatically. Then, 50 mL of BLG solution was gently poured into the glass cylinder, which resulted in a bubble-free liquid height of about 40 mm, and therefore, the total height of the foam and the liquid was about 200 mm. HF and HL were automatically scanned by the device at regular intervals of 1 s (Figure 1a) and with heights’ resolution better than 0.5 mm (see Figure S2 in the Supplementary Materials). The time when the foam height reached a maximum value HFmax (while the liquid height reached a minimum value HLmin) was set as the time origin t0 = 0 from which the foam decay was further monitored. Thus, the following absolute changes in the foam’s and liquid’s heights were calculated as:
ΔHF(ti) = HFmax(t0) − HF(ti),
ΔHL(ti) = HL(ti) − HLmin(t0).
In this way, ΔHF(t) and ΔHL(t) data (Figure 1b) were collected for each of the studied foams.
Foam stability was analyzed in terms of the foam half-lifetime t½, as well as of the characteristic times tdev and ttr [2,6,16]. The latter two parameters define certain stages of the foam decay where different destabilization mechanisms (syneresis and coalescence) predominantly operate. According to this approach, the foam decay can be characterized by three main stages separated by the characteristic time parameters tdev (onset of coalescence) and ttr (maximum rate of coalescence).
The first stage lasts from t0 up to the time of deviation tdev, defined by Equation (3). In this initial period of the foam decay, rapid efflux of liquid from the fresh wet foam takes place; the bubbles get closer, which leads to the formation of foam films, Plateau borders and nodes. The syneresis dynamics is dominated by the gravity-driven drainage in Plateau borders and nodes, but at the same time, films start to drain as well, due to the capillary pressure in the Plateau borders. Indeed, the amount of liquid drained from such freshly formed (thick) films should be a small fraction of all the draining liquid. However, what is most important in this stage of foam decay is that the films presumably do not coalesce yet. Hence, the decrease in foam height is at this point equal to the increase of the height of the solution beneath the foam (see Figure 1b), that is:
ΔHF(t ≤ tdev) = ΔHL(t ≤ tdev)
After tdev, the foam system enters an intermediate period of its decay, where the films can be still draining or have reached their equilibrium thickness. However, in both cases, film rupture can take place that means coalescence of bubbles starts to operate along with the foam syneresis. The liquid volume fraction in the foam progressively decreases with time, that is the syneresis diminishes, and the process of coalescence becomes predominant. Under such conditions, the foam enters the final stage of its decay at the characteristic time ttr, where the rate of bubbles coalescence is at maximum, and the syneresis is negligible, i.e., the state of quasi-static dry foam is reached.
The evaluation of the parameter tdev can be done in several ways, satisfying Equation (3). It is directly obtained from a common plot of the kinetic dependences ΔHF(t) and ΔHL(t), as shown in Figure 1b, the two curves split at tdev. Values for tdev can also be obtained from the initial linear portion of the ΔHF(t) vs. ΔHL(t) dependence (Figure S2 in the Supplementary Materials). Also, the quantity
JFL(t) = ΔHF(t) − ΔHL(t)
can be conveniently plotted as a function of time, which indicates the position of tdev at the onset of JFL(t) > 0 (Figure 1c). The latter dependence is further used to evaluate the other characteristic parameter, time of transition ttr. This characteristic time is defined as coinciding with the inflection point of such function with the limits [2,16]:
JFL(t ≤ tdev) = 0 and JFL(t →∞) → const
To evaluate ttr we fitted each experimental J(t) curve for t > tdev with a cubic spline function and obtained the inflection point at the zero value or null point of its second derivative. The characteristic times tdev and ttr give valuable information about the evolution of the foam decay in terms of recognizing the interplay between syneresis and coalescence, and getting insights into the timescale and the magnitude of their individual and simultaneous actions. 2.3. Foam Films
Thin foam films were studied in a porous plate cell [27,64,65,66] with a radius of the hole in the frit film holder of R = 1 mm, enclosed in a sealed chamber (Figure 2). Further details on the cell and the setup have been reported elsewhere [66]. In such configuration, a foam film is formed in the center of a solution drop entrapped in the hole of the frit that is soaked with the solution under investigation. The disjoining pressure Π of an equilibrium film equals the balance of the acting pressures in the system and the geometry of the measuring cell [65,66]:
Π = Pa − Pr + Pc − Ph = ΔP + 2γ/rc − Δρgzc,
where ΔP is the difference between the pressure applied in the chamber Pa and the reference ambient pressure Pr (usually the atmospheric pressure); Pc and Ph are capillary and hydrostatic pressures, respectively; γ is the surface tension of the solution, from which the film is obtained; rc is the inner radius of the capillary; Δρ is the density difference between the solution and air; g = 9.81 m/s2; and zc is the height of solution in the capillary above the film.
Film stability was determined in terms of the critical disjoining pressure of film rupture Πcr for films in saturated vapor conditions, and with a constant radius equal to that of the film holder (1 mm), i.e., the film occupies the whole area of the hole in the film holder. The experimental pressure limit of the porous frit (which the film holder is made of), i.e., the highest disjoining pressure that could be measured, was about 13 kPa; hence, the measured critical pressures were limited to this value. Due to the fluctuation nature of proposed mechanisms of film rupture [64,65], the reported values for Πcr are averaged from at least five independent measurements, and the standard deviation is given by the error bars in Figure 7. Note that for the most stable films at CBLG = 100 µM and at pH 6 and pH 7, some of the films did not rupture up to the experimental pressure limit.
3. Results 3.1. Foam Films
The pH-dependent electrostatics of BLG covered interfaces [35] play a key role in the properties of the corresponding foam films, such as the foam film thickness [25,61], and the disjoining pressure that keeps the two interfaces of a foam film apart [61]. In the present work, we investigated foam films for two protein concentrations of CBLG = 10 µM and 100 µM but at constant buffer concentration of Cbuff = 10 mM, while the solution pH was adjusted at values of 3, 4, 5, 6 and 7. At 10 mM buffer and in the studied pH range, black films were obtained (Figures 4–6) for all systems studied. These films can be either ‘common black films’ (CBF) or ‘Newton black films’ (NBF) [61]. Foam films of the NBF type are readily obtainable at pH 5 (which is close to the isoelectric point pI) and at low disjoining pressures, while at pH values away from pI, CBF are obtained. These CBF can, however, turn into NBF when the disjoining pressure Π increases, a fact that has been previously revealed by disjoining pressure isotherm measurements [61].
The experiments were performed as follows. In the experimental protocol ‘A’, a foam film was initially formed slowly by a gradual increase of the applied pressure Pa in the sealed chamber. Once the film appeared in the center of the solution drop, Pa was adjusted to yield a disjoining pressure of Π = 100 Pa. For the different solutions in our study, we obtained film radii in the range of 0.2–0.3 mm (Figure 3, Figure 4 and Figure 5). Once the pressure was adjusted, the foam films were left to thin and ideally reach a homogeneous pattern, which would be consistent with films having the same thickness independent of the lateral position within the film. However, this was not always possible due to film dimpling [67], and very slow drainage, typical for protein foam films [19,20,61,67,68,69,70,71], especially at such a low disjoining pressure. Therefore, in many cases, inhomogeneous thick parts remained immobile within the films (Figure 3, Figure 4, Figure 5 and Figure 6). After that, Pa was increased in small steps of 100 Pa at a rate of 100 Pa/min until the film occupied the complete available area. Under the experimental conditions, this was usually the case at about Π = 1 kPa. Note that at pH 3 the foam films for both studied BLG concentrations of 10 µM and 100 µM were the least stable and ruptured already during this stage, without filling the available area in the film holder (Figure 4). For the other more stable systems, Pa was gradually increased in a continuous mode at a rate of 1 kPa/min until the films ruptured. The disjoining pressure at film rupture is defined as Πcr. Several experiments were performed by a second protocol ‘B’, where the initial Pa was set higher, and the films were rapidly formed at Π ≥ 1 kPa (Figure 6). After that, Pa was gradually increased as in protocol ‘A’. In a few experiments, films were left to relax at a given Π in order to observe the film drainage (see in Figure 3 and Figure 6, and Figure S6 in the Supplementary Materials).
Figure 7 shows the variation of the critical pressure of rupture Πcr of BLG black foam films obtained from solutions with various pH and at the two BLG concentrations studied. The reported values for Πcr are averaged from at least five independent measurements done mostly by protocol ‘A’, but also including some results obtained by protocol ‘B’. The standard deviation varied from ±0.5 kPa to ±1.7 kPa, and we could not clearly distinguish certain effects of the different experimental protocols on the film stability. However, rapidly formed films exhibited much more irregular patterns (thickness inhomogeneities) and drained more slowly, which not always led to increased stability. Hence, we conclude that the way of film creation could have a relatively weak impact on the overall film stability.
The lowest film stability was observed for pH 3 and pH 5 at CBLG = 10 µM, where the films ruptured soon after their formation at Π = 100–120 Pa. At pH 3 and CBLG = 100 µM, the film stability was only slightly higher (Πcr < 1 kPa) in full agreement with the disjoining pressure isotherm for such films reported earlier [61]. For pH 4, 6, and 7 at CBLG = 10 µM, the film stability is higher (up to Πcr = 1.6 kPa). At CBLG = 100 µM, Πcr continuously increased with the pH and leveled off at pH 6 and pH 7, for which some of the most stable films did not rupture up to the highest measurable disjoining pressure of about 13 kPa (experimental limit).
At pH 5, where NBF type of foam films are directly formed at low disjoining pressures, we studied film stability also for two intermediate BLG concentrations, namely, 30 µM and 70 µM. The Πcr(CBLG) dependence is presented in Figure 7b, and photos of such films are shown in Figure 5 and Figure S6 in the Supplementary Materials. The increase of CBLG from 10 µM to 30 µM led to much higher stability of the obtained NBFs, which ruptured at disjoining pressures in the range 5.5–7.1 kPa within the current experimental reproducibility. Further increase of CBLG to 70 µM and 100 µM led to a gradual increase of Πcr. It is obvious that the stabilization of the NBF structure (as studied at pH 5 and Cbuff = 10 mM) occurs at protein concentrations in the range 10 µM < CBLG < 30 µM. A critical surfactant concentration, denoted by Ce as the equilibrium concentration for the formation of a stable NBF, has been introduced earlier for foam and emulsion NBFs stabilized by low-molecular-weight surfactant as studied in a tube cell at relatively low disjoining pressures (usually below 140 Pa); for example, Ce = 1.1 mM for SDS solutions at 500 mM NaCl [64,65]. Similar critical concentrations were also used for the characterization of the stability of black foam films of some proteins at pH values close to the isoelectric point; for example 24 µM for lysozyme (pH 10.5–11), 9 µM for human serum albumin (pH 4.9) and 15 µM for bovine serum albumin (pH 4.9) [67]. The present results for BLG show that such critical BLG concentration is of the same order of magnitude. Other characteristic concentrations introduced for surfactant films are, for example, the minimum surfactant concentration required for the formation of black spots in the films, denoted by Cbl, and the concentration, at which black films are always obtained, denoted by Ct [64]. We studied the probability of formation of (NBF) black spots in BLG foam films (pH 5) obtained in a tube cell after aging the solution surfaces for 10 min; the experimental conditions are explained in the Supporting Information, and the results are presented in Figure 7b (also in Figure S6 in the Supplementary Materials). At CBLG = 0.1 µM black spots were not observed at all until film rupture, while at CBLG = 1 µM every film ruptured through the formation of black spots.
In summary, we found the characteristic concentrations of BLG for formation and stability of foam NBFs under the studied solvent conditions (pH 5, Cbuff = 10 mM): 0.1 µM < Cbl < 0.2 µM, Ct = 1 µM and 10 µM < Ce < 30 µM. At the moment, we do not have information in which way small variations in pH and/or in salt concentration would affect these characteristic parameters. However, it was shown earlier that the addition of NaCl can improve the stability of BLG black films (pH 5.3), presumably, due to salting-in effects [71].
3.2. Foams
Figure 8 presents results for the t½(pH), tdev(pH) and ttr(pH) dependencies for foams obtained from BLG solutions at the two studied BLG concentrations of 10 µM and 100 µM. It has been found earlier for surfactant foams that there is a linear relation between t½ and ttr [16]. This is also the case for the present results for BLG foams. Combining the data in Figure 8b and Figure 9b results in a linear t½ ~0.75ttr dependence, which is presented in Figure S3 in the Supplementary Materials. Hence, in our further discussion, we do not consider these parameters separately and will focus on the results for ttr. In contrast, the corresponding values for tdev are two orders of magnitude lower. The data for the two dependences tdev(pH) and ttr(pH) for the two BLG concentrations follow rather similar trends with an offset to longer characteristic times for the higher CBLG; however, some peculiarities were observed, and they are discussed in Section 4. Obviously, the general foam stability is reflected by all of the three used parameters in a similar manner, hence in the following, we use the term ‘foam stability’ as a generalization for all foam results and where appropriate we discuss the individual types of foam (characteristic time) data.
At CBLG = 10 µM, there is a slight increase in foam stability from pH 3 to 4, and then it decreases to a minimum at pH 5 (isoelectric point, pI ≈ 5.1). This comparatively very low foam stability is in correspondence with the foam film data in Figure 7 and with previous results [61], which all show the formation of black films with extremely low stability at CBLG = 10 µM and pH → pI. Further increase of pH to 6 and 7 leads to a significantly increased foam stability. At CBLG = 100 µM, a general increase in foam stability was observed and as mentioned above the trend remains virtually the same. A very similar trend of the t½(pH) dependence was reported for foams obtained from solutions with CBLG ≈ 27 µM (intermediate between those used in our experiments) and comparable buffer concentrations [25]. Interestingly, those authors also found a decrease in the half-lifetime of foams at pH 3 and pH 7 with a further increase of the BLG concentration to about 270 µM, and reported the existence of a maximum in the t½(CBLG) dependence for these pH values, while at pH 5, t½ monotonically increases with the protein concentration as also revealed by our data in Figure 9. The latter trend is in agreement with the increased stability of the NBFs (at pH 5) with increasing CBLG (Figure 7b), while the observed decrease in foam stability with increase of the BLG concentration [25] does not correlate with the increase in the foam film stability found in our work here.
In the special case of foams at pH 5, the effect of increasing BLG concentration from 10 µM to 100 µM is the strongest and changes the shape of the tdev(pH) dependence (Figure 8a). Figure 9 presents measured tdev(CBLG)pH5, t½(CBLG)pH5 and ttr(CBLG)pH5 dependencies. In all cases, a gradual increase in the foam stability was observed with increasing CBLG.
To better illustrate the interrelations between the stabilities of single foam films and macroscopic foams, we plotted in Figure 10 the measured characteristic times tdev(pH, CBLG) and ttr(pH, CBLG) of the foam decay taken from Figure 8 and Figure 9 against the corresponding film critical pressures Πcr(pH, CBLG) taken from Figure 7. We also estimated the approximate values of the capillary pressuresPc,F in the studied foams on the basis of literature data for the Sauter mean bubble radii r32 in foams obtained in BLG solutions at different CBLG and pH [25,37]. The calculations, presented in the SI, yielded a spectrum of estimated capillary pressures between minimumPc,Fmin≈ 170 Pa and maximumPc,Fmax ≈ 1.7 kPa values (indicated as dashed lines in Figure 10) for all studied foams.
4. Discussion
The BLG concentrations of 10 µM and 100 µM used in the experiments were chosen in order to cover a range of foam half-lifetimes from few minutes to hours depending on the solution pH. Half-lifetime of foams, as well as some interfacial characteristics of BLG solutions, were previously studied by Lech et al. [25] as affected by pH at very similar solvent conditions and at several CBLG—therefore, a helpful comparison with our results can be made. In the following, we discuss the decay of a BLG foam column on the basis of the stability of corresponding foam films. Important dynamical, electrical, rheological and structural characteristics of foam interfaces and foam films that affect foam properties are hardly accessible by probing foams in situ, therefore, the usual model systems used are a single isolated interface or film. Although those are compromise options, they are practically the most commonly used. The ‘compromise’ comes from the facts that, on the one hand, the foam is a dynamic (eventually quasi-static at hydrostatic equilibrium) system and interfaces and films can be to a different extent close to equilibrium. On the other hand, a single interface or film is usually investigated under well-defined conditions. In contrast, the lamellae and interfaces in the foam can experience the simultaneous action of several factors (so-called ‘collective effects’ [4,10,18,54]), which cannot be reflected in single film experiments. However, good correlations between the foam film and foam stabilities were found in the present study, and we discuss these in more detail in the following.
Firstly, we briefly consider the pH-dependent adsorption kinetics of BLG at water-air interfaces. It is commonly accepted that fast surfactant adsorption kinetics is the most important prerequisite for foam generation [10,13,31]. Adsorption of surface-active substances and formation of dynamic adsorption layers (DAL) take place during the formation and growth of bubbles, as well as during their rising in solution [10,23,72]. Figure S7 in the Supplementary Materials shows that the ‘short-time’ adsorption kinetics (0.1–80 s as relevant to the time-scale of foam formation) in the studied solutions is relatively fast for pH values in the range 4–7, while for pH 3 it is slower (however significantly improved at the higher CBLG). Rising bubbles in such solutions are supposed to be fully retarded, since the rising bubble velocity technique becomes insensitive for CBLG > 5 µM BLG in the pH range 3–7 [62]. This means that a DAL of BLG is readily formed at the bubble surface in any of the solutions investigated in the present study, however, the exact adsorbed amounts are unknown. Further, longer-time adsorption experiments revealed the tendency that pH-dependent interfacial characteristics of proteins usually exhibit extremum around pI, where an interfacial charge reversal occurs [35]; for example for BLG the surface pressure (Figure S8 in the Supplementary Materials) and the surface excess [35,58,60], as well as the surface shear elastic and viscous moduli [56], are maximal at pH 5 (pH range 3–7). All these can be related to the formation of thick heterogeneous BLG adsorption layers (pH 5) with the occurrence of a loose secondary sublayer of disordered BLG dimers adjacent to the dense primary monolayer as resolved by neutron reflectometry experiments [60].
It is somewhat intriguing that the above considered pH-dependent interfacial properties of BLG contravene the trend in foam stability, which we will discuss in more detail based on the foam film results. Inspecting Figure 10 one can see that the tdev(Πcr) and ttr(Πcr) data can be separated into three groups according to the estimated range of foam capillary pressuresPc,F; a fourth group is assigned to the special case of films and foams at pH 5:
-
Group 1: Πcr <Pc,Fmin . This condition is valid for the foams and films at pH 3 and pH 5 and at CBLG = 10 µM, and can explain the lack of measurable tdev values, where it is possible that the onset of bubble coalescence occurs already in the period of foam formation. These two foams also exhibit the shortest t½ and ttr (in the order of some minutes). For pH 5, this is somehow expected, having in mind that CBLG = 10 µM < Ce and unstable NBF spots emerge in films within seconds (Figure S6 in the Supplementary Materials). For pH 3, the most probable reason is low protein surface excess as suggested by the long induction time in the dynamic surface pressure measurements (Figures S6 and S7 in SI).
-
Group 2:Pc,Fmin< Πcr <Pc,Fmax . This condition is valid for the foams and films at pH 4, 6, 7 and CBLG = 10 µM, but also for those at pH 3 and CBLG = 100 µM, and results in measurable values for tdev. We should mention here that the time scale of tdev is much shorter than the typical timescale of protein foam film drainage [19,20,61,67,68,69,70,71] (Figure 3 and Figure 6). Therefore, bubble coalescence seemingly takes place for films that are far from equilibrium, which is highly probable for bubbles in the top foam layer of a freshly produced foam [73]. Despite the higher protein concentration at pH 3, foam and film stabilities are still relatively low, and the reason for that is not clear yet [25]. Nevertheless, the good film-foam correlations in groups 1 and 2 demonstrate that single foam film stability is an indicator of foam stability. However, this is not the case for the foams in group 3.
-
Group 3: Πcr >Pc,Fmax . This condition is valid for the foams and films at pH 4–7 and CBLG = 100 µM (as well as for intermediate CBLG at pH 5) and suggests higher resistance of the foam bubbles to coalescence. Indeed, tdev and ttr increase, which agrees with the results in ref. [17]. This fact should be related mainly to the increased protein surface excess at the foam interfaces (as known from adsorption studies [58,60]) that leads to at least two important contributions: (i) stronger immobilization of interfaces, which decelerates the liquid efflux [10,28]; and (ii) stabilization of foam films against rupture [61,71], i.e., impact on the processes of syneresis and coalescence, respectively. The longest tdev at pH 5, 6, 7 are attributable to the low calculated rates of liquid efflux vdev (Equations (S1)) and Figure S4 in the Supplementary Materials) in the initial stage of foam decay (t0–tdev); the longest ttr at pH 6 and pH 7 correspond to the most stable foam films. Concerning the latter, we should note that despite the comparable film stability at pH 6 and pH 7, foam stability diminishes at pH 7.
-
Group 4: Πcr >Pc,Fmax ; pH 5. At CBLG > Ce, Πcr continuously increases with increasing CBLG resulting in increasing foam stability [25]. At CBLG = 100 µM, Πcr for pH 5 and pH 4 are comparable (within the experimental error), but at pH 5 the measured tdev is almost two-fold longer, whereas ttr is shorter. While the latter cannot be straightforwardly explained, the former can be attributed to the much lower vdev (Figure S4 in the Supplementary Materials) as a result of the high surface viscosities either in shear [56] or in dilation [59] (Figure S8 in the Supplementary Materials) as measured at early times of adsorption, as well as from a cork effect in the Plateau borders [24], and entrapment of BLG aggregates in the films (Figure 5). Hence, the foam at pH 5 drains slower than that at pH 4, but the rate of bubble coalescence is higher and a possible explanation for that should be due to the influence of aggregates. It has been shown that foams obtained from mixed dispersions of native and aggregated BLG can be either less or more stable than foams from pure native BLG solutions, and the direction of change of foam stability depends on both the aggregate size and the ratio of native/aggregated entities [24].
The pressure condition Πcr >Pc,Fmax (groups 3 and 4) opens the question about the driving force for the observed foam destruction. After tdev, the foam decay is determined not only by syneresis, but also by progressive coalescence of bubbles in the top foam layer [10,73] with a maximum rate at ttr. At the same time, bubble disproportionation (foam coarsening) is expected to take place in the foam. It has been shown that the coarsening rate in protein foams is lower than that in surfactant foams [19], which can be reasonably accounted for with lower permeability coefficients for protein films [71]. Coarsening generally leads to an overall increase in the bubble size. Hence, with time, the newly appeared top foam layers consist of bubbles with progressively increasing sizes, which causes a decrease ofPc,F , but at the same time, the probability for film rupture increases, due to increased film area [4,53]. A possible driving force for rupture of external foam films in the top foam layer is the intensive evaporation of water, since these films are exposed to the ambient atmosphere. Furthermore, the action of “collective effects” in foam has been recognized as a strong destabilizing factor [4,10,18].
To further our discussion, particularly towards the effect of pH, we should once again emphasize the relatively good correlation between the stabilities of foams and single foam films. The major feature of the presented results is the lower film and foam stabilities at pH < pI as compared to those at pH > pI, with the extreme case of pH 3, where foam stability could not be much improved even at a BLG concentration as high as ~1100 µM [25]. A comparison of adsorption and dilational rheology results shows rather similar behaviors for BLG solutions at pH 4, 6, 7 (Figures S6 and S7 in SI) [25,61]. From the viewpoint of surface forces, film thickness measurements did not give any indications to support such a trend in the stability of BLG films [25,61]. Furthermore, disjoining pressure isotherm study revealed that BLG foam films at pH 3.5 and pH 7 have similar (intermediate) stability at a low buffer concentration of Cbuff = 3 mM, but at Cbuff = 10 mM, the stability of the obtained black films increases for pH 7 and decreases for pH 3.5 [25,61]. It is obvious that neither surface properties nor electrostatic interactions could satisfactorily explain the observed pH-dependent film behavior.
The rupture of symmetric thin liquid film stabilized by a low-molecular-weight surfactant can happen through different mechanisms depending on the film thickness, e.g., the fluctuation waves approach for the case of thick films (>30 nm) or the hole-nucleation approach for the case of amphiphile bilayers (NBF and bilayer lipid membrane) [4,64,65]. However, these theories have been developed based on experiments with films at constant low disjoining pressure, and this is not clear yet whether they are applicable to describe the rupture of ‘stable’ films at the critical pressure. So far, to the best of our knowledge, there is not any theory, which considers the case of protein films. Having in mind the almost irreversible adsorption of proteins at liquid interfaces and the pronounced viscoelastic behavior of the corresponding interfacial layers, the phenomenon of rupture of protein black films could obey other mechanisms. Very recently, it was shown for mixed polymer/surfactant aqueous systems that resolving the molecular structure of the adsorbed interfacial layers is a vital step for revealing the surface forces operative in the corresponding foam films, and thus, the film stabilization mechanism [74]. Therefore, looking deeper into the structure of protein layers would give more insight into the stability of protein black films.
In solutions with pH near pI, BLG exists exclusively as a dimer, which dissociates into monomers by moving the pH away from pI; in BLG solutions at pH 3 the monomer is favored, while at pH 4, 6 and 7 the dimer is favored [75,76,77,78,79,80]. In the pH range 3–7, the pH causes relatively weak effects on the secondary and tertiary structures of BLG in solution bulk, and its general conformation is believed to be retained to a certain extent. However, two major structure transitions have been identified [76]: (i) the pH 2.5–4 transition (slight change in secondary and tertiary structures) accompanied by the formation of dimer; and (ii) the pH 4.5–6 transition (slight change in tertiary structure). Adsorption of BLG monomers and dimers at the water/air interface of solutions at pH 3, 5 and 7 were examined via surface excess and layer thickness measurements by neutron reflectometry [60]. At pH 3 and pH 7 the formation of dense monolayers with similar thicknesses (~1.1 nm) was found that consist of flattened BLG globules independent of the fractions of dimers and monomers in the bulk solution. However, any effects of the BLG secondary and/or tertiary structures were not recognized. In another study, SFG spectroscopy experiments revealed that water molecules at the vicinity of BLG covered interfaces form more ordered networks at pH < pI than at pH > pI [35], which could originate from structural specificities of adsorbed BLG entities. It is somewhat intriguing that the trend of the foam film stability with increasing pH (3–7) somehow reflects the pH-induced transitions in the molecular structure of BLG in bulk, as well as the different structuring of interfacial water. At the same time, the “monomer-dimer concept” could not give a satisfactory explanation of these findings.
5. Concluding Comments and Outlook
A good correlation between the stabilities of black foam films and foams obtained for BLG solutions at various pH values (3–7) was found that can be very useful for predicting the stability of protein foams and designing strategies for foam control in applications. This also shows that precisely the critical disjoining pressure of film rupture rather than the film thickness seems to a great extent to be determinant for foam stability. More profound quantitative analysis on the protein foams of interest requires the application of additional methods, for example, the foam pressure drop technique, which allows control over the foam capillary pressure [4,18,52,53].
The comparison of the obtained results to key BLG interfacial properties appeared intricate, thereby, robust interrelations could not be identified. Adsorption and surface viscoelasticity characteristics, which exhibit maximum values around pI [25,35,56,58,59,60], cannot explain the trend of increasing film stability with increasing pH. The most striking example is the remarkably low film and foam stabilities at pH 3, even at very high BLG concentrations [25]. Resolving the molecular structure of interfacial protein layers must be a promising strategy to understand the phenomenon of rupture of protein black films at the critical disjoining pressure. Since the thickness of BLG monolayers is almost independent of pH ([60] and references therein), one should investigate in greater detail eventual pH-induced changes in the molecular (tertiary) structure of adsorbed BLG (dimer/monomer) entities in order to reveal the mechanism of film rupture. Some insights into possible mechanisms could be gained by attempting to apply the hole-nucleation theory developed for the description of the mechanism of rupture of low-molecular-weight amphiphilic bilayers [4,64,65] to protein bilayers [61,81] obtained at constant low disjoining pressure (such data are shown in Figure 7b). However, the comparison of experimental data on protein bilayers with that theory can be complicated, due to the pronounced viscoelastic behavior of interfacial protein layers.
Supplementary Materials
The following are available online at https://www.mdpi.com/2075-163X/10/7/636/s1, Figure S1: titration curves and zeta potential, Figure S2: ΔHF vs. ΔHL plot for an exemplary foam, Figure S3: t½ vs ttr plot of the data in Figure 8b and Figure 9b, Figure S4: vdev vs. pH, Figure S5: scheme of the Scheludko-Exerowa tube cell, Figure S6: Wbs vs. CBLG (pH 5), Figure S7: kinetic dependencies of the surface pressure, Figure S8: kinetic dependencies of the surface pressure and dilational viscoelasticity modulus.
Author Contributions
Conceptualization, G.G.G., V.U. and R.M.; methodology, G.G.G., V.U. and I.R.; software, G.G.G. and V.U.; investigation, G.G.G., V.U., I.R. and R.M.; resources, U.K. and B.B.; writing-original draft preparation, G.G.G. and V.U.; writing-review and editing, G.G.G., V.U., C.G.-D., D.Z.G., M.L., U.K., R.M. and B.B.; visualization, G.G.G. and V.U.; supervision, G.G.G. and R.M.; funding acquisition, C.G.-D., D.Z.G., M.L., U.K., R.M. and B.B. All authors have read and agreed to the published version of the manuscript.
Funding
This work was financially supported by Nestlé Research Center, Lausanne and by a DFG-AiF cluster project on "Protein Foams" Mi418/20-1 (DFG-199448917). G.G.G. and B.B. gratefully acknowledge financial support by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement No 638278).
Acknowledgments
We sincerely thank Daniel Baier and Cornelia Rauh from Technical University of Berlin, Germany for allowing us to use the lab facilities and the Krüss®, DFA-100 instrument.
Conflicts of Interest
The authors declare no conflict of interest.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The complexity and high sensitivity of proteins to environmental factors give rise to a multitude of variables, which affect the stabilization mechanisms in protein foams. Interfacial and foaming properties of proteins have been widely studied, but the reported unique effect of pH, which can be of great interest to applications, has been investigated to a lesser extent. In this paper, we focus on the impact of pH on the stability of black foam films and corresponding foams obtained from solutions of a model globular protein—the whey β-lactoglobulin (BLG). Foam stability was analyzed utilizing three characteristic parameters (deviation time, transition time and half-lifetime) for monitoring the foam decay, while foam film stability was measured in terms of the critical disjoining pressure of film rupture. We attempt to explain correlations between the macroscopic properties of a foam system and those of its major building blocks (foam films and interfaces), and thus, to identify structure-property relationships in foam. Good correlations were found between the stabilities of black foam films and foams, while relations to the properties of adsorption layers appeared to be intricate. That is because pH-dependent interfacial properties of proteins usually exhibit an extremum around the isoelectric point (pI), but the stability of BLG foam films increases with increasing pH (3–7), which is well reflected in the foam stability. We discuss the possible reasons behind these intriguingly different behaviors on the basis of pH-induced changes in the molecular properties of BLG, which seem to be determining the mechanism of film rupture at the critical disjoining pressure.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer