Abstract
Central composite experimental design was used for fast, simple, and accurate high-performance liquid chromatography (HPLC) determination of hydrochlorothiazide, amlodipine and valsartan in combined dosage forms. This method avoids the disadvantages of the traditional analytical approach, which is time-consuming, involves a large number of runs, and does not allow establishing the multiple interacting parameters. On the basis of preliminary experiments and physicochemical characteristic of analyzed substances, three independent variables (methanol content, pH of the mobile phase, and column temperature) were selected as input, while as dependent variables, six responses (retention time of hydrochlorothiazide, retention time of amlodipine, retention time of valsartan, asymmetry of hydrochlorothiazide peak, asymmetry of amlodipine peak, and asymmetry of valsartan peak) were chosen. Face centered central composite design enables an estimation of factors which have the most importance. After optimizing experimental conditions, a separation was conducted on a Zorbax C8 (150 mm×4.6 mm, 5 μm) column with a mobile phase consisting of methanol-acetonitrile-acetate buffer 40:20:40 (v/v/v), pH adjusted to 3.5 with acetic acid, flow rate of 1 mL/min and column temperature of 40 degC. The method was successfully applied to the simultaneous separations of these active drug compounds in their commercial dosage forms.
Introduction
Recent guidelines for the treatment of hypertension have focused on the need for multiple medications to get most patients to goal blood pressure (BP). Two to three different classes of antihypertensive agents are frequently required, increasing the risk of poor compliance with therapy. Hence, the guidelines have recommended starting with combination therapy in patients with BP that is over 20 mm Hg systolic or 10 mm Hg diastolic above goal. The latest advance in treatment regimen has been the development of triple-therapy combinations of an angiotensin receptor blocker, amlodipine and hydrochlorothiazide [1, 2]. Hydrochlorothiazide [HCT, 6-chloro-1,1-dioxo-3,4-dihydro-2H-1,2,4-benzothiadiazine-7-sulfonamide, Fig. 1 (a)] is a diuretic drug of the thiazide class, frequently used for treatment of hypertension and congestive heart failure. It reduces blood volume by acting on the kidneys to reduce sodium resorption in the distal convulated tubule, resulting in a decrease of blood pressure. Amlodipine [AML, (RS)-3-ethyl 5-methyl 2-[(2-aminoethoxy)methyl]-4-(2-chlorophenyl)-6-methyl-1,4-dihydropyridine-3,5-dicarboxylate, Fig. 1 (b)] is a long-acting calcium channel blocker which belongs to dihydropyridine (DHP) class, used as an antihypertensive and in the treatment of angina pectoris. Amlodipine acts by relaxing the smooth muscle in the arterial wall, decreasing total peripheral resistance and hence reducing blood pressure. Valsartan [VAL, (S)-3-methyl-2-(N-{[2’-(2H-1,2,3,4-tetrazol-5-yl)biphenyl-4-yl]methyl}pentanamido) butanoic acid, Fig. 1 (c)], is an angiotensin II receptor (AT1) antagonist. Angiotensin II receptor antagonists represent a relatively new pharmacological class [3] which acts mainly by selective blockade of AT1 receptors and reduces the effects of angiotensin II. By blocking the action of angiotensin, valsartan dilates blood vessels and reduces blood pressure.

Structural formulae of HCT (a), AML (b) and VAL (c).
There are many manuscripts focusing on the determination of HCT, AML and VAL only or in combination with other drugs [4–21]. Several published papers describe the simultaneous determination of binary combinations HCT, AML and VAL [22–25]. Galande et al. [26] describe estimation of all three drugs in combination by UV Spectrophotometry. The literature revealed only one paper dealing with the determination of HCT, AML and VAL by HPLC in human plasma [27]. Neither of the already published methods reports HPLC method for simultaneously determination all three components in a solid dosage form.
Central composite design (CCD) and response surface methodology (RSM) have been applied to obtain precise, simple and rapid HPLC method. Face centered CCD enables an estimation of factors which have the most importance. RSM 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 [28–33]. CCD requires a minimum number of experiments to be performed and has been widely used for fitting a second-order model. The independent variables including methanol content, pH of mobile phase and column temperature were defined as factors and retention times and asymmetries peaks of HCT, AML and VAL (tR HCT , tR AML , tR VAL , As HCT , As AML and As VAL ) were used separately as the responses in the CCD.
A market tablet formulation (ExforgeHCT, Novartis GmbH, Wien, Austria) containing 25 mg HCT, 10 mg AML and 160 mg VAL was analyzed. The significant feature of this combination lies in the fact that hydrochlorthiazide and amlodipine are present in minute amounts compared to valsartan, 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.
Chemicals and reagents
The hydrochlorothiazide, amlodipine and valsartan standards (HCT = 101.12%, AML = 100.97%, VAL = 100.13%), ExforgeHCT tablets (Noavartis GmbH, Wien, Austria), consisting of 25 mg HCT, 10 mg AML and 160 mg VAL, were used. All solvents: methanol, acetonitrile and sodium-acetate (purchased by Chromosol, Sigma-Aldrich, Munich, Germany), and acetic acid (Fluka, Eindhoven, The Nederlands) were of a grade suitable for high-performance liquid chromatography analysis.
Preparation of standard solutions
Stock solutions were prepared by dissolving standard substances in mobile phase to obtain concentrations of 0.25 mg/mL for HCT, 0.1 mg/mL for AML and 1.6 mg/mL for VAL. Solutions for method optimization were prepared by dissolving stock solutions in mobile phase in order to obtain a concentration of 0.025 mg/mL for HCT, 0.01 mg/mL for AML and 0.16 mg/mL for VAL. Standard solutions for linearity testing for the calibration curves were a series of ten solutions prepared in the concentration range from 0.0125 to 0.125 mg/mL for HCT, 0.005 to 0.05 mg/mL for AML and 0.08 to 0.8 mg/mL for VAL. Solutions for accuracy testing containing placebo were prepared as the laboratory mixture of solutions of HCT, AML and VAL in the ratio related to the investigated tablets. For quantitative analysis, three solutions corresponding to 80, 100 and 120% of the label amount were prepared. To estimate precision, three series (0.025, 0.05, and 0.075 mg/mL for HCT; 0.01, 0.02, and 0.03 mg/mL for AML, 0.08, 0.16, and 0.24 mg/mL for VAL) were prepared with ten solutions for each of the concentrations.
Sample solutions
A tablet mass corresponds to 25 mg HCT, 10 mg AML and 160 mg VAL (ten tablets were accurately weighed and ground to fine powder) was dissolved in a 100-mL volumetric flask with mobile phase, placed into an ultrasonic bath, and filtrated. Filtrate (5 mL) was diluted with mobile phase to 50 mL.
Apparatus and software
The HPLC analyses were done by using a Shimadzu chromatographic system, consisting of a Shimadzu HPLC system (Shimadzu Technologies, Kyoto, Japan), a pump (LC-20AD VP), UV-Visible detector (SPD-20AV VP), column oven (CTO-20AC VP), column Zorbax C8 (150×4.6 mm, 5 μm), degasser system (DGU-20AS) and injector with a 10 μL sample loop. Experimental design, structural analysis, statistical analysis and calculation of desirability function were performed by using Marvin Sketch 5.8.2 (Chem Axon Ltd., Somerville, MA, USA, and Budapest, Hungary) and Design-Expert 8.0.7.1 (Stat-Ease Inc.).
Procedure
Separations were performed on a Zorbax C8 column (150 mm×4.6 mm, 5 μm particle size) with detection at 254 nm. Mobile phases were prepared according to the experimental plan given in Table 1. The resulting mobile phases were degassed and vacuum filtered through a 0.45-μm membranes filter (Alltech Associates, Lokeren, Belgium).
CCD design and results of experiments
CCD design and results of experiments
Structural analysis
The presence of several functional groups in the molecular structures, such as biphenyl, imidazole, dihydropyridine, and benzene (Fig. 1), makes a reversed-phased (RP)-HPLC method with SPDA detection suitable for the determination. Since the RP-HPLC method is based on use of a polar mobile phase, a complete description of the ionization profile of the examined substances was 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).
As it can be seen, HCT has three pKa values (9.09, sulfonamide group; 9.83 cyclic sulfonamide; 11.31, secondary amine), AML has two pKa values of nitrogen based groups (4.16, dihydropyrimidine nitrogen; 9.45, amine nitrogen) and VAL has two pKa values (4.37, carboxyl group; 7.40, tetrazole nitrogen) (Fig. 2). Considering the chemical structures and pKa values, it is possible to assume ionized structures of HCT, AML and VAL. Having in mind the acid-base properties of HCT, AML and VAL, the pH interval from 2.5 to 4.5 was chosen for further investigation. At the pH values from 2.5 to 4.5 molecules are present in its ionizing or non-ionizing various forms. Major dominant species of HCT and VAL at pH 2.5 and 3.5 is non-ionized structure form, while two nitrogens in AML are ionized (Fig. 3).

pKa values of functional groups of HCT, AML and VAL.

Major dominant molecular species at pH 2.5 and pH 3.5.
At maximum pH value 4.5 HCT has still been in non ionizing form, AML has only protonated amine nitrogen, and VAL carboxyl group is ionized (Fig. 4).

Major dominant molecular species at pH 4.5.
To predict the order of peak appearance on chromatogram, logP was analyzed for HCT (logP = –0.58), AML (logP = 1.64) and VAL (logP = 5.27). Several papers report the effect of logP on the retention time of substance in RP-LC. Papers report that an increase of logP leads to an increase of logtR [34, 35].
The objective of the optimization 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 column temperature. CCD has been used for fitting a second-order model. Since a good separation is characterized by good resolution, run time is very important (from a practical point of view) and since the shape of the peak is determined by its asymmetry, eight responses were chosen as dependent variables: retention time of hydrochlorothiazide (tRHCT), retention time of amlodipine (tRAML), retention time of valsartan (tRVAL), asymmetry of the hydrochlorothiazide peak (As HCT ), asymmetry of the amlodipine peak (As AML ), asymmetry of the valsartan peak (As VAL ) and resolution (Rs AML and Rs VAL ). Retention parameters were measured using Class VP software with the USP computing option selected.
On the basis of structures, physicochemical characteristics of these three molecules and composition of mobile and stationary phase it is possible to assume the following order of retention: HCT, AML and VAL. 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 30% to 50% was selected for investigation.
The analysis of HCT, AML and VAL was started on a non-polar stationary phase (Zorbax C8 column, 150 mm×4.6 mm, 5 μm particle size) with the mobile phase consisting of methanol-acetonitrile- acetate buffer (pH of mobile phase was adjusted to 4.0 with acetic acid). The column temperature was set at 25 °C and the flow rate at 1 mL/min. Separation was acceptable with methanol-acetonitrile-acetate buffer mixtures ranging from 30:20:50 to 50:20:30 (v/v/v), but peak shape and run time needed to be improved (Fig. 5).

Chromatogram before optimization.
In order to evaluate the effect of the most important factor, a face centered central composite design (CCD) which includes six replicates at the zero level was chosen. A central composite design runs needed 20 runs to complete a whole design. As dependent variables, eight responses were chosen: retention time of hydrochlorothiazide (tRHCT), retention time of amlodipine (tRAML), retention time of valsartan (tRVAL), asymmetry of the hydrochlorothiazide peak (As HCT ), asymmetry of the amlodipine peak (As AML ), asymmetry of the valsartan peak (As VAL ) and resolution (Rs AML and Rs VAL ).
Since the peak resolution of AML and VAL was greater than 2.0 on preliminary chromatogram in all further runs this parameter has been excluded (all runs meet the criteria for peak resolution (Rs > 1.5) and do not need further optimization of this dependent variable). Retention time of non-retained peak (t0) was determined by injection of methanol and the k value of HCT was calculated. The obtained result was greater than 1 (not shown in chromatogram). At t0 mobile phase has a small response and cannot be seen on chromatograms. The matrix of experiments and results obtained as an average value of three runs are presented in Table 1.
Appropriate calculations were done with Design-Expert 8.0.7.1 software (Stat-Ease Inc. Minneapolis, MN, USA). A quadratic interaction model was suggested for the relationship between input and output except for As AML where linear model was suggested. Quadratic model interaction is presented by Equation (1):
The estimates of the coefficients for CCD regression models
CCD provides a precise estimation of the experimental errors and the measure of the adequacy of the models (lack of fit). The results were analyzed by the analysis of variance (ANOVA) method, and the results are presented in Table 3.
Statistical parameters of model obtained by ANOVA
The lack-of-fit test was determined by performing a Fischer F test. The high value of F with a very low probability (only model terms with corresponding p-value less than 0.05 are significant at 95% confidence level) implies that there was no evidence of the models lack of fit (lack of fit value is not significant) and the models could be accepted as an adequate representation of the data. In addition, the values of r 2 and adjusted r 2 taking into account the degrees of freedom indicated that the regression model fits the data well. The data collected from the performed CCD design led to the following conclusions: it was noticed that methanol contents in the mobile phase has the largest influence on retention time of hydrochlorothiazide (tRHCT), amlodipine (tRAML) and retention time of valsartan (tRVAL). This effect has a negative sign, which means that the increase of methanol contents leads to a decrease in retention time. Influence of methanol on asymmetries of peaks has positive sign, which implies that increasing of methanol content lead to increase asymmetry of peaks. These results indicate that change of pH does not significantly affect the retention time and peak asymmetry of HCT, while on the retention time and asymmetry of AML and VAL has varying impact (positive sign for tRAML and As VAL and a negative sign for tRVAL and As AML ). Effect of temperature is presented by a negative sign, which implies that increasing in temperature shortens retention time of valsartan and amlodipine and reduced asymmetry of peaks.
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, a multicriterion methodology was employed by means of Derringer’s desirability function [28, 36]. It is based on constructing desirable ranges for each response (individual desirable function, di) 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) as
The relative importance pi is a comparative scale for weighting each of the resulting d i in the overall desirabi-lity product and it varies from the least important (pi = 1) to the most important (pi = 5). It is noteworthy that the outcome of the overall desirability D depends on the pi value that offers users flexibility in the definition of desirability functions.
The 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 were selected. Individual desirability functions range from 0 (undesired response) to 1 (a fully desired response). D value close to 1 means that the combination of different criteria was globally optimal (Figs. 6 and 7). If any of the responses or factors fall outside their desirability range, the overall function becomes zero. The goals of multi-criteria optimization for each response in this paper are presented in Table 4.

Graphical representation of the constraints accepted fot the determination of global desirabilty and obtained optimal conditions (Solution 1).

Graphical representation of the constraints accepted fot the determination of global desirabilty and obtained optimal conditions (Solution 6).
Criteria for multivariate optimization of the individual responses
Software has released 42 solutions that meet the set criteria with resulting desirability (D = 1.000). Review of all the solutions, it was observed that two of the 42 are those that belong to the already existing points of the experiments (Table 5).
Two solutions that meet criteria for analysis
These two solutions are chosen for solutions on the basis of which should be validated method for analysis of samples (Figs. 8 and 9).

Overlay plot highlighted with yellow shows solutions obtained by optimization keeping the temperature at 30°C.

Overlay plot highlighted with yellow shows solutions obtained by optimization keeping the temperature at 40°C.
Figure 10 and Fig. 11 shows graphically and numerically the standard errors of design, depending on the pH and methanol mobile phase. Comparing numerical values of standard errors of design (0.344 and 0.701) it can be concluded that the temperature should be set at 40 °C. The representative chromatogram taken under these conditions are represented in Figs. 12 and 13.

Relation of standard error of design to methanol content and pH, temperature is set as constant value at 30°C.

Relation of standard error of design to methanol content and pH, temperature is set as constant value at 40°C.

Chromatogram after optimization keeping temperature at 30°C.

Chromatogram after optimization keeping temperature at 40°C.
After setting the optimal conditions, the proposed method was validated [37]. Interfering peaks were not detected at the retention time of HCT, AML and VAL, indicating good selectivity of the method. Linear dependence of the peak areas versus concentration was determined for the proposed ranges. Parameters of the linear regression equations were calculated and are presented in Table 6. 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.
Statistical parameters for individual calibration curves
Statistical parameters for individual calibration curves
The precision of the procedure was assessed by analyzing nine solutions containing known quantities of the investigated compounds. Low values of relative standard deviation for repeatability, RSD < 2.5%, and high recovery (Table 7) indicate very good precision of the proposed method.
Precision of the RP-HPLC method
The applicability of the proposed method was examined by analyzing commercially available ExforgeHCT tablets.
The experimental design methodology was used for simultaneous HPLC determination of hydrochlorothiazide, amlodipine and valsartan in combined dosage forms. The significant feature of these combinations lies in the fact that hydrochlorothiazide and amlodipine are present in minute amounts compared to valsartan which makes for a more complicated and tedious analysis. The chemometric approach for optimization of chromatographic separation of hydrochlorothiazide, amlodipine and valsartan has been demonstrated. Central composite design and 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. Central composite design requires a minimum number of experiments to be performed and has been widely used for fitting a second-order model. 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 ten 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.
