The response surface methodology using Box-Behnken design was applied to optimize the fibrin scaffold. Four variables including Ca
2+ concentrations, cell numbers, various ratio of plasma/RPMI 1640 and thickness of fibrin scaffold were considered to evaluate the conditions for fibrin scaffold optimization. Urea secretion in culture media was opted as the response for different run series of the runs.
Table 1 shows the design matrix for four independent variables for urea secretion in 27 experiments, which were run in three times. Then mean of each response (urea secretion) was computed. The obtained data from the experiments were analyzed by linear multiple regression using Minitab 15 and presented in
Tables 2 and
3. The corresponding second-order response model founded after the regression analysis was (Day 5 of incubation):
Where Y is the predicted response (urea secretion) and X
1, X
2, X
3 and X
4 are the coded values of Ca
2+ concentration, cell number, plasma/RPMI1640 ratio and thickness, respectively. As illustrated in
Tables 4 and
5, relatively high F value and very low P value indicated that the experimental model was in excellent conformity with the experiment for both responses (day 5 and day 10). ANOVA showed that linear term of the polynomial model was much significant. Data in
Tables 4 and
5 show low value of F test and high P value (P > 0.05) for lack of fit, which was non-significant indicating that the model is fit. The linear regression coefficients R
2 = 0.9795 and 0.9351 and the adjusted determination coefficients R
2 (Adj) were 0.9556 and 0.8594 for the model (for day 5 and day 10, respectively), which showed accuracy of the model for the Box-Behnken design. R
2 values close to one proposed the robust model resulting in actual values of responses which were closer to the predicted one (
16). Among interactive terms for day 5 of incubation, only X
2X
3 did not have a significant influence (P > 0.05) on responses (urea secretion), while X
1X
3 (Ca
2+ concentration vs. plasma/RPMI1640 ratio) and X
2X
4 (cell number vs. thickness) coefficient presented a negative effect on HepG2 cell viability (urea secretion for day 5); therefore, they were not suitable for viability of HepG2 cell line. These negative impacts may attribute to the fact that high Ca
2+ provided turbid gels that were more viscous than transparent gels. This higher viscosity may hinder the formation of more densely cross-linked networks leading to failure of withstanding high mechanical loading needed for the cell viability (
28). Among the four variables, the plasma/RPMI 1640 ratio and thickness had a negative impact on cell viability. However, P > 0.05 for the thickness means that no significant correlation is present. The detrimental impact of plasma/RPMI 1640 may attribute to the high fibrin concentration, which was proportional to the plasma volume, resulted in a more rigid and dense fibrin gel that halted cell migration and proliferation (
11).
Due to the high number of variables, one by one comparison is time consuming and may lead to misinterpretation. Therefore, the final optimum levels of four components were calculated by means of Minitab Response Surface Optimizer function. Optimized values of the factors (for day 5) were found to be: Ca
2+ concentration 0.15 mol, cell number 105, plasma/RPMI 1640 -1:4 and thickness 2.3 mm. A similar pattern was also found for day 10 results. In addition, this model forecasted up to 21.19 mg/dL of urea secretion proportional to the cell viability and activity. By comparing the forecasted and observed values of the Box-Behnken design (
Table 1), their good correlation supported precision of the response model and existence of an optimal point (
16,
29). However, most researchers investigated optimal conditions for manufacturing fibrin scaffold by means of one-factor-at-a-time methods; there was no report of using statistical design. For the first time, we used the Box-Behnken Design to determine the optimum components for fabricating fibrin scaffold. In this regard, Eyrich et al. reported the optimum fibrin scaffold components as fibrinogen concentration of 25 mg/mL, 3 million cells per construct and Ca
2+ concentration of 20 mM in pH between 6.8 and 9 using one-factor-at-a-time method (
23). In a similar study, Willerth et al. found optimal fibrin scaffold constituents of 10 mg/mL for fibrinogen and 2 NIH units/mL of thrombin (
9). In another study, Ferreira et al. evaluated various scaffolds for blood-hematopoietic stem cell expansion and found that fibrin-based scaffold provided the best 3D support, which gave the highest numbers of engraftment and multilineage differentiation in comparison to other 3D biomaterial scaffolds. These features may be due to the efficient cell adhesion to the substrate, which is known to be part of the natural process happening in the liver niche that controls cell proliferation and differentiation. Moreover, this efficient adhesion is sufficient enough for cell migration and homing abilities (
30). In comparison to our results, it was concluded that moderate concentration of fibrinogen plus Ca
2+ concentration would result optimum fibrin scaffold, though at this fibrinogen concentration, the proliferation and migration of HepG2 cell line would be occurred optimally.
Additional studies are needed to elucidate other components, which exist in human plasma for better fibrin scaffold manufacturing. In addition, it is necessary to evaluate the immunological reactions to ensure lack or minimum allergenicity. Using human plasma instead of purified fibrinogen provides an easy method for fibrin scaffold fabrication and is considered an economic advantage due to elimination of the costs of purification steps in the preparation of pure fibrinogen with no need for using thrombin. Moreover, in human plasma, some growth factors or other components may exist in human plasma, which are essential for cell proliferation or differentiation.