Application of response surface methodology for optimization of metal–organic framework based pipette-tip solid phase extraction of organic dyes from seawater and their determination with HPLC

This paper describes the application of response surface methodology (RSM) to develop a miniaturized metal organic framework based pipette-tip solid phase extraction for the extraction of malachite green (MG), rhodamine B (RB), methyl orange (MO) and acid red 18 (AR) dyes from seawater samples and their determination by high performance liquid chromatography. The effects of various parameters such as pH of the sample solution, type and amount of added salt, type and volume of eluent solvent, concentration of surfactant (triton X-114), sample volume, and number of cycles of extraction and desorption were investigated and optimized by two methods of one-variable-at-a-time and RSM based on Box–Behnken design. Under optimum conditions, the linear range of the method was 0.5–200.0 µg/L for RB and MG and 1.0–150.0 µg/L for AR and MO. Limits of detection of the analytes were obtained in the range of 0.09–0.38 µg/L. Reproducibility of the method (as RSD %) was better than 6.4%. The method has been successfully used for analysis of four dyes in seawater of Chabahar Bay. Electronic supplementary material The online version of this article (10.1186/s13065-019-0572-0) contains supplementary material, which is available to authorized users.


Introduction
Rhodamine B (RB) (Fig. 1a), is among the oldest and most commonly used synthetic dyes that have been recently identified as possible illegal additives in foods exported from European Union and China [1]. It belongs to the class of xanthenes dyes, a basic red cationic dye that is highly soluble in water, methanol and ethanol. This dye is used widely as a colorant in textiles and plastic industries. RB is harmful if swallowed with human beings and cause irritation to the skin, eyes and respiratory tract. Also, it has been shown to have carcinogenicity, reproductive and developmental toxicity, neurotoxicity and chronic toxicity towards human and animals [2]. Malachite green (MG, Fig. 1b), although a forbidden dye, has been widely applied illegally as a fungicide and parasitical and in the fish industry as an antimicrobial, antiseptic and ectoparasitic agent, because of its high efficiency and low cost [3][4][5]. Acid red 18 (AR, Fig. 1c), is a popular food color, not toxic but can be harmful if used in excess [6,7]. Methyl orange (MO, Fig. 1d), have many application as textile dyeing stuff and staining agents in laboratories [8]. These dyes are of the most abundant applied dying agents throughout the world and therefore can find their way to the environmental sources such as seawater as hazardous pollutants [9,10].
Metal-organic frameworks (MOFs) are three dimensional crystalline porous materials having different geometries and functional groups within the channels/cavities, which are synthesized using mixing organic linkers and metal salts, often under hydrothermal or solvothermal conditions. The unique characteristic of the hybrid solids are adjustable pore-sizes and controllable structural properties, extra ordinarily large porosity, low density and their very high surface areas. MOFs have been considered as promising candidate materials for different applications including adsorption, removal, separation, selective extraction and pre-concentration of various analytes [27,28].
Pipette-tip solid phase extraction (PT-SPE) is a representative SPE technique because of its miniature device and use of reduced amount of reagents and less time consumption [29]. For PT-SPE, an ordinary pipette tip acts as the extracting column, packed with sorbent. This technique has been successfully used in many applications [30,31]. Response surface methodology (RSM) can be summarized as a compilation of statistical tools and method for constructing and exploring estimated function relationship between a response variable and set of design variable. It is the collection of mathematical and numerical methods that are suitable for modeling and analysis of the problems having numerous variables influencing the response, and objective is to optimize the response. The most extensive application of RSM can be found in industrial world, where a number of input variables affect some performance measures, called the response, in ways which are not easy or unfeasible to depict by a rigorous mathematical formulation [32,33].
In the present work, we synthesized a novel Co-MOF and used it for simple, fast and sensitive PT-SPE of RB, MG, AR and MO organic dyes in seawater samples and their determination with HPLC. Parameters affecting PT-SPE were optimized by two methods of one variable-at-atime and RSM, based on Box-Behnken design. This is the first report on using Co-MOF for pipette-tip solid phase extraction of dyes in Chabahar Bay (Oman Sea).

Apparatus
A Knauer HPLC (Germany) equipped with a EA4300F smart line pump and a smart line auto sampler 3950, was used for all analyzes. Detection system was a diode array spectrophotometer, used at wavelengths of 448 nm for MO, 510 nm for AR, 555 nm for RB and 618 nm for MG. Analytical column was a 250 × 4.6 mm Eurospher 100-5 C 18 utilizing the same pre-column. ChromGate V3.1.7 software was used for chromatographic data handling. The injection loop volume was 20 µL. A model 630 Metrohm (Switzerland) pH meter was employed for pH determination.

Reagents
All dyes and chemical reagents were of analytical grade and were purchased from Merck KGaA (Darmstadt, Germany). The HPLC grade methanol, acetonitrile and water were also obtained from the same company. Milli-Q ® water (18.3 MʹΩ/cm) was used throughout the run after filtering through 0.22 mm Nylon membrane. Triton X-114 (5% v/v) solutions was prepared at 70:30 (v/v) water/methanol and used as the surfactant. Stock solution of each dye with a concentration of 500 mg/L was prepared with dissolving of 0.0500 g of each dye in distilled water in 100 mL flasks. Working solutions were prepared daily by proper dilution of stock solutions.

Synthesize and characterization of Co-MOF adsorbent
Synthesize of Co-MOF was according to the work of Sargazi et al. [34]. Briefly, 5.62 mmol of cobalt nitrate and 1.85 mmol of pyridine 2, 6-dicarboxylic acid were dissolved in 14 mL of ethanol. Obtained solution was transferred into a Teflon reactor with a tight cap and kept for 7 h at 85 °C. The product was washed with dimethylformamide. After mixing and dissolving the reactants, the clear solution radiated in the ultrasound bath for 13 min at working condition of 160 W, 1 kJ, and 21 kHz. Synthesized adsorbent was stored in 4 °C. Scanning electron microscopy ( Fig. 2) showed an average size of 17 µm for synthesized MOF. By BET (Brunauer, Emmett and Teller), the specific surface area of Co-MOF was determined 3000 m 2 /g.

Pipette-tip solid phase extraction procedure
PT-SPE of dyes was performed using an Extra GENE tip mounted on a variable 150 µL volume pipettor (Dragon Labs, USA). 8 mL of an aliquot of sample solution containing appropriate amounts of dye was transferred into a 10 mL flask and proper amount of triton X-114 (0.15% v/v for MG and AR and 0.20% v/v for RB and MO) and 150 mg of KCl was added. Then pH of solutions were adjusted to the desired value (pH = 3.0 for MG and RB, 6.0 for AR and 6.9 for MO) with drop-wise addition of either 1 mol/L of HCl or 1 mol/L of NaOH. PT-SPE carried out by loading the sample solution into the cartridge and washing out with 0.5 mL of methanol-water (1:1). After the analyte retained on the MOF sorbent, it was eluted using 300 µL (for RB, MO and AR) and 250 µL (for MG), of methanol contain 5% acetic acid. Finally eluted solvent was filtered through a 0.45 µm filter and was injected into HPLC for analysis.

Chromatographic conditions
Various mobile phases were investigated consisting of methanol, acetonitrile and water in different combinations and pH settings. Finally a gradient of 85% B at 0-3.5 min and 100% B at 3.5-10 min was selected; in which eluent A was water and eluent B was acetonitrile which was adjusted to the pH 5.25 using acetic acid at a flow rate of 0.8 mL/min. The column oven temperature was maintained at room temperature and the mobile phase was degassed using a stream of helium prior to use.

Optimization of MOF-PT-SPE
In order to achieve the best efficiency of the MOF-PT-SPE, different factors affecting extraction efficiency were optimized using two methods of one-variable-at-a-time and RSM based on Box-Behnken design. A standard aqueous solution at concentration of 150.0 µg/L for AR and MO and 250 µg/L for RB and MG was used for optimization experiments. Each experiment repeated at least three times.

Effect of amount of sorbent
The effect of amount of MOF for preconcentration and determination of selected dyes in pipette tip was investigated in the range 1.0-2.5 mg. The results showed that the percent of extraction increases to 2.0 mg of MOF and then the recovery decreases. So, 2.0 mg of sorbent in pipette tip was used for further experiments.

Effect of type and amount of salt
To investigate effect of type and amount of salt on extraction efficiency of dyes, NaCl, KCl and Na 2 SO 4 as common salts were selected and used for MOF-PT-SPE of dyes. Among them, KCl improved the extraction better than the other salts and hence selected as spiked salt in further works. To study effect of amount of KCl on extraction efficiency, various brine sample solutions containing different quantity of KCl in the range of 25-200 mg were prepared. The results indicated that the extraction efficiency of dyes is quantitative for amount of KCl greater than 200 mg. Hence, next runs were performed with saturation of the samples using 200 mg of KCl.

Effect of concentration of triton X-114
The concentration of triton X-114 as surfactant can effect on the extraction efficiency of dyes by MOF-PT-SPE; so, we tried to optimize its concentration. We found that by increasing the concentration of the surfactant, the extraction efficiency was also increases, but in the amounts more than 0.20 (for RB and MO) and 0.15 (for MG and AR) %v/v of triton X-114, a decrease in the extraction efficiency of dyes was observed. This is probably due to the dilution of the analytes in larger volumes of the surfactant.

Box-Behnken design
Four factors in three levels were utilized to consider and optimize the process factors which potentially have an effect on the extraction efficiency of the analytes by MOF-PT-SPE. The investigated factors and input variable for four dyes were pH (X 1 or A), eluent volume (µL) (X 2 or B), number of extraction cycles (X 3 or C), and number of eluent cycles (X 4 or D). Table 1 shows the levels of these variable which were coded as − 1 (low), 0 (central point) and 1 (high). The design of real runs is given in Additional file 1: Table S1.
The following quadratic equation (Eq. 1) can be used to explain the behavior of the system: In Eq. 1, Y is output; i.e. is the response of HPLC, which is the dependent variable; i and j are the index numbers of the model; β 0 is the free or offset term, called intercept term; X 1 , X 2 , …, X k are coded independent variables; B i is the first-order (linear) main effect,; B ii is the quadratic (squared) effect; β ij is the interaction effect; and ε is the random error which allows for description or uncertainties between predicted and determined values [35].
For 4 selected dyes, subsequent equations explain the relationship between the four variables and response of HPLC (output, Y): (1)  .21 for RB, MG, AR and MO, respectively) and a very low probability value (p), the ANOVA of the regression models shows that quadratic models are also significant. In Fig. 3 two dimensional response surfaces as the function of other variable are shown.

Analytical performance Linear range, limit of detection and enrichment factor
The linearity of the proposed method was examined under the optimized conditions. Over a concentration range of 0.5-200.0 µg/L for RB and MG; and 1.0-50.0 µg/L for AR and MO, the calibration curve was linear. The least square equations over the dynamic linear range are indicated in Table 2. The limit of detection of the method for all target analytes was calculated using 3S b /m equation (where S b is the standard deviation of 7 consecutive measurements of the blank and m is shop of the calibration curve) and was 0.09, 0.17, 0.33 and 0.38 for RB, MG, AR and MO, respectively.
To achieve a high enrichment factor (EF), the effect of the sample volume on the recovery of dyes was investigated in the range of 2 to 10 for all of the analytes. The results showed that the extraction efficiency of selected dyes were very efficient (> 97%) in a sample volume of 8 mL and at the eluent solvent of 300 µL (for RB, MO and AR), 250 µL (for MG).
By mathematical calculation from the volume ratio of the sample to extracting phase, and a recovery of 97%, it is expected to have a pre-concentration factor of 26.6 (for RB, MO and AR), 32.0 (for MG). The real enrichment factors were experimentally achieved were 25.8 (for RB, MO and AR), 31.0 (for MG), and 28.2 (for AR). Table 3 compares the characteristic data of present method with those reported in the literature.

Determination of dyes in seawater samples
The performance of proposed method was investigated by extraction and determination of dyes in five seawater samples taken from different spots of Oman Sea, close to Chabahar Bay (southern-east part of Iran). No salt was added to the real samples since they are fully salt saturated by themselves. Since no dyes could be detected in them, to evaluate the effect of sample media on recovery, they were spiked at the concentration of 10 µg/L with dyes. Results are presented in Table 4. Figure 4 shows sample chromatograms obtained for the analysis of seawater sample, taken from station 3. Significant raise of signal can be observed. Reproducibility of the method (as RSD%) was found to be in the range of 0.7-4.6% for RB, 0.6-4.0 for MG, 1.9-6.4 for MO and 0.7-6.3 for AR.  These results show that the proposed technique can be used for determination of selected dyes in very complicated matrices such as seawater.

Conclusion
In this paper, the combination of pipette tip solid phase microextraction by means of a novel metal organic framework with HPLC was successfully used for the analysis of dyes in seawater. This technique has enough simplicity and sensitivity to be employed for routine analysis of dyes in such complicated media. An additional advantage of the suggested technique is its easy operation. Besides, the technique is feasible for high number of samples due to its short processing time.

Additional file
Additional file 1: