Optimization of Phycobiliprotein Solubilization from a Thermotolerant Oscillatoria sp.

: The present study evaluated the effect of multiple variables (drying time, drying temperature, biomass/solvent ratio, glass beads/biomass ratio, extraction time, and extraction speed) in the solubilization of three different phycobiliproteins (C-PC, APC, and PE) from a thermotolerant Oscilla-toria sp. The strain was grown in BG11 media (28 ◦ C, light: dark cycle of 12:12 h at 100 µ mol · m − 2 · s − 1 , 20 days) and the experiments were conducted according to a two-level randomized factorial design with six center points (38 runs). Results show that biomass/solvent ratio, glass beads/biomass ratio, and extraction time, are the most signiﬁcant variables in the extraction of all three proteins, whereas the glass beads/biomass ratio and extraction time signiﬁcantly affect their purity. The optimized conditions allow a statistical increase in the concentration of C-PC, APC, and PE extracted from the biomass; however, the purity was lower in comparison with the expected value. The latter occurs due to a larger biomass/solvent ratio and longer extraction times, which enhanced the solubility of other hydrophilic metabolites (proteins and carbohydrates, etc.).


Introduction
The global demand for organic and sustainable food products has risen over the past decade.Furthermore, as a result of the expansion in the world's population, it is vital to investigate and identify alternative food sources [1] that comply with national and international regulations [2].Microalgae and cyanobacteria are now considered promising sources of natural molecules (such as carbohydrates, proteins, lipids, etc.) suitable for food, feed, cosmetic, and pharmaceutic applications [3][4][5][6][7][8][9][10][11].Among those components, phycobiliproteins (PBPs) have gained popularity in the food and pharmaceutical sectors [12].PBPs are a group of brilliant-colored proteins responsible for photosynthetic activity in different microbial groups such as cyanobacteria, cryptophytes, red algae, and glaucocystophytes [2].PBPs are composed of several classes called phycocyanin (C-PC), allophycocyanin (APC), and phycoerythrin (PE).Each of these classes possesses a unique spectrum referred to as blue (610−620 nm), blue-green (650−655 nm), and pink (540−570 nm) [3,7,13].Due to its antioxidant and anti-inflammatory properties, and its flexibility to be used by various industries such as the medical, pharmaceutical, cosmetic, and food industries, phycocyanin (C-PC) is one of the most expensive dye proteins on the market, with an estimated value of 5000-33,000 USD per gram [8].C-PC has shown great success in the food sector as a replacement for its synthetic counterparts such as Patent Blue V and Brilliant Blue FCF (E133) [12].
The cyanobacterial and algal biomass is known for its high humidity content (up to 80%) biomass.This high humidity can significantly reduce its shelf life and thus reduce the quality of the metabolites of interest [16].Therefore, it is critical to remove excess water totally (or partially) [24].In the case of phycobiliproteins, temperatures higher than 45 • C degrees can reduce the concentration and quality of these pigments by up to 50% (w/w) [25].However, as pointed out by Pez Jaeschke et al. [1], the information available on the use of biomass form (fresh, frozen, oven-dried, or freeze-dried) is still unclear.Another factor that influences the extraction efficiency is the destruction time and speed.Since the extraction solvent employed is an aqueous buffer, other hydrophilic substances such as carbohydrates and total proteins will be diluted into the buffer, affecting PBPs' final purity [26].Moreover, longer extraction times will eventually increase the system's temperature, which will degrade the PBPs [12].
Since most of the works published in recent years focused on the analysis of one or two factors, the objective of this study was to evaluate for the first time the effect of multiple factors (drying, biomass/solvent ratio, biomass/glass beads ratio, mixing speed and time) on the extraction yield and purity of phycobiliproteins (C-CP, APC, and PE) from a thermotolerant strain of Oscillatoria sp. with a high concentration of those colorant proteins.

Experimental Design
The evaluation of six factors that affect the extraction yield and purity of PBPs were studied using a 2-level randomized factorial design with six center points (38 runs), using Design-Expert ® software (version 13.0, Stat-Ease Inc., Minneapolis, MN, USA).Table 1 shows the levels studied for each factor.

Culture Conditions
For each experiment, Oscillatoria sp. was cultured (in triplicate) in 500 mL flasks with a working volume of 250 mL of BG-11 culture media.Each flask was mixed through the injection of filtered air with 1% (v/v) CO 2 at a flow rate of 0.15 L air •min −1 , with a light:dark cycle of 12:12 h at 100 µmol•m −2 •s −1 , 27 ± 1 • C, for 20 days.

Biomass and PBPs Quantification
After 20 days, each flask was removed from the air inlet and allowed to self-flocculate, due to the filamentous morphology of Oscillatoria sp., the cyanobacteria precipitate in 20 min.The concentrated biomass was transferred to a 500 mL plastic flask and centrifuged at 3500 rpm (20 • C, 20 min), and the supernatant was withdrawn.The pellet was poured into non-stick food-grade silicone molds.The different samples were dried using a foodgrade dehydrator (Hamilton Beach ® 32100a).The mass obtained was recorded using a digital balance.The dehydrated biomass was subjected to extraction (one sample plus 5 replicates) using the method described by Zuorro et al. [28].Briefly, a known amount of dried biomass was mixed with a volume of cold phosphate buffer solution (0.05 M, pH 6.8) and a known amount of glass beads (0.5 mm diameter) according to the experimental design.The solution was mixed using an automatic vortex (Multi Reax, Heidolph, Germany).The mixture was stored in a refrigerator to promote the solubilization of the phycobiliproteins (4 • C, 24 h).PBPs were separated from cell debris by centrifugation (3400 rpm, 30 min, 20 • C).The deep blue supernatant was collected and measured in a spectrophotometer at different wavelengths (620, 652, 562, and 280 nm).The concentration of C-PC, APC, and PE was calculated using the Equations ( 1)-( 3) described by Bennett and Bogorad [29].The purity of C-PC, APC, and PE was determined using the Equations ( 4)-( 6) proposed by Patil [30] and Antello et al. [31].The average of the results obtained for each experiment was used to perform the ANOVA analysis according to the Design-Expert ® software.

Effect of Multiple Parameters on the Concentration of PBPs
The results of the ANOVA analysis (Table 2) on the effect of multiple variables (drying time (A), drying temperature (B), biomass/solvent ratio (C), glass beads/biomass ratio (D), extraction time (E), and extraction speed (F)) in the extraction of phycobiliproteins (C-PC, APC, and PE) from Oscillatoria sp.biomass shows that in the case of C-PC, the F-value obtained (85.78) implies the model is significant.There is only a 0.01% chance that an F-value this large could occur due to noise; therefore, the results obtained did not happen by chance.
The p-values obtained for the different variables analyzed show that drying time, drying temperature, biomass/solvent ratio, glass beads/biomass ratio, extraction time, and extraction speed (A, B, C, D, E, and F, respectively) are statistically significant (p-value < 0.05) variables that affect the extraction efficiency.Another interesting result obtained is that several interactions between the analyzed variables (AB, BC, BE, BF, CE, DF, EF, ABD, ABE, ABF, ACD, ADE, ADF, and AEF) are significant model terms.The Lack of Fit F-value of 1.32 implies the Lack of Fit is not substantial relative to the pure error, which means that the model fits.
In the case of allophycocyanin (APC), an F-value of 13.38 implies the model is significant, with only a 0.01% chance that an F-value this large could occur due to noise.In the extraction of APC, the statistically significant variables (p-value < 0.05) were biomass/solvent ratio, glass beads/biomass ratio, and extraction time (C, D, and E, respectively).Moreover, different interactions between the variables were found to be statistically significant (AB, BC, CE, DF, ABD, ABF, ACD, ADE, and ADF).The Lack of Fit F-value of 0.39 implies the Lack of Fit is not significant relative to the pure error, which means that the model fits.
Finally, the analysis of the extraction of phycoerythrin (PE) shows that an F-value of 18.54 implies the model is significant.There is only a 0.01% chance that an F-value this large could occur due to noise.Drying time, biomass/solvent ratio, glass beads/biomass ratio, extraction time, and extraction speed (A, C, D, E, and F, respectively) and several interactions (BF and ADE) are the most important terms affecting the extraction of the colorant protein PE.The Lack of Fit F-value of 3.73 implies the Lack of Fit is not significant relative to the pure error, which means that the model fits.

Effect of Multiple Parameters on the Purity of PBPs
One of the most critical responses in the extraction of phycobiliproteins is their purity.Since PBPs are proteins, the selected extraction method will inevitably allow the solubilization of other groups of proteins that will reduce the quality of the extract.The results of the ANOVA analysis (Table 3) on the purity of PBPs shows that, in the case of C-PC, the F-value obtained (16.17) for this protein implies the model is significant.There is only a 0.04% chance that an F-value this large could occur due to noise.In this case, the purity of C-CP is affected only by glass beads/biomass ratio, extraction time (D, E, respectively), and several interactions (AB, AE, BC, BE, BF, EF, ABC, ABD, ABE, ABF, ACD, ADE, and ADF); in this case, these terms are considered to be significant model terms.The Lack of Fit F-value of 0.09 implies the Lack of Fit is not significant relative to the pure error, which means that the model fits.The different surface responses (three-dimensional plots) obtained from the extraction of C-PC, APC, and PE are shown in Figure 1.According to the results, a high biomass/solvent ratio of up to 0.5% w/v coupled with 30 h of drying using a food-grade dehydrator maximizes the extraction of the studied PBPs.In the case of APC, an F-value of 5.37 implies the model is significant, with only a 0.21% chance that an F-value this large could occur due to noise.According to the p-values obtained for the different model terms, only glass beads/biomass ratio, extraction time (D and E, respectively), and several interactions between the analyzed variables are significant model terms (AB, AE, BC, BF, ABC, ABD, ABF, ACD, and ADE).The Lack of Fit F-value of 0.05 implies the Lack of Fit is not significant relative to the pure error, which means that the model fits.

Effect of Multiple Parameters on the Purity of PBPs
Finally, the analysis of the extraction of PE shows that an F-value of 13.75 implies the model is significant.There is only a 0.01% chance that an F-value this large could occur due to noise.Drying time, glass beads/biomass ratio, extraction time (A, D, and E, respectively), and several interactions (AB, AE, BC, BE, BF, DE, ABC, ABD, ABE, ABF, ACD, ADE, and ADF) are the most significant terms affecting the extraction of the colorant protein PE.The Lack of Fit F-value of 0.36 implies the Lack of Fit is not significant relative to the pure error, which means that the model fits.
The surface responses (three-dimensional plots) of the model equation fitted to the data on the purity index of the extracted C-PC, APC, and PE are shown in Figure 2. Unlike the extraction results, the purity of PBPs is favored with a lower biomass/solvent ratio of 0.2 % w/v (or lower), and 30 h of drying using a food-grade dehydrator can effectively enhance the purity of the extracted PBPs.
Table 4 represents the highest scenario that maximizes C-PC, APC, and PE concentration and purity using all the different variables analyzed.To test the viability of the proposed method, Oscillatoria sp. was grown in three 2 L flasks (1.2 L of working volume, 0.72 L•min −1 of filtered air, 12:12 h light:dark cycle at 100 µmol•m −2 •s −1 , 27 ± 1 • C, for 20 days).The biomass obtained was prepared according to Table 3, with 12 samples (one experiment with 11 replicates).
The surface responses (three-dimensional plots) of the model equation fitted to the data on the purity index of the extracted C-PC, APC, and PE are shown in Figure 2. Unlike the extraction results, the purity of PBPs is favored with a lower biomass/solvent ratio of 0.2 % w/v (or lower), and 30 h of drying using a food-grade dehydrator can effectively enhance the purity of the extracted PBPs.The results obtained for the concentration and purity of C-PC, APC, and PE were analyzed using an unpaired t-test (Figure 3).The analysis shows that the concentration of C-PC, APC, and PE under the conditions presented in Table 4 were statistically higher than expected (<0.0001).However, the purity of the obtained extracts was lower than expected.The latter indicates that the proposed factors can effectively enhance the solubility of C-PC, APC, and PE; however, they also increase the solubility of other proteins (globular proteins, etc.), which in turn reduce the purity of PBPs.
Typically, most studies in this area evaluate the effect of two or a maximum o variables at a time.Moreover, most of these works employ One Factor at a Time (O
Typically, most studies in this area evaluate the effect of two or a maximum of four variables at a time.Moreover, most of these works employ One Factor at a Time (OFAT) experimental processes, increasing the experimentation time required.Design of Experiments (DoE) is a series of statistical tools that improve the number of the experiments needed to obtain a statistically significant response [32][33][34][35]; in recent years, its application in biotechnology has proven to be a valuable tool to reduce experimental costs and improve the quality of the results [36][37][38][39][40].One example of the application of DoE in the optimization of PBP extraction is the work reported by Sintra et al. [11], where the researchers evaluated the effect of solid-liquid ratio, time of extraction, and temperature in the extraction of C-PC and chlorophylls in Anabaena cylindrica using the response surface methodology (RSM).In another study, Ferreira-Santo et al. [10] found that the application of ohmic heating (OH) at shorter times (44 • C and 30 min, respectively) allowed a significant increase in C-phycocyanin extraction yield compared to conventional heating and freeze-thawing.Another example is the experiment reported by Tavanandi et al. [18], where they analyzed the effect of ultrasonication with thawing and maceration.Other more advanced research, such as the work of Ilter et al. [14], examined the impact of up to five variables (homogenization rate (rpm), amplitude (%), microwave power (W), biomass/solvent ratio (%), and extraction time of classical, ultrasound and microwave extraction) in the extraction of C-PC using central composite rotatable design (CCRD).
In the case of the PBP extraction, the common factors affecting the efficiency of C-PC, APC, and PE were biomass/solvent ratio (C), glass beads/biomass ratio (D), and extraction time (E).By comparison, the most repetitive interaction found was drying time-glass beads/biomass ratio-extraction time (ADE).In the case of the PBP purity, the common factors affecting the efficiency of C-PC, APC, and PE were glass beads/biomass ratio (D) and extraction time (E).In addition, the most repetitive interactions found were AB, AE, BC, BF, ABC, ABD, ABF, ACD, and ADE, where drying time is the most repetitive factor (A).One exciting result obtained from the ANOVA analysis for the extraction of PBPs is that only C-PC and PE showed a fit for the curvature model.The latter may occur due to the high content of both proteins in comparison with APC in the thermotolerant strain of Oscillatoria sp.The same result was obtained in analyzing the purity of all three PBPs.
In each of the papers mentioned previously, the authors employed freeze-dried or spray-dried biomass.Since PBPs are sensitive to high temperatures (>45 • C) [25], these drying processes are preferred.To the best of the author's knowledge, this is the first time a food dehydrator has been used for the drying process of cyanobacterial biomass for PBP extraction.A food dehydrator is an inexpensive piece of equipment that can be acquired in almost every country.This research showed excellent results in maintaining the quantity and quality of PBPs while removing the moisture within the cyanobacterial biomass using longer times and mild temperatures (30 h and 40 • C, respectively).Another poorly evaluated variable is the glass beads/biomass ratio; the latter occurs because when glass beads are employed, they are commonly added in excess to better destroy cyanobacterial or algal biomass.
This research shows that a glass beads/biomass ratio of 15% w/w enhances biomass contact with the beads, allowing better PBP solubilization.Another variable worth mentioning is the destruction time; according to the results, longer times increase the content of PBPs extracted from the cells but reduce their purity.These results are congruent with those reported by Pott et al. [26], who found that longer times will eventually solubilize other proteins and hydrophilic metabolites such as carbohydrates.Another variable that affects the extraction time is the biomass/solvent ratio.According to the results, larger proportions of biomass in the buffer will increase the final content of PBPs but decrease their purity.One possible explanation for this behavior between PBP concentration and purity can be found in the drying time-glass beads/biomass ratio-extraction time (ADE) interaction; however, it is necessary to explore this interaction using other experimental Processes 2022, 10, 836 10 of 12 design methods, such as central composite design, to determine if both the concentration and purity of PBPs can be optimized in parallel.

Conclusions
This study indicates that the Design of Experiments is a robust methodology for identifying critical variables that enhance the extraction and purity of PBPs in Oscillatoria sp.The obtained results prove that, in the case of the PBP extraction, the common factors are biomass/solvent ratio, glass beads/biomass ratio, and extraction time, whereas the factors affecting purity are glass beads/biomass ratio and extraction time.The optimized conditions allowed the concentration of extracted C-PC, APC, and PE to be significantly increase; however, the purity of the proteins was lower than expected.The latter occurs due to a larger biomass/solvent ratio and longer extraction times, which enhance the solubility of other hydrophilic metabolites (proteins and carbohydrates, etc.).Future studies should provide further insight into the optimization of the interaction between the drying time, glass beads/biomass ratio, and extraction time variables.

Figure 1 .
Figure 1.Surface response (three-dimensional plots) of the model equation fitted to the data on extraction concentration (g/L) of C-PC (a), APC (b), and PE (c).

Figure 1 .
Figure 1.Surface response (three-dimensional plots) of the model equation fitted to the data on extraction concentration (g/L) of C-PC (a), APC (b), and PE (c).

Figure 2 .
Figure 2. Surface response (three-dimensional plots) of the model equation fitted to the data on purity of C-PC (a), APC (b), and PE (c).

Figure 3 .
Figure 3. Unpaired t-test between the expected and expected results from extraction (a-c) and (d-f) of C-PC, APC, and PE.

Table 1 .
Variables evaluated with their levels for the extraction of PBPs.

Table 2 .
Analysis of variance (ANOVA) of the model obtained for C-PC, APC, and PE extraction.

Table 3 .
Analysis of variance (ANOVA) of the model obtained for phycobiliprotein purity.

Table 4 .
Variables for optimal biomass concentration on both strains were studied.