Bioprocessing optimization for efficient simultaneous removal of methylene blue and nickel by Gracilaria seaweed biomass

The pollution of water by heavy metal ions and dyes, particularly from industrial effluents, has become a global environmental issue. Therefore, the treatment of wastewater generated from different industrial wastes is essential to restore environmental quality. The efficiency of Gracilaria seaweed biomass as a sustainable biosorbent for simultaneous bioremoval of Ni2+ and methylene blue from aqueous solution was studied. Optimization of the biosorption process parameters was performed using face-centered central composite design (FCCCD). The highest bioremoval percentages of Ni2+ and methylene blue were 97.53% and 94.86%; respectively, obtained under optimum experimental conditions: 6 g/L Gracilaria biomass, initial pH 8, 20 mg/L of methylene blue, 150 mg/L of Ni2+ and 180 min of contact time. Fourier Transform Infrared Spectroscopy (FTIR) spectra demonstrated the presence of methyl, alkynes, amide, phenolic, carbonyl, nitrile and phosphate groups which are important binding sites involved in Ni2+ and methylene blue biosorption process. SEM analysis reveals the appearance of shiny large particles and layers on the biosorbent surface after biosorption that are absent before the biosorption process. In conclusion, it is demonstrated that the Gracilaria seaweed biomass is a promising, biodegradable, ecofriendly, cost-effective and efficient biosorbent for simultaneous bioremoval of Ni2+ and methylene blue from wastewater effluents.


Results and discussion
Several studies have been performed on the adsorption of only one contaminant, despite the fact that both dyes and metals are present together in large quantities of industrial wastewaters 25,26 . Therefore, high priority was given to the simultaneous removal of heavy metal ions and dyes in wastewater treatment processes 27 .
In the current study, the dried biomass of marine red algae Gracilaria was used as biosorbent for simultaneous bioremoval of methylene blue and nickel from aqueous solutions. Gracilaria is a genus of red algae (Rhodophyta). In the current study, Gracilaria verrucosa had been used because it is a cost effective biosorbent, very abundant macroalgae in marine areas and easily collected. Very scarce studies have been done on the biosorption by G. verrucosa 23 . The red algae have a complex cell walls consist of mannan fibrils, xylan, cellulose, as well as a large polysaccharide matrix which contain gelans, agar, carrageen and protein, which are of economic importance along with minerals. Calcium can be absorbed by some red seawater algae and stored in their bodies as calcium carbonate 28,29 . These cell wall components and calcium carbonate can provide various functional groups, increasing the metal ions and dye binding sites on the algal surface.
Statistical optimization of simultaneous bioremoval of methylene blue and nickel by Gracilaria seaweed biomass using FCCCD. Different variables can affect the biosorption processes such as biosorbent dosage, contact time, temperature, pH of the solution and initial concentration of the metal ions 30,31 . Optimizing the conditions for the biosorption process is crucial to achieving maximum biosorption. The conventional single variable optimization method has been used for optimization of the process, but this method is www.nature.com/scientificreports/ not only laborious, expensive and time consuming, but also it does not reflect the impact of interactions between the independent process factors 32 . The optimization of process parameters by statistical experimental designs has been used for many decades. Mathematical and statistical techniques such as response surface methodology (RSM) are very advantageous as they help to determine the optimum level of each variable and help to explore the interactions between different variables for maximum response 33 . Optimization of biosorption process parameters was carried out using FCCCD to study the individual, quadratic and interaction impacts among the process parameters including algal biomass concentration, initial pH level, nickel concentration, contact time and methylene blue concentration on simultaneous bioremoval of Ni 2+ ions and methylene blue from aqueous solution by dry Gracilaria seaweed biomass. FCCCD was also used to predict the best working conditions for the biosorption process for simultaneous bioremoval of methylene blue and Ni 2+ . The FCCCD had eight center points, 32 factorial and 10 axial points which resulted in a total of 50 experimental trials that were used for optimization of the selected variables. In order to determine experimental errors, the eight replicates at the center points were performed. Table 1 shows a design matrix containing coded and actual values of the five process parameters, the experimental and predicted removal percentages of methylene blue and Ni 2+ , and the residuals. Based on the variations in the levels of the five process parameters, the experimental results show substantial differences in the removal percentages of both methylene blue and Ni 2+ . Removal percentage of methylene blue ranged from 52.78 to 94.86% and the removal percentage of Ni 2+ ions ranged from 58.33 to 97.53%.
The highest removal percentages of Ni 2+ and the methylene blue were achieved in the run no. 31 with percent of 97.53% and 94.86% for Ni 2+ ions and the methylene blue; respectively, where algal biomass concentration 6 g/L, initial pH was 8, Ni 2+ concentration was 150 mg/L, the methylene blue concentration was 20 mg/L and contact time was 180 min at room temperature. On the other hand, the minimum removal percentages of Ni 2+ ions and methylene blue were observed in run number 47 where algal biomass concentration 2 g/L, initial pH was 6, Ni 2+ concentration was 100 mg/L, the methylene blue concentration was 10 mg/L for 60 min of contact time at 30 °C. The adsorption efficiency decreased in the run number 47 may be as a result of low level of algal biomass or contact time.
Multiple regression analysis and ANOVA. Both nickel removal (%) and methylene blue removal (%) data were statistically analyzed by using multiple regression analysis and the results were presented in Tables 2, 3, 4 and 5. Linear, interactions and quadratic effects of the selected five process parameters were estimated. The analysis includes the coefficient of determination (R 2 ) which determines the efficiency of the polynomial regression model, the predicted R 2 value, the adjusted R 2 value, the effect of each variable, Fisher test (F test) and probability P value. The Ni 2+ regression model has R 2 value = 0.9902, meaning that 99.02% (Table 2) of variations in the removal percentages of Ni 2+ were due to the independent factors and the model cannot describe just 0.98% of the total variations. In addition, the methylene blue regression model has R 2 value = 0.9944 (Table 3), meaning that 99.44% of the variations in the removal percentages of methylene blue were due to the independent factors and the model cannot describe just 0.56% of the total variations. A regression model with a high R-squared value greater than 0.9 is known to has a great, positive correlation 34 . On the other hand, the values of the adjusted determination coefficient (Adj R 2 ) are 0.9835 and 0.9905 for Ni 2+ and the methylene blue; respectively. The very high value of Adj R 2 verified that the model was very significant (Tables 2, 3). The predicted R 2 values of 0.9746 and 0.9840 for Ni 2+ and the methylene blue removal %; respectively, were in an excellent agreement with the adjusted R 2 values. That revealed a strong agreement between the predicted and observed values of Ni 2+ and methylene blue removal %. The predicted and adjusted values of R 2 must be within 20% of each other, so that we can say that there is a reasonable agreement between them 35 . Therefore, the models used in this study are optimal in the range of experimental variables in the prediction of methylene blue and Ni 2+ removal %.
The negative coefficient values suggest an antagonistic correlation among the variable(s) and the removal percentage while the positive coefficient implies a synergistic relationship between the variable(s) and the removal percentage. Consequently, the negative coefficients values of linear effects, mutual interactions effects between two variables and quadratic effects of the selected process parameters means that they exert a negative effect on methylene blue and Ni 2+ removal % by the dry biomass of Gracilaria seaweed, whereas the positive coefficient values mean that they increase methylene blue and Ni 2+ removal % by the dry biomass of Gracilaria seaweed in the tested range of the selected five process parameters. It can be seen from the coefficient values ( Table 2) that the algal biomass concentration, initial pH level, nickel concentration and methylene blue concentrations and contact time had positive effects on Ni 2+ removal %. However, all the variables (algal biomass concentration, initial pH level, nickel concentration, contact time and methylene blue concentration) had positive effects on the methylene blue removal % (Table 3).
In order to assess the significance of each parameter, P-values were used as a tool. In this case, the variables showing P-values below 0.05 were assumed to have significantly effects 36 . Analysis of variances (ANOVA) for a quadratic regression model of nickel removal (%) indicates that the model is highly important as is apparent from a very small probability value [P-value ˂ 0.0001] with the Fisher's F test (F-value = 146.74) ( Table 2). The P-values showed that linear coefficients for algal biomass concentration, initial pH level, Ni 2+ concentration, methylene blue concentration and contact time were significant for nickel removal with probability values of < 0.0001 for all ( Table 2). The interaction effects between X 1 X 2 (Gracilaria seaweed biomass, initial pH), X 1 X 3 (Gracilaria seaweed biomass, initial Ni 2+ concentrations), X 2 X 3 (initial pH, initial Ni 2+ concentrations), X 2 X 4 (initial pH, methylene blue concentrations), X 2 X 5 (initial pH, contact time), X 3 X 5 (initial Ni 2+ concentrations, contact time) are significant [P-value ˂ 0.05]. While, the interaction effects between X 1 X 4 (Gracilaria seaweed biomass, methylene blue concentrations), X 1 X 5 (Gracilaria seaweed biomass, contact time), X 3 X 4 (initial Ni 2+ concentrations, methylene blue concentrations), X 4 X 5 (methylene blue concentrations, contact time) are not  Table 2). On the other hand, the ANOVA of the quadratic regression model of methylene blue removal (%) indicates that the model is highly important as is apparent from a very small probability value [P-value ˂ 0.0001] with the Fisher's F test (F-value = 256.92) ( Table 3). All the variables (algal biomass concentration, initial pH level, nickel concentration, methylene blue concentration and contact time) were significant for removal of methylene blue with probability values of < 0.0001 ( Table 3). The lower P-values denote that the factors are more significant on the methylene blue removal. The interaction effects between all the variables are significant except interactions between X 3 X 4 (nickel concentration and methylene blue concentration) are not significant [P-value ˃ 0.05]. In addition, the P-values of the coefficients show that the quadratic effects of all the five variables are significant (Table 3).
Statistically analyzed data for nickel removal (%) indicates that the coefficient of variation percentage (C.V. = 1.43%) ( Table 2) is substantially small, suggesting that the experiments carried out have a high degree of precision and reliability 37 . Adequate precision determines the level of noise; the level higher than 4 is preferable and implies the model reliability. The current adequate precision ratio of the Ni 2+ removal model is 53.44, which means that the model was reliable. The predicted residual error sum of squares (PRESS) is an indication of how well every point is suitable for the design. The regression model PRESS value is 107.32. The model's mean and standard deviation are 83.63 and 1.19; respectively (Table 2). Meanwhile, statistically analyzed data for methylene blue removal (%) shows a comparatively low value of the coefficient of variation percentage (C.V. = 1.07%) ( Table 3) confirms a high degree of experimental accuracy and reliability. The present methylene blue removal model had an acceptable, adequate precision ratio of 75.00 and this means the reliability of the model. The model's mean and standard deviation are 82.12 and 0.88; respectively. Value of PRESS for statistically analyzed data of methylene blue removal (%) is 63.64 (Table 3).   Tables 4 and 5 displays the fit summary results applied to the selection of the highest order polynomial model, which has the significant model terms and an insignificant lack of fit test. In addition, the model summary statistics indicate what model has higher adj. and pred. R 2 and lower standard deviation. The fit summary results (Tables 4, 5) demonstrated that, the quadratic models of both nickel and methylene blue removal by the dry biomass of Gracilaria seaweed are very significant with a very low P-value < 0.0001. Lack of Fit Test for Ni 2+ removal is non-significant (P-value = 0.9590) with F-value = 0.39. The model summary statistics for Ni 2+ removal quadratic model (Table 4) recorded the lowest standard deviation of 1.19 and the highest adjusted and predicted R 2 of 0.9835 and 0.9746; respectively. Furthermore, Lack of Fit Test for methylene blue removal is non-significant (P-value = 0.5069) with F-value = 1.06. The model summary statistics for methylene blue removal quadratic model (Table 5) recorded the smallest standard deviation of 0.88 and the highest adjusted and predicted R 2 of 0.9905 and 0.9840; respectively.
The mathematical relationships between the chosen independent factors and the responses (nickel removal and methylene blue removal percentages) are given by the following polynomial equations of the second order: where X 1 -X 5 are the coded values of Gracilaria seaweed concentration, initial pH value, nickel concentration, methylene blue concentration and contact time; respectively.
Above the second-order polynomial Eqs. (1), (2) indicate how the individual variables, their quadratic or their interactions among them affected the removal (%) of nickel or methylene blue from aqueous solution by using the biosorbent dry biomass of the Gracilaria marine algae. The negative coefficient values imply that the linear effects, quadratic effects or interaction effects between the variables have negative effects on the removal (%) of nickel or methylene blue (i.e. decreases the biosorption percentage), whereas positive coefficient values imply that the linear effects, quadratic effects or interaction effects between different factors increase the biosorption percentage of nickel or methylene blue in the tested ranges from the aqueous solutions.
The normal probability plot (NPP) of the residuals. The NPP of the residuals is an essential statistical tool for verifying the suitability of the model 38 . Figure 2A shows the NPP of the residuals for Ni 2+ removal (%) data analysis that indicating that the residuals distributed around the diagonal line of the normal distribution. This means that the model is adequate. Figure 2B  Effects of process variables on nickel and methylene blue biosorption (three dimensional surface plots). The three-dimensional (3D) surface plots were generated in order to identify the optimum conditions for the bioprocess (removal of Ni 2+ and methylene blue from aqueous solution) and to visualize the rela- The predicted value of nickel removal (%) = + 94.77 The predicted value of methylene blue removal (%) = + 92.09 + 2.36X 1 + 2.95X 2 + 2.14X 3 www.nature.com/scientificreports/ tionships between the responses and the interactions between the selected process parameters. 3D graphs for the five variables combined in pairs (algal biomass concentration, initial pH level, nickel concentration, methylene blue concentration and contact time) were created by plotting removal of Ni 2+ (percent) on Z-axis versus two process parameters while other parameters maintained fixed at their midpoints. The three dimensional surface plots (3D) (Figs. 3A, 4A) shows the simultaneous effect of algal biomass concentration, (X 1 ) and initial pH level (X 2 ) on nickel removal (%) and methylene blue removal (%); respectively, while nickel concentration, methylene blue concentration and contact time (X 3 -X 5 ) were kept at their zero levels. The removal percentages of both nickel (Fig. 3A) and methylene blue (Fig. 4A) increases with the increase of initial concentration of algal biomass from 2 to 6 g/L. The increased removal percentages of nickel and methylene blue can be mainly due to the increased surface area and the availability of unsaturated adsorption binding sites on the surface of alga needed for the biosorption process reaction by increasing the algal biomass concentration. The percentage removal of nickel and methylene blue then decreased with increasing the algal biomass concentration from 6.5 to 10 g/L. At high algal biomass concentrations, the agglomeration of the biomass can be a reason for the decrease in the removal efficacy. An increase in the initial pH led to an increase in the percentage of removal up to pH 8 and then further increase in pH decreased the removal percentage. This means that the maximum removal percentages of nickel and methylene blue could be obtained by 6.4 g of the algal biomass and initial pH 8. This may be explained by the slightly alkaline pH 8 improved the access of nickel ions and methylene blue to the adsorption active sites of the algal biomass.
The three dimensional surface plots (3D) (Figs. 3B, 4B) shows the simultaneous effect of algal biomass concentration (X 1 ) and nickel concentration (X 3 ) on nickel removal (%) and methylene blue removal (%); respectively, while initial pH level, methylene blue concentration and contact time (X 2 , X 4 , X 5 ) were kept at their zero levels. The removal percentages of both nickel (Fig. 3B) and methylene blue (Fig. 4B) increases with the increase of initial concentration of algal biomass. The percentage removal of nickel and methylene blue then decreased with increasing the algal biomass concentration from 6.5 to 10 g/L. An increase in the nickel concentration resulted in an increase in the nickel removal (%) up to 171 mg/L (by using 6.4 g/L of the algal biomass) and then further increase in the nickel concentration decreased the removal percentage (Fig. 3B). While, an increase in the nickel www.nature.com/scientificreports/ concentration up to 160 mg/L (by using 6.4 g/L of the algal biomass) resulted in an increase in the percentage removal of methylene blue and then further increase in the nickel concentration decreased the removal percentage (Fig. 4B). The surface plots (Figs. 3C, 4C) shows removal (%) of nickel and methylene blue as a function of algal biomass concentration (X 1 ) and methylene blue concentration (X 4 ) while other parameters maintained fixed at their midpoints. By using high and low algal biomass concentrations, removal (%) of nickel and methylene blue decreased, suggesting that the removal was highly dependent on the concentration of algal biomass. The interaction effects between algal biomass concentration (X 1 ) and methylene blue concentration (X 4 ) are not significant (P-value = 0.2703) for removal (%) of nickel (Table 2). While, the interaction effects between algal biomass concentration (X 1 ) and methylene blue concentration (X 4 ) are significant for removal (%) of methylene blue (P-value = 0.0393) (Table 3).
Similarly, Figs. 3D and 4D shows that the percentages of removal of both nickel and methylene blue increased with an increase in the time of contact. Figures 3E and 4E with the increased concentration of nickel above 164 mg/L, the percentage of removals decreases. This may attributed to the saturation of binding sites on the surface of alga. Figures 3F and 4F shows that the percentages of removal of both nickel and methylene blue have increased with an increase in the initial pH level to pH 8 and then further increase in pH decreased the removal percentage. The effect of both initial pH level and contact time are also reflected in Figs. 3G and 4G. Similarly, the effect of both nickel and methylene blue concentrations are also reflected in Figs. 3H and 4H. The percentage removal of nickel increased with the increased concentration of nickel from 100 to 171 mg/L (Figs. 3H, 4H). In addition, a slight decrease in the removal percentage resulted from a further increase in the concentration of nickel. The high percentage removal of methylene blue was obtained using about 24.2 mg/L of methylene blue (Fig. 4H). In addition, the effects of contact time, nickel concentrations and methylene blue concentrations are also presented in Figs. 3I,J and 4I,J. Initial pH. Different variables can influence biosorption processes such as biosorbent dosage, contact time, initial metal concentration and pH. Among these variables, the initial pH level has been identified to be the most important variable that controls the pollutant biosorption process. The pH levels influences the activity of the functional groups of the biomass, the net charge on the biosorbent cell surface as well as affect the chemistry of metals solutions 39 . The biosorption process of heavy metals is highly pH-dependent, may be minimized or increased by competition between heavy metal ions and H + at the same cellular active sites on the biosorbent cell surface 40 .
Experimental results in this study clearly indicate that the rise in pH enhances the biosorption of nickel by the red algae Gracilaria biomass, and reaches the maximum at pH 8, however; a very high pH increase can lead to nickel precipitation that must be avoided 41 . pH affects the stability and state of nickel species [Ni 2+ , insoluble nickel(II) hydroxides at different pH values] in aqueous environments, as well as the properties of the adsorbent surface 5 . In acid, neutral, and fairly basic solutions, Ni 2+ is the predominant form. When pH is > 8, precipitation of insoluble nickel(II) hydroxides occurs 5 . In this study, the strong pH dependence of nickel biosorption could be attributed to the strong electrostatic attraction at higher pH taking place between the Gracilaria seaweed biomass and nickel ions.
In addition, the influence of pH on the biosorption process of metal ions could be clarified by the association-dissociation of many functional groups that have existed on the algal biomass surface, like hydroxyl and carboxyl groups. Since there is an excess of H 3 O + and H + at low pH, carboxylic groups are acidic and are not dissociated in the solutions (existing in a protonated state). As a result, the repulsive forces between these protonated carboxylic groups and heavy metal ions (positively charged) result in lower biosorption capacity at low pH levels 31,42 . At low pH, high concentrations of H + ions in solution compete with nickel ions for available functional groups on the adsorbent, consequently reducing the effective exchange of metal ions (Ni(II) ions) on the sorbent surface as the sorbent surface takes up more H +5 . However, the concentration of the proton in solution is lower than that of the metal ions with a further rise in the pH level, and the surface of the sorbent accumulates more negative charges that attract more metal ions. This facilitates the electrostatic attraction of positively charged Ni(II) ions toward fairly negatively charged adsorbent surface 43 . At higher pH levels (pH 8) deprotonation of functional groups with negative charges such as hydroxyl, amine and carboxyl on the surface of the sorbent increases the electrostatic attraction of positively charged heavy metal ions 31,44 and improves the metal biosorption 5 .
At higher pH levels (more than pH 8), the decrease in the biosorption efficiency may be correlated with the repulsion between anionic solution negative charges and the negative charges of the biosorbents' surface 45 . Precipitation of insoluble anionic hydroxide complexes (Ni(OH) 2 ) occurs at the high increase in pH (pH > 8.3) that reduces the free metal ions concentration and decrease the biosorption capacity for metals 31,42,44 . Vakili et al. 5 reported that acidic pH is not favorable for the adsorption of nickel ions and maximum sorption capacities were obtained at pH ranging from 6 to 8.
Biosorption of methylene blue was done by using the red algae Gracilaria as an biosorbent. Experimental results in this study clearly indicate that the rise in pH enhances the biosorption of methylene blue dye by the red algae Gracilaria biomass, and reaches the maximum biosorption at around pH 8. The result can be explained by the protonation and de-protonating of the methylene blue dye as well as the red algae Gracilaria surface as biosorbent. The ionic attraction regulates the adsorption of the dye on the surface of biomass 46 . At acidic pH values, the surface of the biosorbent carries a positive charge due to the presence of H + in the solution. This leads to electrostatic repulsion between the cationic methylene blue dye and positively charged Gracilaria surface as biosorbent. On the contrary, at alkaline pH values, due to the existing of OH − ions in excess, the surface of the biosorbent carries a negative charge. This results in electrostatic attraction between cationic methylene blue dye and negatively charged Gracilaria surface leading to enhanced adsorption at higher pH values. Hammud et al. 47 reported that the maximum MB removal by Carolina sp. was attained at pH 6.8. While, Khaled et al. 48 49 also found maximum removal of MB at pH 7 using Sargassum muticum. However, Vilar et al. 50 recorded that MB removal using Gelidium biomass based material was unaffected in the pH range from 4 to10 but was impaired at low pH. The biosorption of the colored dye cations on the surface of the adsorbent has been mainly affected by the charges of the algal cell walls, which in turn were affected by the solution pH 51 . The algal cell walls, mainly made up of carbohydrates and protein, have functional groups such as sulfate, carboxyl, hydroxyl, amino, etc. that can serve as metal binding sites 52 . The maximum adsorption percentage of methylene blue removal value was 92% from solutions using marine algae Ulva lactuca was found at around pH 8 53 . A fewer anionic adsorption sites were generated on the surface of the biosorbent (the dried algal biomass) at low pH values, and the sorption process decreased probably due to electrostatic repulsion between excess H + in the solution and the cationic dye molecules for sorption sites on the dried algal biomass 51 .
Contact time. The time of contact influences the efficiency of the biosorption process. As the time of contact increases, the removal percentage increases as well. In the present study, experimental results clearly indicate that the percentage of simultaneous removal of Ni(II) ions and methylene blue by the red algae Gracilaria as biosorbent increases as the contact time increase up to the optimum, and reaches the maximum at 177.88 min. The increase in the rate of removal might be because of the availability of vacant functional groups on the biosorbent surface. Further increase in contact time leads to the occupation of active sites which results in the saturation of the algal cell surface causing an equilibrium state and no additional adsorption 54 .
In the case of different forms of biosorbents, the optimal contact time for maximum removal percentage is varied. Ibrahim 42 reported that the optimum contact time for red macroalgae was 60 min. While, Ahmad et al. 31 reported that the optimum contact time for the immobilized algal mass was 300 min (240 min for free suspended mass). Akhtar et al. 4 found that adsorption of nickel onto the biomass of Chlorella sorokiniana immobilized in loofa sponge was rapid and reached equilibrium after only 15 min. On the contrary, the biosorption capacity increased with increasing contact time up to 110 min, after which the equilibrium was attained and the biosorption capacity is less constant 53 . Pahlavanzadeh et al. 39 stated that the biosorption capacity of nickel by dried algal biomass increased with increased time of contact and the equilibrium was attained within 120 min. Xu et al. 55 reported that adsorption attained equilibrium after 168 h (7 days). Obviously, such variability of the equilibrium time depends on the used adsorbent. This difference could be related to their structural and physicochemical characteristics diversity (e.g. surface area, pore volume, pore size, permeability, and density and availability of adsorption sites) 56 . Materials with higher adsorption capacity needs a longer time to reach equilibrium, implying that the advantage of a larger capacity is counter balanced by a longer treatment time, and vice versa 5 .

Initial Ni(II) and methylene blue concentrations. The initial concentrations of metal ions and methyl-
ene blue also affect the biosorption process. As a general, increasing the initial concentrations of metal and dye usually leads to an increase in the percentage of removal as it offers a driving force that causes metal ions and dye to pass onto the surface of the biosorbent particles. An increase in the initial concentration of nickel ions leads to an increase in the sorption capacities of nickel ions on the surface of adsorbents 5 .
Algal biomass concentration. The availability and low cost are the key reasons for selecting biomass for largescale industrial uses 57 . Algae are known as among the most promising biosorbents types. Algal biomass has a wide range of advantages of removing heavy metals from wastewater, such as cost-effective renewable natural biomass, high metal removal efficiency, high absorption capacity, biomass regeneration and re-uses 58 . Large-scale biosorption using dead biomass can be easier and more applicable than living microorganisms, since the latter require a supply of nutrients 59 . The effect of metabolic processes of living microorganisms on biosorption is often not appreciated.
The increase in simultaneous biosorption capacity of methylene blue and nickel with an increase of the biosorbent concentration (red algae Gracilaria) is mostly attributed to the accessibility of many more unoccupied active sites on the surface of biosorbent 60 . The biosorption capacity of the surface of the biosorbent is specified by the availability of the active sites, for example, phosphate, imidazole, sulfate, sulphuryl, phosphoryl, amine, hydroxyl, carboxyl, etc. 61 . However, decreased biosorption capacity at higher concentrations of Gracilaria biomass could be attributed to the agglomeration of Gracilaria biomass, which in turn could reduce the intercellular spacing leading to a reduction of the total effective surface area available for the biosorption process and consequently decrease the number of active binding sites on the algal biomass surface available to simultaneous removal of dye and metal ions 51,60,62 .
However, Kavitha et al. 63 stated that the decrease in biosorption process efficiency with the higher biosorbent dose can be attributed to blocking the binding sites from metal ions. Hannachi et al. 23 reported the biosorption potential of G. verrucosa for Zn 2+ ions removal. They found that the biosorption percentage of Zn (II) increases with an increase in the initial concentration of G. verrucosa biomass from 1 to 4 g/L and there is no substantial increase in the biosorption percentage of Zn(II) when the concentration of G. verrucosa biomass increases beyond 4 g/L. This has been explained by the increased number of unsaturated active adsorption sites on the biosorbent surface at a high concentration of algal biomass, and the metal ions available in the solution are insufficient to bind to all available binding sites, resulting in a low percentage of metal biosorption 23 .
On the other hand, Vijayaraghavan et al. 64 observed that the biosorption percentage of methylene blue (MB) from aqueous solution increased with increase in Gracilaria corticata biomass concentration up to 5 g/L and further increase in the concentration of Gracilaria corticata biomass had no impact on the removal of percentages. The biosorption percentage of the dye decreases by increasing the concentration of the algal biomass and this has been explained by the lack of sufficient dye molecules at high sorbent dosages to cover all the active adsorption sites on the algal biomass, resulting in low dye biosorption percentage 64 . On the contrary, Arumugam et al. 65 stated that the metal ion removal is proportional to the increasing seaweed dosage from 2 to 10 g. They Scientific Reports | (2020) 10:17439 | https://doi.org/10.1038/s41598-020-74389-y www.nature.com/scientificreports/ suggested that a higher concentration of the algal biomass provides higher surface area and more adsorption active sites to bind to a constant concentration of metal ions. Therefore, more metal ions bind to the surface of seaweed and increase the removal percentage. Numerous studies have been conducted on the adsorption of single-pollutant, despite the fact that large quantities of both dyes and metals are present together in the industrial wastewaters. High priority was therefore given to the simultaneous removal of dyes and heavy metal ions during wastewater treatment processes. Biosorption is a method used to treat wastewater efficiently by removing toxic metal ions and dyes from contaminated effluents using cost-effective and efficient natural biosorbents. Marine algae (seaweeds) are sustainable natural biomass, have been shown to be highly effective biosorbents with a high binding ability to a wide range of metals and dyes contaminants in aqueous effluents. Dried biomass of marine algae have proven to be more favorable than living cells as biosorbents because it does not require nutrient supply, it can be preserved and used prolonged periods at room temperature and can be regenerated and reused for multiple treatment cycles. Biosorption process has been affected by many environmental factors. Optimization of the operating process variables is crucial to achieving maximum biosorption. Response surface methodology is a mathematical and statistical technique used over several decades to optimize different process variables. It is very advantageous than conventional single variable optimization method as it helps to determine the optimal levels of the variables and helps to explore the interaction effects between various process factors in order to achieve maximum response.
The model adequacy. In the application of regression models, the main goal is to assess the adequacy of the model to predict the response variable. Figure 5A shows the actual versus predicted percentages for a simultaneous removal percentage of methylene blue from aqueous solution by Gracilaria seaweed biomass. Figure Figure 5B shows Box-Cox plot of model transformation has generated for a simultaneous removal percentage of methylene blue from aqueous solution by Gracilaria seaweed biomass. Box-Cox plot of model transformation method can help checks data that are not normally distributed by transforming the data for normalization. As can be seen in Fig. 5B, the optimal value of Lambda (λ) of 1 is lies between the two vertical red lines so that data transformation is not required.

Desirability function (DF).
The main goal of the experimental design is to determine the optimal predicted conditions for maximizing the responses. The desirability function (DF) was used to find the optimal predicted conditions for maximum response 66 . The values for DF varied from zero (undesirable) to one (desirable). The numerical optimization finds the points that maximize the desirability function. The DF option in the Design Expert Software (version 7.0.0) was used for the optimization process. Figure 6 shows the optimization plot displays the desirability function and the optimum predicted values for the maximum simultaneous removal percentage of methylene blue and nickel ions. The optimal predicted conditions attained using the desirability function for the maximum simultaneous removal percentages of the nickel ions and methylene blue were the algal biomass concentration of 6.94, initial pH level of 7.76, initial nickel ions concentration of 154.67 mg/L, initial concentration of methylene blue of 24.29 mg/L and contact time of 177.88 min. The conditions could be resulted in the removal percentages of 98.54% and 96.25% (with DF of 1) for nickel ions and methylene blue; respectively. In order to verify the removal percentages of the nickel ions and methylene blue under the optimal predicted conditions, the experiments have been done in triplicate and compared with the predicted values. The experimental removal percentages averages of nickel ions and methylene blue were 98.9% and 97.31; respectively. The verification showed that the observed and predicted values of a strong agreement imply that the DF effectively determines the optimal predicted conditions for the simultaneous removal of nickel ions and methylene blue.
Surface morphology. The surface morphology of Gracilaria seaweed biomass was analyzed by scanning electron microscopy before and after the biosorption process. The micrograph before the biosorption process revealed a regular surface of the biosorbent (Fig. 7A). After biosorption of methylene blue and nickel ions, the surface of the biosorbent was filled with adsorbate ions. Figure 7B clearly shows the presence of new glossy massive particles and layers on the surface of the biosorbent which are absent on the surface of the biosorbent before the biosorption process. Obvious morphological changes were seen in the cell surfaces (Fig. 7B), such as shrinking of surface. These changes may be due to strong crosslinking between nickel ions and negatively charged functional groups on the surface of Gracilaria seaweed biomass.  Table 6) to detect any variation in the morphological features and surface characteristics due to the interaction of methylene blue and nickel ions with the functional groups. Biosorption of metal ions occurs via an ion exchange process on the cell surface. FTIR spectroscopy was often used to detect changes in the frequency of vibration in seaweeds 67 . The extent of band shift reflects the degree to which functional groups interact with the adsorbed ions 67 . The spectra of Gracilaria seaweed biomass were measured within wave number range from 400 to 4000 cm −1 . FTIR spectrum of dry Gracilaria seaweed biomass sample before simultaneous bioremoval of methylene blue and nickel ions showed the characteristic absorption peaks at 3748, 3431, 2929, 2269, 2199, 2144, 1648, 1544, 1441, 1277, 1053, 759, 705 and 616 cm −1 . FTIR spectrum of dry Gracilaria seaweed biomass sample after simultaneous bioremoval of methylene blue and nickel ions showed the previously mentioned characteristic signals were shifted to 3767, 3446, 2932, 2272, 1646, 1436, 1256, 1053 and 604 cm −1 . Significant shifts in the absorption peaks after the biosorption confirm the significant role of the functional groups on Gracilaria biomass surface in the biosorption of Ni 2+ and methylene blue dye from aqueous solution.
In the FTIR spectrum of the dry Gracilaria seaweed biomass sample before biosorption of methylene blue and nickel ions, the characteristic absorption peak at 3748 cm −1 is attributed to vibration of free Si-OH groups 68 shifted after the biosorption process to peak at 3767 cm −1 which assigned to -OH, N-H stretching. The peak at 3431 cm −1 (before biosorption) represents -OH stretching vibrations 69 shifted to 3446 cm −1 (after biosorption) which assigned to -OH or N-H stretching vibrations 70 . The peak at 2929 cm −1 is assigned to C-H stretching 71 shifted after the biosorption process to peak at 2932 cm −1 which assigned to C-H stretching 72 . Whereas, Júnior et al. 73 recorded that the signal at 2269 cm -1 corresponding to the N=C=O (isocyanate) groups. It is found that the signal at 2199 cm −1 was assigned to the nitrile group (C≡N) 74 . However, peak at 2144 cm −1 , assigned to the vibrations of hydrocarbons, alkynes (C≡C) non terminal groups 75 . Significantly, the signal at 2144 cm -1 can be assigned to the C-H stretching in alkanes or triple bond stretch of alkyne C≡C 76 . The three peaks at 2269, 2199 and 2144 cm −1 are transformed into a single peak at 2272 cm −1 (after biosorption process) which is belonging to -O-C≡N groups (free isocyanate groups) 77 . The major peak shift suggests that the nitrile, alkanes, alkynes and isocyanate groups participated in the process of biosorption. The peak at 1648 cm −1 in the spectral range (before   80 . In addition to these peaks, one peak at 1544 cm −1 in the spectral range may be attributed to stretching vibration of carboxyl (-CO), C=C stretching or amide-II (C-N stretch, protein N-H bend) groups 23,81,82 , this peak disappeared in the FTIR spectrum of the Gracilaria biomass sample after biosorption process. The peaks at 1530-1560 cm -1 are attributed to stretching vibration of amino groups (NH stretching) 67 .
In parallel, the CH 2 -bending vibrations of the methyl group in lipids and proteins peak at 1441 cm −183 is shifted to lower wave number at 1436 cm −1 of the methyl peak (O-CH 3 methyl esters) 57 , that indicated the involvement of the methyl group in the biosorption process via ion exchange. The peak at 1277 cm −1 was attributed to stretching of C-O stretching vibrations shifted to peak at 1256 cm −1 (-21) which attributed to stretching vibration of the PO 2− asymmetric vibrations (phosphate I) 79,84 . The signal at 1053 cm −1 was attributed to the PO 4 3− group vibration 85 and not shifted after the biosorption process. The presence of intense peaks in the region of 1000-1100 cm −1 is a characteristic of polysaccharides 80 . Additionally, two peaks have appeared in the FTIR spectrum of the Gracilaria biomass sample before the biosorption process of methylene blue and nickel ions at 759 and 705 cm −1 which attributed to the stretching vibration of C-S 86 . These peaks disappeared in the FTIR spectrum of the Gracilaria biomass sample after biosorption process. Gracilaria species are known to produce   80 . Whereas, the signal at 616 cm −1 (before biosorption) shifted to peak at 604 cm −1 (after biosorption) which represents the phosphate groups 87,88 . These shifts in the absorption peaks confirm the association of these functional groups on Gracilaria biomass in the biosorption process of methylene blue and nickel ions. Finally, the existing FITR spectra indicated that the key groups participating in the biosorption of methylene blue and nickel ions were methyl, alkynes, amide, phenolic, carbonyl, nitrile and phosphate groups.

Materials and methods
Preparation of the biosorbent. The red alga, G. verrucosa, used in this study, which classified under Rhodophyta, was harvested from the Mediterranean Sea coast of Abu-Qir, Alexandria, Egypt during summer season, 2019. Algal biomass was thoroughly washed with running tap water to remove salts, sand and any other external particles that may be attached to the surface. Algal biomass was finally washed with distilled water, immediately after the running tap water. The washed algal biomass was then air dried to a constant weight was obtained. The dried algal biomass was crushed and the particles with an average size of 1-1.5 mm were selected and stored until used for biosorption experiments.
Preparation of nickel and methylene blue dye solutions. Nickel standard solution was prepared by the precise dissolving of weighted nickel (II) sulfate (NiSO 4 ) in distilled water at a concentration of 1000 mg/L. The desired work concentrations were then prepared by dilution of nickel stock solution.
The methylene blue stock solution was prepared by the precise dissolving of weighed dye in distilled water at a concentration of 1000 mg/L. The desired work concentrations were then prepared by dilution of the methylene blue stock solution.
The biosorption experiments. The biosorption experiments were conducted in in batch mode in 250 mL flasks with 100 mL working volume to evaluate the efficiency of Gracilaria seaweed biomass for simultaneous bioremoval of methylene blue and nickel ions from aqueous solution (schematically illustrated in Fig. 9). Dry Gracilaria seaweed biomass was thoroughly mixed with the previously prepared solutions of methylene blue and nickel with different concentrations as demonstrated in the 50 trials (Table 1). The initial pH has been adjusted with the addition of 0.1 N NaOH or 0.1 N HCl to each solution. The suspensions were incubated and agitated at 150 rpm in a shaker incubator at 30 °C for period contact time.

Optimization of biosorption experiments for simultaneous bioremoval of methylene blue and nickel by Face-centered central composite design (FCCCD).
Five-variables FCCCD with eight center runs was generated in this study with Design Expert Software (version 7.0.0). FCCCD was used to determine the optimum level of each variable to achieve the maximum efficiency for simultaneous bioremoval of methylene blue and nickel ions from aqueous solution. In addition, help to explore the linear, quadratic and interactions effects between the selected process parameters that have a significant impact on the bioremoval of methylene blue and nickel ions. The five variables used were as follows: initial Ni 2+ concentrations, initial pH, Gracilaria seaweed biomass (biosorbent), methylene blue concentrations and contact time. As shown in Table 1, each variable was evaluated at three coded levels (− 1, 0 and + 1). All variables, actual and coded levels and the Table 6. FTIR spectral analysis of Gracilaria seaweed biomass before and after biosorption of methylene blue and nickel ions. for methylene blue removal (%) and Ni 2+ removal (%) in each trial were the average of the triplicate. All experiments were performed at 30 °C to save energy. The correlations between the selected independent variables and the responses (methylene blue and nickel biosorption percentages) were determined using the equation of second-degree polynomial as follows: In which Y is the predicted methylene blue or nickel biosorption, X i is the coded levels of the process parameters, β i (linear coefficient), β 0 (regression coefficients), β ij (interaction coefficients) and β ii (quadratic coefficients). Statistical analysis. "Design Expert version 7 for Windows and STATISTICA softwares were used for the experimental designs, statistical analysis and to plot the three-dimensional surface plots". Analytical methods. The content of each flask for FCCCD experiment was centrifuged at 6000 × g and analyzed using Atomic Absorption spectroscopy according to "standard methods for the examination of water and wastewater 23rd edition 2017" 89 . The capacity of Gracilaria seaweed biomass as biosorbent for Ni 2+ removal was determined as follows: where, C i , C f are the initial and final concentrations of nickel ions (mg/L); respectively.
The residual methylene blue concentrations were determined spectrophotometrically with a UV/Vis spectrophotometer at the wavelength of the highest absorbance (λ max ) that was 670. The methylene blue removal % was determined as follows: where, C i and C f are the initial and residual concentrations of methylene blue (mg/L); respectively.
Fourier-transform infrared (FTIR) spectroscopy. FTIR analyses were performed in order to elucidate the distinctive surface functional groups of Gracilaria seaweed biomass that may be responsible for binding of methylene blue and nickel ions. The samples of Gracilaria seaweed dry biomass were analyzed before and after (5) Removal efficiency (%) = C i − C f C i × 100 Figure 9. Schematic illustration of batch experimental biosorption steps to assess the efficiency of Gracilaria seaweed biomass for simultaneous bioremoval of methylene blue and nickel ions from aqueous solution.