1. Introduction
In recent decades, biopolymers have been a subject of a serious scientific interest for the development of nanosized drug delivery systems, due to their specific characteris-tics 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 nanoparticu-late 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 ag-glomerates 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 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 over-come 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 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 its 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 of 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 ca-sein-dextran 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, which 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, 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 protein-polysaccharide bonds, which largely depends on the nature of the biopolymers, their quantitative ratio, pH of the medium and ionic strength [
20]. Regarding physical protein-polysaccharide complexes, non-covalent interactions can be due to hydrogen bonding, 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 ca-sein-polysaccharide 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 acidic pH (2.5), while at pH 7 up to 90% of the incorporated nutraceutical was re-leased, demonstrating the potential use of the carrier as a drug delivery system for targeted delivery in the small intestine. In another study, drug delivery system for epigal-locatechin-3-gallate was obtained via polyelectrolyte complexation between casein and chitosan [
24]. The obtained nanostructures demonstrated average particle size of 150 nm and 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, increase of their average particle size from 264 nm to 294 nm was observed. Furthermore, the drug delivery system showed 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 ca-sein-polysaccharide 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).
Type I fucoidan is composed of an α-(1→3)-α-L-fucopyranose backbone linked to sulfate radicals at the C2 and C4 positions, while type II fucoidan is composed of α-(1→3)- and α -( 1→4)-α-L-fucopyranose backbone, with sulfate radicals attached at the C2, C3, and C4 positions [
28,
29].
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 crosslinking 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.
Figure 1.
Fucoidan structure: α-(1→3)-α-L-fucopyranose skeleton linked with sulfate radicals at C2 and C4 position (Type I) and α-(1→3)- and α-(1→4) )-α-L-fucopyranose skeleton with sulfate radicals attached to C2, C3 and C4 positions (Type II) [
27].
Figure 1.
Fucoidan structure: α-(1→3)-α-L-fucopyranose skeleton linked with sulfate radicals at C2 and C4 position (Type I) and α-(1→3)- and α-(1→4) )-α-L-fucopyranose skeleton with sulfate radicals attached to C2, C3 and C4 positions (Type II) [
27].
Figure 2.
Schematic representation of the preparation of casein/fucoidan composite nanoparticles (Generated using BioRender™).
Figure 2.
Schematic representation of the preparation of casein/fucoidan composite nanoparticles (Generated using BioRender™).
Figure 3.
Pareto chart reflecting the standardized effect of independents on yield.
Figure 3.
Pareto chart reflecting the standardized effect of independents on yield.
Figure 4.
3D plot showing the effect of casein to fucoidan ratio and glutaraldehyde concentration on particle yield.
Figure 4.
3D plot showing the effect of casein to fucoidan ratio and glutaraldehyde concentration on particle yield.
Figure 5.
Normal probability plot of residuals with respect to yield (a); Plot of residuals versus predicted values (b).
Figure 5.
Normal probability plot of residuals with respect to yield (a); Plot of residuals versus predicted values (b).
Figure 6.
Pareto chart showing the standardized effect of independents on particle size.
Figure 6.
Pareto chart showing the standardized effect of independents on particle size.
Figure 7.
3D surface plot showing the effect of casein to fucoidan ratio and glutaraldehyde concentration on particle size.
Figure 7.
3D surface plot showing the effect of casein to fucoidan ratio and glutaraldehyde concentration on particle size.
Figure 8.
SEM micrographs of casein/fucoidan composite nanostructures of different models obtained at 10 000x magnification, accelerating voltage 20 kV using an Everhart–Thornley detector.
Figure 8.
SEM micrographs of casein/fucoidan composite nanostructures of different models obtained at 10 000x magnification, accelerating voltage 20 kV using an Everhart–Thornley detector.
Figure 9.
Normal probability plot of residuals with respect to size (a); Plot of residuals versus predicted size values (b).
Figure 9.
Normal probability plot of residuals with respect to size (a); Plot of residuals versus predicted size values (b).
Figure 10.
Pareto chart showing the standardized effect of independents on size.
Figure 10.
Pareto chart showing the standardized effect of independents on size.
Figure 11.
3D surface plot showing the influence of casein to fucoidan ratio and glutaraldehyde concentration on the zeta potential of the particles.
Figure 11.
3D surface plot showing the influence of casein to fucoidan ratio and glutaraldehyde concentration on the zeta potential of the particles.
Figure 12.
Normal probability plot of residuals versus zeta potential (a); plot of residuals versus predicted zeta values (b).
Figure 12.
Normal probability plot of residuals versus zeta potential (a); plot of residuals versus predicted zeta values (b).
Figure 13.
DSC thermograms of casein (A), fucoidan (B) and casein-fucoidan composite nanostructures (C).
Figure 13.
DSC thermograms of casein (A), fucoidan (B) and casein-fucoidan composite nanostructures (C).
Figure 14.
FTIR spectra of casein, fucoidan and casein-fucoidan composite nanostructures.
Figure 14.
FTIR spectra of casein, fucoidan and casein-fucoidan composite nanostructures.
Table 1.
Variables in the 3 k-p experimental design.
Table 1.
Variables in the 3 k-p experimental design.
Variables |
Levels |
Independent variables |
-1 |
0 |
+1 |
X1: Casein:Fucoidan ratio |
1:1 |
2:1 |
3:1 |
X2: Glutaraldehyde concentration (%) |
3 |
2 |
1 |
X3: Spray intensity (%) |
30 |
40 |
50 |
Dependent variables |
|
|
|
Y1: Yield (%) |
|
|
|
Y2: Particle size (µm) Y3: Zeta potential (mV) |
|
|
|
Table 2.
Compositions and main characteristics of the nine nanoparticle models. Cas/FN - casein to fucoidan ratio; GLD - glutaraldehyde concentration, Sp.int. – Spray intensity.
Table 2.
Compositions and main characteristics of the nine nanoparticle models. Cas/FN - casein to fucoidan ratio; GLD - glutaraldehyde concentration, Sp.int. – Spray intensity.
Models |
Independent variables |
Dependent variables |
Cas/FN |
GLD (%) |
Sp.int.(%)
|
Size ± SD (µm) |
ζ ± SD (mV) |
Yield ± SD (%) |
C1F1G1Sp.30 |
-1 |
-1 |
-1 |
264 ± 32 |
-25.71 ± 1.0 |
66.41 ± 2.2 |
C1F1G2Sp.40 |
-1 |
0 |
+1 |
356 ± 11 |
-20.12 ± 0.9 |
48.91 ± 2.9 |
C1F1G3Sp.50 |
-1 |
+1 |
0 |
282 ± 25 |
-23.72 ± 1.8 |
53.72 ± 1.9 |
C2F1G1Sp.50 |
0 |
-1 |
+1 |
742 ± 252 |
-10.21 ± 1.5 |
20.66 ± 1.5 |
C2F1G2Sp.40 |
0 |
0 |
0 |
512 ± 360 |
-13.43 ± 1.6 |
27.81 ± 2.6 |
C2F1G3Sp.30 |
0 |
+1 |
-1 |
484 ± 121 |
-15.51 ± 1.4 |
38.03 ± 1.8 |
C3F1G1Sp.40 |
+1 |
-1 |
0 |
1043 ± 986 |
-3.23 ± 0.7 |
7.31 ± 3.4 |
C3F1G2Sp.30 |
+1 |
0 |
-1 |
842 ± 263 |
-5.65 ± 0.9 |
18.54 ± 2.7 |
C3F1G3Sp.50 |
+1 |
+1 |
+1 |
1515 ± 353 |
-1.03 ± 1.9 |
5.21 ± 2.2 |
Table 3.
Estimation of main effects of factors and coefficients for predictive mathematical models.
Table 3.
Estimation of main effects of factors and coefficients for predictive mathematical models.
Factor |
Yield (%) |
Particle size (µm) |
ζ (mV) |
Effect |
SE |
p |
β |
Effect |
SE |
p |
β |
Effect |
SE |
p |
β |
Mean/Int. |
34.59 |
0.38 |
0.00 |
207.54 |
0.62 |
0.06 |
0.01 |
0.34 |
-14.35 |
0.66 |
0.00 |
-83.330 |
Cas/Fn (1) |
-45.86 |
0.92 |
0.00 |
-2.68 |
0.80 |
0.15 |
0.03 |
0.01 |
19.64 |
1.59 |
0.00 |
1.244 |
GLD (2) |
0.11 |
0.96 |
0.91 |
-6.40 |
0.09 |
0.16 |
0.63 |
-0.68 |
-0.15 |
1.67 |
0.78 |
-7.010 |
Sp.Int. (3) |
-15.32 |
0.96 |
0.00 |
-0.49 |
0.32 |
0.16 |
0.18 |
-0.04 |
5.31 |
1.67 |
0.08 |
0.374 |
1 by 2 |
7.86 |
1.36 |
0.02 |
0.30 |
0.00 |
0.23 |
0.97 |
0.00 |
-1.08 |
2.35 |
0.69 |
-0.040 |
1 by 3 |
1.48 |
1.36 |
0.38 |
0.00 |
0.10 |
0.23 |
0.69 |
0.00 |
-2.32 |
2.35 |
0.42 |
-0.010 |
2 by 3 |
-6.37 |
1.41 |
0.04 |
-0.31 |
0.37 |
0.24 |
0.26 |
0.01 |
4.72 |
2.44 |
0.19 |
0.236 |
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.
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.
Factor |
SS |
df |
MS |
F |
p |
Yield: R2=0.99; R2Adj.=0.99; MS=1.30 |
Cas/Fn (1) |
3188.93 |
1 |
3188.93 |
2451.84 |
0.00 |
GLD (2) |
0.01 |
1 |
0.01 |
0.01 |
0.91 |
Spray.Int (3) |
326.21 |
1 |
326.21 |
250.81 |
0.00 |
1 by 2 |
41.43 |
1 |
41.43 |
31.86 |
0.02 |
1 by 3 |
1.54 |
1 |
1.54 |
1.18 |
0.38 |
2 by 3 |
26.25 |
1 |
26.25 |
20.18 |
0.04 |
Error |
2.60 |
2 |
|
|
|
Total SS |
3632.52 |
8 |
|
|
|
Particle size: R2=0.94; R2Adj.=0.77; MS=0.03 |
Cas/Fn (1) |
0.99 |
1 |
0.99 |
26.40 |
0.03 |
GLD (2) |
0.01 |
1 |
0.01 |
0.30 |
0.63 |
Spray.Int (3) |
0.14 |
1 |
0.14 |
3.95 |
0.18 |
1 by 2 |
0.00 |
1 |
0.00 |
0.00 |
0.97 |
1 by 3 |
0.00 |
1 |
0.00 |
0.20 |
0.69 |
2 by 3 |
0.09 |
1 |
0.09 |
2.40 |
0.26 |
Error |
0.07 |
2 |
|
|
|
Total SS |
1.36 |
8 |
|
|
|
Zeta potential: R2=0.98; R2Adj.=0.95; MS=3.87 |
Cas/Fn (1) |
585.15 |
1 |
585.15 |
150.98 |
0.00 |
GLD (2) |
0.37 |
1 |
0.37 |
0.09 |
0.78 |
Spray.Int (3) |
39.31 |
1 |
39.31 |
10.14 |
0.08 |
1 by 2 |
0.82 |
1 |
0.81 |
0.21 |
0.69 |
1 by 3 |
3.79 |
1 |
3.79 |
0.97 |
0.42 |
2 by 3 |
14.44 |
1 |
14.44 |
3.72 |
0.19 |
Error |
7.75 |
2 |
3.87 |
|
|
Total SS |
633.89 |
8 |
|
|
|
Table 5.
Experimental and model-predicted values of relative yields of particles (%).
Table 5.
Experimental and model-predicted values of relative yields of particles (%).
Batch |
Observed |
Predicted |
Residuals |
Error (%) |
C1F1G1Sp.30 |
66.41 |
66.52 |
-0.11 |
-0.16 |
C1F1G2Sp.40 |
48.91 |
49.12 |
-0.21 |
-0.42 |
C1F1G3Sp.50 |
53.72 |
53.74 |
-0.02 |
-0.03 |
C2F1G1Sp.50 |
20.66 |
20.69 |
-0.03 |
-0.14 |
C2F1G2Sp.40 |
27.81 |
26.34 |
1.46 |
5.28 |
C2F1G3Sp.30 |
38.03 |
38.36 |
-0.33 |
-0.86 |
C3F1G1Sp.40 |
7.31 |
7.74 |
-0.45 |
-5.88 |
C3F1G2Sp.30 |
18.54 |
18.58 |
-0.04 |
-0.21 |
C3F1G3Sp.50 |
5.21 |
5.49 |
-0.24 |
-5.37 |
Table 7.
Experimental and model-predicted values of relative particle size (µm).
Table 7.
Experimental and model-predicted values of relative particle size (µm).
Batch |
Observed |
Predicted |
Residuals |
Error (%) |
C1F1G1Sp.30 |
0.264 |
0.243 |
0.02 |
-0.21 |
C1F1G2Sp.40 |
0.356 |
0.329 |
0.02 |
-0.23 |
C1F1G3Sp.50 |
0.282 |
0.267 |
0.01 |
-0.17 |
C2F1G1Sp.50 |
0.742 |
0.719 |
0.02 |
-0.28 |
C2F1G2Sp.40 |
0.512 |
0.768 |
-0.25 |
2.60 |
C2F1G3Sp.30 |
0.484 |
0.442 |
0.04 |
-0.37 |
C3F1G1Sp.40 |
1.043 |
0.986 |
0.05 |
-0.50 |
C3F1G2Sp.30 |
0.842 |
0.811 |
0.03 |
-0.38 |
C3F1G3Sp.50 |
1.515 |
1.471 |
0.04 |
-0.44 |
Table 8.
ANOVA significance testing of the proposed model.
Table 8.
ANOVA significance testing of the proposed model.
R2
|
R2adj
|
SS |
df |
MS |
SS |
df |
MS |
F |
p |
0.810 |
0.696 |
1.102 |
3 |
0.367 |
0.258 |
5 |
0.0517 |
7.105 |
0.0297 |
Table 9.
Experimental and model-predicted values of particle zeta potential (mV).
Table 9.
Experimental and model-predicted values of particle zeta potential (mV).
Batch |
Observed |
Predicted |
Residuals |
Error (%) |
C1F1G1Sp.30 |
-25.71 |
-25.91 |
0.20 |
-0.81 |
C1F1G2Sp.40 |
-20.12 |
-20.35 |
0.23 |
-1.18 |
C1F1G3Sp.50 |
-23.72 |
-23.89 |
0.17 |
-0.75 |
C2F1G1Sp.50 |
-10.21 |
-10.49 |
0.28 |
-2.75 |
C2F1G2Sp.40 |
-13.43 |
-10.82 |
2.60 |
-19.4 |
C2F1G3Sp.30 |
-15.51 |
-15.88 |
0.37 |
-2.39 |
C3F1G1Sp.40 |
-3.23 |
-3.73 |
0.50 |
-15.63 |
C3F1G2Sp.30 |
-5.65 |
-6.03 |
0.38 |
-6.75 |
C3F1G3Sp.50 |
-1.03 |
-1.47 |
0.44 |
-43.05 |
Table 10.
ANOVA significance testing of the proposed model.
Table 10.
ANOVA significance testing of the proposed model.
R2
|
R2adj.
|
SS |
df |
MS |
SS |
df |
MS |
F |
p |
0.9642 |
0.9642 |
611.22 |
3 |
203.74 |
22.66 |
5 |
4.533 |
44.94 |
0.00048 |