A Novel Method for the Preparation of Casein–Fucoidan Composite Nanostructures

The aim of the study was to develop casein–fucoidan composite nanostructures through the method of polyelectrolyte complexation and subsequent spray drying. To determine the optimal parameters for the preparation of the composite structures and to investigate the influence of the production and technological parameters on the main structural and morphological characteristics of the obtained structures, 3(k-p) fractional factorial design was applied. The independent variables (casein to fucoidan ratio, glutaraldehyde concentration, and spray intensity) were varied at three levels (low, medium, and high) and their effect on the yield, the average particle size, and the zeta potential were evaluated statistically. Based on the obtained results, models C1F1G1Sp.30, C1F1G2Sp.40, and C1F1G3Sp.50, which have an average particle size ranging from (0.265 ± 0.03) µm to (0.357 ± 0.02) µm, a production yield in the range (48.9 ± 2.9) % to (66.4 ± 2.2) %, and a zeta potential varying from (−20.12 ± 0.9) mV to (−25.71 ± 1.0) mV, were selected as optimal for further use as drug delivery systems.


Introduction
In recent decades, biopolymers have been the subject of a serious scientific interest for the development of nanosized drug delivery systems due to their specific characteristics and numerous advantages over synthetic polymers, such as low immunogenicity, high biocompatibility, biodegradability, and low toxicity [1].Biopolymers such as pep-tides and polysaccharides have been widely used for the preparation of nanoparticulate drug delivery systems that can provide controlled and targeted release, improving biodistribution of the encapsulated drug and minimizing its systemic toxicity [2].Compared to polysaccharides, proteins are preferred as structure-forming biopolymers for the production of drug delivery systems due to their ability to form smaller-in-size particles [3].Among the proteins used for the formulation of nanosized carriers, casein has shown optimal structural and physicochemical properties [4].Casein is a collective term for a family of calcium (phosphate)-binding phosphoproteins commonly found in mammalian milk.Bovine milk casein consists of four peptides, namely αs1, αs2, β, and k, which differ in amino acid, phosphorus, and carbohydrate content, but all of them possess amphiphilic properties [5].Cysteine amino acid residues, which are involved in the formation of disulfide bonds, are found only in the polypeptide chains of k-casein.The lack of secondary structures of the casein molecules, due to the proline-rich amino acid residues and the ability to bind to calcium phosphate, leads to the formation of electrostatic, hydrogen, and hydrophobic interactions [6].As a result, self-association of the casein peptides occurs, resulting in the formation of stable agglomerates known as casein micelles [6].The inner part of the micelle is composed of αs1, αs2, and β caseins, while the outer layer, which stabilizes the Polymers 2024, 16, 1818 2 of 23 micelle, contains glycosylated k-casein [7].Casein micelles exhibit pH-dependent behavior.A decrease in the surface negative charge of the casein micelles leads to their contraction, while an increase in the charge leads to an electrostatic repulsion of the casein molecules [8].
A major problem in the development of casein-based nanostructures is the low drug encapsulation efficiency and drug loading index [9].This is related to the use of large amounts of drug delivery systems to ensure the desired dose of the incorporated therapeutic agent [9].Thus, the development of suitable technological approaches to overcome these drawbacks is necessary for the preparation of stable and effective drug delivery systems.One of the main approaches is the formation of casein-polysaccharide complexes.
Due to its specific physicochemical properties, casein can participate in the formation of composite nanostructures because of polyelectrolyte or chemical interactions with one or more polymers.Compared to chemically and physically cross-linked casein nanoparticles, composite nanoparticles are characterized by higher stability and the possibility of modified drug release [10].
Protein-based carriers for drug delivery, including casein, are characterized by physical instability, resulting in protein aggregation and precipitation at a pH close to the isoelectric point of the protein.The incorporation of polysaccharides into the carrier structure induces steric hindrance, which limits protein aggregation [11].As a result, protein-polysaccharide composite nanoparticles have increased stability and more efficiently encapsulate hydrophilic and hydrophobic therapeutic agents [12].These complexes also provide protection to the included substance from oxidation and hydrolysis under the influence of environmental factors such as pH, light, and temperature [13].In addition, the colloidal network composing the nanoparticles, built of protein-polysaccharide complexes, provides a delayed release of the incorporated therapeutic agent [14].Polysaccharide complexation of proteins can be due to both covalent and non-covalent binding of biopolymers.Covalent protein-polysaccharide conjugation is due to enzymatic cross-linking (oxidases and transglutaminases), chemical cross-linking (genipin, glutaraldehyde, and polyethylene glycol), or through the Maillard reaction [15].
In a study conducted by Pan et al., the Maillard reaction was used to prepare caseindextran copolymer [16].The obtained structures showed pH-dependent behavior, with micelle formation occurring at a pH equal to the isoelectric point of casein, and demonstrated successful incorporation of the fluorescent agent pyrene [16].
Covalent interaction between proteins and polysaccharides also occurs because of cross-linking with transglutaminase, which catalyzes acyl-transfer reactions between the γ-carboxamide groups of glutamine residues and free amino groups [17].
Glutaraldehyde is among the most widely used chemical cross-linking agents, and provides a high degree of cross-linking of protein-polysaccharide complexes [18].Nevertheless, the use of high concentrations of the cross-linking agent is associated with the occurrence of adverse drug reactions, such as breathing disorders, development of asthma, and dermatitis, as well as irritation of the nasal mucosa and ocular conjunctiva [19].
Physical mixing of two biopolymers leads to the formation of non-covalent proteinpolysaccharide bonds, which largely depends on the nature of the biopolymers, their quantitative ratio, pH of the medium, and ionic strength [20].Regarding physical proteinpolysaccharide complexes, non-covalent interactions can be due to hydrogen bonding and electrostatic and hydrophobic interactions [21].Most polysaccharides, except for chitosan, are negatively charged, and at pH values below the isoelectric point of casein, the formation of a casein-polysaccharide complex occurs [22].
Several studies on the preparation of casein-polysaccharide complexes have been reported in the literature, but only a few are related to the preparation of composite caseinpolysaccharide drug delivery nanocarriers.Markman et al. [23] reported the preparation of vitamin D2-loaded casein-maltodextrin nano capsules.The obtained conjugates show a high percentage of therapeutic agent incorporation efficiency (up to 90%).Furthermore, the carrier demonstrated high colloidal stability and prevention of vitamin D2 oxidation in acidic pH conditions.Regarding the release, the carrier showed a pH-dependent behavior.Lack of release of the incorporated nutraceutical was observed at an acidic pH (2.5), while at pH 7 up to 90% of the incorporated nutraceutical was released, demonstrating the potential use of the carrier as a drug delivery system for targeted delivery in the small intestine.In another study, a drug delivery system for epigal-locatechin-3-gallate was obtained via polyelectrolyte complexation between casein and chitosan [24].The obtained nanostructures demonstrated an average particle size of 150 nm, a high degree of biocompatibility, and improved intestinal permeability of epi-gallocatechin-3-gallate.As a result of polyelectrolyte complexation, the preparation of alginate-coated casein nanoparticles loaded with doxorubicin has also been reported [25].The results showed that, after coating the casein particles with alginate, an increase in their average particle size from 264 nm to 294 nm was observed.Furthermore, the drug delivery system showed a high drug loading index and encapsulation efficacy (95%) of the cytostatic agent-doxorubicin.In vitro release was performed in two different media (pH 7.4 and pH 5), simulating systemic circulation and tumor cell pH.The results demonstrated a delayed release in alkaline pH and a more pronounced release in acidic pH [25].
One of the polysaccharides with potential application for the preparation of caseinpolysaccharide complexes is the biopolymer fucoidan.Fucoidan belongs to a family of sulfated, L-fucose-rich polysaccharides found in the cell wall of various species of brown algae (Phaeophyta: Laminariaceae, Fucaceae, Chordariaceae, and Alariaceae) [26].It can also be obtained from sea cucumbers (Holothuroidea: Stichopodidae, Holothuri-idae), sea urchin eggs (Echinoidea: Strongylocentrotidae, Arbaciidae), and sea grasses (Cymodoceaceae) [27].Over 90% of the total sugar content in fucoidan consists of deoxyhexose-L-fucose [28].Along with fucose, other monosaccharides can be found in the fucoidan structure: galactose, glucose, xylose, mannose, rhamnose, arabinose, uronic acid, and acetyl groups [29].Although the structure and composition of fucoidans differ depending on the origin and extraction method, fucoidans can be categorized into two main types according to their monomeric subunits (Figure 1).
the carrier demonstrated high colloidal stability and prevention of vitamin D2 oxida in acidic pH conditions.Regarding the release, the carrier showed a pH-dependent havior.Lack of release of the incorporated nutraceutical was observed at an acidic (2.5), while at pH 7 up to 90% of the incorporated nutraceutical was released, demons ing the potential use of the carrier as a drug delivery system for targeted delivery in small intestine.In another study, a drug delivery system for epigal-locatechin-3-ga was obtained via polyelectrolyte complexation between casein and chitosan [24].The tained nanostructures demonstrated an average particle size of 150 nm, a high degre biocompatibility, and improved intestinal permeability of epi-gallocatechin-3-gallate a result of polyelectrolyte complexation, the preparation of alginate-coated casein n particles loaded with doxorubicin has also been reported [25].The results showed t after coating the casein particles with alginate, an increase in their average particle from 264 nm to 294 nm was observed.Furthermore, the drug delivery system show high drug loading index and encapsulation efficacy (95%) of the cytostatic agent-d rubicin.In vitro release was performed in two different media (pH 7.4 and pH 5), si lating systemic circulation and tumor cell pH.The results demonstrated a delayed rel in alkaline pH and a more pronounced release in acidic pH [25].
Depending on the sulfate content and molecular weight, fucoidans show different biological effects (antioxidant, antibacterial, antiviral, anticoagulant, hypolipidemic, hypoglycemic, antitumor, anti-inflammatory, and immunomodulatory).In addition to its biological activity, fucoidan also possesses suitable physicochemical properties as a drug carrier, which is confirmed by numerous reports in the literature about micro-and nanoparticles developed on this basis, liposomes, films, hydrogels, etc. [30,31] The polymer backbone has anionic properties due to the presence of negatively charged sulfate ester groups.Depending on the specific chemical characteristics and charge density, fucoidan can interact with various biomolecules, such as proteins and other polysaccharides [32].It is soluble in water and in acidic solutions, and its solubility depends on the molecular weight and the number of sulfate groups.The highly water-soluble fucoidans form solutions with low viscosity and are therefore not generally used as gelling agents.On the other hand, the viscoelastic properties of fucoidan are highly dependent on the origin of the polysaccharide, its concentration, molecular weight, sulfate content, pH, and temperature [27].The interaction of fucoidan with oppositely charged biopolymers, such as chitosan [33], poly (2-hydroxyethyl methacrylate) [34], soy protein isolate (SPI) [35], and lactoglobulin [36], can lead to the formation of polyelectrolyte complexes with improved characteristics for biomedical application.Polyelectrolyte complexation (complex coacervation) is one of the most widely used methods for obtaining fucoidan micro-and nanoparticles.Other common techniques for formulating drug systems from the polysaccharide are simple coacervation, ionotropic gelation, spray drying, emulsification, etc. [37].
With the possibilities of nanotechnology, the development of innovative carriers capable of providing selective and specific therapy for several diseases has long been the subject of fundamental pharmaceutical science.Polymers of natural origin such as casein and fucoidan enable the development of biodegradable nanocarriers providing modified release and targeted delivery, which significantly improves the distribution of the incorporated drugs in the body and reduces the risk of adverse drug reactions.Their complex biological effects and the ability to form electrostatic interactions provide an opportunity for the development of casein-fucoidan composite nanoscale drug delivery systems.
The aim of the present work is to develop casein-fucoidan composite nanostructures by polyelectrolyte complexation and cross-linking of the complex with glutaraldehyde.By further applying 3(k-p) fractional factorial design, the work aims to evaluate the effect of dependent variables on independent variables and outline optimal technological and process parameters for the development of effective drug delivery systems.

Preparation of Casein/Fucoidan Composite Nanoparticles
Casein/fucoidan composite nanoparticles were obtained by polyelectrolyte complexation, cross-linking of the complex with glutaraldehyde, and subsequent spray drying (Figure 2).
Sodium caseinate was dissolved in deionized water pre-acidified with 0.1 M HCl solution to pH 2 to ensure protonation of the amino groups in the protein structure.For 15 min, the resulting solution was added dropwise under continuous stirring at 25,000 rpm (Miccra MiniBatch D-9, MICCRA GmbH, Heitersheim, Germany) to an aqueous fucoidan solution (pH 7), negatively charged due to the presence of sulfate esters groups in its structure.As a result of the interaction of the oppositely charged biopolymers, a polyelectrolyte complex was formed, which was further cross-linked with glutaraldehyde (0.004 µL of glutaraldehyde solution for 1 mL of casein solution).The obtained nanosuspension was spray dried under the following conditions: mesh size of the spray membrane 4.0 µm; inlet temperature 30 • C; solution feed rate 30%; sputtering intensity 40%; drying gas rate 120 L/min; pressure 30 nbar.To derive optimal production conditions and study the influence of the production parameters (casein/fucoidan ratio, glutaraldehyde concentration) and the technological parameter sputtering intensity on the size, yield, and zeta potential, a 3(k-p) partial experimental design was applied.The advantage of the 3(k-p) design over the full factorial design is that not all possible combinations between independent variables are considered, which greatly reduces the number of models to be developed and therefore shortens the run time [18].Each of the independent variables was varied at three levels: low (−1), medium (0), and high (+1), resulting in nine possible combinations.The three levels of varying the ratio of casein to fucoidan were 1:1, 2:1, and 3:1, using 0.25% (w/v), 0.33% (w/v), and 0.37% (w/v) casein solutions and 0.13% (w/v), 0.17% (w/v), and 0.25% (w/v) fucoidan solutions.Glutaraldehyde concentration was 1% (w/v), 2% (w/v), and 3% (w/v) and the sputtering intensity was 30%, 40%, and 50%.concentration) and the technological parameter sputtering intensity on the size, yield, a zeta potential, a 3(k-p) partial experimental design was applied.The advantage of the 3 p) design over the full factorial design is that not all possible combinations between in pendent variables are considered, which greatly reduces the number of models to be veloped and therefore shortens the run time [18].Each of the independent variables w varied at three levels: low (−1), medium (0), and high (+1), resulting in nine possible co binations.The three levels of varying the ratio of casein to fucoidan were 1:1, 2:1, and using 0.25% (w/v), 0.33% (w/v), and 0.37% (w/v) casein solutions and 0.13% (w/v), 0.1 (w/v), and 0.25% (w/v) fucoidan solutions.Glutaraldehyde concentration was 1% (w/v), (w/v), and 3% (w/v) and the sputtering intensity was 30%, 40%, and 50%.

Production Yield
The production yields of the obtained nanoparticles from different batches were c culated according to Equation (1): Production yield % Spray dried nanoparticles mg Casein mg Fucoidan mg 100

Particle Size Analysis, Size Distribution, and Zeta Potential
The particle size of the obtained casein/fucoidan composite structures was analyz by dynamic light scattering (Microtrac, York, PA, USA).The system measures the part size in the range of 0.8 nm to 6.5 µm with non-invasive backscattering technology usin 3-mW helium/neon laser at a 780 nm wavelength.The system allows measurements o potential in the range from −200 mV to +200 mV.All analyses were performed in triplica at 25 °C.

Scanning Electron Microscopy
Visualization of the obtained structures was performed using scanning electron croscopy (Prisma E SEM, Thermo Scientific, Waltham, MA, USA).The samples w loaded on a copper sample holder and sputter coated with carbon followed by gold us a vacuum evaporator (Quorum Q150T Plus, Quorum Technologies, Lewes, UK).T

Particle Size Analysis, Size Distribution, and Zeta Potential
The particle size of the obtained casein/fucoidan composite structures was analyzed by dynamic light scattering (Microtrac, York, PA, USA).The system measures the particle size in the range of 0.8 nm to 6.5 µm with non-invasive backscattering technology using a 3-mW helium/neon laser at a 780 nm wavelength.The system allows measurements of ζ-potential in the range from −200 mV to +200 mV.All analyses were performed in triplicates at 25 • C.

Scanning Electron Microscopy
Visualization of the obtained structures was performed using scanning electron microscopy (Prisma E SEM, Thermo Scientific, Waltham, MA, USA).The samples were loaded on a copper sample holder and sputter coated with carbon followed by gold using a vacuum evaporator (Quorum Q150T Plus, Quorum Technologies, Lewes, UK).The images were recorded at 20 kV acceleration voltage, 10,000× magnifications using an ETD (Everhart-Thornley) detector.

Transmission Electron Microscopy
The particle size and morphology of the structures were investigated by TEM (Talos F200X, Thermo Fisher Scientific, Waltham, MA, USA).Nanocomposite suspension was added dropwise onto a formvar/carbon-coated copper grid, and then the TEM observation of the samples was performed at an operating voltage of 200 kV.

Fourier Transform Infrared Spectroscopy (FTIR)
The structural characteristics of casein, fucoidan, and the Cas/FN composite nanoparticles were characterized by FTIR.The spectra were collected in the range from 600 cm −1 to 4000 cm −1 with a resolution of 4 nm and 16 scans, using a Nicolet iS 10 FTIR spectrometer (Thermo Fisher Scientific, Pittsburg, PA, USA).The instrument is equipped with a diamond attenuated total reflection (ATR) accessory and the spectra were analyzed with OMNIC ® software package (Version 7.3, Thermo Electron Corporation, Madison, WI, USA).

Differential scanning calorimetry (DSC)
The thermal stability and phase state of the substances in the composition of the final structures were investigated using a differential scanning calorimeter DSC 204F1 Phoenix (Netzsch Gerätebau GmbH, Selb, Germany).An indium standard (Tm = 156.6 • C, ∆Hm = 28.5 J/g) was used for temperature and heat flux calibration.The sample for analysis was placed in an aluminum crucible and hermetically sealed, then heated in a temperature range of 25 • C to 300 • C at a heating rate of 10 • C/min.The results were processed by Proteus 8.0 software (Netzsch Gerätebau GmbH, Selb, Germany), and the obtained thermograms reflect the variation of the heat flux relative to the mass of the sample as a function (mW/mg) of the temperature.

Experimental Design and Statistical Analysis
To determine the optimal conditions for the preparation of the composite structures and to study the influence of the production and technological parameters over the main structural and morphological characteristics of the particles, a 3(k-p) fractional factorial design was applied.
Three optimal experimental responses were studied: production yield (Y1), average particle size (Y2), and zeta potential (Y3).They were the results of the individual influence and the linear interactions of the three independent variables.Hence, the responses were modeled by the following linear polynomial model: The results of these experiments were compared using analysis of variance (ANOVA), which was able to determine if the factors were significant.To test whether the terms were significant in regression model, t-tests were performed using a 95% (α = 0.05) level of significance.Furthermore, F-test was used to determine if there was an overall regression relationship between the response variable Y and independent variables X1, X2, X3 and their linear interactions at a 95% (α = 0.05) level of significance.The statistical analysis was performed using Statistica 8 (TIBCO Software Inc., Palo Alto, CA, USA).

Preliminary Experiments
Since the proposed method for the preparation of the structures is unique, it was necessary to define the technological parameters for the preparation of the composites by the method of spray drying.

Optimal Spray Head Mesh Size
To determine the optimal spray head mesh size, a suspension containing caseinfucoidan composite nanostructures obtained at a casein:fucoidan ratio of 1:1 and cross-linked with 1% glutaraldehyde solution was sprayed through three spray mesh heads with a pore size of 4 µm, 5.5 µm, and 7 µm, at inlet temperature 30 • C, solution feed rate 30%, spray intensity 30%, and drying gas rate 120 L/min.The results showed that the increase in the spray mesh head leads to an increase in the particle size of the obtained structures.Since the aim was to obtain particles with sizes between 100 nm and 500 nm, 4 µm was chosen as the optimal spray head mesh size, providing the formation of composite particles with an average size of 137.22 ± 21 nm.

Inlet Temperature and Solution Feed Rate
Since casein denatures at temperatures above 60 • C, to maintain the stability of the protein a temperature of 30 • C was chosen as optimal for the spray drying.
To determine the optimum solution feed rate of the process, a suspension containing casein-fucoidan composite nanostructures prepared at a casein:fucoidan ratio of 1:1, crosslinked with 1% glutaraldehyde, and spray dried at an inlet temperature of 30 • C, spray mesh head 4 µm, spray intensity of 30%, gas flow rate 120 L/min, and solution feed rate varying from 30% to 60%.The results showed that spraying at 20% solution feed rate led to the formation of smaller particles, 122 ± 27 nm, but the spray drying time was prolonged.Spray drying at a solution feed rate of 30% and 40% resulted in the production of structures with an average particle size of 147 ± 28 nm and 388 ± 21 nm.Further increase in the solution feed rate up to 50% resulted in the production of larger particles with an average size of 602 ± 35 nm.Based on the obtained results, the optimum solution feed rate to was set to 30%.

Gas Flow Rate
The optimum gas flow rate was found in the atomization of a suspension containing casein-fucoidan composite nanostructures prepared at a casein:fucoidan ratio of 1:1, crosslinked with 1% glutaraldehyde, and spray dried at an inlet temperature of 30 • C, spray mesh head 4 µm, spray intensity of 30%, solution feed rate 30%, and varying gas flow rate from 90 L/min to 120 L/min.The results showed that spray drying at a gas flow rate of 90 L/min was inefficient and resulted in deposition of the polymer to the walls of the glass chamber.Increasing the gas flow rate to 100 L/min led to effective drying of a small part of the feeding suspension, but at the same time a part of the polymer was deposited on the walls of the glass chamber.Further increase of the gas flow rate to 120 L/min resulted in efficient drying of the suspension and the formation of composite structures demonstrating a significant yield (64.8 ± 3.0) %.As a result of these preliminary tests, a gas flow rate of 120 L/min was selected as optimal for the spray drying of the composites.

Spray Intensity
Based on the literature review regarding the technological parameters affecting the process of nano spray drying, it was found that the spray intensity largely affects the main characteristics of the obtained particles.Therefore, the influence of this technological parameter on the structural-morphological characteristics of the obtained particles was determined by applying appropriate factorial design.

Experimental Design
To determine the optimal conditions and to study the influence of the production parameters (casein/fucoidan ratio, glutaraldehyde concentration) and the technological parameter spray intensity on the size, yield, and zeta potential, a 3(k-p) fractional factorial design was applied.The dependent variables were production yield (Y1), average particle size (Y2), and zeta potential (Y3), and the independent variables were casein:fucoidan ratio (X1), glutaraldehyde concentration (X2), and spray intensity (X3) (Table 1).Each of the independent variables was varied at three levels: low (-1), medium (0), and high (+1), resulting in nine possible combinations (Table 2).To estimate the effects of the independents on the dependent variables and to generate polynomial equations for predicting the response, linear regression analysis and analysis of variance (one way ANOVA) were applied using Statistica 8 (TIBCO Software Inc., CA, USA).The main effects of the independent variables, which represent the mean values reflecting the change of an independent variable from its lowest to its highest level on the dependent variables, are presented in Table 3.The main regression coefficients presented in the table, which apply to the linear equation, can be used to predict the responses at different levels of variation.
β0 represents the intercept with the y-axis; X1, X2, X3 are the observed effects of the independent variables; X1 X2, X1 X3, and X2 X3 are the observed effects of the interaction, while β1 to β6 demonstrate the linear regression coefficients.After the evaluation of the Polymers 2024, 16, 1818 9 of 23 main effects, the determination of the significant independent variables for each dependent variable is performed by analysis of variation (one way ANOVA) (Table 4).In the table, the sum of squares (SS) provides information for estimating the main effects of the factors, and the F-values represent the ratios of the corresponding mean squared effects and mean squared errors [38].The table also presents the p-values that show the statistical significance of each factor (p < 0.05).
Based on the data obtained from the analysis of variations, in terms of F-and p-values, the statistical significance of each of the independent variables, as well as their combined interaction on the investigated dependent variables, can be determined.The independent variables that demonstrate the highest F-values and p < 0.05 are the main factors with a statistically significant impact on the studied dependent variables [38].
The results show that the linear increase in casein/fucoidan ratio and spray intensity, and the linear interaction between casein/fucoidan and glutaraldehyde, as well as the interaction between glutaraldehyde and spray intensity, have a statistically significant effect on production yield (p < 0.05).Regarding particle size and zeta potential, only a linear increase in the casein/fucoidan ratio has a statistically significant effect (p < 0.05).The standardized effects of the independent variables are represented by Pareto charts, which rank the effects of the independent variables on the dependent variables in terms of their significance.

Effect of Independent Variables on Particle Yield
The standardized effects of the independent variables are represented by the Pareto chart (Figure 3), which confirms the results of the analysis of variance (one way ANOVA) and shows that casein/fucoidan ratio (X1), spray intensity (X3), linear interaction between casein/fucoidan ratio, and glutaraldehyde concentration (X1X2), as well as the linear interaction between glutaraldehyde concentration and spray intensity (X2X3), have a statistically significant effect on the yield (p > 0.05) (Table 4).
increase in the casein/fucoidan ratio has a statistically significant effect (p < 0. standardized effects of the independent variables are represented by Pareto chart rank the effects of the independent variables on the dependent variables in terms significance.

Effect of Independent Variables on Particle Yield
The standardized effects of the independent variables are represented by th chart (Figure 3), which confirms the results of the analysis of variance (one way A and shows that casein/fucoidan ratio (X1), spray intensity (X3), linear interaction casein/fucoidan ratio, and glutaraldehyde concentration (X1X2), as well as the li teraction between glutaraldehyde concentration and spray intensity (X2X3), have tically significant effect on the yield (p > 0.05) (Table 4).Particle yield varied from (5.21 ± 2.2) % to (66.41 ± 2.2) % (Table 2).The mo gested that increasing the casein/fucoidan ratio and increasing the spray intensity decrease in the yield.A similar antagonistic effect on yield was also observed w simultaneous increase in the casein/fucoidan ratio and spray intensity.Furtherm model suggested that an increase in glutaraldehyde concentration and spray inten to an increase in the particle yield, as shown in Table 4.The two-way linear int between X1 and X2, as well as X2 and X3, demonstrated the change in the respo occurred when two independent variables were changed simultaneously.Positiv cients indicated a synergistic effect and negative coefficients indicated an antagon fect on the response [39].
The influence of the statistically significant independent variables, as well complex influence on the yield, is further presented through 3D surface plots (F The results prove that there is a tendency to increase the yield from (53.72 ± 1.9) % ± 2.2) %, when decreasing casein/fucoidan from 3:1 to 1:1 and decreasing glutara concentration from 3% (w/v) to 1% (w/v) (Figure 4a).Probably, at higher concentra casein and glutaraldehyde, the viscosity of the feed suspension increases, which an inefficient spray drying [40].Particle yield varied from (5.21 ± 2.2) % to (66.41 ± 2.2) % (Table 2).The model suggested that increasing the casein/fucoidan ratio and increasing the spray intensity led to a decrease in the yield.A similar antagonistic effect on yield was also observed with the simultaneous increase in the casein/fucoidan ratio and spray intensity.Furthermore, the model suggested that an increase in glutaraldehyde concentration and spray intensity led to an increase in the particle yield, as shown in Table 4.The two-way linear interaction between X1 and X2, as well as X2 and X3, demonstrated the change in the response that occurred when two independent variables were changed simultaneously.Positive coefficients indicated a synergistic effect and negative coefficients indicated an antagonistic effect on the response [39].
The influence of the statistically significant independent variables, as well as their complex influence on the yield, is further presented through 3D surface plots (Figure 4).The results prove that there is a tendency to increase the yield from (53.72 ± 1.9) % to (66.4 ± 2.2) %, when decreasing casein/fucoidan from 3:1 to 1:1 and decreasing glutaraldehyde concentration from 3% (w/v) to 1% (w/v) (Figure 4a).Probably, at higher concentrations of casein and glutaraldehyde, the viscosity of the feed suspension increases, which leads to an inefficient spray drying [40].
On the other hand, as the spray intensity decreases from 40% to 30% and the glutaraldehyde concentration increases from 1% (w/v) to 3% (w/v), there is a tendency toward an increase in the yield from (38.03 ± 1.8) % to (48.91 ± 2.9) % (Figure 4b).Probably, under these conditions, an effective cross-linking of the polyelectrolyte complex occurs, reducing the viscosity of the feed suspension, which in turn leads to its effective atomization at 30%.
Increasing the spray intensity leads to the generation of a larger number of droplets per unit time, which coalesce and form irregularly shaped structures [40,41].On the other hand, as the spray intensity decreases from 40% to 30% and the glutaraldehyde concentration increases from 1% (w/v) to 3% (w/v), there is a tendency toward an increase in the yield from (38.03 ± 1.8) % to (48.91 ± 2.9) % (Figure 4b).Probably, under these conditions, an effective cross-linking of the polyelectrolyte complex occurs, reducing the viscosity of the feed suspension, which in turn leads to its effective atomization at 30%.Increasing the spray intensity leads to the generation of a larger number of droplets per unit time, which coalesce and form irregularly shaped structures [40,41].
A linear regression analysis is applied to examine the fit of experimental to predicted data obtained in terms of particle yield [42].A linear polynomial equation is generated based on the yield regression coefficients: Yield = 207.54-2.686.X1 -6.4087.X2 − 0.4995.X3 + 0.3075X1.X2 + 0.0059X1.X3 − 0.3188.X2.X3 A correlation is observed between the experimental and the predicted data, which is confirmed by the coefficient of determination of the linear polynomial equation (R 2 = 0.99).The value for the adjusted coefficient of determination (R 2 adj = 0.99) confirms that 99% of the obtained data can be explained by the model.The observed experimental data as well as the data predicted by the model are presented in Table 5.One way ANOVA (one way ANOVA) was applied to estimate the significance of the model and its application for predicting the response at different levels of variation of the independent variables (Table 6).The obtained data present a significant model describing the influence of the independent variables, values Prob > F value and p < 0.05 (p = 0.002).A linear regression analysis is applied to examine the fit of experimental to predicted data obtained in terms of particle yield [42].A linear polynomial equation is generated based on the yield regression coefficients: Yield = 207.54− 2.686.X1 − 6.4087.X2 − 0.4995.X3 + 0.3075X1.X2 + 0.0059X1.X3 − 0.3188.X2.X3 A correlation is observed between the experimental and the predicted data, which is confirmed by the coefficient of determination of the linear polynomial equation (R 2 = 0.99).The value for the adjusted coefficient of determination (R 2 adj = 0.99) confirms that 99% of the obtained data can be explained by the model.The observed experimental data as well as the data predicted by the model are presented in Table 5.One way ANOVA (one way ANOVA) was applied to estimate the significance of the model and its application for predicting the response at different levels of variation of the independent variables (Table 6).The obtained data present a significant model describing the influence of the independent variables, values Prob > F value and p < 0.05 (p = 0.002).For the developed model, the possibility of observing constant variation and normal distribution of errors is assumed.To verify this claim, a normal probability plot of the residuals is generated in Figure 5a.The residuals are the differences in values between the experimental and predicted values of the response.The closer the residuals are to the theoretically expected, the more reliable the model is for predicting the response.From Figure 5a it can be clearly observed that all the experimental values are close to the theoretically expected (Figure 5a). Figure 5b is a plot representing the residuals versus the predicted response values.The results presented in Figure 5b show that the residuals are symmetrically located and tend to cluster in the low values of the Y-axis.Based on the obtained data, it can be assumed that the proposed model is adequate in terms of its ability to predict the response at different levels of variation of the independents.For the developed model, the possibility of observing constant variation and normal distribution of errors is assumed.To verify this claim, a normal probability plot of the residuals is generated in Figure 5a.The residuals are the differences in values between the experimental and predicted values of the response.The closer the residuals are to the theoretically expected, the more reliable the model is for predicting the response.From Figure 5a it can be clearly observed that all the experimental values are close to the theoretically expected (Figure 5a). Figure 5b is a plot representing the residuals versus the predicted response values.The results presented in Figure 5b show that the residuals are symmetrically located and tend to cluster in the low values of the Y-axis.Based on the obtained data, it can be assumed that the proposed model is adequate in terms of its ability to predict the response at different levels of variation of the independents.

Effect of Independent Variables on Particle Size
The standardized effects of the independent variables on particle size are represented by the Pareto chart (Figure 6), which confirms the results from the analysis of variance (Table 4) and shows that only the casein/fucoidan ratio has a statistically significant effect on the yield of the nanoparticles (p > 0.05), while the remaining variables (X1X2) and their combined effects (X1X2, X1X3, X2X3) do not have a statistically significant effect on particle size.

Effect of Independent Variables on Particle Size
The standardized effects of the independent variables on particle size are represented by the Pareto chart (Figure 6), which confirms the results from the analysis of variance (Table 4) and shows that only the casein/fucoidan ratio has a statistically significant effect on the yield of the nanoparticles (p > 0.05), while the remaining variables (X1X2) and their combined effects (X1X2, X1X3, X2X3) do not have a statistically significant effect on particle size.
The relationship between the independent and dependent variables is further represented by 3D surface plots (Figure 7).The figure visualizes the trend of the decrease in particle size from (1.515 ± 0.35) µm to (0.265 ± 0.03) µm with the decrease in the casein/fucoidan ratio from 3:1 to 1:1 and decrease in the glutaraldehyde concentration from 3% (w/v) to 1% (w/v).On one hand, probably, when equivalent amounts of the two biopolymers are used, an effective polyelectrolyte complexation occurs leading to the formation of more dense particles [43].On the other hand, low concentrations of glutaraldehyde are sufficient to saturate the free functional groups and ensure effective covalent cross-linking of the complex and the formation of smaller composite nanostructures [44].The relationship between the independent and dependent variables is furthe sented by 3D surface plots (Figure 7).The figure visualizes the trend of the dec particle size from (1.515 ± 0.35) µm to (0.265 ± 0.03) µm with the decrease in the ca coidan ratio from 3:1 to 1:1 and decrease in the glutaraldehyde concentration f (w/v) to 1% (w/v).On one hand, probably, when equivalent amounts of the two mers are used, an effective polyelectrolyte complexation occurs leading to the fo of more dense particles [43].On the other hand, low concentrations of glutaraldeh sufficient to saturate the free functional groups and ensure effective covalent cross of the complex and the formation of smaller composite nanostructures [44].The relationship between the independent and dependent variables is furth sented by 3D surface plots (Figure 7).The figure visualizes the trend of the dec particle size from (1.515 ± 0.35) µm to (0.265 ± 0.03) µm with the decrease in the ca coidan ratio from 3:1 to 1:1 and decrease in the glutaraldehyde concentration f (w/v) to 1% (w/v).On one hand, probably, when equivalent amounts of the two mers are used, an effective polyelectrolyte complexation occurs leading to the fo of more dense particles [43].On the other hand, low concentrations of glutaraldeh sufficient to saturate the free functional groups and ensure effective covalent cross of the complex and the formation of smaller composite nanostructures [44].The results were confirmed by the performed scanning and transmission microscopy.On the SEM micrographs (Figure 8A) it can be noticed that the part tained at 1:1 casein/fucoidan ratio have smooth surface and a spherical shape, a size ranged from (0.264 ± 0.03) µm (C1F1G1Sp.30)to (0.356 ± 0.01) µm (C1F1G2Sp The results were confirmed by the performed scanning and transmission electron microscopy.On the SEM micrographs (Figure 8A) it can be noticed that the particles obtained at 1:1 casein/fucoidan ratio have smooth surface and a spherical shape, and their size ranged from (0.264 ± 0.03) µm (C1F1G1Sp.30)to (0.356 ± 0.01) µm (C1F1G2Sp.40).The observations were also confirmed by transmission electron microscopy (Figure 8B), where spherical particles of a size around 200 nm were found.Increasing the casein/fucoidan ratio resulted in the formation of particle agglomerates with irregular shape and size ranging from 0.842 ± 0.26 µm (C3F1G2Sp.30)up to 1.04 ±0.98 µm (C3F1G1Sp.40)and 1.51 ± 0.35 µm (C3F1G3Sp.50).Probably, increasing the concentration of casein leads to inefficient polyelectrolyte complexation between the polymers, and an inefficient amount of glutaraldehyde needed for the formation of denser particles.This leads to an increase in the viscosity of the feeding suspension and an inefficient spray drying process.
Polymers 2024, 16, x FOR PEER REVIEW 14 of 24 observations were also confirmed by transmission electron microscopy (Figure 8B), where spherical particles of a size around 200 nm were found.Increasing the casein/fucoidan ratio resulted in the formation of particle agglomerates with irregular shape and size ranging from 0.842 ± 0.26 µm (C3F1G2Sp.30)up to 1.04 ±0.98 µm (C3F1G1Sp.40)and 1.51 ± 0.35 µm (C3F1G3Sp.50).Probably, increasing the concentration of casein leads to inefficient polyelectrolyte complexation between the polymers, and an inefficient amount of glutaraldehyde needed for the formation of denser particles.This leads to an increase in the viscosity of the feeding suspension and an inefficient spray drying process.
(A) Regarding the ability of the model to predict the response at different levels of variation of the independent variables, one way ANOVA is performed.Prob > F-value and p < 0.05 (p = 0.029) prove the ability of the model to predict the response (Table 8).To examine the variance and normal distribution of errors, a normal probability plot of residuals and a plot of residuals versus predicted values are generated Figure 9.The results show that all the experimental values are close to the theoretically expected ones (Figure 9a).Furthermore, the residuals are symmetrically located and tend to cluster in the low Y-axis values (Figure 9b).The obtained data confirms that the developed model can be used to predict the response at different levels of variation of the independent variables.Regarding the ability of the model to predict the response at different levels of variation of the independent variables, one way ANOVA is performed.Prob > F-value and p < 0.05 (p = 0.029) prove the ability of the model to predict the response (Table 8).To examine the variance and normal distribution of errors, a normal probability plot of residuals and a plot of residuals versus predicted values are generated Figure 9.The results show that all the experimental values are close to the theoretically expected ones (Figure 9a).Furthermore, the residuals are symmetrically located and tend to cluster in the low Y-axis values (Figure 9b).The obtained data confirms that the developed model can be used to predict the response at different levels of variation of the independent variables.

Effect of Independent Variables on Particle Zeta Potential
The standardized effects of the independent variables on the zeta potential of the particles are represented by Pareto chart (Figure 10), which shows that only the casein/fucoidan ratio has a statistically significant effect on the zeta potential of the nanoparticles (p > 0.05), while the other variables (X2 and X3) and their combined effect (X1X2, X1X3, X2X3) do not have a statistically significant effect.

Effect of Independent Variables on Particle Zeta Potential
The standardized effects of the independent variables on the zeta potential of the particles are represented by Pareto chart (Figure 10), which shows that only the casein/fucoidan ratio has a statistically significant effect on the zeta potential of the nanoparticles (p > 0.05), while the other variables (X2 and X3) and their combined effect (X1X2, X1X3, X2X3) do not have a statistically significant effect.
The relationship between the independent and the dependent variables is further represented by 3D surface plots (Figure 11).The figure shows a trend of an increase of the zeta potential from (−25.71 ± 1.0) mV to (−1.03 ± 1.9) mV with an increasing casein/fucoidan ratio from 1:1 to 3:1, which is probably due to the ineffective coating of the casein particles with fucoidan [37].The hypothesis is also confirmed by the larger size of these particles, which can be explained by the insufficient amount of fucoidan required for the stabilization of the obtained structures [37].The relationship between the independent and the dependent variables is fu represented by 3D surface plots (Figure 11).The figure shows a trend of an increase o zeta potential from (−25.71 ± 1.0) mV to (−1.03 ± 1.9) mV with an increasing casein/fuco ratio from 1:1 to 3:1, which is probably due to the ineffective coating of the casein par with fucoidan [37].The hypothesis is also confirmed by the larger size of these part which can be explained by the insufficient amount of fucoidan required for the stab tion of the obtained structures [37].Based on the obtained regression coefficients regarding the zeta potential, a l polynomial equation is derived as follows:  The relationship between the independent and the dependent variables is fur represented by 3D surface plots (Figure 11).The figure shows a trend of an increase o zeta potential from (−25.71 ± 1.0) mV to (−1.03 ± 1.9) mV with an increasing casein/fuco ratio from 1:1 to 3:1, which is probably due to the ineffective coating of the casein part with fucoidan [37].The hypothesis is also confirmed by the larger size of these parti which can be explained by the insufficient amount of fucoidan required for the stabi tion of the obtained structures [37].Based on the obtained regression coefficients regarding the zeta potential, a linear polynomial equation is derived as follows: Zeta potential = −83.33+ 1.24X1 −7.02X2 + 0.37X3 -0.04X1.X2 -0.01X1.X3 + 0.24X2.X3 A good correlation between experimental and predicted data is established (R 2 = 0.96), while the value of R 2 adj = 0.94 confirms that 94% of the obtained data can be explained by the model (Table 9).

Determining the Compatibility and Physical State of Casein/Fucoidan Composite Nanoparticles
To study the physical state and to understand the protein denaturation before and after functionalization with fucoidan, differential scanning calorimetry (DSC) was applied (Figure 13).After heating, casein exhibited a single endothermic denaturation transition peak (Tp) at 102.7 • C due to its dehydration and changes in the molecular integrity of the protein [45].
The first endothermic peak occurs due to the removal of physically absorbed water.The second exothermic peak is in the temperature range 200 °C-275 °C, which represents the second and the third decomposition steps occurring due to devolatilization processes.The last decomposition step usually occurs in the range 375 °C-1000 °C (outside the temperature range of the current study) and corresponds to oxidative degradation of the polymer [46].
The decomposition path of the synthesized casein-fucoidan composites followed a different pattern compared to the parent components.The degradation of the composites occurred as a two-stage process with two endothermic peaks observed in the DSC thermograms.The peak at 81.2 °C corresponds to moisture loss.The second peak at 213.3 °C is in the temperature range of the second and the third decomposition stages of fucoidan.However, transition from exothermic into endothermic process occurred, which could be related to an interaction between the components and formation of a conjugate of casein with fucoidan [47].To determine the interaction between casein and fucoidan molecules and the changes in the chemical structure during the formation of the Cas/FN composite nanostructures, FTIR analysis was conducted (Figure 14).Casein shows peaks at 1646 cm −1 in the amide I region and 1530 cm −1 in the amide II region, which could be assigned to the stretching of the carbonyl group (C=O) and to the symmetric stretching of N-C=O bonds, respectively [48].The decomposition path of the synthesized casein-fucoidan composites followed a different pattern compared to the parent components.The degradation of the composites occurred as a two-stage process with two endothermic peaks observed in the DSC thermograms.The peak at 81.2 • C corresponds to moisture loss.The second peak at 213.3 • C is in the temperature range of the second and the third decomposition stages of fucoidan.However, transition from exothermic into endothermic process occurred, which could be related to an interaction between the components and formation of a conjugate of casein with fucoidan [47].
To determine the interaction between casein and fucoidan molecules and the changes in the chemical structure during the formation of the Cas/FN composite nanostructures, FTIR analysis was conducted (Figure 14).Casein shows peaks at 1646 cm −1 in the amide I region and 1530 cm −1 in the amide II region, which could be assigned to the stretching of the carbonyl group (C=O) and to the symmetric stretching of N-C=O bonds, respectively [48].To determine the interaction between casein and fucoidan molecules and the changes in the chemical structure during the formation of the Cas/FN composite nanostructures, FTIR analysis was conducted (Figure 14).Casein shows peaks at 1646 cm −1 in the amide I region and 1530 cm −1 in the amide II region, which could be assigned to the stretching of the carbonyl group (C=O) and to the symmetric stretching of N-C=O bonds, respectively [48].The FTIR spectrum of fucoidan showed a strong absorption band at 1224 cm −1 , possibly related to the O=S=O vibration of sulfuric acid (Figure 14), as the band at 849 cm −1 was the characteristic band of the sulfate bond of fucoidan.Furthermore, the band at 1058 cm −1 was attributed to the C-O-C stretching vibration of the glycosidic bond [49].For the Cas/FN composite nanoparticles, a clear shift in the amide I region from 1646 to 1642 cm −1 and a shift in the amide II region from1530 to 1542 cm −1 were observed, due to the electrostatic interactions and hydrogen bonding between protein and polysaccharide [18] (Figure The FTIR spectrum of fucoidan showed a strong absorption band at 1224 cm −1 , possibly related to the O=S=O vibration of sulfuric acid (Figure 14), as the band at 849 cm −1 was the characteristic band of the sulfate bond of fucoidan.Furthermore, the band at 1058 cm −1 was attributed to the C-O-C stretching vibration of the glycosidic bond [49].For the Cas/FN composite nanoparticles, a clear shift in the amide I region from 1646 to 1642 cm −1 and a shift in the amide II region from1530 to 1542 cm −1 were observed, due to the electrostatic interactions and hydrogen bonding between protein and polysaccharide [18] (Figure 15).The FTIR spectrum of fucoidan showed a strong absorption band at 1224 cm −1 , possibly related to the O=S=O vibration of sulfuric acid (Figure 14), as the band at 849 cm −1 was the characteristic band of the sulfate bond of fucoidan.Furthermore, the band at 1058 cm −1 was attributed to the C-O-C stretching vibration of the glycosidic bond [49].For the Cas/FN composite nanoparticles, a clear shift in the amide I region from 1646 to 1642 cm −1 and a shift in the amide II region from1530 to 1542 cm −1 were observed, due to the electrostatic interactions and hydrogen bonding between protein and polysaccharide [18] (Figure 15).

Conclusions
Within the goal of this study, a novel, one-of-its-kind technological process for the preparation of casein-fucoidan composite nanostructures combining polyelectrolyte complexation and spray drying was developed.This research demonstrated that the main physicochemical properties of the obtained structures can be manipulated by various process parameters.The influence of the independent productional and technological factors was studied and optimized with the application of suitable surface response experimental design (3(k-p)).The results from the mathematical analysis showed that the casein to fucoidan ratio and spray intensity seemed to be crucial for the production yield of the obtained structure.Regarding the particle size and zeta potential, they were mainly affected by the casein to fucoidan ratio used for the preparation of the complexes.Therefore, the

Conclusions
Within the goal of this study, a novel, one-of-its-kind technological process for the preparation of casein-fucoidan composite nanostructures combining polyelectrolyte complexation and spray drying was developed.This research demonstrated that the main physicochemical properties of the obtained structures can be manipulated by various process parameters.The influence of the independent productional and technological factors was studied and optimized with the application of suitable surface response experimental design (3(k-p)).The results from the mathematical analysis showed that the casein to fucoidan ratio and spray intensity seemed to be crucial for the production yield of the obtained structure.Regarding the particle size and zeta potential, they were mainly affected by the casein to fucoidan ratio used for the preparation of the complexes.Therefore, the statistical experimental design methodology has clearly shown its usefulness for the optimization of the productional and technological parameters for the preparation of this innovative complex.Based on the results from the experimental design, the models obtained at casein/fucoidan ratio 1:1, cross-linked with 1% (w/v), 2% (w/v), and 3% (w/v) glutaraldehyde solution, and spray dried at 30%, 40%, and 50% spray intensity (C1F1G1Sp.30,C1F1G2Sp.40and C1F1G3Sp.50),demonstrating production yields varying from (48.91 ± 2.9) % to 66.41 ± 2.2, particle size ranging from (0.264 ± 0.03) µm to (0.356 ± 0.01) µm, and zeta potential from (−20.12 ± 0.9) mV to (−25.71 ± 1.0 mV) were selected as optimal models for further use as drug delivery systems for various bioactive substances.

Figure 3 .
Figure 3. Pareto chart reflecting the standardized effect of independents on yield.

Figure 3 .
Figure 3. Pareto chart reflecting the standardized effect of independents on yield.

Figure 4 .
Figure 4. 3D plot showing the effect of the casein to fucoidan ratio and glutaraldehyde concentration (a) and glutaraldehyde concentration and spray intensity (b) on particle yield.

Figure 4 .
Figure 4. 3D plot showing the effect of the casein to fucoidan ratio and glutaraldehyde concentration (a) and glutaraldehyde concentration and spray intensity (b) on particle yield.

Figure 5 .
Figure 5. Normal probability plot of residuals with respect to yield (a); plot of residuals versus predicted values (b).

Figure 5 .
Figure 5. Normal probability plot of residuals with respect to yield (a); plot of residuals versus predicted values (b).

Figure 6 .
Figure 6.Pareto chart showing the standardized effect of independents on particle size.

Figure 7 .
Figure 7. 3D surface plot showing the effect of casein to fucoidan ratio and glutaraldehyd tration on particle size.

Figure 6 .
Figure 6.Pareto chart showing the standardized effect of independents on particle size.

Figure 6 .
Figure 6.Pareto chart showing the standardized effect of independents on particle size.

Figure 7 .
Figure 7. 3D surface plot showing the effect of casein to fucoidan ratio and glutaraldehyd tration on particle size.

Figure 7 .
Figure 7. 3D surface plot showing the effect of casein to fucoidan ratio and glutaraldehyde concentration on particle size.

Figure 9 .
Figure 9. Normal probability plot of residuals with respect to size (a); plot of residuals versus predicted size values (b).

Figure 9 .
Figure 9. Normal probability plot of residuals with respect to size (a); plot of residuals versus predicted size values (b).

Figure 10 .
Figure 10.Pareto chart showing the standardized effect of independents on size.

Figure 11 .
Figure 11.3D surface plot showing the influence of casein to fucoidan ratio and glutaralde concentration on the zeta potential of the particles.

Figure 10 .
Figure 10.Pareto chart showing the standardized effect of independents on size.

Figure 10 .
Figure 10.Pareto chart showing the standardized effect of independents on size.

Figure 11 .
Figure 11.3D surface plot showing the influence of casein to fucoidan ratio and glutaralde concentration on the zeta potential of the particles.

Figure 11 .
Figure 11.3D surface plot showing the influence of casein to fucoidan ratio and glutaraldehyde concentration on the zeta potential of the particles.

Figure 13 .
Figure 13.DSC thermograms of casein (A), fucoidan (B), and casein-fucoidan composite nanostructures (C).Unlike casein, fucoidan shows two peaks-an endothermic peak at 78.8 • C and an exothermic peak at 208.9 • C. The thermal decomposition of fucoidan is a four-step process.The first endothermic peak occurs due to the removal of physically absorbed water.The second exothermic peak is in the temperature range 200 • C-275 • C, which represents the second and the third decomposition steps occurring due to devolatilization processes.The last decomposition step usually occurs in the range 375 • C-1000 • C (outside the temperature range of the current study) and corresponds to oxidative degradation of the polymer [46].

Figure 15 .
Figure 15.Schematic representation of the mechanism of formation of casein/fucoidan nanocomposites.

Figure 15 .
Figure 15.Schematic representation of the mechanism of formation of casein/fucoidan nanocomposites.

Table 1 .
Variables in the 3(k-p) experimental design.

Table 2 .
Compositions and main characteristics of the nine nanoparticle models.Cas

Table 3 .
Estimation of main effects of factors and coefficients for predictive mathematical models.

Table 4 .
Summary results of the performed analysis of variance (one way ANOVA) regarding the influence of the independent variables and their combined influence on the dependent.
SS-sum of squares; df-degree of freedom; MS-mean sum of squares; F-value; p-value.

Table 5 .
Experimental and model-predicted values of relative yields of particles (%).

Table 5 .
Experimental and model-predicted values of relative yields of particles (%).

Table 6 .
ANOVA significance testing of the proposed model.

Table 6 .
ANOVA significance testing of the proposed model.

Table 8 .
ANOVA significance testing of the proposed model.

Table 8 .
ANOVA significance testing of the proposed model.