Abstract

Silk fibroin (SF) and sodium alginate (SA) are natural polymers with interesting properties to produce biomaterials. Blends of these polymers form a complex protein–polysaccharide system where phase separation can be observed. Therefore, the thermodynamic analysis of this system is important to understand the interaction between the polymers and predict the final state and composition of the phases found in these blends. This study explored blends with a different initial composition of SF, SA, and water at 25°C and neutral pH. The influence of the proportion between mannuronic and guluronic acids on SA composition was investigated. After phase separation, two phases were identified, and the equilibrium data were fitted on three different thermodynamic models: Flory–Huggins, non-random two-liquids, and universal quasichemical. Cohn equation was also used to investigate the potential of SA to precipitate the SF in solution. The results show that the proportion of mannuronic and guluronic acids on the SA can significantly influence on equilibrium data and on the SF/SA interaction parameter, hence, becoming a variability factor if this parameter is not under control in formulations.

1. Introduction

The systems formed by protein and polysaccharides are complex, and phase separation is often observed. Silk fibroin (SF, a fibrous protein extracted from silkworm cocoons [1]) and sodium alginate (SA, a polysaccharide extracted from brown algae [25]) blends have been studied over the years due to their intrinsic characteristics, such as biodegradability, biocompatibility, and low toxicity [69]. However, it forms a complex system where partial miscibility is observed [6, 10]. When in solution, phase separation is observed for SF–SA blends, and two phases are formed, one which is liquid and rich in SA and another which is solid-like, rich in SF, as shown in Figure S1 of the supplementary materials [10, 11].

The thermodynamic understanding of SF–SA interactions is important because their blends have been used extensively in tissue engineering, as dressing, scaffolds, and controlled release devices with a differentiated microstructure as a consequence of the interaction between the polymers [1215]. The biocomposites of SF–SA also presented interesting mechanical properties, such as high flexibility and resistance to tear [15], due to the polymer interaction. Therefore, fundamental studies about the polymer interactions could stimulate the development of new products in this area. Moreover, the SA content of β-D mannuronic acid (M block) and α-L guluronic acid (G block) seemed to influence the physicochemical, rheological, and mechanical properties of SA [16] and, consequently, could interfere on the properties of SF–SA blends. Therefore, it is extremely important to relate the thermodynamic properties of this system, to the influence of guluronic and mannuronic acids of SA, in order to fine-tune the characteristics of the blends, such as morphology, mechanical, and thermal properties.

The interaction between SF and SA occurs by intra and intermolecular hydrogen bonds [17, 18] between the hydroxyl/carboxyl and amine groups (Figure 1); however, the thermodynamic properties are not completely clear. Lopes et al. [19] studied SA and SF blends, observing that there is a small area in the phase diagram where the formation of two phases is observed and that SA, as other hydrophilic polymers, seems to interfere in the solvation of the SF in solution, having a salting-out effect [19, 20]. It is claimed that this effect probably accelerates the formation of β-sheet structures and the solid-like phase formation [10]. Therefore, fundamental studies, like the present one, are important to predict the system behavior and allow the blend production with specific properties.

In this study, equilibrium data of SF and SA blends were acquired and adjusted in three different thermodynamic models: Flory–Huggins, non-random two-liquids (NRTL), and universal quasichemical (UNIQUAC). The influence of the SA composition (proportion of M and G blocks) on the phase equilibrium and on the interaction parameters was investigated. Two alginates extracted from different seaweeds were used. Cohn equation was used to analyze the potential of SA to precipitate the SF and form the solid-like phase. Since the properties of the system can be directly influenced by blend composition, it is important to understand the parameters that can influence on the polymer interaction and consequently on blend miscibility.

2. Materials and Methods

2.1. Natural Polymers

The SF cocoons produced by Bombyx mori silkworm were donated by Bratac (Brazil). Two different SAs were used, both from FMC Polymers (USA), one extracted from Macrocystis pyrifera seaweed, which will be called alginate-M, and the other extracted from Laminaria hyberboria seaweed, which will be called alginate-G.

2.2. Sodium Alginate Characterization
2.2.1. Determination of Viscosimetric Molar Mass

The viscosimetric molar mass () was estimated using the intrinsic viscosity data and the Mark–Houvink equation (Equation (1)).where K and α are constants equal to 1.228 × 10−4 and 0.963, respectively, and is the intrinsic viscosity.

The intrinsic viscosity was measured by using SA solution with different concentrations (0.05, 0.10, 0.15, 0.20, 0.25, and 0.3 g/dL) diluted on 0.1 mol/L NaCl solution.

The tests were carried out using an Ostwald–Cannon–Fenske capillary viscometer and a thermostatic recirculation bath (Nova Ética, Brazil) at 25 ± 0.1°C. The volume of SA solution used was equal to 8.0 mL, and the flow time was marked with a digital stopwatch. For each solution, five time measurements were taken, and the average was used to calculate the viscosity. The relative () and inherent () viscosities were calculated using Equations (2) and (3), respectively, as follows:where t is the alginate solution flow time in the viscometer (s); t0 is the pure solvent flow time in the viscometer (s); C is the concentration of the alginate solution (g/dL).

2.2.2. Determination of Mannuronic and Guluronic Acid Residues

The molar ratio of β-D-mannuronic (M) and α-L-guluronic (G) acid residues, the molar fractions of GG, MM, and GM/MG and the sequences FGG, FMM, and F(GM/MG) were determined by 1H nuclear magnetic resonance (1H-NMR) spectroscopy in deuterated water (D2O) at 70°C using a Bruker 300 MHz spectrometer. The data obtained were analyzed using the MestReNova 14.1.0 software (Mestrelab Research S.L., Spain) and calculated using Equations (5)–(9) [21].where I5.05, I4.67, and I4.45 are the peak intensities of 5.05, 4.67, and 4.45 ppm, respectively; FG is the molar fraction of the G blocks; FM is the molar fraction of the M blocks; FGG, FMG, FMM, and FGM are the molar fractions of groups GG, MG, MM, and GM, respectively.

2.3. SF and SA Solutions Preparation

SF solution was prepared following the method previously published [6]. Briefly, silkworm cocoons were degummed using 1 g/L Na2CO3 solution in a thermostatic bath at 85°C (Nova Analitica, Brazil) for 30 min, and the process was repeated two times. The obtained fibers were washed with deionized water and maintained at room temperature until completely dry. They were milled and passed through a 10 mesh sieve and then solubilized in a ternary solution of CaCl2 : ethanol : water (1 : 2 : 8 molar rate). To remove the salts, the solution was dialyzed using a membrane (MWCO 3.5 kDa, Thermo Fisher Scientific, EUA) and ultrapure water (1 : 10 rate). To remove the insoluble parts that are formed during dialysis, the solution was centrifuged for 30 min and 2,300 relative centrifugal force.

SA powder extracted from M. pyrifera and Laminaria hyperboria was dissolved in deionized water with a concentration of 2% (w/v). The solution was kept for 3 days, so the water could solvate the powder, and then the solution was homogenized by stirring. To perform the SA quantification, it was necessary to label it with 4,6 diamidino-2-phenylindole (DAPI). First, 1-ethyl-3 (3-dimethylaminopropyl) carbodiimide (EDC) (0.2 mg/mL) was added to the SA solution under stirring at 4°C. After EDC dissolution, DAPI was added in a concentration equal to 0.1 mg/mL. The alginate-DAPI system was kept in an ice bath and under agitation for 48 hr. The resulting solution was dialyzed for 72 hr using 1 : 100 distilled water, changed every 24 hr.

2.4. Blend Preparation

To prepare the SF/SA blends, SF solution was added to the SA solution slowly and under magnetic stirring in different mass proportions. There was no pH adjustment after the preparation of the solutions, which were around 7.0 for both polymeric solutions. The mass proportions that were used were based on previous work [10] and are shown in Table 1. The blends were kept at room temperature, around 25°C, in a thermostatic bath (Tecnal, Brazil). The nomenclature used in the work to refer to the blends of SF and SA will be SFxSAyWAz, where x, y, and z are the mass fractions of the components (in percentage), SF is silk fibroin, SA is sodium alginate, and WA is water.

2.5. Optical and Fluorescence Microscopy

The optical microscopy was used to analyze the morphology of the mixture between SF solution and SA after mixing the two solutions. For this, the SA solution was placed on a glass slide; then, the SF solution was added to the SA solution and gently mixed with a spatula. Nikon optical microscope (Tokyo, Japan) model Eclipse E200 was used to take images of the blend.

The fluorescence microscopy was used to investigate the composition of the globular structures observed after blend formation. The SA labeled with DAPI, as described before, was used, and the blends were analyzed with a Zeiss model LSM 710 fluorescence microscope (Germany).

2.6. Components Phase Quantification

The amount of water was quantified by evaporation in both phases. To perform the quantification of SA in the liquid phase, it was necessary to label it with a fluorescent molecule, DAPI, as described before. After labeling it, the alginate was quantified by fluorescence spectroscopy in a Varioskan LUX Plate Reader (ThermoFisher Scientific, USA) at 358 nm (emission) and 461 nm (excitement).

SF was quantified, in liquid phase, by Bradford’s method of protein quantification. This method was chosen because it does not suffer interference from polysaccharides. The standard curve was prepared using SF solution instead of bovine serum albumin to obtain more reliable results. The analysis was performed using 96-well plates, in which 10 µL of the sample was mixed with 190 µL of Bradford reagent (Scienco, Brazil). After 30 min, the plate was read at 595 nm in a Varioskan LUX Plate Reader (ThermoFisher Scientific, USA).

In the solid phase, both SF and SA were quantified by mass balance, using Equation (10).where , and are the total blend mass, the mass of the component in the solid phase and the mass of the component in the liquid phase, respectively.

3. Thermodynamic Model Calculation

3.1. Flory–Huggins Modeling

Equation (11) was used to estimate Flory–Huggins parameters. For each component, in each phase, Equation (11) can be written as shown in Equations (12)–(14), obtaining a linear equations system, which can be easily solved.where ai is the activity, , , and are the volumes fractions of component i, j, and k, Vi is the volume of component i, Vs is the volume of the solvent, T is the temperature, is the interaction parameter of i and j components, and is the interaction parameter of k and j components.

The volume fraction was calculated from the mass fraction using Equation (15). The density and specific volume values for SA [22], SF [23], and water are shown in Table 2.

3.2. NRTL and UNIQUAC Models

The modeling of experimental data for the NRTL and UNIQUAC models was performed using the TME-LLE version 2.0 software written in Fortran language [24]. The software was developed to estimate the binary interaction parameters for Gibbs free energy models with any number of components and different temperatures through the minimization of the objective function (FO) shown in Equation (16) using the modified simplex method.where n is the number of tie lines of each system, P is the number of components, and are the experimental and calculated molar fraction of component i in the liquid phase, respectively, and and are experimental and calculated molar fraction of component i in the solid phase, respectively.

The r and q structure parameters used in the UNIQUAC model were estimated using the subgroup contribution method, according to Equations (17) and (18), and are shown in Table 3. The parameter values of Rk and Qk were taken from the literature [25]. This was necessary, as there is no information in the literature about the r and q parameters for SA and SF. For fibroin, the subgroups considered were those present in the chemical structure of alanine, serine, and glycine, which are the main amino acids present in the structure of SF. As for SA, the subgroups of α-L-guluronic acid and β-D-mannuronic acid were considered.where is the numbers of k type subgroups in a molecule i, Rk and Qk are the relative volume and the relative area of each subgroup in the molecule, respectively.

The root mean square deviation (RMSD) was used to evaluate the correlation between experimental and theoretical data obtained by the thermodynamical models, and it is shown in Equation (19).where n is the number of tie lines of each system, P is the number of components, and the exp and cal refer to the calculated and experimental values of mass fraction, respectively.

3.3. Cohn Equation

For protein systems with or without electrolytes that present liquid–solid separation with the formation of precipitates, such as the fibroin-alginate blend, a simple equation that adjusts the equilibrium data is the Cohn [26] equation, that relates the protein solubility with the salting-out constant (Ks).

The study of SF solubility in different SA solutions was carried out using samples with constant volume of silk solutions with different concentrations of SA solutions. Samples containing 4.5 mL of SA solution and 0.5 mL of 2% SF solution were prepared. The concentration of SA solutions varied from 0.4% to 1.4% (m/m). The blends were kept in a thermostatic bath at 25°C for 8 days. The phases were separated by centrifugation, and the polymers were quantified in the liquid phase as previously described. After quantification, the data were adapted to the Cohn equation (Equation (20)) and the solubility constants β (Ks) were determined.where S is the solubility of the protein, Ks is the protein solubility constant, m is the concentration of the precipitating agent, and β is the constant related to the net charge of the protein; for nonelectrolytes, this parameter is related to the solubility of the protein in the absence of precipitating agents.

4. Results and Discussion

4.1. Sodium Alginate Characterization

The results of viscosity and viscosimetric molar mass for alginate-M and alginate-G are shown in Table 4. The viscosity and molar mass of the alginates are significantly different; the alginate extracted from L. hyberboria seaweed has a higher viscosity and molar mass than the alginate extracted from M. pyrifera seaweed.

The hydrogen NMR spectra (1H-NMR) for the two types of alginate used are shown in Figure S2 of the supplementary materials. The spectra showed characteristic signals between 5.05 and 4.55 ppm, referring to the hydrogen bond of the G and M blocks, respectively [21, 27]. The fractions of the M and G blocks and their proportion for the studied alginates are shown in Table 5. The spectra of alginates showed chemical shifts close to 5.3 ppm, which can be attributed to the reducing end hydrogen [27], which is found in higher values when compared to the hydrogen in the inner chain. The presence of this sign suggests alginates with lower molar mass [21].

The ratio between M and G blocks reported in the literature for SAs extracted from brown seaweed is between 0.4 and 1.94, being higher in SAs from M. pyrifera seaweed [28, 29]. The amount of G and M blocks present in the polymer chain is directly related to its configuration and, consequently, to its interaction with other polymers and ions. As shown in Table 4, alginate-M, extracted from M. pyrifera, has a higher proportion of mannuronic groups (M) in the chain, which are predominantly alternated (MG/GM), while alginate-G, extracted from Laminaria hyperborea, has a higher proportion of guluronic groups (G and GG).

4.2. Optical and Fluorescence Microscopy

Figure 2 shows the optical microscopy image obtained after mixing the SA and SF solutions. It is observed the formation of micelle-like structures, with a mean diameter of 0.247 ± 0.120 mm, indicating that phase separation occurs immediately after the polymers mixing. During the preparation of the SF and SA blends, it is possible to notice that when SF solution is added to SA solution, the blend becomes cloudy, with a milk-like aspect. The turbidity could be associated with the formation of these structures, which could interact and coalesce over time, trapping water molecules and forming a solid-like phase observed after some days in these blends, as described in previous work [10].

The fluorescence microscopy was performed to investigate the composition of these micelle-like structures. The alginate-M and alginate-G were labeled with DAPI as described previously; the images obtained are shown in Figure 2. It is possible to observe the SA particles labeled with DAPI before the addition of fibroin, as they show homogeneous fluorescence. After the addition of fibroin, the same micelle-like structures observed with the optical microscope were detected with the fluorescence microscope. It can be seen that they contain SA but are heterogeneous, indicating the presence of SF. It is possible to observe SA outside the structures, indicating that part of the polymer remains in the second phase. Similar structures were observed in SF and SA membranes [6, 11]. De Moraes et al. [6] found that the structures observed were formed by both polymers and that after the formation of the membrane, it is possible to extract the SF and obtain porous membranes, providing a very interesting microstructure for this material. So, it is possible that the globules observed in the membranes are formed during the blend preparation. The formation of micelle-like structures has also been observed in other polysaccharide and protein systems, such as alginate-casein and alginate-soy protein [30, 31].

4.3. Phase Equilibrium

After 8 days of the blend preparation, two distinct phases are observed, one liquid and another solid-like, as described in previous work [10]. The quantification was carried out after reaching thermodynamic equilibrium, that is, on the 8 days after the preparation of the blend [11]. The mass fraction values of SF, SA, and water in each phase are shown in Table S1 of supplementary materials, and the ternary diagrams for the blends prepared with the different SAs are shown in Figure 3.

After the phase separation, SA is expected to remain in the liquid phase, as it is a hydrophilic polymer and has greater affinity with the water. However, blends with a higher mass fraction of SA also presented a greater amount of this polymer in the solid phase. The analysis was performed by comparing blends with the same initial mass fraction of SF and different mass fractions of SA, as shown in Figure 4. It is possible that the micelles-like structures formed after the preparation of the blends agglomerate and precipitate, forming the solid-like phase. Thus, one hypothesis is that blends with a higher initial mass fraction of SA form micelles with a higher concentration of this polymer, trapping it and enriching the solid phase. This behavior was observed by other authors for the fibroin-alginate system and for other protein–polysaccharide systems such as alginate-pea protein and κ-carrageenan-soy protein [3133].

4.4. Thermodynamic Model Fitting
4.4.1. Flory–Huggins

The values of the Flory–Huggins parameters are shown in Table 5, and the calculated mass fraction values from the model are shown in Tables S2 and S3 in supplementary materials, which have been plotted in the ternary diagrams shown in Figure 5. In Figure 5, each point was represented by a color to facilitate the reader’s visualization and interpretation. Fibroin–alginate interaction parameter values (χSFSA) were close to those found in previous works in which glass transition temperature data were used to calculate them [10]. The values of χSAWA and χSFWA varied in relation to the type of SA used but remained in the same order of magnitude. The parameters were adjusted from the experimental data, and since they are different for systems with different types of alginates, it is expected that the values will vary a little. For SA, the value of the alginate–water interaction parameter reported in the literature is around 0.025 [34], and for SF, this value is close to 0.95 [35].

The binodal curve can be experimentally determined by the cloud point method, in which the clear solutions of the polymers are mixed in a controlled method until the blend becomes cloudy, and the mixture is considered two-phase. For the fibroin–alginate system, this experiment is difficult to carry out in practice due to the micelles-like formation and the rapid turbidity of the blend, even in low mass fractions of SF and SA. However, it was possible to estimate the binodal curve from the Flory–Huggins model, shown in Figure 5. The RMSD value of the model was about 2%. Despite the values are small, the poor adequacy of the experimental data to the Flory–Huggins model is easily noticed in the ternary diagram. As it is a small region of the phase diagram, small variations in the mass fraction can cause significant changes in the diagram’s tie lines.

4.4.2. NRTL and UNIQUAC Models

The experimental data of mass fraction were also used to calculate the interaction parameters of the NRTL and UNIQUAC models, shown in Table 6. The mass fraction values calculated from the thermodynamic models are shown in Tables S4S7 in supplementary materials and were plotted in the ternary diagrams shown in Figures 6 and 7. For the NRTL model, the values of the adjustable (Ai,j and Aj,i) and nonrandom (αi,j) fibroin–alginate and alginate–water parameters presented significantly different values when changing the type of SA used in the blend. The same was not observed for the adjustable parameters related to the fibroin–water pair, which have different values but of the same order of magnitude.

The UNIQUAC model presented a good fit to the experimental data, shown in Figure 7. Unlike the NRTL model, the adjustable parameters (Ai,j and Aj,i) did not show significant differences with the change in the type of SA. The values had the same order of magnitude for both systems, except for the parameter A1,2, referring to the fibroin-alginate adjustable parameter. The differences between the observed parameters may be a consequence of the change in the configuration of the SA chain, since no significant changes were observed for the fibroin–water pair. The fraction of groups M and G in the SA chain causes a change in the arrangement of the functional groups. However, it is important to recognize that the molar mass can have an effect on the entropy of the mixture, since the greater its value, the greater the polymer chain and the smaller the possibilities of molecular arrangement. These two factors (disposition of M/G blocks and molar mass) can cause changes in the interaction between SA and the other components. The NRTL and UNIQUAC models had lower RMSD values than the Flory–Huggins model, and the one that best suited the experimental data was the UNIQUAC. This may happen because the UNIQUAC model uses the shape and volume data of the molecules involved in the system to determine the adjustable parameters, which can generate better results.

4.4.3. Cohn Equation

Another thermodynamic approach for the system is based on the study of the salting-out effect of SA in the ternary system, which was investigated by the Cohn equation, which can be used to analyze liquid–solid phase separations. The parameters β and Ks obtained from the adjustment of the experimental data (Figure 8) for each system are shown in Table 7. The parameter Ks is called the salting out constant and expresses the ability of SA to induce this effect. The value obtained for alginate-G was higher than alginate-M, which suggests that it has a greater capacity to precipitate SF.

The salting-out phenomenon occurs when water-soluble molecules are excluded from the aqueous phase due to the high concentration of salt [36]. Despite not being, by definition, a salt, SA promotes a similar effect. Due to its hydrophilicity, the polymer remains solvated in water, stimulating fibroin–fibroin interactions and its subsequent precipitation. A similar phenomenon is observed in fibroin-PEG systems [37]. Badra [32] investigated SF and SA systems and modeled the equilibrium data to the Cohn equation, obtaining Ks values of the same order of magnitude. The author used the same SA at different pHs and found that the increase in pH did not cause significant effects on the salting-out phenomenon. In addition to the difference in the proportion of M and G groups, the SAs used differ significantly in relation to their viscosity and molar mass, which are higher in G-alginate. It is possible that, as it has a longer chain, alginate-G has a greater molar volume, requiring more water molecules to solvate it, and, therefore, the salting out effect is more pronounced when alginate-G is used in the SF/SA blend.

5. Conclusion

It was possible to identify differences between the compositions of the SAs used. The proportion between M and G groups has a significant influence on the thermodynamic equilibrium of SF and SA blends. By using the equilibrium data, the phase diagram data were fitted to thermodynamic models of Flory–Huggins, NRTL, and UNIQUAC, in addition to verifying the salting-out effect of each SA from the Cohn equation. The thermodynamic model that presented the best fit to the equilibrium experimental data was UNIQUAC, presenting the lowest mean square deviation value. For NRTL and UNIQUAC models, it was possible to observe differences in the interaction parameters between the systems that used alginate-M and alginate-G, indicating that these thermodynamic models were more sensitive to variations in the SAs used in each system. Regarding the solubility parameters, alginate-G showed a higher constant value of Ks, being more efficient in precipitating SF.

These results represent an advance in the thermodynamic study of blends using SF and SA, natural polymers that have high potential for application as biomaterials. The results emphasize the importance of the source of natural polymers and the precise characterization of their composition, since these factors may have a significant influence on blend characteristics and, consequently, on blend final applications. It is expected that the data presented will be used as a basis and guide for further studies, aiming the application of fibroin and alginate in biomaterials or theoretical studies of the interaction of this protein with other macromolecules.

Data Availability

No data are used in this study.

Additional Points

Highlights. Influence of sodium alginate mannuronic and guluronic groups on the interaction with silk fibroin. Thermodynamic modeling of silk fibroin and sodium alginate blends using Flory–Huggins, NRTL, and UNIQUAC models. Use of the Cohn equation to analyze the potential of sodium alginate to precipitate the silk fibroin in solution.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

The authors acknowledge the financial support of the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil (CAPES)—Finance Code 001.

Supplementary Materials

Figure S1: graphical abstract illustrating a silk fibroin and sodium alginate blend after phase separation in liquid and solid phases. Figure S2: hydrogen RMN spectra (1H-RMN) for alginate-M (a) and alginate-G (b) using D2O as solvent at 70°C. Table S1: average of molar fraction (n = 3) of silk fibroin (SF), sodium alginate (alginate-M and alginate-G), and water in liquid and solid phases. Table S2: values of mass fraction for silk fibroin, alginate-M, and water calculated with Flory–Huggins model. Table S3: values of mass fraction for silk fibroin, alginate-G, and water calculated with Flory–Huggins model. Table S4: values of mass fraction for silk fibroin, alginate-M, and water calculated with NRTL model. Table S5: values of mass fraction for silk fibroin, alginate-G, and water calculated with NRTL model. Table S6: values of mass fraction for silk fibroin, alginate-M, and water calculated with UNIQUAC model. Table S7: values of mass fraction for silk fibroin, alginate-G, and water calculated with UNIQUAC model. (Supplementary Materials)