Simultaneous Analysis of Irbesartan and Hydrochlorothiazide: An Improved HPLC Method with the Aid of a Chemometric Protocol

Experimental design method was used for HPLC determination of irbesartan and hydrochlorothiazide in combined dosage forms. The traditional approach for optimization of experiments is time-consuming, involves a large number of runs and does not allow establishing the multiple interacting parameters. The main advantages of the experimental design method include the simultaneous screening of a larger number of factors affecting response and the estimation of possible interactions. On the basis of preliminary experiments, three factors-independent variables were selected as inputs (methanol content, pH of the mobile phase and temperature) and as dependent variables, five responses (resolution, symmetry of irbesartan peak, symmetry of hydrochlorothiazide peak, retention factor of irbesartan and retention factor of hydrochlorothiazide) were chosen. A full 23 factorial design, where factors were examined at two different levels (“low” and “high”) was used to determine which factors had an effect on the studied response. Afterwards, experimental design was used to optimize these influent parameters in the previously selected experimental domain. The novelty of our method lies in the optimization step accomplished by Derringer′s desirability function. After optimizing the experimental conditions a separation was conducted on a Supelcosil C18 (150 mm × 4.6 mm, 5 μm particle size) column with a mobile phase consisting of methanol-tetrahydrofuran-acetate buffer 47:10:43 v/v/v, pH 6.5 and a column temperature of 25 °C. The developed method was successfully applied to the simultaneous separation of these drug-active compounds in their commercial pharmaceutical dosage forms.

Abstract: Experimental design method was used for HPLC determination of irbesartan and hydrochlorothiazide in combined dosage forms. The traditional approach for optimization of experiments is time-consuming, involves a large number of runs and does not allow establishing the multiple interacting parameters. The main advantages of the experimental design method include the simultaneous screening of a larger number of factors affecting response and the estimation of possible interactions. On the basis of preliminary experiments, three factors-independent variables were selected as inputs (methanol content, pH of the mobile phase and temperature) and as dependent variables, five responses (resolution, symmetry of irbesartan peak, symmetry of hydrochlorothiazide peak, retention factor of irbesartan and retention factor of hydrochlorothiazide) were chosen. A full 2 3 factorial design, where factors were examined at two different levels ("low" and "high") was used to determine which factors had an effect on the studied response. Afterwards, experimental design was used to optimize these influent parameters in the previously selected experimental domain. The novelty of our method lies in the optimization step accomplished by Derringers desirability function. After optimizing the OPEN ACCESS experimental conditions a separation was conducted on a Supelcosil C 18 (150 mm × 4.6 mm, 5 m particle size) column with a mobile phase consisting of methanol-tetrahydrofuranacetate buffer 47:10:43 v/v/v, pH 6.5 and a column temperature of 25 °C. The developed method was successfully applied to the simultaneous separation of these drug-active compounds in their commercial pharmaceutical dosage forms.
Some published papers have dealt with the optimization procedure for determination either of IRB (and other angiotensin-II-receptor antagonist) 15-18 or HCT 19,20, but to the best of our knowledge neither of the already published methods has been related to HPLC method optimization for simultaneous determination of IRB and HCT. The aim of this work was the evaluation of the chromatographic behavior of IRB and HCT using an appropriate experimental design (ED). The main advantage of such approach is simultaneous optimization of influencing factors and response variables which enables prediction of chromatographic retention and postulation of optimum conditions for separation. Several statistical approaches could be used for estimating the retention behavior as a function of chromatographic conditions. In this work, full factorial design (FFD) and response surface methodology (RSM) have been applied. Factorial design enables an estimation of investigated factors which have the most importance. Response surface methodology is generally employed in order to provide a description of the response pattern in the region of the studied observations and to assist in finding the region where the optimal response occurs 21.
Marketed tablet formulation (Co-Irda tablets, Nobel Ilac, Turkey) containing 150 mg of irbesartan and 12.5 mg of hydrochlorthiazide have been analyzed. The significant feature of this combination lies in the fact that hydrochlorthiazide is present in minute amounts compared to irbesartan, which makes an analysis more complicated and tedious. In the present paper, a fast, simple and accurate HPLC method has been proposed without the tedious extraction procedure.

Results and Discussion
When applying experimental design methodologies, it is advisable to keep the number of variables as low as possible in order to avoid very complex response models and large variability. When the number of influencing factors is up to four, full factorial design (FFD) is recommended.
The investigation was carried out in several steps. The objective of the first step in the investigation was to perform a screening of the factors that could potentially influence chromatographic retention, thus the independent variables were defined during the preliminary study. Some chromatographic parameters, such as flow rate were excluded as its influence can usually be predicted by common chromatographic theory knowledge. The factors generally selected to optimize the chromatographic separation of ionisable compounds are pH and the content of organic solvent of the mobile phase. The variations of these parameters induce a variation of the degree of ionization and thus affect chromatographic behavior. In addition, column temperature affects retention behavior, thus three factors-independent variables were selected as inputs: methanol content, pH of the mobile phase and temperature. Since a good separation is characterized by good resolution and since run time is very important (from a practical point of view), five responses were chosen as dependent variables: resolution (k R ), symmetry of the irbesartan peak (Sym IRB ), symmetry of the hydrochlorothiazide peak (Sym HCT ), retention factor of irbesartan (Rt IRB ) and retention factor of hydrochlorothiazide (Rt HCT ).
The presence of several functional groups in the molecular structures, such as biphenyl, imidazole and benzene (Figure 1), makes a RP-HPLC method with PDA detection suitable for the determination. Since the RP-HPLC method is based on using a polar mobile phase, a complete description of the ionization profile of the examined substances has been used for the evaluation of retention behavior and also for the separation. The degree of ionization of the drug strongly affects solubility and retention. Additionally, the knowledge of dissociation constant of ionisable compounds at different pH values and the solvent composition is also significant to determine the optimal separation conditions in reversed phase liquid chromatography (RP-LC). Considering the chemical structures, it is possible to establish a number of proton acceptor and donor groups (Figure 2), to assume ionized structures of IRB and HCT (Figures 3 and 4) and degrees of ionization depending on the pH (Tables 1 and 2).     As it can be seen, IRB has two pKa values (4.12, imidazole nitrogen; 7.40, tetrazole nitrogen) and at pH 4 it is present as a mixture of molecular (43.34%) and monoprotonated (56.64%) forms. The structure of HCT has three pKa values (9.09, sulphonamide group; 9.83 cyclic sulphonamide; 11.31, secondary amine). At higher pH values, HCT is present as a mixture of different deprotonated species. Having in mind the acid-base properties of IRB and HCT, the pH interval from 4.0 to 6.5 was chosen for further investigation. In this pH interval, IRB is partially ionized and HCT is completely unionized and under these conditions the following order of retention could be expected: HCT, then IRB. The retention of a substance is a function of the volume fraction of the organic modifier in the mobile phase. Taking into account the variation of the retention factors of compounds with polarity of the mobile phase, a range of methanol concentrations from 35% to 55% was selected for investigation.
The analysis of IRB and HCT was started on a non-polar stationary phase (Supelcosil C 18 column, 150 mm × 4.6 mm, 5 m particle size) with the mobile phase consisting of methanol-tetrahydrofuranacetate buffer (pH of mobile phase was adjusted to 4 with acetic acid). The column temperature was set at 35 °C and the flow rate at 0.75 mL/min. Acceptable separation was achieved with methanoltetrahydrofuran-acetate buffer mixtures ranging from 35:10:55 to 55:10:35 (v/v/v), but peak shape and run time needed to be improved. In order to evaluate the effect of the most important factor, a 2 3 full factorial design (FFD) with three replicates at the zero level was chosen. The experimental data was coded in order to follow the significance of factors in an easier way. Factors and their "low" (−1), "high" (+1) and "zero" (0) values are presented in Table 3. Table 3. Factors and levels.

Factors
Factor levels A two-level factorial design runs and three replicates of the central point needed 11 runs to complete a whole factorial design. As dependent variables, five responses were chosen: Resolution (k R ), symmetry of irbesartan peak (Sym IRB ), symmetry of hydrochlorothiazide peak (Sym HCT ), retention factor of irbesartan (Rt IRB ) and retention factor of hydrochlorothiazide (Rt HCT ). The matrix of experiments and results obtained as an average value of three runs are presented in Table 4. Appropriate calculations were done with the Design-Expert 7.0 software (Stat-Ease Inc. Minneapolis, MN, USA). A second-order interaction model was suggested as a model of relationship between input and output and is presented by Equation (1): 13 , b 23 ) and b ijk represent the regression coefficient and A,B,C represent independent variables. The calculated coefficient and model of polynomial regression is presented in Table 5. The repetition of the central experimental point provided a precise estimation of the experimental errors and the measure of the adequacy of the models (lack of fit). The results were analysed by ANOVA method and the results are presented in Table 6.
The lack of fit test was determined by performing Fischer-F test. The high value of F with a very low probability (only model terms with corresponding p-value lesser than 0.05 are significant at 95% confindence level) implies that there was no evidence of the models lack-of-fit and the models could be accepted as an adequate representation of the data. In addition, the values of R 2 and R 2 adjusted taking into account the degrees of freedom indicated that the regression model fits the data well. The exception is a model which considers the retention time of HCT. Since p is greater than 0.05, the model needs reduction in order to improve the relationship between parameters. The data collected from the performed FFD design led to the following conclusions: it was noticed that methanol content in the mobile phase and pH have the largest influence on k R , Rt IRB and less influence on the other responses. This influence had a minus sign, which means that the higher pH values and percentage of methanol in the mobile phase will reduce resolution and retention time of IRB. At the same time, the temperature of the column has a positive effect on Sym IRB and a negative effect on Sym HCT . Since the selected responses were not affected in the same manner an additional optimization procedure was needed.
In order to get the best chromatographic performance, the multicriteria methodology was employed by means of Derringers desirability function 22. It is based on constructing desirable ranges for each response (individual desirable function, d i ) and establishing an overall desirability function (the Derringer desirability function). The Derringers desirability function is defined as the geometric mean of individual desirability functions and can be expressed by Equation (2): where n is the number of responses, and p n is the weight of the responses. Weight of the response is the relative importance of each individual functions d i and may range from 0.1 to 10. With a weight of 1, d i varies in a linear way. In this study, weights equal to 1 was selected. Individual desirability functions range from 0 (undesired response) to 1 (a fully desired response). A value of D close to 1 means that the combination of different criteria is globally optimal. If any of the responses or factors falls outside their desirability range, the overall function becomes zero. There are several ways for calculating the desirability function depending on the goal desired. The meanings of the goal parameters are: The goals of multicriteria optimization for each response in this paper are shown in Table 7. Table 7. Criteria for multivariate optimization of the individual responses.

Goal
Lower Desirability function calculations were performed using Design-Expert  7.0. Obtained results are graphically presented ( Figure 5). For better visualization of the results, the global desirability function D was presented in a form of a three-dimensional plot and presented in Figure 6.  After setting the optimal conditions, the proposed method was validated. Interfering peaks were not detected at the retention time of IRB and HCT, indicating the good selectivity of the method. Linear dependence of the peak areas versus concentration was determined for the proposed ranges.
(a) (c) (b) Parameters of the linear regression equations were calculated and are presented in Table 8. The statistical significance of the intercept was tested using Student's t-test. The limit of detection (LOD) and limit of quantification (LOQ) were calculated as LOD = 3/S and LOQ = 10/S, where  is the standard deviation of the response and S the intercept determined from the corresponding calibration curve. Precision of the procedure was assessed by analyzing nine solutions containing known quantities of the investigated compounds. Law values of relative standard deviation for repeatability, RSD < 2.5%, and high recovery (Table 9) indicate very good precision of the proposed method. The applicability of the proposed method was examined by analysing commercially available Co-Irda tablets.

Drugs and Reagents
The irbesartan and hydrochlorothiazide standards and Co-Irda tablets (Nobel Ilac, Istanbul, Turkey) consisting of 150 mg IRB and HCT were obtained from Zada Pharmaceuticals d.o.o. (Tuzla, Bosnia and Herzegovina).
All solvents: methanol, tetrahydrofuran and sodium-acetate (purchased by Chromosol, Sigma-Aldrich, Munich, Germany), and acetic acid (Fluka, Eindhoven, The Netherlands) were of a grade suitable for high-performance liquid chromatography analysis. The HPLC analyses were done by using a Thermo Finnigan Surveyor chromatographic system equipped with a PDA detector and sample injections were made through an injector valve with a 5 L sample loop. Separations were performed on a Supelcosil C 18 column (150 mm × 4.6 mm, 5 m particle size) with detection at 271 nm. Mobile phases were prepared according to the plan of the experiments given in Table 4. The resulting mobile phases were degassed and vacuum filtered through a 0.45 m membranes filter (Alltech Associates, Lokeren, Belgium). The flow rate was 0.75 mL/min.

Software
Experimental design, statistical analysis and desirabiity function calculation were performed by using MarvinSketch 5.8.2 (Chem Axon Ltd., Somerville, MA, USA, and Budapest, Hungary) and Design-Expert  7.0 (Stat-Ease Inc.).

Solutions
Stock solutions were prepared by dissolving standard substances in methanol to obtain concentrations of 0.8 mg/mL for IRB and 0.2 mg/mL for IRB.

Solutions for Method Optimization
Stock solutions were diluted with methanol to obtain a concentration of 0.08 mg/mL of IRB and 0.02 mg/mL of HCT.

Standard Solutions for Linearity Testing
For the calibration curves, a series of eight solutions were prepared from stock solution in the concentration range from 0.08 to 0.4 mg/mL for IRB and from 0.02 to 0.1 mg/mL for HCT.

Solutions for Accuracy Testing
The laboratory mixture containing placebo, IRB and HCT was prepared in the ratio related to the investigated tablets. For the quantiative analysis, three solutions corresponding to 80%, 100% and 120% to those in tablets were prepared.

Solutions for Estimating Precision
In order to estimate precision, three series (0.16, 0.24 and 0.32 mg/mL for IRB; 0.04, 0.06 and 0.08 mg/mL for HCT) were prepare with ten solutions for each of the concentrations.

Sample Solutions
A tablet mass which corresponds to 150 mg IRB and 12.5 mg HCT was dissolved in 100 mL volumetric flask with methanol, placed into a ultrasonic bath and filtrated. 2 mL of filtrate were diluted with methanol to 10 mL.

Conclusions
Experimental design methodology was used for simultaneous HPLC determination of irbesartan and hydrochlorothiazide in combined dosage forms. The significant feature of these combinations lies in the fact that hydrochlorothiazide is present in minute amounts compared to irbesartan which makes for a more complicated and tedious analysis.
The chemometric approach for optimization of chromatographic separation of irbesartan and hydrochlorothiazide has been demonstrated. The chemometric methodology chosen for the particular objectives was very successful in the retention behavior exploration. Since there was a mix of linear responses with different targets, Derringers desirability function was applied. After defining a global desirability according to the accepted constraints, optimal chromatographic conditions were established.
The proposed HPLC method was validated according to ICH guidelines. From the study of validation parameters, it was observed that the method is specific, accurate, precise, reproducible and is not time-consuming (run time is less than six minutes). Since there was no interference from other components present in the dosage forms, complicated procedures for extraction were not required. The results obtained in this study corroborate that the proposed HPLC method can be used for routine quantitative analyses of the investigated compounds in a mixture or for their individual determination in pharmaceutical dosage forms.