Quality by design-based capillary electrophoresis method development for the quantitative analysis of four iridoid compounds in Gentiana macrophylla Radix

In this study, the capillary electrophoresis-photodiode array detector was employed for the analysis of four iridoid compounds in Gentiana macrophylla Radix (RGM), and the method was optimized based on the concept of analytical quality by design (AQbD). The peak areas relative standard deviation (n = 3) and resolution of the four analytes were selected as critical method attributes. Fractional factorial design test combined with Pareto analysis were employed to screen critical method parameters (buffer concentration, pH, sodium dodecyl sulfate [SDS] micelle concentration, temperature, and voltage). Subsequently, three main factors (buffer concentration, buffer pH, and SDS concentration) were selected by central composite design test for constructing the design space. The optimal separation conditions as follows: capillary column (50.2 cm × 50 µm, detection length 40 cm). Working background electrolyte consisted of 51 mmol/L borax solution (pH = 9.47) and 40 mmol/L SDS. The samples were injected by pressure (5 s at 0.5 psi) and the detection was performed at 254 nm. Applied voltage was 20 kV and column temperature was 23°C. The developed method is rapid and reliable for the quantitative analysis of four iridoid compounds in RGM, providing a reference for the application of AQbD concept in the analysis of natural products.     

Unique green chromatography method for the determination of serotonin receptor antagonist (Ondansetron hydrochloride) related substances in a liquid formulation,robustness by quality by design-based design of experiments approach

Serotonin receptor antagonist drug Ondansetron hydrochloride injectable formulation containing all related substances was identified and quantified by a single, simple, sensitive, eco-friendly, and green high-performance liquid chromatography method. The disseverment of all impurities was achieved with the Discovery Cyano (250 × 4.6) mm, 5 μm column. The gradient program was composed of pH 5.7 phosphate buffer as mobile phase A and acetonitrile as mobile phase B. The flow rate, column compartment temperature, and detection wavelengths were 0.9 mL/min, 30°C, and 216 nm, respectively. The method was validated as per current regulatory guidelines. The obtained %relative standard deviation for the precision results was between 0.55 and 2.72% for all impurities. The correlation coefficient values from the linearity experiment for impurities and analyte were more than 0.995. The accuracy results were obtained between 88.4 and 113.0% for all impurities. Both sample and standard solutions showed 24 h stability at benchtop and refrigerator conditions. All impurities and analytes met the specificity and mass balance for all forced degradation conditions. Quality-by-design-based design of experiments was utilized to establish the method's robustness. Method greenness was assessed by using the current advanced tool green analytical procedure index, National Environmental Methods Index, and analytical eco-scale.     

Applications of Monte-Carlo simulation and chemometric techniques for development of bioanalytical liquid chromatography method for estimation of rosuvastatin calcium

In the present work, we investigated the development of a bioanalytical HPLC method of rosuvastatin (RSV) calcium as per the Quality by Design (QbD)-based systematic chemometric tools. At first, the method objectives were framed and critical analytical attributes (CAAs) were chosen. Risk assessment and factor screening was performed using Hybrid Risk Matrix and Plackett–Burman design for identifying vital factors influencing the critical method parameters (CMPs). Monte-Carlo simulation analysis was conducted which confirmed excellent process robustness (Ppk >1.33) for the studied ranges of CMPs. Furthermore, systematic method development was carried out using custom experimental design, where mobile phase ratio, pH, and injection volume were taken as CMPs at three levels. The obtained trials were evaluated for peak area, retention time, theoretical plates, and peak tailing as CAAs. Mathematical response surface modeling was carried out and optimal chromatographic solution was identified using response optimizer plots. Method transfer was made to bioanalytical scale for estimation of the analyte in rat plasma samples. Extensive method validation was performed as per the ICH Q2 guideline, which indicated validation parameters within the acceptable limits. Overall, the studies construed successful development of QbD compliant HPLC method of rosuvastatin with potential utility bioanalytical testing.     

Analytical QbD-based systematic bioanalytical HPLC method development for estimation of quercetin dihydrate

The current work entails development of rapid, sensitive, and inexpensive high-performance liquid chromatographic method of quercetin dihydrate using the quality by design approach. Quality target method profile was defined and critical analytical attributes (CAAs) were earmarked. Chromatographic separation was accomplished on a C₁₈ column using acetonitrile and ammonium acetate buffer (35:65) %v/v (containing 0.1% acetic acid, pH 3.5) as mobile phase at 0.7?mL/min flow rate with UV detector at 237?nm. Screening studies using fractional factorial design revealed that organic modifier, injection volume, column temperature, and buffer strength have significant influence on method CAAs, namely, peak area, retention time, and peak tailing. The critical method parameters were systematically optimized using Box–Behnken design. Response surface mapping was used along with numerical optimization and desirability function for identifying the optimal chromatographic conditions. Linearity was observed in the drug concentration ranging between 2 and 50?µg/mL. Accuracy analysis revealed mean % recovery between 93.6 and 96.2%, while precision study revealed mean % recovery between 93.7 and 96.5%. Limits of detection and quantification of the developed method were found to be 12.1 and 36.6?ng/mL. Overall, the studies construed successful development of chromatographic method of quercetin with enhanced method performance.     

QbD‐oriented development and validation of a bioanalytical method for nevirapine with enhanced liquid–liquid extraction and chromatographic separation

The present studies describe the systematic quality by design (QbD)‐oriented development and validation of a simple, rapid, sensitive and cost‐effective reversed‐phase HPLC bioanalytical method for nevirapine in rat plasma. Chromatographic separation was carried out on a C₁₈ column using isocratic 68:9:23% v/v elution of methanol, acetonitrile and water (pH 3, adjusted by orthophosphoric acid) at a flow rate of 1.0 mL/min using UV detection at 230 nm. A Box–Behnken design was applied for chromatographic method optimization taking mobile phase ratio, pH and flow rate as the critical method parameters (CMPs) from screening studies. Peak area, retention time, theoretical plates and peak tailing were measured as the critical analytical attributes (CAAs). Further, the bioanalytical liquid–liquid extraction process was optimized using an optimal design by selecting extraction time, centrifugation speed and temperature as the CMPs for percentage recovery of nevirapine as the CAA. The search for an optimum chromatographic solution was conducted through numerical desirability function. Validation studies performed as per the US Food and Drug Administration requirements revealed results within the acceptance limit. In a nutshell, the studies successfully demonstrate the utility of analytical QbD approach for the rational development of a bioanalytical method with enhanced chromatographic separation and recovery of nevirapine in rat plasma. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of chemometric approach for development and validation of high performance liquid chromatography method for estimation of ropinirole hydrochloride

A systematic Quality by Design approach was employed for developing an isocratic reversed‐phase liquid chromatographic technique for the estimation of ropinirole hydrochloride in bulk drug and pharmaceutical formulations. LiChrospher RP 18‐5 Endcapped column (25 cm × 4.6 mm id) at ambient temperature (25 ± 2°C) was used for the chromatographic separation of the drug. The screening of factors influencing chromatographic separation of the active pharmaceutical ingredient was performed employing fractional factorial design to identify the influential factors. Optimization of the selected factors was carried out using central composite design for selecting the optimum chomatographic conditions. The mobile phase employed was constituted of Solvent A/Solvent B (65:35 v/v) (Solvent A [methanol/0.05 M ammonium acetate buffer, pH 7, 80:20 v/v] and Solvent B [high performance liquid chromatography grade water]) and used at 0.6 mL/min flow rate, while UV detection was performed at 250 nm. Linearity was achieved in the drug concentration range 5–100 µg/mL (R² = 0.9998) with limits of detection and quantification of 1.02 and 3.09 µg/mL, respectively. Method validation was performed as per ICH guidelines followed by forced degradation studies, which indicated good specificity of the developed method for detecting ropinirole hydrochloride and its possible degradation products in the bulk drug and pharmaceutical formulations.     

Chaotropic chromatography method development for the determination of aripiprazole and its impurities following analytical quality by design principles

In this paper, development of robust and reliable chaotropic chromatography method for the determination of aripiprazole and its impurities, following Analytical Quality by Design principles is presented. The efficient baseline separation and accurate determination of aripiprazole and its four impurities from tablets were set as Analytical Target Profile. In line with it, the influence of Critical Method Parameters (acetonitrile content, concentration of perchloric acid in water phase, and column temperature) on predefined Critical Method Attributes (separation of the critical pair of peaks, retention of the first and last eluting peak) was investigated with aid of the Central Composite Design. Further on Design Space, where Critical Method Parameters meet predefined acceptance limits with a high level of probability (π ≥ 85%), was computed as a result of performed Monte Carlo simulations. A normal operating conditions corresponding to 34% of acetonitrile, 66% of 42.5 mM perchloric acid, and column temperature at 35°C were selected from created Design Space. Robustness testing of the quantitative performances of the developed method was conducted combining Plackett–Burman design with alias matrix approach. Through the additional validation testing, reliability of the developed method for the use in the routine practice was completely confirmed.     

A Quality by Design and green LC technique for the determination of mast cell stabilizer and histamine receptor antagonist (Olopatadine HCl) in multiple formulations

Mast cell stabilizer and histamine receptor antagonist olopatadine hydrochloride (OPT) assay method predicated on LC have been established for the analysis in multiple formulations. The current method dealt with ophthalmic solution, nasal spray, and tablet formulation products. The isocratic chromatography method was optimized and validated with a Boston green C8 column (150 × 4.6 mm, 5 μm i.d.). Sodium dihydrogen phosphate buffer (pH 3.5) with acetonitrile in the ratio of 75:25 (v/v) was used as a mobile phase at a flow rate of 1.0 mL min⁻¹ and at the column temperature of 30°C, and the detection was done at 299 nm. The method was validated as per International Council for Harmonisation (ICH) guidelines and United States Pharmacopoeia (USP). The accuracy results ranged from 99.9 to 100.7%, % relative standard deviation (RSD) from the precision was 0.5, and correlation coefficient from the linearity experiment was > 0.999. Solution stability was established for 24 h at room temperature and refrigerator conditions, and it was found that the solutions were stable. Using quality by design-based experiment designs, critical quality attributes were studied and it was found that the method was robust. In all the forced degradation studies peak purity was passed, and no interference was found at the retention time of the active component. The method validation data demonstrated that the developed method is linear, precise, accurate, specific, robust, and stable for the determination of OPT from multiple formulations. Analytical eco-scale tool, Green Analytical Procedure Index (GAPI) tool, and the National Environmental Method Index (NEMI) were used to evaluate the greenness of the method, and the analytical eco-score of 77 for the presented method was found to be excellent.     

Quality by Design Approach for Development of W/O Type Microemulsion-Based Transdermal Systems for Atenolol

The objective of the present investigation was to develop microemulsion-based transdermal systems of highly water soluble drug, Atenolol, by quality by design technique. Atenolol-loaded W/O microemulsions were optimized using D-optimal design with concentrations of oil, surfactants mixture, and water as independent variables, which was converted into microemulsion-based gel (MBG). The results of in vitro permeation of the optimized batch of Atenolol-loaded MBG revealed significant increase in permeability parameters as compared to its convention gel. All results suggested suitability of W/O type MEs as carriers for transdermal delivery of highly water soluble drug, Atenolol.     

Quality by design-guided development of a capillary electrophoresis method for the simultaneous chiral purity determination and impurity profiling of tamsulosin

Analytical Quality by Design principles using the design of experiments were applied for the development of a capillary electrophoresis method for the determination of enantiomeric purity and chemically related impurities of tamsulosin. From initial scouting experiments, a dual cyclodextrin (CD) system composed of sulfated β-CD and carboxymethyl-α-CD was selected as the chiral selector. A fractional factorial resolution V+ design was used for the identification of the critical process parameters, while a face-centered central composite design and Monte Carlo simulations were employed for final optimization and defining the design space of the method. The experimental conditions of the working point were: 30 mM sodium phosphate buffer, pH 3.0, containing 40 mg/mL sulfated β-CD and 7 mg/mL carboxymethyl-α-CD, capillary temperature 18°C, applied voltage -23 kV. Following the assessment of robustness by applying a Plackett-Burman design, the method was validated according to the International Council for Harmonization of Technical Requirements for Pharmaceuticals for Human Use guideline Q2(R1). The method allowed the quantification of the chiral impurity and three other related impurities at the 0.1 % level with acceptable accuracy and precision.     

Application of Central Composite Design for the Determination of Exfoliating Agents in Cosmetics by Capillary Electrophoresis with Electroosmotic Flow Modulation

Multivariate analysis within central composite design is applied to simplify an optimization procedure and explore the interactions among experimental parameters in analytical chemistry. In this study, central composite design was used to identify the optimal capillary electrophoresis conditions with electroosmotic flow modulation to determine seven exfoliating agents in cosmetics. The influence of phosphate concentration, cetyltrimethylammonium bromide concentration, and methanol percentage on the response was evaluated by the use of the chromatographic exponential function to simultaneously investigate the resolution and separation under sixteen sets of capillary electrophoresis conditions. The optimized conditions were 150 mM phosphate solution (pH = 7) containing 0.5 mM cetyltrimethylammonium bromide, 3 mM γ-cyclodextrin, and 25% methanol as the running buffer. To shorten the analysis time, an electroosmotic flow modulating agent (cetyltrimethylammonium bromide) was added to the separation buffer. Calibration plots were linear (r ≧ 0.998) with high precision and accuracy in the homemade cosmetic matrix. The exfoliating agents in two commercial cosmetic products were determined using the optimized conditions, and the results correlated well with results obtained by high-performance liquid chromatography–mass spectrometry.     

An Experimental Design Approach to Quantitative Expression for Quality Control of a Multicomponent Antidiabetic Formulation by the HILIC Method

A rapid and reproducible hydrophilic liquid chromatography (HILIC) process was established for concomitant determination of remogliflozin etabonate (RE), vildagliptin (VD), and metformin (MF) in a formulation. A face-centered central composite experimental design was employed to optimize and predict the chromatographic condition by statistically studying the surface response model and design space with desirability close to one. A HILIC column with a simple mobile phase of acetonitrile (65% v/v) and 20 mM phosphate buffer (35% v/v, pH 6, controlled with orthophosphoric acid) was used to separate RE, VD, and MF. RE, VD, and MF were separated in 3.6 min using an isocratic mode mobile phase flow at a flow rate of 1.4 mL at room temperature, and the analytes were examined by recording the absorption at 210 nm. The developed HILIC method was thoroughly validated for all parameters recommended by ICH, and linearity was observed in the ranges 20–150 µg/mL, 10–75 µg/mL, and 50–750 µg/mL for RE, VD, and MF, respectively, along with excellent regression coefficients (r² > 0.999). The calculated percentage relative deviation and relative error ascertained the precision and accuracy of the method. The selectivity and accuracy were further confirmed by the high percentage recovery of added standard drugs to the formulation using the standard addition technique. The robustness of the HILIC processes was confirmed by developing a half-normal probability plot and Pareto chart, as the slight variation of a single factor had no significant influence on the assay outcomes. Utilization of the optimized HILIC procedure for concurrent quantification of RE, VD, and MF in solid dosage forms showed accurate and reproducible results. Hence, the fast HILIC method can be regularly employed for the quality assurance of pharmaceutical preparations comprising RE, VD, and MF.     

Integrated quality by design (QbD) and design of experiments (DoE) approach for UFLC determination of telaprevir in rat serum

An ultrafast liquid chromatographic bioanalytical method was developed and validated for the determination of telaprevir in Wistar albino rat serum. Principles of quality by design (QbD) were implemented for enhancing the bioanalytical liquid–liquid extraction of telaprevir from rat serum. A Box–Behnken design was utilized in the studies by selecting extraction time, centrifugation speed, and vortex time as the critical method variables for evaluating their effect on the critical analytical attribute, i.e., %recovery of telaprevir. Chromatographic separation was achieved within a run time of 10?min using a C-18 column and mobile phase comprising of methanol:borate buffer of pH 9 (90:10 v/v) flowing at 1.2?mL/min. Photodiode array detection was performed at 270?nm. Results of validation studies were satisfactory. The method was linear over a concentration of 25–10,000?ng/mL. Limit of detection for the developed method was 10?ng/mL. Further, design of experiments (DoE) used for inter-day accuracy and precision study suggested superior method reliability. This integrated QbD- and DoE-based approach ensured the development of a validated and reliable analytical method for optimum bioanalysis of telaprevir in biological matrix.     

Application of Cu-based metal-organic framework (Cu-BDC) as a sorbent for dispersive solid-phase extraction of gallic acid from orange juice samples using HPLC-UV method

In the current work, the Cu-based metal-organic framework (MOF) (Cu-BDC) was synthesized through a simple and one-pot solvothermal method and used as a sorbent for dispersive solid-phase extraction (DSPE) of gallic acid (GA) from orange juice samples followed by high-performance liquid chromatography-ultraviolet (HPLC-UV) determination. The prepared Cu-BDC was fully characterized using different analysis including Fourier-transform infrared spectroscopy (FT-IR), X-ray Powder Diffraction (XRD), Scanning Electron Microscopy (SEM), Energy-dispersive X-ray spectroscopy (EDS) and Brunauer-Emmett-Teller (BET). The prepared Cu-BDC had exceptional properties in terms of high specific surface area (400 m² g⁻¹) and high adsorption capacity (300 mg g⁻¹) toward GA. The influence of effective factors on extraction recovery of GA such as sorbent dosage, ultrasonic time, washing solvent volume and eluent solvent volume was evaluated using a central composite design (CCD) based response surface methodology (RSM). Under the optimized conditions, the presented method showed a wide linear range of 0.5–15.0 mg L⁻¹ with a detection limit of (0.145 mg L⁻¹) and acceptable repeatability (RSD < 6.0%) which shows the applicability of the proposed method for accurate determination of GA from complicated orange juice samples.     

Quality by Design enabled enhanced bioanalytical extraction and UFLC determination of vilazodone from rat serum

An ultrafast liquid chromatographic bioanalytical method was developed and validated for the determination of vilazodone in Wistar rat serum. Principles of quality by design were implemented for enhancing the bioanalytical liquid–liquid extraction of vilazodone from rat serum. A Box–Behnken design was utilized in the studies by selecting extraction time, centrifugation speed, and vortex time as the critical method variables for evaluating their effect on the analytical attribute, i.e., %recovery of vilazodone. Chromatographic separation was achieved within a run time of 10?min using a C-18 column and mobile phase comprising of methanol:phosphate buffer of pH 7 (85:15 v/v) flowing at 1.5?mL/min. Photodiode array detection was performed at 242?nm. Results of validation studies were satisfactory. The method was linear over a concentration of 100–2,000?ng/mL with acceptable accuracy and precision. Limits of detection and quantitation for the developed method were 50 and 100?ng/mL, respectively. This QbD-based approach was found suitable for routine bioanalysis of vilazodone in the biological matrix.     

Development of a validated liquid chromatographic method for quantification of sorafenib tosylate in the presence of stress‐induced degradation products and in biological matrix employing analytical quality by design approach

The current research work envisages an analytical quality by design‐enabled development of a simple, rapid, sensitive, specific, robust and cost‐effective stability‐indicating reversed‐phase high‐performance liquid chromatographic method for determining stress‐induced forced‐degradation products of sorafenib tosylate (SFN). An Ishikawa fishbone diagram was constructed to embark upon analytical target profile and critical analytical attributes, i.e. peak area, theoretical plates, retention time and peak tailing. Factor screening using Taguchi orthogonal arrays and quality risk assessment studies carried out using failure mode effect analysis aided the selection of critical method parameters, i.e. mobile phase ratio and flow rate potentially affecting the chosen critical analytical attributes. Systematic optimization using response surface methodology of the chosen critical method parameters was carried out employing a two‐factor–three‐level–13‐run, face‐centered cubic design. A method operable design region was earmarked providing optimum method performance using numerical and graphical optimization. The optimum method employed a mobile phase composition consisting of acetonitrile and water (containing orthophosphoric acid, pH 4.1) at 65:35 v/v at a flow rate of 0.8 mL/min with UV detection at 265 nm using a C₁₈ column. Response surface methodology validation studies confirmed good efficiency and sensitivity of the developed method for analysis of SFN in mobile phase as well as in human plasma matrix. The forced degradation studies were conducted under different recommended stress conditions as per ICH Q1A (R2). Mass spectroscopy studies showed that SFN degrades in strongly acidic, alkaline and oxidative hydrolytic conditions at elevated temperature, while the drug was per se found to be photostable. Oxidative hydrolysis using 30% H₂O₂ showed maximum degradation with products at retention times of 3.35, 3.65, 4.20 and 5.67 min. The absence of any significant change in the retention time of SFN and degradation products, formed under different stress conditions, ratified selectivity and specificity of the systematically developed method.     

Analytical quality by design in the development of a cyclodextrin‐modified capillary electrophoresis method for the assay of metformin and its related substances

Quality by Design (QbD) is a new paradigm of quality to be applied to pharmaceutical products and processes, recently encouraged by International Conference on Harmonisation guidelines. In this paper QbD approach was applied to the development of a CE method for the simultaneous assay of metformin hydrochloride (MET) and its main impurities. QbD strategy was focused on electrophoretic process understanding, and the analytical method was thoroughly evaluated by applying risk assessment and chemometric tools. Method scouting allowed CD‐CZE based on the addition of carboxymethyl‐β‐CD to Britton‐Robinson acidic buffer to be chosen as operative mode. Seven critical process parameters (CPPs) were selected, related to capillary, injection, BGE and instrumental settings. The effect of the different levels of the CPPs on critical quality attributes (CQAs), e.g. critical resolution values and analysis time, was evaluated in a screening study. Response surface methodology led to draw contour plots and sweet spot plots. The definition of design space was accomplished by applying Monte‐Carlo simulations, thus identifying by risk of failure maps a multivariate zone where the CQAs fulfilled the requirements with a selected probability. Finally, a control strategy was designed and the method was applied to a real sample of MET tablets.     

Validation of a high-performance liquid chromatographic method for determination of isoflavonoids in soybeans. Study of the extraction procedure by experimental design

Summary An improved analytical method, reversed-phase high-performance liquid chromatography on a narrow-bore C₁₈ column, has been developed for the simultaneous determination of genistein, daidzein, formononetin, and biochanin A. The method was validated in terms of detection limits, quantitation limits (LOQ), linearity, and precision.LOQ in the 0.04–0.1 μg mL⁻¹ range were calculated, enabling determination of these compounds of nutritional concern at trace levels. Good linearity was demonstrated over three orders of magnitude of concentration for each analyte (r ²=0.998–1.000). The intra-day repeatability was evaluated in terms ofRSD (%) at two concentration levels for each analyte (RSD (%) <1.8%). Good inter-day reproducibility of data was proved by performing homoscedasticity and ANOVA tests (P>0.05 at the 95% confidence level). The method was applied to the determination of genistein and daidzein in yellow soybeans, after optimization of the method for extraction of isoflavonoid aglycones from soybeans by experimental design, i.e. central composite design. Extraction recoveries up to 87±4% were obtained when the corresponding glycosidic forms (genistin and daidzin) were added to soybean samples.     

Application of a multivariate approach for analyte focusing by micelle collapse‐micellar electrokinetic chromatography for analyzing sunscreen agents in cosmetics

The operating parameters that affect the performance of the online preconcentration technique “analyte focusing by micelle collapse‐MEKC (AFMC‐MEKC)” were examined using a multivariate approach involving experimental design to determine the sunscreen agents in cosmetics. Compared to the single‐variable approach, the advantage of the multivariate approach was that many factors could be investigated simultaneously to obtain the best separation condition. A fractional factorial design was used to identify the fewest significant factors in the central composite design (cCD). The cCD was adopted for evaluating the location of the minimum or maximum response in this study. The influences of the experimental variables on the response were investigated by applying a chromatographic exponential function. The optimized condition and the relationship between the experimental variables were acquired using the JMP software. The ANOVA analysis indicated that the Tris pH value, SDS concentration, and ethanol percentage influenced the separation quality and significantly contributed to the model. The optimized condition of the running buffer was 10 mM Tris buffer (pH 9.5) containing 60 mM SDS, 7 mM γ‐CD, and 20% v/v ethanol. The sample was prepared in 100 mM Tris buffer (pH 9.0) containing 7.5 mM SDS and 20% v/v ethanol. The SDS concentration in the sample matrix was slightly greater than the CMC value that makes the micelle be easily collapsed and the analytes be accumulated in the capillary. In addition, sunscreen agents in cosmetics after 1000‐fold dilution were successfully determined by AFMC‐MEKC.     

Optimization of a novel method based on solidification of floating organic droplet by high-performance liquid chromatography for evaluation of antifungal drugs in biological samples

In this study, a simple, rapid, and highly efficient liquid-phase microextraction method based on solidification of floating organic droplet was coupled with high performance liquid chromatography-photo diode array detection (HPLC-PDA) for determination of ketoconazole, clotrimazole, and miconazole as antifungal drugs. Central composite design (CCD) was used for optimization of several factors affecting the extraction efficiency. The optimized conditions were established to be 550 rpm for stirring rate, 35 min for extraction time, 57 °C for extraction temperature, 8.5 for solution pH, 10 μl for organic solvent volume, and 7% (w/v) of NaCl for ionic strength. Limit of detections (LODs) of the extraction method ranged from 0.01 to 0.1 μg L⁻¹ and the linear dynamic ranges (LDRs) ranged from 0.1 to 300 μg L⁻¹ for the three antifungal drugs. Relative standard deviations (RSDs) of the proposed method were 5-11%. Preconcentration factors in the range of 306-1350 were obtained at extraction time of 35 min. Finally, performance of the proposed method was evaluated for the extraction and determination of the drugs’ levels in microgram per liter in samples and satisfactory results were obtained.