Simultaneous Determination of Amlodipine, Hydrochlorothiazide, and Valsartan in Pharmaceutical Products by a Combination of Full Spectrum Measurement and Kalman Filter Algorithm

In this work, a new solution has been found for selecting the approximate initial value of concentration (by means of the classical least squares) and variance (calculated by the Horwitz equation) for the Kalman filter algorithm. With this solution, the Kalman method is less error-prone and has a better repeatability than the least squares method when using the full spectrum. A protocol for simultaneous determination of amlodipine (AML), hydrochlorothiazide (HYD), and valsartan (VAL) in pharmaceutical products was developed based on the spectrophotometry-chemometric method using full spectrummeasurement in combination with the Kalman filter algorithm written in Microsoft Excel 2016 and Visual Basic for Applications (VBA). e method was validated on the Exforge HCT tablets with good repeatability (RSD) (varied from 2.2% to 2.3% (n� 3) for all the three studied compounds) and good recovery (90.0%–94.0% for AML, 90.3%–94.5% for HYD, and 98.5%–103.1% for VAL (n� 3)). e results were in good agreement with the measurements achieved from the high-performance liquid chromatography (HPLC) method.


Introduction
In the development history of analytical chemistry, chemometrics has been used and applied successfully for identification and quantification of a mixture in different matrices, particularly in pharmaceutical samples. e commonly used chemometric methods could be named as partial least squares (PLS), classical least squares (CLS), principal component regression (PCR), artificial neural networks (ANNs), derivative spectrum, and Kalman filter.
Generally, each method has its own advantages and disadvantages. CLS [1] uses all the spectrum data to form a system of m equations and n unknowns (m > n), in which the transformation matrix basing on least squares principle yields acceptable results in terms of relative errors. However, if the input data contain many noises (or errors) and/or there are reactions between analytical compounds causing photometric effects to the absorbance, then this CLS method fails to reduce the noises and consequently large number of errors is expected. Meanwhile, the ANN method [1] requires long time for network formation and a number of algorithms. erefore, to build a suitable model, each potential model developed needs to be tested until the optimal network structure is defined. e derivative spectrum method [1,2], in the meantime, is inapplicable when the sample contains components of overlapping or similar absorbance spectrum since it is difficult to find the suitable wavelength and their derivative spectrum still has maximum absorbance overlapping. e Kalman filter method is able to remove maximally the noise effects, therefore, limiting the measurement errors. is method in combination with spectrophotometry was applied to analyze mixtures of metal ions and components in multicomponent pharmaceutical dosages [3,4]. In the previous works [5][6][7][8], equations for calculating employing Kalman filter algorithm were given. However, the method for selection of initial values including initial concentrations and initial variances for the Kalman filter still has not been dealt with. Based on the pretests (data not shown), it could be addressed that the calculated results strongly depend on the method to select the initial guesses for the Kalman filter. If the selected initial guesses of concentration (C est (0)) and variance (P est (0)) are not suitable or selected concentration is too different from the real value, then it would cause divergent result or big errors. In general, until now, there has not been a comprehensive solution for selection of initial guesses for the Kalman filter.
From the mentioned gaps above, the three main objectives of this study are (i) to find a suitable method to define the initial guesses for the Kalman filter algorithm, (ii) to build a software based on the Kalman filter algorithm and spectrophotometry (called the Kalman method from now on) to determine simultaneously components in a mixture, (iii) to apply the developed Kalman method for simultaneous determination of the compounds with overlapping absorbance spectra in a multicomponent pharmaceutical product, specifically in this study were amlodipine, hydrochlorothiazide, and valsartan in Exforge HCT tablet.

Sample Preparation.
e drug sample was prepared following the method of Galande [9]: 20 tablets were balanced to calculate the average tablet weight (M), which were then finely milled and homogenized. Take m grams, in accordance with one tablet weight (theoretically, mass of AML is 10 mg, of HYD is 12.5 mg, and of VAL is 160 mg), into a 250 mL glass bottle with cap, 50 mL of methanol was added, and the bottle was shaked well before applying ultrasonic extraction for 30 min. e solution was then filtered into a 100 mL volumetric flask and filled up with methanol. e obtained solution was diluted 100 times before measurement by using a V-630 UV/Vis Spectrometer JASCO (Japan).
Accordingly, an average mass of one tablet was M � 4.1252 g.
Mass of an individual compound in a tablet was calculated as follows: where M is the average mass of one tablet (g), m is the mass of the weighted sample (g), C is the concentration of the target compound (μg/mL), V is the initial volume (100 mL), and K is the dilution factor (K � 100). Accordingly, when m � M, we have x � C · 10(mg/tablet). (2)

Kalman Filter Method.
e Kalman filter is a linear parameter estimation technique. In analytical chemistry, it is used to estimate the concentrations of components in a mixture from the absorbance spectra. e initial state of concentration (at the first wavelength) is required. e next concentration state will be estimated based on the initial one. Basically, the model consists of two equations as follows: e equation to describe the chemical system: And the other to describe the measurement process: where C (k) is a vector of state concentrations at point k (which is the wavelength), w (k) is the vector of noise contribution to the system model at point k, A (k) is the measurement at point k, ε (k) is the state transition matrix, and v (k) is the corresponding measurement noise. e Kalman filter algorithm applied in this study for multicomponent spectrophotometry analysis consists of the following equations: (i) State (concentration) estimate extrapolation: (ii) Error covariance extrapolation: (iii) Kalman gain: (iv) State estimate updation: (v) Error covariance updation: e above calculation steps are performed from the first wavelength to the last wavelength. Finally, the calculation program will produce the result: the concentration of each 2 Advances in Materials Science and Engineering constituent in the system and the covariance of the error. is variance is usually the smallest at the last wavelength [5][6][7][8].

Relative Error (RE %).
e relative error was calculated as follows [10,11]: where C is the measured concentration (µg/mL) and C 0 is the known concentration (standard solution) (µg/mL).

Repeatability.
Repeatability was assessed via the relative standard deviation (RSD) value [10,11]: where S is the standard deviation and is the mean concentration after n times of measurement (µg/mL).
For internal laboratory quality control, the method repeatability was approved if the obtained RSDs were lower than a half of the RSD value calculated from the Horwitz function [9][10][11]:

Method Recovery.
Method recovery was calculated based on the spiked samples as follows [10,11]: where a is the spike concentration (µg/mL), C T is the measured concentration after spiking (µg/mL), and C a is the measured concentration before spiking (µg/mL).

Analytical
Procedure. e analytical procedure is shown in Figure 1. e three main analytical steps are shown as follows: Step 1. Prepare the standard solutions or samples Step 2. Measure molecular absorbance spectra, data were recorded as .txt or .dat files Step 3. Extract the files to the computer and run the developed Kalman-Excel program to calculate the specific concentration

Absorbance Spectra of Standard Solutions.
Absorbance spectra of four solutions, namely, AML 5 μg/mL, HYD 5 μg/mL, VAL 5 μg/mL, and mixture of AML 5 μg/mL, HYD 5 μg/mL, and VAL 5 μg/mL in methanol at the wavelength range 230-340 nm were scanned and are shown in Figure 2. As shown in Figure 2, absorbance spectra of AML, HYD, and VAL in ethanol overlapped between 230 nm and 340 nm wavelength, causing difficulty in simultaneously determination of these compounds in mixture. is problem, however, could be solved smoothly using a combination of spectrophotometry and chemometrics.
Within the wavelengths from 230 nm to 340 nm at 0.5 nm intervals, the measured absorbance spectrum of the standard mixture was almost fit with the theory spectrum (estimated from the additive property of absorbance); therefore, absorbance of mixture containing AML, HYD, and VAL had additive property. In other words, it was able to use full spectrum for simultaneous determination of AML, HYD, and VAL by using the combination of spectrophotometry and chemometric method.

Selection of Initial Guesses to Start the Kalman Filter.
As mentioned above, the main challenge to use the Kalman filter is to choose a proper method to identify the initial guesses. A wrong selection would cause an improper calculation. In a mixture containing different substances, the initial guesses are the estimated concentrations in accordance with specific variances of individual substances. Based on the previous studies, there have been two solutions to select the initial guesses: Group 1. Random selection of initial guesses, which means the values of concentration (C) could be randomly assigned, such as 0 or 0.5 µg/mL; and of variance (P) could be 1, etc. [5,8].
Group 2. Assumption of initial guesses, which means either (i) C and P values are subjectively selected based on personal experience and the properties of the samples, (ii) some preliminary experiments are conducted to define the initial C and P values [1,5], or (iii) Beer-Lambert's law is applied at some selected wavelengths to calculate the initial C value (for individual substances in the mixture) by solving linear equation systems, while the variance P is calculated based on a specific guideline for statistical errors in analytical chemistry (i.e., applying Horwitz equation to calculate relative standard deviation (RSD), accordingly the standard deviation and variance at that concentration are defined) [5,7,10,12].
In general, so far there have been no comprehensive method to select suitable initial guesses (C and P values) for the Kalman filter algorithm, raising a challenge for the analytical chemists who want to apply the Kalman filter in their studies. In this study, we investigated three different methods to select the initial C and P values, specifically based on the selection of group 1, group 2, and the proposed selection method of this study.

Random Selection of Initial Guesses.
In this method, C est(0) and P est(0) values could be randomly selected [11]; therefore, this study chose 0.3 µg/mL for C est(0) and 1 for P est(0) . e C est(0) value of 0.3 µg/mL was delivered from the Advances in Materials Science and Engineering common limit of detection (LOD) value of spectrophotometry of about 0.1 µg/mL, of which the limit of quantification would be around 0.3 µg/mL.
Apply the Kalman method for monospectral data and a mixture (AML, HYD, and VAL) of three substances (in the range of 220 nm-340 nm), and the results are shown in Table 1.
e results in Table 1 demonstrated that, in the cases of AML and HYD of mixture H1, when C 0 was close to the selected concentration (0.3 µg/mL), this method provided acceptable relative error (RE) values (20% for AML and −6% for HYD). Otherwise it brought big measurement errors with RE values fluctuated from 40% to 98% for all the target compounds in all the rest investigated mixtures. e cause of this result could be interpreted that the applied Kalman filter considered the random selected concentration of 0.3 µg/mL belonged to another distribution with a certain real value, which was different with this current distribution in accordance with the real values of the studied compounds (1.00 to 9.00 µg/mL). In other words, when the initial concentration is too different from its real value, the Kalman filter using the random selection method for initial guesses will not result in convergence but divergence, causing unacceptable statistical errors.
In short, the initial guesses (C est(0) and P est(0) ) selected in random is not suitable if the selected concentration is too different from the real value, which would cause divergent result or big errors.

Assumption of Initial Guesses.
is assumption was applied in some previous studies [5,7], in which a series of .txt (a) AML 5μg/mL (b) HYD 5 μg/mL (c) VAL 5μg/mL (d) AML 5μg/mL, HYD 5μg/mL, and VAL 5 μg/mL (e) AML 5 μg/mL, HYD 5 μg/mL, and VAL 5μg/mL (theory) preliminary experiments were conducted with different initial guesses (C est(0) and P est(0) ) to identify the suitable initial concentration and variance for the Kalman filter. Moreover, there were studies suggested automatically carrying out the preliminary experiments (employing self-written computer programs); however, the calculation speed was generally slow, which required some repeated calculations to select proper C est(0) and P est(0) values [1,5].
In this study, two different solutions were tested for an assumption of C est(0) and P est(0) values (in mixtures containing 2 or 3 compounds): Solution 1. Based on the equation systems of 2 (or 3) unknowns (concentrations) at 2 (or 3) adjacent wavelengths (this is the equation showing the relationship between absorbance and concentration in a mixture with absorption coefficient (α) calculated from the individual spectrum of the studied compound), concentration of each compound will be calculated which then used as C est(0) for the Kalman filter. Meanwhile, the P est(0) value was randomly selected, for instance by 1. e results are shown in Table 2.
Solution 2. Select 0.3 µg/mL as C est(0) for each compound in a mixture. e initial variance P est(0) was calculated based on the Horwitz equation, which resulted in a value of 0.003 in accordance with a concentration C � 0.3 µg/mL � 3 · 10 −7 [13]. e results are shown in Table 2.
Apply the Kalman method for monospectral data and a mixture of two substances (in the range of 220 nm-340 nm), and the results are shown in Tables 2 and 3.
Solution 1 provided big relative errors (RE % fluctuating from 14% to 82%, except for the case of AML in mixtures H2 and H4, Table 2). is Solution 1 required complicated steps and depended on the two initial wavelengths selected to solve the equation system, which helped us to define the initial concentration. On the contrary, the application on real samples was strongly affected by the matrix, causing big errors.
In Solution 2, although a different approach was applied to define the initial variance (employing the Horwitz equation), the Kalman filter in this stage brought big errors (RE % varying from 7% to 97%).
In short, the above two tested methods to identify initial guesses by either random selection or assumption failed to bring an acceptable result (judged via relative errors). e only exception resulted when the initial concentration selected was close to the real value of concentration in the mixture. It was, therefore, necessary to find a new approach in identifying the initial concentration approximately to the actual value of the analyzed compound in the mixture.

Selection of Approximate Initial Guesses.
Based on the literature review and laboratory experiments for a mixture containing 3 substances, a new method to select initial guesses was proposed, specifically: (i) Apply the classical least squares method to solve the system of m linear equations with n unknowns (in which m was the number of wavelengths selected to scan the absorption spectrum of the mixture and n was the number of compounds in the mixture), followed by the Gaussian elimination method to solve a system of n linear equations with n unknowns for concentrations of compounds in the mixture. e resulting concentrations of the compounds in the mixture were selected as initial concentrations. (ii) Apply the Horwitz equation to estimate the variance corresponding to the concentration of each compound in the mixture, which was considered as the initial variance P est(0) . e calculation of initial variance P est(0) corresponding to the initial concentration C est(0) was as follows: From equation (11), ⇒ S � RSD Horwitz · C est(0) 100 . (15)

Advances in Materials Science and Engineering
In which RSD Horwitz was calculated as (12) and C est(0) represented by fractions. Finally, Noticeably, when conducting repeated measures in a laboratory, if the repeatability (represented via RSD value) smaller or equal to a half of the theory RSD value delivered from the Horwitz equation (RSD ≤ 1/2 RSD Horwitz ), then the method repeatability is acceptable [10]. Accordingly, the initial variance corresponding to the initial concentration C est(0) would be a quarter of the P est(0) value calculated from (13).
e overall experiment was conducted as follows: apply the CLS method to define C est(0) values of AML, HYD, and VAL from the spectrum data-which were the absorbance values of individual compound and mixture solutions measured from 230 nm-340 nm wavelength. From the result of C est(0) , calculate the specific P est(0) values of AML, HYD, and VAL. Finally, provide C est(0) and P est(0) to the designed computer program using the Kalman filter to get the results, as shown in Table 4. e output data provided the concentrations close to the actual values. In other words, this developed method showed good results with small relative errors (RE ≤ 4 %), compared to the two previous methods-using random selection and assumption of initial guesses.
To ensure the applicability of the developed method using selection of approximate initial guesses to identify the initial guesses for the Kalman filter, it is necessary to validate the method in both standard solutions and real samples (pharmaceutical products).

Relative Errors of the Method.
To check the performance of the method, four different mixtures (AML/HYD/VAL, µg/mL) of the three target compounds were prepared, including H1: 0.250/0.325/4.000; H2: 0.50/0.65/8.00; H3: 1.00/1.30/16.00; and H4: 5.00/5.00/5.00. e absorbance spectrum of the prepared mixtures was scanned from 230 nm to 340 nm wavelength. e Kalman-Excel program was then applied to calculate the concentration of each compound in order to identify RE values. e results are shown in Table 5.
Under different mixtures, the obtained RE values of AML measurements varied from 0.4 to 2.2%, of HYD from −1.5 to 1.3%, and of VAL from −3.6 to 0.8% (Table 5). ese low RE values demonstrated the high similarity of the standard concentrations and the measured concentrations of the three studied compounds. In other words, the developed method has good trueness.

Method Repeatability for Laboratory-Prepared Samples.
e similar experiment as described in Section 3.1 was conducted, in which each mixture was prepared and analyzed in triplicate. Method repeatability was assessed based on the comparison between the calculated RSD values and 1/2 RSD Horwitz . e results are shown in Table 6.
e results in Table 6 show that the RSD value for both AML and VAL measurements (n � 3) in 4 different mixtures was 0.4%, for HYD fluctuated from 0.3 to 0.5%. For internal laboratory quality control, the method repeatability was approved if the obtained RSD was lower than a half of the RSD value calculated from the Horwitz function [10]. Accordingly, this developed method has good repeatability.

Method Repeatability and Trueness for Pharmaceutical Samples
Repeatability. e repeatability of the procedure to simultaneously determine AML, HYD, and VAL in the pharmaceutical sample (Exforge HCT tablets, n � 3) is described in Section 2.4.2. Accordingly, the final masses of AML, HYD, and VAL per tablet after preparation were 10.00 mg, 12.50 mg, and 160.00 mg, which were considered as the expected contents. e developed Kalman-Excel program was used for calculating the concentrations of each target compounds. e results of content of target substances and repeatability are shown in Table 7.
Average mass of AML per tablet was 9.61 mg with the RSD value of 2.3% (the expected mass was 10 mg, RSD Horwitz was 8%), of HYD was 11.63 mg with RSD value of 2.2% (the expected mass was 12.50 mg, RSD Horwitz was 5.5%), and of VAL was 169.17 mg with RSD value of 2.2% (the expected mass was 160.00 mg, RSD Horwitz was 5.3%). Apparently, all of the RSD values were lower then the corresponded RSD Horwitz , implying that the developed method was successfully applied to analyze simultaneously AML, HYD, and VAL in pharmaceutical product.

Trueness.
To assess the trueness, in this study, two different approaches were considered, which were (i) define method recovery and (ii) compare the results with the ones analyzed by a validated method: high-performance liquid chromatography (HPLC). (1) Method Recovery. ree replicate samples (B1, B2, and B3) were prepared from Exforge HCT tablets. Different spiked levels of AML, HYD, and VAL were added (Table 4). e absorbance spectrum was scanned from 230 to 340 nm wavelength (0.5 nm step), followed by data computing and concentration calculating by the developed Kalman-Excel program.
e results of method recovery are shown in Table 8.
e average recoveries of AML, HYD, and VAL were 92.9%, 93.3%, and 101.8%, respectively (Table 8). According to AOAC, for the measured concentration from 1 ppm to 10 ppm, the required recovery should fluctuate from 80% to 110% [1]. Based on this, the developed method performed good trueness for all the analyzed substances, suggesting that the excipients caused almost no effects to the analytical results.
(2) Method Trueness Assessment Based on Comparison with HPLC Analytical Measurement. Exforge HCT tablets were sent to the Drug, Cosmetic and Food Quality control Center of ua ien Hue Province for analysis using the HPLC method. AML was analyzed following Vietnam Pharmacopoeia IV guideline, USP 38 was used for VAL and HYD analysis.
Student's t-test [13,14] was used to compare the analytical results of the two methods. e result of comparing mean values of two methods is shown in Table 9. Table 9 shows that all the t exp values were smaller than the corresponding t (0.05; f ) of the three target compounds, demonstrating that the analytical results obtained from the developed Kalman-Excel method were in agreement with the ones obtained from HPLC measurements.

Conclusions
In this work, a new solution has been found for the first time, selecting the approximate initial value of the concentration (by means of the classical least squares) and variance (calculated by using the Horwitz equation) for the Kalman filter algorithm. is new solution allows convenient application of the chemometric-spectrophotometric method using the Kalman filter algorithm (Kalman method) to simultaneously determine two or three substances in their mixture with an UV-Vis absorption spectrophotometer. e Kalman method is less error-prone and has a better repeatability than the least squares method when using the full spectrum. A computer program that uses the Visual Basic for Applications programming language written on the basis of Microsoft software Excel 2016 based on the Kalman filter algorithm has been written, which allows quick and convenient calculation when applied on practical testing of pharmaceutical products in laboratories.
First, the process of simultaneous analysis of three active ingredients, i.e., amlodipine, hydrochlorothiazide, and valsartan, was established in multicomponent pharmaceutical formulation by the Kalman method using full spectrum without any separation technique. e process exhibited good repeatability and trueness for all the three analyzed compounds with RSD <2.5% (n � 3), recovery varied from 93 to 102%, and the received analytical results were identical with ones of HPLC method. e process was not only simple to implement but also reduced the cost of analysis compared to the standard method of high-performance liquid chromatography (HPLC).

Data Availability
e data used to support the finding of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.