Int J App Pharm, Vol 11, Issue 3, 2019, 211-215Original Article


CENTRAL COMPOSITE DESIGN APPLIED IN HPLC OPTIMIZATION FOR ANALYSIS OF TARTRAZINE AND AURAMINE O IN POWDER DRINKS

ANNISA DYAH LESTARI1,2, ABDUL ROHMAN1,3*, SUDIBYO MARTONO1

1Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Universitas Gadjah Mada, Yogyakarta, 2The National Agency of Drug and Food Control, District of Banjarmasin 70124, South Kalimantan Republic of Indonesia, 3Institute of Halal Industry and System , Universitas Gadjah Mada, Yogyakarta.
Email: abdul_kimfar@ugm.ac.id

Received: 25 Feb 2019, Revised and Accepted: 06 Apr 2019


ABSTRACT

Objective: This study was intended to optimize reversed-phase high-performance liquid chromatography (RP-HPLC) method for the determination of Tartrazine (TAR) and Auramin O (AUO) in powder drinks using experimental design of central composite design (CCD) approach.

Methods: TAR and AUO in powder drink product has same properties, therefore both analytes were analysed using C18 column (XBridge Shield RP 18 250 mm x 4.6 mm i.d., 5 µm) using Shimadzu LC 20AD chromatograph equipped with photo-diode array (PDA) detector at 300-650 nm. Some factors responsible for RP-HPLC separation of TAR and AUO including the concentration of buffer, the ratio of mobile phase and flow rate were optimized using CCD. The responses evaluated were peak area, retention time, and tailing factor. The mobile phase used was acetonitrile and ammonium acetate buffer, and acetonitrile composition was optimized at 84-86% for separation of TAR and AUO, delivered at a flow rate of 0.8–1.2 ml/min, using ammonium acetate buffer at 19-21 mmol.

Results: CCD showed that separation of TAR and AUO was influenced by flow rate, the ratio of acetonitrile and ammonium acetate concentration. These factors affected significantly to retention time, peak area, and tailing factor. The optimal condition obtained based on CCD was flow rate of 1.2 ml/min, the ratio of acetonitrile 86%, and ammonium acetate concentration of 19 mmol.

Conclusion: CCD can be used to get optimum condition for analysis of TAR and AUO in powder drink product.

Keywords: Tartrazine, Auramin O, HPLC, CCD, Powder drink


INTRODUCTION

Tartrazine (TAR) is one of the synthetic dyes food additives (BTP) that are permitted to be used on food products to improve the appearance, colour, and texture of foods [1]. Tartrazine (TAR) is azo dyes and auramin O (AUO) is diphenyl methane dyes. The chemical structures of TAR and AUO were shown in fig. 1. TAR is allowed in a food product with a certain maximum value limit. AUO is one of the synthetic dyes that are prohibited to be used in food products. Several studies in various countries have shown that there are cases of counterfeiting of traded food products, including counterfeiting of added dyes because illegal synthetic dyes are cheaper than legal food colouring. Because of its similarity colour, TAR can be replaced by AUO [2].

The synthetic colorants including TAR and AUO are suspected to be unhealthy and unsafe substances for humans, and as a consequence, synthetic colorants became perceptible as undesirable or harmful by consumers [3, 4]. Therefore, the synthetic colours have been the subject of numerous toxicological investigations and their values are established by national and international legislation, especially for their use in food, drinks, drugs and cosmetics [5, 6]. TAR is dangerous if used over safety limits and AUO can cause toxicity and even death for consumers. Therefore, analytical methods capable of detecting and quantifying TAR and AUO must be developed in order to ensure the food safety. Among analytical methods, reversed phase HPLC with a variety of detectors was available for analysis of TAR and AUO.

Fig. 1: The chemical structures of Tartrazine (TAR) and Auramine O (AUO)

Simultaneous determination of TAR and AUO has been carried out by Tonogai et al. [7] using HPLC, however, there is no further studies related to simultaneous analysis of two synthetic dyes. Most quantitative analysis of synthetic dyes based on the types of legal or non-illegal colouring agents. The similar solubility of TAR and AUO poses a challenge in the separation of these two dyes simultaneously due to the proximity of polarity and the possibility of large counterfeiting which has the same colour when used in food products. Some methods have been reported for determination of TAR and AUO individually which included thin layer chromatography (TLC) [8], HPLC with the detector of photo-diode array [9-11], liquid chromatography-mass spectrometry (LCMS) [8] and FTIR spectroscopy [12]. The optimization method of HPLC for the simultaneous analysis of TAR and AUO using experimental design is very interesting. Experimental design is a tool having the ability to reveal possible interactions between variables, while saving time and simplifying work [13, 14]. Experimental designs have been widely used to determine the optimum conditions for chromatographic separation in the field of food and pharmaceuticals. Purba et al. [15] have optimized HPLC conditions to determine Acid Orange 7 and Sudan II in blusher product based on response surface methodology using box behnken design (BBD) approach. In this study, reversed-phase high-performance liquid chromatography (RP-HPLC) using experimental design of central composite design (CCD) approach was optimized for the separation of Tartrazine (TAR) and Auramin O (AUO) in powder drinks.

MATERIALS AND METHODS

Powder drink products were obtained from local markets in Yogyakarta. Reference standards of Tartrazine (CI 19140, Control Number: 110397), Auramin O (CI 41000, Control Number: B0114315) were acquired from the national agency of drug and food control (NADFC) of Republic of Indonesia. All solvents used for the mobile phase were of HPLC grade and obtained from E. Merck (Darmstadt, Germany). Aquabidest was obtained from Ikapharmindo (Indonesia).

Preparation of reference standards

An approximately of 5.00 mg of each TAR and AUO was accurately weighed using analytical balance (Metler Toledo MX5) with the sensitivity of 0.01 mg and was added into volumetric flask 5 ml. TAR and AUO were dissolved in 3 ml aquabidest, sonicated using sonicator (Elma ultrasonic, Germany) for 5 minute, and made to volume with aquabidest (5 ml) to get solution with a concentration of 1000 µg/ml.

Preparation of samples

An approximately of 100.0 mg of powder drink products was accurately weighed using analytical balance (Metler Toledo MX5) with the sensitivity of 0.1 mg, added with 0.5 ml of each standard solutions (TAR and AUO), added with 3 ml aquabidest, sonicated for 5 min, and added with aquabidest to volume 10 ml. The solution was filtered with PTFE 0.45 µm. The solution was injected into HPLC system.

HPLC instrumentation

TAR and AUO were analysed using chromatograph of Shimadzu LC 20AD chromatograph equipped with photo-diode array (PDA) (Shimadzu LC 20AD, M20A PDA Detector) at the wavelength of 300-650 nm. Separation of analytes was performed using C18 column (XBridge Shield RP 18 250 mm x 4.6 mm i.d., 5 µm). The mobile phase was modified from the method of determining Tartrazine [10], used water as solvent, the composition of acetonitrile was optimized at 84-86% for separation TAR and AUO, delivered at a flow rate of 0.8–1.2 ml/min, using ammonium acetate buffer at 19-21 mmol.

Experimental design using CCD

The most relevant multivariate techniques is response surface methodology (RSM), an optimization based on fit of the polynomial equation to data experiment. Symmetrical design of RSM, namely central composite design (CCD) and box behnken design (BBD) are frequently used in HPLC method optimization because they can resolve HPLC separation-related problems which the number of factors is higher than 2 [16,17]. CCD and BBD differs in the selection of experimental point, variables number, as well as number of run and block. Central composite design (CCD), which is a widely used the form of RSM, encompasses the advantages of factorial design [16]. The factors (independent variables) evaluated were flow rate of mobile phase (X1), the ratio of acetonitrile for separation of TAR and AUO (X2), and ammonium acetate buffer concentration (X3). While, the responses or dependent variables evaluated included retention time TAR and AUO (Y1 and Y2), peak area TAR and AUO (Y3 and Y4), tailing factor TAR and AUO (Y5 and Y6).

Table 1: Central-composite design using dependent variables of flow rate (ml/min) (X1), concentration of acetonitrile (%) (X2), and ammonium acetate concentration (X3) with response variables of retention time of Tartrazine (TAR) (Y1), retention time Auramine O (AUO) (Y2), peak area TAR (Y3), peak area AUO (Y4), tailing factor TAR (Y5) and tailing factor AUO (Y6) used in HPLC optimization for separation of TAR and AUO

Std Run Dependent variables Responses
Flow Rate (mL/min) (X1) Conc. ACN (%) (X2) Conc. Ammonium Acetate (mM) (X3) Retention time TAR (Y1) Retention time AUO (Y2) Area TAR (Y3) Area AUO (Y4) TF TAR (Y5) TF AUO (Y6)
15 1 1 85 20 2.074 3.843 2607452 7382238 1.481 0.951
13 2 1 85 18.32 2.063 3.871 2607739 7379260 1.464 0.983
7 3 0.8 86 21 2.585 4.777 3253099 9244923 1.541 0.935
8 4 1.2 86 21 1.74 3.182 2160304 6139588 1.54 0.966
6 5 1.2 84 21 1.744 3.175 2162932 6115172 1.48 0.943
19 6 1 85 20 2.08 3.82 2601642 7378946 1.485 0.945
20 7 1 85 20 2.078 3.857 3083444 7368462 1.374 1.133
4 8 1.2 86 19 1.729 3.207 2169451 6173501 1.52 0.989
1 9 0.8 84 19 2.588 4.785 3250445 9176912 1.5 0.93
9 10 0.66 85 20 3.133 5.762 3952517 11122714 1.539 0.919
17 11 1 85 20 2.08 3.861 3079782 7341409 1.382 1.138
5 12 0.8 84 21 2.603 4.75 3246682 9133512 1.502 0.915
3 13 0.8 86 19 2.577 4.8 3261900 9184790 1.514 0.95
12 14 1 86.68 20 2.073 3.838 2604962 7358650 1.541 0.96
2 15 1.2 84 19 1.735 3.202 2169033 6160740 1.474 0.96
18 16 1 85 20 2.078 3.86 3075765 7314584 1.386 1.144
11 17 1 83.32 20 2.093 3.837 2606311 7365544 1.482 0.92
10 18 1.34 85 20 1.557 2.857 1944879 5537886 1.507 0.965
16 19 1 85 20 2.077 3.862 3069520 7273512 1.391 1.149
14 20 1 85 21.68 2.1 3.757 2595196 7359074 1.525 0.934

Conc. = concentration; TAR = tartrazine; AUO = Auramine O; ACN = acetonitrile.

Data analysis

All experiments using CCD along with statistical parameters were performed using Design-Expert version 8.0.4.1. The responses evaluated were retention time, peak area, tailing factor of TAR and AUO. Factors (independent variables) significantly affected the responses (dependent variables) if R2≥ 0.8 and Adjusted R²>0.8. The difference between predicted R² with the adjusted R² must be less than 0.2. The confirmation of optimal method was performed using six injection replicates. The statistical test of independent t-test used for comparing results obtained from CCD and from actual experiments was carried out using Minitab software version 17 (Minitab Corp., USA).

RESULTS AND DISCUSSION

Quantitative analysis of different dyes is most often performed in reversed phase (RP) or ion pair (IP) systems, while is usually based on measurements via UV–VIS detection, especially using diode array detector (DAD). The similar polarity between TAR and AUO might be copied by isocratic elution method, therefore the separation of TAR and AUO is in one condition. Reversed phase (C18) column retained TAR and AUO in high concentration of non-polar solvent, therefore an experimental design approach was used. Central-composite design (CCD) was used for HPLC separation of TAR and AUO. CCD was performed using 20 runs, applying 3 independent variables (factors) namely flow rate (X1), ratio of acetonitrile for separation of TAR and AUO (X2), and ammonium acetate buffer concentration (X3) along with response variables of retention time TAR (Y1), retention time AUO (Y2), peak area TAR (Y3), peak area AUO (Y4), tailing factor TAR (Y5) and tailing factor AUO (Y6). CCD using these factors and responses resulted during optimization were compiled in table 1.

Based on analysis of variance (ANOVA) results, the equation obtained using X1, X2, and X3 as independent variables with the retention time of TAR (Y1) as the response was:

Y1=–0.065364–7.77684X1+0.143565 X2+0.118803 X3+0.011875X1X2–0.001875X1X3–0.000625 X2X3+2.29397X120.00872 X22–0.001402 X32 (Adj. R20.9954) (Eq.1).

The statistic results revealed that adjusted R2 obtained was>0.8 (acceptable) [18], exhibiting that the experimental model was a good fit using polynomial equation. The difference between predicted R² with the adjusted R² in all responses was less than 0.2. Eq.1 informed that variables of flow rate (mL/min) (X1), have a negative effect on the retention time of TAR, while ratio of ACN (%) (X2) and ammonium acetate concentration (X3) have a positive effect on retention time of TAR.

The variables of X1, X2, and X3, quadratic form of X1 contributed significantly for response of Y1 (P<0.05) based on one way ANOVA results. The variables of X1 affected negatively, meaning that the increased levels of flow rate (X1) would decrease the retention time of TAR (decreased sensitivity), while the increased ratio of acetonitrile (X2) and ammonium acetate concentration (X3) could increase the retention time of TAR. Contour plot of retention time TAR along with 3D surface graph was shown in fig. 2.

Similarly, the equation for retention time AUO (Y2) using multiple linear regression were:

Y2=–41.79701–10.59627X1+1.14694X2+0.430479X3–0.018750X1X2+0.03750X1X3+0.01750X2X3+4.00147X12-0.006818X22-0.015127X32(Adj. R²of 0.9961) (Eq.2).

The contour plot along with along with 3D surface graph of retention time of AUO was shown in fig. 3. Statistic parameter of Y2 revealed adjusted R2 (Adj. R2) was>0.8 (acceptable) [17, 18] exhibiting that the experimental model was a good fit using the polynomial equation. Based on ANOVA results from variables of X1, X2, and X3, the quadratic form of X1 contributed significantly for the response of Y2 (P<0.05). The variables of X1 affected negatively, meaning that the increased levels of flow rate (X1) would decrease the retention time of AUO (decreased sensitivity), while the increased ratio of acetonitrile (X2) and ammonium acetate concentration (X3) could increase the retention time of AUO.

Fig. 2: The contour plot of retention time (in minute) of tartrazine (TAR) [A] and 3D surface graph of retention time of TAR [B] as a results of variables of flow rate (ml/min), concentration of acetonitrile (%), and ammonium acetate concentration

Fig. 3: The contour plot of retention time of Auramine O (AUO) [A] and 3D surface graph of retention time of AUO [B] as a results of variables of flow rate (ml/min), the concentration of acetonitrile (%), and ammonium acetate concentration

Eq. 3 revealed the response of peak area TAR (Y3). The statistic results for Y3 informed that adj. R2was>0.8. The variables of X1, X2, X3 and X4, linear form of X1 and quadratic form of X2 and X3 contributed significantly for the response of Y3 (P<0.05). The variables of X1, X2 and X3 affected positively, meaning that the increased levels of flow rate (ml/min) (X1), acetonitrile (%) (X2) and ammonium acetate concentration (mM) (X3) would increase peak area of TAR (increased sensitivity).

Y3=–8.39256 x 108–2.27267 x 106X1+1.88030 x 107X2+4.55925 x 106X3–12551.25000X1X2–1677.50000X1X3–1010.50000X2X3+2.7207 x 105X12–1.10407 x 105X22–1.1881 x 105 X32 (Adj. R²of 0.8836) (Eq. 3).

Eq. 4 showed the correlation between response of peak area of Auramine O (AUO) and independent variables of X1, X2, and X3 along with its interaction. The statistic results for Y4 showed that Adj. R2 obtained was in the acceptable limits [18]. The ANOVA results revealed that variables of X1, X2, X3 and X4, quadratic form X1 contributed significantly for the response of Y4 (P<0.05).

Fig. 4: The contour plot of tailing factor of tartrazine (TAR) [A] and 3D surface graph [B] as a results of variables of flow rate (ml/min), the concentration of acetonitrile (%), and ammonium acetate concentration

Y4=2.64137 x 107–1.93302 x 107X1+1.41257 x 105X2–1.18438 x 106X3–51320.00000X1X2–60133.75000X1X3+14398.50000X2X3+8.50377 x 106X12–2160.58476X22+339.03771X32(Adj. R² of 0.9956)(Eq. 4)

Eq. 5 and 6 corresponded to the response of tailing factor of TAR (Y5) and AUO (Y6). The statistic results for Y5 revealed that Adj. R2 obtained was<0.8, which was not acceptable. Based on ANOVA results variables of X1, X2, X3 and X4, quadratic form of X1, X2 and X3 contributed significantly for the response of Y5 (P<0.05). Based on ANOVA results in variables of X1, X2, X3 and X4, quadratic form of X1 and X2 contributed significantly for the response of Y6 (P<0.05).

Y5=254.92078–4.62826X1–5.59046X2–1.44966X3+0.033X1X2–0.001875X1X3+0.004875X2X3+0.907344X12+0.032228X22+0.026218X32 (Eq. 5).

(Adj. R² of 0.5453)

Y6=–352.13095+1.84260X1+7.89106X2+1.64942X3+0.007500X1X2–0.006250X1X3–0.000750X2X3–1.13995X120.046305X220.039765X32

(Adj. R² of 0.2471) (Eq. 6)

Fig. 4 and fig. 5 showed the contour plot along with along with 3D surface graph of tailing factor of TAR and AUO.

Fig. 5: The contour plot of tailing factor of Auramine O (AUO) [A] and 3D surface graph [B] as a results of variables of flow rate (ml/min), the concentration of acetonitrile (%), and ammonium acetate concentration

Central composite design using three optimum factors namely flow rate (1.2 ml/min), ratio of acetonitrile (86 %) and ammonium acetate concentration (19 mmol) was successfully used to get an optimum condition of HPLC method for analysis of TAR and AUO in powder drink samples. The HPLC chromatogram obtained using this condition was shown in Fig.6. It is clear that both TAR and AUO were clearly separated using optimum condition suggested by CCD.

Fig. 6: Separation of tartrazine (TAR) and Auramine O (AUO) using HPLC condition as suggested by central-composite design. See text for HPLC condition

CONCLUSION

CCD design can be used to get optimum condition for analysis of TAR and AUO in powder drink product. The optimum conditions suggested for separation TAR and AUO based on CCD was the mobile phase containing ACN 86% with a flow rate of 1.2 ml/min, with ammonium acetate buffer concentration of 19 mmol.

ACKNOWLEDGMENT

The author acknowledges to Indonesian National Agency of Drug and Food Control, Indonesian National Agency of Drug and Food Control in Banjarmasin, South Kalimantan and Indonesian National Agency of Drug and Food Control in Yogyakarta for financial support and instrument facilities.

AUTHORS CONTRIBUTIONS

ADL conducted research activity, compiled data, and prepared manuscript. AR and SM designed research activities, prepared manuscript and made critical thinking on the manuscript.

CONFLICT OF INTERESTS

The authors have declared “no conflicts of interest with respect to the research, authorship, and/or publication of this article”.

REFERENCES

  1. Kucharska M, Grabka J. A review of chromatographic methods for determination of synthetic food dyes. Talanta 2010; 80:1045–51. 

  2. Bachalla N. Identification of synthetic food colors adulteration by paper chromatography and spectrophotometric methods. IAIM 2016;3:182-91.

  3. Kamel MM, El-lethey SH. The potential health hazard of tartrazine and levels of hyperactivity, anxiety-like symptoms, depression and anti-social behaviour in rats. J Am Sci 2011;7:1211–8.

  4. Sha O, Zhu X, Feng Y, Ma W. Determination of sunset yellow and tartrazine in food samples by combining the ionic liquid-based aqueous two-phase system with high-performance liquid chromatography. J Anal Methods Chem 2014;2014:964273.

  5. Khera KS, Munro IC. A review of the specifications and toxicity of synthetic food colors permitted in Canada. CRC Critical Rev Toxicol 1979;6:81–133.

  6. Tonogai Y, Kingkate A, dan Halilamian C. Quantitative determination of colorants in dried shrimp and shrimp paste using ion-exchange extraction and high performance liquid chromatography. J Food Protec 1983;46:592–5.

  7. Tatebe C, Zhong X, Ohtsuki T, Kubota H, Sato K, Akiyama H. A simple and rapid chromatographic method to determine unauthorized basic colorants (rhodamine B, auramine O, and pararosaniline) in processed foods. Food Sci Nutr 2014; 2:547–56.

  8. Li J, Ding X, Zheng J, Liu D, Guo F, Liu H. Determination of synthetic dyes in bean and meat products by liquid chromatography with tandem mass spectrometry: other techniques. J Sep Sci 2014;37:2439–45.

  9. Rejczak T, Tuzimski T. Application of high-performance liquid chromatography with diode array detector for simultaneous determination of 11 synthetic dyes in selected beverages and foodstuffs. Food Anal Methods 2017;10:3572–88.

  10. Minioti KS, Sakellariou CF, Thomaidis NS. Determination of 13 synthetic food colorants in water-soluble foods by reversed-phase high-performance liquid chromatography coupled with diode-array detector. Anal Chim Acta 2007;583:103–10.

  11. Sulekova M, Hudak A, Smrcova M. The determination of food dyes in vitamins by RP-HPLC. Molecules 2016;21:1368-72.

  12. Karimi S, Feizy J, Mehrjo F, Farrokhnia M. Detection and quantification of food colorant adulteration in saffron sample using chemometric analysis of FT-IR spectra. RSC Adv 2016;6:23085–93.

  13. Hibbert DB. Experimental design in chromatography: a tutorial review. J Chromatogr B 2012;910:2–13.

  14. Bezerra MA, Santelli RE, Oliveira EP, Villar LS, Escaleira LA. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta 2008;76:965–77.

  15. Purba N, Rohman A, Martono S. The optimization of HPLC for quantitative analysis of acid orange 7 and Sudan II in cosmetic products using box behnken design. Int J Appl Pharm 2019;11:130-7.

  16. Khanam N, Alam MI, Iqbal QMA, Ali MY, Siddiqui AR. A review on the optimization of drug delivery system with experimental designs. Int J Appl Pharm 2018;10:12-7.

  17. Siregar C, Prabaningdyah NK, Choiri S, Riyanto S, Rohman A. Optimization of HPLC using central composite design for determination of curcumin and demethoxycurcumin in the tablet dosage form. Dhaka Univ J Pharm Sci 2018;16:137-45.

  18. Setyawan EI, Setyowati EP, Rohman A, Nugroho AK. Central composite design for optimizing extraction of EGCG from green tea leaf (Camellia sinensis L.). Int J Appl Pharm 2018;10:211–6.