Abbreviations: ANOVA: Analyzes of Variance; DPPH: 2, 2-Diphenyl-1-picrylhydrazyl; TPTZ: 2, 4, 6-Tris(2-pyridyl)-s-triazine; DRSC: DPPH Radical Scavenging Capacity; FRAP: Ferric Reducing Antioxidant Power; RSM: Response Surface Methodology; TPC: Total Phenolic Content; UAE: Ultrasound-assisted extraction; UV: Ultraviolet Radiation; HCl: Hydrochloric Acid
1.
Introduction
Seaweeds (Marine macroalgae) are a group of marine multicellular algae, rich in minerals, vitamins and polysaccharides. They are known as potential sources of bioactive substances with strong biological activities such as antibacterial, anticancer, antioxidant, anti-fungal and antiviral properties [1]. They are divided into three main groups based on their colors, including Rhodophyceae (red seaweed), Chlorophyceae (green seaweed) and Phaeophyceae (brown seaweed) [2]. They are abundant on the coastline of Vietnam and are used for food, medicinal and cosmeceutical purposes [3]. Additionally, recent studies on seaweeds have proven their usefulness as biofuels, fertilisers, fish feed, food ingredients, cosmeceuticals and their application in bioremediation and anti-biofilm activity [4,5,6,7,8,9,10,11]. Their pharmacological and biomedical application including antibiotics, antiviral, antifungal, anticancer, anticoagulant, anti-inflammatory and antioxidant properties have also been reported [12,13,14,15,16,17].
Red seaweeds considered as the largest group of marine macroalgae have been found to contain a variety of bioactive compounds including polysaccharides (alginate, agar, and carrageenan), lipids, phenolics, steroids, glycosides, saponins, alkaloids and triterpenoids [18]. Recent studies have reported that isolated compounds from red seaweeds possessed a range of biological activities such as antioxidant, antimicrobial, antidiabetic, anti-inflammatory and anticancer properties [12,19,20,21,22]. These findings offer red seaweeds as a promising source of bioactive compounds which can be used in the cosmetic, pharmaceutical and food industries.
Ultrasound-assisted extraction (UAE) is one of recent novel extraction techniques. The method applies the cavitation phenomenon which occurs when an extraction solvent in contact with a sample, that is subjected to ultrasounds at high frequency pulses to generate local hotspots at macroscopic scale with high shear stress and temperature by producing cavitation bubbles [23,24]. The collision of cavitation bubbles created in the bath, results in pressure and temperature changes thereby facilitating the rate of mass transfer of analytes to the solvent. Although the method can successfully be applied at low temperature, short time and less solvent, high extract yield can be obtained.
Response surface methodology (RSM) is an effective tool for optimisation of the extraction process. It is a compilation of mathematical and statistical methods, helpful for developing the models as well as analysing the effects of parameters and their interactions [25,26]. RSM has been successfully used to optimise the extraction conditions of antioxidants and phenolic compounds from different seaweed species [27,28,29].
This study aimed to investigate total phenolic contents and antioxidant activities of different red seaweeds (Galaxaura arborea, Mastophora rocea and Laurencia intermedia Yamada). L. intermedia showed the most potential one and was selected to investigate the effects of ultrasound-assisted extraction conditions including temperature, time, solvent concentration and ratio of sample to solvent on the phenolics and antioxidant power. In addition, the optimum UAE conditions would be established for maximising the yields of TPC and antioxidants produced from L. intermedia. The crude extract and derived fractions (n-hexane, ethyl acetate and aqueous fractions) from L. intermedia were then prepared for examination of phenolic content and antioxidant activities.
2.
Materials and methods
2.1. Sample collection and preparation
The red seaweed samples (Galaxaura arborea, Mastophora rocea and Laurencia intermedia Yamada) were harvested during March to July 2017 at the coast of Nha Trang Bay, Khanh Hoa Province, Vietnam. The samples were authenticated by Mr. Do Anh Duy (Research Institute for Marine Fisheries, Hai Phong City, Vietnam). The voucher specimens (NT-02) of L. intermedia Yamada was deposited at the Department of Pharmacognosy, Faculty of Pharmacy, University of Medicine and Pharmacy in Ho Chi Minh City, Vietnam. Samples were individually rinsed to remove particulates and epiphytes, and then air-dried in the shade. Dried seaweed samples were milled into a fine powder to fifty (50) micron. Dried powder was vacuum packaged and stored at −20 ℃ for phytochemical and biological analysis.
2.2. Chemicals
Folin-Ciocalteu's phenol reagent, 2, 4, 6-Tris(2-pyridyl)-s-triazine (TPTZ) were obtained from Sigma- Aldrich Company (USA); 2, 2-diphenyl-1-picrylhydrazyl (DPPH) was purchased from Alfa Aesar (UK). The rest of the chemicals and reagents were of analytical grade.
2.3. Screening test to select the starting sample
A conventional extraction method was used in screening TPC and antioxidant activity of three seaweed species (G. arborea, M. rocea and L. intermedia Yamada). Briefly, 1 g of the dried seaweed sample was mixed with 50 mL of methanol 75%. The mixture were placed in the Orbital Shaking Water Bath (VS—1205 SW2) at 60 ℃ for 2 hours. Post-conventional extraction, the mixture was cool in a cold-water bath and centrifuged at 8500 rpm for 20 min at 4 ℃ (MEGA 17R Small High Speed refrigerated centrifuge) to obtain the extract. The total phenolic content and the antioxidant capacity of the extract were determined and the species with the highest phenolic content and antioxidant activity was selected for subsequent experiments.
2.4. Experimental design
Based on single factor experiments (data not shown), ethanol was found to be the suitable solvent for extraction of TPC and antioxidants from L. intermedia, and four key variables with optimal ranges were selected for RSM optimisation, including ethanol concentration (0–75%), temperature (40–60 ℃), time (20–60 min) and sample to solvent ratio (1–6 g/100 mL). The Box-Behnken design was employed to investigate the interactions of these variables on the extraction efficiency of TPC and antioxidant capacity. The design consisted of 27 experimental runs with three center points as shown in Table 3. The data were fitted into a second-order polynomial model as follows (Equation 1):
where k is the number of variables; Xi represents independent variables affecting the responses Y (TPC, DPPH and FRAP); and βo, βi, βii and βij are the regression coefficients for intercept, linear, quadratic and interaction terms, respectively.
2.5. Ultrasound-assisted extraction (UAE)
Ultrasound-assisted extraction was performed in an ultrasound bath (Branson 2510, 42 KHZ, 100 W, Branson Ultrasonic Corp., Danbury, USA). Seaweed samples were mixed with ethanol with appropriate concentrations, placed into an extraction flask and sonicated for different times and temperatures. After the extraction, the flasks were immediately removed from the ultrasound bath and cooled to room temperature by cooling water. The seaweed extracts were centrifuged using a MEGA 17R small high-speed refrigerated centrifuge at 8500 rpm, 4 ℃ for 20 min. Extracts were collected and used for the determination of the TPC and antioxidant capacity.
2.6. Fractionation
The crude extract was obtained using the optimum extraction conditions for L. intermedia determined as above. The extract was condensed using a vacuum rotary evaporator (Laborota 4001, Heidolph, Germany). The crude extract was then fractionated using n-hexane and ethyl acetate to generate n-hexane, ethyl acetate and aqueous fractions, and their phenolic content and antioxidant activity were further evaluated.
2.7. Total phenolic content (TPC) determination
The TPC of the extracts was determined as described by Pham et al. [30]. To 0.5 mL of extract, 2.5 mL of 10% (v/v) Folin-Ciocalteu reagent was added and left at room temperature for 8 min, then thoroughly mixed with 2 mL of 7.5% (w/v) Na2CO3. The resulting mixture was incubated in the dark at room temperature for 1 hour. The absorbance was then measured at 765 nm using a UV–VIS spectrophotometer (Biochrom Libra S50, Biochemical Ltd., Cambridge, UK). Gallic acid was used as a standard with the results expressed as mg of gallic acid equivalents per 100 g of dried material (mg GAE/100 g).
2.8. DPPH radical scavenging capacity (DRSC) assay
DRSC was assessed as described by Pham et al. [30]. DPPH stock solution was prepared by dissolving 24 mg of DPPH in 100 mL of methanol. The stock solution was kept in the dark at −20 ℃ for further use. To 10 mL of the stock solution, 45 mL methanol was added to prepare the working solution at absorbance value 1.1 ± 0.02 at 515 nm. To 2.85 mL of the working solution, 0.15 mL of the sample was added and left in the dark at room temperature for 3 hours. The absorbance was measured at 515 nm using a UV–VIS spectrophotometer (Biochrom Libra S50, Biochemical Ltd., Cambridge, UK). Using trolox as a standard, and the results were expressed as mg trolox equivalents per 100 g of dried weight (mg TE/100 g).
2.9. Ferric reducing antioxidant power (FRAP) assay
FRAP was measured according to the method of Pham et al. [31] with some modifications. FRAP working solution was prepared by mixing 300 mM acetate buffer, 10 mM TPTZ in 40 mM HCl, and 20 mM FeCl3 in the ratio of 10: 1: 1. To 0.15 mL of sample, 2.85 mL of the working FRAP solution was added and the absorbance was read at 593 nm after incubation in the dark for 30 min at room temperature. Results were expressed as mg trolox equivalents per 100 g of dried material (mg TE/100 g).
2.10. Statistical analysis
All experiments were conducted in triplicate. JMP 15 was used to design optimisation experiments, to generate the model equations, to graph the 3D and 2D contour plots of the responses and to predict the optimum values for the independent variables. Minitab 16 and excel statistical tools were used for data analysis.
3.
Results
3.1. Screening antioxidant capacity and phenolic content of different red seaweeds
The data in Table 1 compare the antioxidant capacity and total phenolic content of three screened red seaweed's species (G. arborea, M. rocea and L. intermedia Yamada). According to the results obtained, L. intermedia Yamada demonstrated the highest total phenolic content and antioxidant activity. It also shows no significant difference of the DRSC and TPC of G. arborea and M. rocea. L. intermedia Yamada was therefore selected for subsequent investigation.
3.2. Fitting the models
Table 2 shows the Box-Behnken design (BBD) with the experimental and estimated data. The actual values strongly correlated with the estimated values, that was revealed through high R2 values (0.95, 0.97 and 0.99 for TPC, DRSC and FRAP, respectively) (Table 3). Additionally, the p-values for lack of fit of TPC, DRSC and FRAP were 0.0909, 0.0768 and 0.1951, respectively, greater than 0.05, revealing the lack of fit was not significant. In Table 3, the estimated regression coefficients for the quadratic polynomial model and analyses of variance for the experimental results of TPC, DRSC and FRAP, and p-values of models are also displayed. The models developed to predict the TPC and antioxidant activities were found to be highly reliable for predicting the responses. These models can be fitted into the following second-order polynomial equations (2–4):
3.3. Influence of extraction parameters on TPC and antioxidant capacity
The influence of extraction temperature, extraction time, ethanol concentration and sample to solvent ratio on TPC and antioxidant capacity was shown in Table 3 and Figures 1–3. The results showed that in the tested narrow ranges of temperature (40–60 ℃), time (20–60 min), ethanol concentration (0–75%), TPC was not significantly affected when these parameters were changed (p > 0.05) with the exception of sample to solvent ratio. It was seen that temperature, ethanol concentration and sample to solvent ratio had significant effect on DPPH scavenging ability, while temperature, time and solvent concentration significantly affected FRAP (p < 0.05, Table 3). The data in Table 3 show minor interaction effect of the variables. The only interaction recorded are temperature and ethanol concentration for FRAP (p < 0.05, Table 3), and ethanol concentration and sample to solvent ratio for DRSC (p < 0.05, Table 3).
3.4. Optimisation and verification of the models
The optimal UAE conditions for maximum recovery of both phenolic compounds and antioxidant properties from L. intermedia were determined at fixed ultrasound power (100 W) and frequency (42 KHZ) to be temperature of 50 ℃, time of 60 min, ethanol concentration of 30%, and sample to solvent ratio of 2 g/100 mL. Verification experiments were conducted under these optimal UAE conditions to verify and confirm the adequacy of the models, and the results were presented in Table 4. Under the optimum conditions, the experimental values of TPC, DRSC and FRAP were 161.79 mg GAE/100 g, 32.30 mg TE/100 g and 87.77 mg TE/100 g, respectively, while the predicted values of respective assays were 162.34 mg GAE/100 g, 33.78 mg TE/100 g and 90.09 mg TE/100 g (Table 4). The results of mean comparison show no significant difference between experimental and estimated values of all assays (p > 0.05), revealing the reliability and adequacy of the developed models.
3.5. TPC and antioxidant activities of the fractions
The antioxidant activities and TPC of fractions derived from L. intermedia were assessed (Table 5). Ethyl acetate fraction demonstrated to exhibit the highest yield of phenolic content, over eighteen times higher than that of the aqueous fraction, and ten times higher than that of n-hexane fraction. Similar trend was found in the DRSC and FRAP assessments. Particularly, DRSC of the ethyl acetate fraction was over 53-fold and 81-fold greater than those of the aqueous and n-hexane fractions, respectively. The ethyl acetate was also found to possess the highest FRAP value compared to the aqueous and n-hexane fractions.
4.
Discussion
Determining the fit of the models strongly demonstrates the reliability of the developed mathematical models and the predicted conditions. It further confirms the strong correlation between the predicted and experimental values of the total phenolic content and antioxidant activity of L. intermedia.
The R2 (determination coefficient) is the ratio of the explained variation to the total variation and it measures the degree of fit [32]. The variance analysis for the TPC showed R2 (0.95), lack of fit test (0.091) and p-value < 0.0001, strongly indicating the model adequacy. The R2 (0.97) for DRSC also demonstrates a good fit between the actual and predicted values. Moreover, there was no significance in the lack of fit (0.077) indicating the model accuracy. The p-value was found to be lower than 0.0001, further providing firm evidence of the reliability of DRSC model. Regarding FRAP, an insignificant lack of fit (0.1951) and the low p-value of the model (<0.0001) showed a good fit and strong reliability of the model. The high R2 value (0.99) show a very strong correlation between the actual and predicted values of FRAP (99% matching).
The influence of four extraction parameters, including temperature, time, ethanol concentration and sample to solvent ratio on the phenolic content of L. intermedia was analysed. Temperature, time and solvent concentration had no significant impact on the total phenolic content (p > 0.05). However, the linear term of sample to solvent ratio, and the quadratic term of sample to solvent ratio were noticed to have statistical effect on TPC (p < 0.05). These terms were considered as the most influential variable for the extraction of phenolic compounds from L. intermedia. In particular, increasing the ratio of sample to solvent led to a significant decrease in TPC. The rest of the terms show no significant effect (p > 0.05). The result was agreed by previous studies done by Ahmed et al. [33] and Bamba et al. [34], who reported that an increase in the total phenolic content of koreeb seeds flour and blueberry pomace could be obtained when the solid-solvent ratio was decreased. These findings are in line with the mass transfer principle, the higher the volume of solvent is used, the greater the concentration gradient will be, creating a driving force for the transfer of the solutes from the sample matrix to the external solvent [35].
The linear (X1, X3, X4), quadratic (X12, X22, X32) and interaction (X3X4) effects had significant impact on the DRSC (p < 0.05). They demonstrated to be the most influential variables for the DPPH radical scavenging activity from L. intermedia. The effects of the rest of the factors were insignificant (p > 0.05). The DRSC increased with a decrease in the interaction of ethanol concentration and sample to solvent ratio (X3X4) and the quadratic effects of temperature and ethanol concentration, favoring the scavenging capacity of L. intermedia. This is indicated by negative coefficients. This finding was supported by the results of Ahmed et al. [33].
According to the results for FRAP, the linear (X1, X2, X3), quadratic (X12, X22, X32) and interaction (X1X3) effects had significant impact on the ferric reducing antioxidant ability. The rest of the terms were insignificant with p-values higher than 0.05. Negative coefficients for the linear effect of ethanol concentration (X3, ), the interaction between temperature and ethanol concentration (X1X3), and between time and ethanol concentration (X2X3), as well as the quadratic effect of temperature, time, and ethanol concentration denotes that decreases in these variables resulted in an increase in FRAP value. Mokrani and Madani [36] reported the response of FRAP under these interactions similar to our findings.
The 3D response surface plots were constructed using Equations (2-4) in order to provide a better understanding of the interaction between factors. The graphs were generated in JMP 15 statistical software by plotting the responses using the z-axis against two independent variables, while keeping the other independent variables constant. The interaction between sample to solvent ratio and ethanol concentration for DRSC, and temperature and ethanol concentration for FRAP, were the only significant effects recorded.
Solvent fractionation is an important step to separate the compounds based on the polarity. Here, the optimised conditions for the extraction of phenolic compounds and antioxidants were used to obtain the crude extract for the fractionation. The highest yield of TPC and the strongest antioxidant activity were found in the ethyl acetate fraction, while the values of TPC, DRSC and FRAP in the n-hexane fraction were substantially lower. Ethyl acetate is a polar solvent though less polar than water, but the high phenolic content and antioxidant activity contain in the ethyl acetate fraction would suggest that most of the bioactive compounds contained in this species are polar. This better explains the reason why low solvent concentration favor high phenolic content and ferric reducing antioxidant power in this current study. Similar results were reported by Hacke et al. [37] who found that ethyl acetate fractions derived from Cymbopogon citratus possessed significant quantities of phenolic acids (caffeic, chlorogenic, gallic and rosmarinic acids) and flavonoids (catechin, epicatechin, rutin, luteolin and apigenin). In the study done on Anthemis praecox Link aerial parts, Belhaoues et al. [38] similarly found ethyl acetate fraction possessed the highest amount of phenolic compounds. Mariem et al. [39] also reported that the highest antioxidant activity was found in the ethyl acetate fraction derived from medicinal halophyte Retama raetam. Hence, ethyl acetate fraction from L. intermedia could be further employed to isolate individual phenolic compounds and antioxidants.
Phenolic compounds exhibit free-radical scavenging properties, protecting algae thalli from negative effects of UV radiation [40]. High level of phenolic content is always coupled with antioxidant and antibacterial activities. In this study, there was a strong correlation found between the total phenolic content and the DPPH radical scavenging capacity with high R2 value (0.865). TPC and FRAP also had a close correlation with R2 value = 0.956. Similarly, phenolic compounds of Helicteres hirsuta leaf and stem were found to be well correlated with antioxidant capacity [30,31]. However, the study of Rahiman et al. (2013) [41] found that no correlation between the phenolic content and antioxidant activity of Ocimum sanctum (Lamiaceae), Cucumis sativus. (Cucurbitaceae), Capsicum frutescens (Solanaceae) and Coriandrum sativum. Terpinc and co-workers [42] reported similar findings when they investigated the correlation between total phenolic content and antioxidant capacity of oil cake extract. Variations in the correlation of phenolic contents and antioxidant properties can be explained due to the abundance and complexity of compounds contain in different materials. Also, it is difficult to predict this correlation as the result of the interactions that these compounds undergo under various extraction conditions [43].
5.
Conclusions
In conclusion, the models developed to predict the TPC and antioxidant activities were found to be highly reliable for predicting the responses. The optimisation process of the total phenolic content (TPC), DPPH radical scavenging capacity (DRSC) and ferric reducing antioxidant power (FRAP) from L. intermedia using ultrasound-assisted extraction (UAE) methods was fully investigated in this study. The optimum conditions using UAE were determined to be temperature of 50 ℃, time of 60 min, ethanol concentration of 30%, and sample to solvent ratio of 2 g/100 mL at a fixed UAE power and frequency at 100 W and 42 KHZ, respectively. It was found that only sample to solvent ratio had significant impact on the TPC, whereas temperature, ethanol concentration and sample to solvent ratio significantly affected DRSC. FRAP value was also influenced by extraction time and ethanol concentration. This is the first study reporting the optimisation of L. intermedia using UAE, which is very useful for further investigation on the species.
Acknowledgements
The authors gratefully thank the Project of Khanh Hoa Province (No. ĐT-2017-20902-ĐL to Dr. The Han Nguyen) for the financial support. The authors would like to thank all members of project KC.09.05/16-20 for their help in seaweed collection and preparation.
Conflict of interest
The authors declare no conflict of interest.