Abstract

This work used dispersive Raman spectroscopy to evaluate acetaminophen in commercially available formulations as an analytical methodology for quality control in the pharmaceutical industry. Raman spectra were collected using a near-infrared dispersive Raman spectrometer (830 nm, 50 mW, 20 s exposure time) coupled to a fiber optic probe. Solutions of acetaminophen diluted in excipient (70 to 120% of the commercial concentration of 200 mg/mL) were used to develop a calibration model based on partial least squares (PLSs) applied to Raman spectra of solutions and, subsequently, obtain linearity, accuracy, precision (repeatability), and sensitivity of the method using the near-infrared spectroscopy (NIRS) as a gold standard method. This model was used to predict the acetaminophen concentration in commercial samples from different lots of acetaminophen formulations (200 mg/mL) with a PLS-prediction error of about 0.6%. Commercial medicines had PLS predicted concentrations errors below 2.5%, whereas NIRS had an error of about 3.7% compared to the label concentration. It has been demonstrated the applicability of Raman spectroscopy with fiber probe for quality control in pharmaceutical industry of commercial formulations.

1. Introduction

The acetaminophen is a well-known analgesic and antipyretic that has been widely used because of the low side effects and still being well tolerated when administered in therapeutic quantities. Routine techniques for measuring acetaminophen include chromatographic methods (high performance liquid chromatography—HPLC), spectrophotometric methods (ultraviolet/visible, near-infrared and infrared spectroscopy, including Fourier-transform near-infrared spectroscopy (NIRS), and fluorescence spectroscopy), flow-injection spectroscopy, chemiluminescence methods, capillary electrophoretic methods, and electroanalytical methods [1].

The classical techniques for pharmaceutical drugs evaluation, such as chromatography, are sensitive and very specific, but these techniques are time consuming and of high costs, needing samples to undergo a preprocessing prior analysis, which may ultimately destroy the sample. To address this problem, optical methods have been proposed for the quantitative determination of active pharmaceutical ingredients (APIs) in drugs and medicines, most of them based on absorption and reflectance spectrometric techniques in the ultraviolet-visible and infrared regions [2, 3], with advantages such as rapid API evaluation with minimum or even without sample destruction [4].

Vibrational spectroscopy such as NIRS and Raman spectroscopy have advantages which lie on spectral selectivity and representativity and have been proposed as a quality control tool. The NIRS has grown substantially as a quality control technique in the pharmaceutical industry, including analysis of acetaminophen [5]. It allows chemical and physical information about samples to be derived from spectra. However, due to the nature of NIRS data (dependent of changes in the molecular dipole moment) and high sensitivity of water bonds, the low sensitivity of the technique restricts its scope to major components. Moreover, near-infrared spectra exhibit strong overlap [6]. Raman spectroscopy provides information about the chemical structure, that is, the molecular composition of the material under study, as well as its quantification [7, 8]. Dispersive Raman spectrometer makes possible a fast scanning/imaging through the use of an imaging spectrograph and multichannel detector such as CCD (charge-coupled device) camera [8], with lower intrinsic fluorescence by using near-infrared excitation (785 and 830 nm wavelength) [810].

Raman spectroscopy has been used in the quantitative analysis of API, noninvasively and nondestructively, with advantages such as the need for small amount of sample, virtually no sample preparation, no sample destruction, and short time to perform the analysis [4]. Raman spectra can be obtained from different pharmaceutical presentations, even directly in the container, aiming the quantification of the API and other formulation constituents with low error [6, 7, 11]. Raman spectra could be obtained direct in biological tissues, quantifying injected drugs in situ [4].

Chemometric techniques such as partial least squares (PLSs) have been used in the analytical chemistry in order to predict the amount of a specific compound in a sample using the spectral information [12]. Chemometrics have a great significance in the pharmaceutical area, since it deals with the analysis of chemical data of multivariate nature, collected through multichannel detectors, where the intensities of the signal in different wavelengths are recorded in one single spectrum [13]. Chemometrics uses mathematical and statistical methods to provide as much chemical information as possible with the analysis of the obtained data, even in the presence of experimental noise, interferences of unknown nature, collinearities, and nonlinearities, besides it is sufficiently robust to allow no need for changes in the input parameters with the inclusion of new samples to the model [1214].

The process analytical technology (PAT) aims at the development of strategies to the control of the production processes in the pharmaceutical industry [1517]. The PAT Guide attempted for the need for research into new and faster analytical methodologies for quality control of drugs to be used on-line and at-line. Optical techniques and multivariate statistics such as PLS have been claimed for a consistent, reliable, rapid, and low-cost analysis of chemical species in the industrial field for analytical instrumentation [1517], which could successfully attend the PAT requirements. In conjunction, the Expert Working Group under the International Conference on Harmonisation of Technical Requirements for Registration of Pharmaceuticals for Human Use (ICH) developed guidelines and set parameters to be determined in order to validate the analytical techniques [18, 19], according to PAT directions.

The purpose of this study was to assess the feasibility of a quantitative evaluation of acetaminophen (commercial drop’s presentation) based on dispersive Raman spectroscopy using fiber optic probe and multivariate statistics PLS (RS/PLS quantitative model) and to verify selectivity, linearity, range, precision, and accuracy, in order to check the ability of Raman spectra and PLS to be used as an analytical technique to attend industries requirements and compare with the NIRS technique already used for this. The developed model was then employed to predict the concentration of acetaminophen in samples of commercial medicines containing this API with 200 mg/mL (drop’s), aiming at its use as an analytical technique of choice. The results were compared to the ones obtained by NIRS, assigned by the manufacturer of the drugs.

2. Material and Methods

2.1. Solutions Preparation

Solutions of nine different concentrations of acetaminophen were prepared by the laboratory of a leading manufacturer, representing the range from 70 to 130% of the concentration found in the commercially available drug, diluted in the same excipient used in the commercial drug (Table 1). Also, samples of the excipient (placebo) and of the API diluted in polyethylene glycol (the main component of the excipient) were obtained. The placebo is composed basically by polyethylene glycol, sodium cyclamate, sodium saccharin, sodium benzoate, sodium metabisulfite, sodium hydroxide, citric acid, yellow dye, flavors, and deionized water. The amount of each excipient component could not be revealed, because it is copyrighted. Solutions were prepared in 6 replicates (one used to develop the model and the remaining in the prediction), and the nominal concentration of acetaminophen in each solution was confirmed by the manufacturer using NIRS. The dilution of 260 mg/mL was withdrawn due to the low solubilization of the API in the excipient. Also, it was obtained six flasks of commercial medicines containing 200 mg/mL of acetaminophen (drop’s presentation) from different lots in local drugstores.

2.2. Raman Spectroscopy

The solution samples and commercial medicines were submitted to the near-infrared dispersive Raman spectroscopy technique using the Raman spectrometer detailed elsewhere [21], connected to a fiber optic probe as shown in Figure 1. It is composed by an 830 nm diode laser (Micro Laser Systems Inc., Model L4830S, CA, USA), coupled to a Raman probe made up of low-OH fused-silica optical fibers with “6 collections around 1 excitation” configuration, as described elsewhere [20], using a 𝑓 = 2 5  mm lens. The measured power at excitation fiber tip was about 50 mW. Solutions and medicines were placed in a quartz cuvette of 4 mL for Raman measurements. The fiber optic cable distal excitation/collection end was immersed 1 mm inside the solution, and Raman signal was collected in about 20 s. The signal emitted from the samples and collected through the optical cable was coupled to the entrance slit of the spectrograph (Chromex, Model 250IS, MA, USA), which had a notch filter at 830 nm before the entrance (Spectral Iridian Technologies, model PN-ZX 000080, ON, Canada) to reject the scattered light from the laser. The spectrograph slit was set to 150 μm, providing a Raman spectral resolution of about 10 cm−1 [22]. The signal dispersed by the spectrograph was detected by a liquid-nitrogen cooled CCD camera and controller (Princeton Instruments, model LN/CCD-1024-EHR1 and ST130, NJ, USA) connected to a computer for acquisition and storage of the spectra. Each sample was scanned at the same day and experimental conditions.

The procedures of spectra calibration, fluorescence background subtraction, cosmic rays removal, and normalization were performed using the software OriginPro 7.5 (OriginLab Corp., MA, USA). The Raman shift calibration was done by collecting the spectrum of naphthalene and correlating the spectrum band positions (pixel) with the known band shift (cm−1) with a 3rd-order polynomial fitting. To remove the fluorescence background, it was performed the baseline correction, in which straight lines were fitted in the spectra valleys and subtracted from the gross spectrum, remaining the high frequency Raman bands. Spectra were then normalized using the band of polyethylene glycol (main constituent of the excipient) at 1450 cm−1 and plotted in the spectral range of 600 to 1800 cm−1.

2.3. PLS Quantitative Model

The PLS uses the nominal concentration information to obtain factors or latent variables, which are correlated to the principal components vectors, calculated by the PCA (principal component analysis) method [13, 23]. The first factors describe the direction of the maximum variance (the first principal components) which, in the majority of the cases, is correlated to the concentrations of the API [13, 23, 24]. A multivariate calibration model using PLS was developed with a spectral dataset of API dilutions in the concentration range of 140 to 240 mg/mL. Then, this model was validated using a dataset composed of samples randomly chosen from the replicate dilutions in the same range, and the standard error of validations was calculated.

In the PLS model, the Raman spectra and the nominal concentrations were considered as independent (ordinates) and dependent (abscissas) variables, respectively, and, under Matlab 4.2 and PLS Toolbox 1.5, it was used the command [ 𝑝 , 𝑞 , 𝑤 , 𝑡 , 𝑢 , 𝑏 , s s q ] = p l s ( 𝑥 , 𝑦 , l v , 1 ) to build the model, where 𝑥 = spectra, 𝑦 = concentrations, and lv = number of factors (latent variables). The PLS outputs are 𝑝 = loading of independent variables, 𝑞 = loading of dependent variables, 𝑤 = importance of the independent variables, 𝑡 = scores of independent variables, 𝑢 = importance of the dependent variables, 𝑏 = vector of PLS internal coefficient, and s s q d i f = variances between 𝑥 and 𝑦 obtained by the PLS. Since there are only few variables that are responsible to the spectrum of the drug (acetaminophen and excipient), one expect that only a few lv are needed to build the model.

In order to define the number of lv to be used in the model and to predict the concentrations of the dilutions and drug samples with the least error of prediction, it was established a cross-validation procedure, using the Raman spectra (independent variables) from the replicate dataset, with the number of lv changing from 1 to 5, through the function 𝑦 p r e d = p l s p r e d ( 𝑥 , 𝑏 , 𝑝 , 𝑞 , 𝑤 , l v ) . The results of the nominal concentrations versus the ones estimated by the PLS model ( 𝑦 p r e d ) were plotted, and the standard error of prediction (SEP) was calculated for each lv. The purpose of the SEP is to determine the error of the PLS model when applied to a different (in our case the replicate) dataset in the same range. The SEP can be calculated as follows (using Excel function STEYX): S E P ( % ) = 1 ( 𝑛 2 ) 𝑛 𝑖 = 1 𝐶 𝑎 𝑖 𝐶 𝑎 𝑖 2 𝑛 𝑖 = 1 𝐶 𝑖 𝐶 𝑖 𝐶 𝑎 𝑖 𝐶 𝑎 𝑖 2 𝑛 𝑖 = 1 𝐶 𝑖 𝐶 𝑖 2 , ( 2 . 1 ) where 𝐶 𝑎 𝑖 is the nominal concentration of the API, 𝐶 𝑖 is the concentration predicted by the PLS model, and 𝑛 is the number of samples. Figure 2 shows the SEP values depending on the number of lv used in the cross-validation. It was found that the SEP was 1.19 mg/mL for lv = 1, increasing for lv = 2 and decreasing for higher lv values.

The spectra of dilutions and commercial drugs were then used in the PLS model as new independent variables using the number of lv = 1 and the predicted concentrations were then calculated.

2.4. Quantitative Evaluation: Figures of Merit

The predicted concentration of the replicate samples were used to estimate the parameters [18, 19]: range, linearity, selectivity, precision, and accuracy of the RS/PLS method.

Range
In this study, it was defined the range of 70% to 120% (140 to 240 mg/mL), due to problems of solubility of the API at 130%.

Selectivity
It is the ability of the method to measure a compound in the presence of other components, such as adjuvants. For this parameter, it was considered a solution of the acetaminophen with 200 mg/mL in pure excipient and in polyethylene glycol (6 replicates of each).

Linearity
It is the ability of the method to demonstrate whether there is proportionality between the predicted and the nominal concentrations of the API in the sample, within a specific range. This parameter was determined by constructing a curve with solutions of five concentrations: 160, 180, 200, 220, and 240 mg/mL. We calculated the curve slope (m) which determines the proportionality of the analytical function, the linear coefficient (h), and the correlation coefficient (r). The criterion for acceptance is 𝑟 = 0 . 9 9 .

Precision (Repeatability)
It is referred to the agreement between the results within a short period of time with the same analyst and same instrumentation. It is done by using at least nine determinations under the linear range of the method, that is, three concentrations (low, intermediate or 100%, and high) with three replications each. For the study, it was used the concentrations 160, 200, and 240 mg·mL−1 in triplicate, totaling nine solutions. The precision was expressed as relative standard deviation (RSD), according to the formula R S D = S D D A C × 1 0 0 , ( 2 . 2 ) where SD is standard deviation and DAC is determined average concentration. The maximum RSD for the method should not exceed 5%.

Accuracy
It is the closeness of the results obtained by the method under study in relation to the true value (nominal). It is obtained by adding known amounts of the drug to the components of the excipient. It was determined from nine determinations in the range of the procedure, that is, concentrations of 160, 200, and 240 mg/mL in triplicates. The equation used to calculate the accuracy (AC) is A C = D A C V A × 1 0 0 , ( 2 . 3 ) where VA is value added.

3. Results and Discussion

The Raman spectra of the acetaminophen solutions in the range of 600 to 1800 cm−1, with concentrations ranging from 140 to 240 mg/mL and normalized by the band at 1450 cm−1, are shown in Figure 3. These spectra are dominated by peaks from acetaminophen and in some extent by the excipient, mainly ethylene glycol. Figure 4 shows the Raman spectra of acetaminophen: Figure 4(a) is the spectrum obtained indirectly by subtracting solution G (200 mg/mL in excipient) from solution J (only excipient); Figure 4(b) is the Raman spectra of acetaminophen, obtained from Sigma-Aldrich web page and used as a reference of band positions. It can be seen that the reference spectra present bands in the positions closer to the spectra obtained indirectly, despite a few differences in the resolution and intensity of particular bands. Therefore, the main bands of acetaminophen are, according to the literature, amide I (C=O) at 1667 cm−1, N–H at 1622 cm−1, amide II (C–N stretching, N–H bending) at 1556 cm−1, C–H at 1329 cm−1, and phenyl ring at 799 cm−1 [2527].

By applying the PLS model in the Raman spectra of replicate samples (cross-validation) with selected number of lv, one could calculate the predicted concentration and the SEP that the model is subjected to. Figure 5 shows the real versus predicted concentration for the cross-validation. It can be seen an excellent agreement of the data points to the zero-error diagonal, revealing an excellent model fit. The SEP obtained in the validation procedure (using 1 lv) was of 1.19 mg/mL or 0.6%.

The results of the parameters obtained by RS/PLS are as follows.

Selectivity
The higher relative error of acetaminophen in standard excipient compared to acetaminophen in pure PEG was 2.4%, showing that the technique allows detection of the major component of the sample in different excipient (Table 2).

Linearity
It was determined the concentrations of five solutions (160, 180, 200, 220, and 240 mg/mL) and compared with the nominal concentration values. From the curve, it was obtained 𝑚 = 1 . 0 2 5 and = 6 . 2 0 4 . The linear regression coefficient obtained was 𝑟 = 0 . 9 9 , indicating adequate correlation between the data predicted concentration and the dilution effectively performed.

Precision (Repeatability)
It was used three concentrations with six replicate samples of the following concentrations: 160, 200, and 240 mg/mL. It was calculated the RSD for each concentration of RS/PLS compared to the nominal dilution (Table 3). The acceptance criteria established was RSD ≤ 5%. For solutions with low and intermediate concentrations, the obtained RSDs were found to be below the maximum limit specified. For high concentrations, the resulting RSD was found to be 6.0%, which may be due to an outlier. When this outlier is withdrawn (Table 3), the error is minimized, being a RSD value comparable to the lower and intermediate concentrations.

Accuracy
It was calculated by using the predicted concentrations of the 160, 200, and 240 mg/mL, in triplicate, being determined by means of RSD and recovery, as shown in Table 3. The acceptance criteria established for relative deviation were ≤2%, being that satisfactory recovery values are found between 98 and 102%. The solutions with intermediate and high concentrations showed relative deviation of 1.0 and 1.5%, respectively, being within the limits of acceptance, unlike the solution of low concentration, which did not provide a satisfactory level. The same occurred with the recovery values, where the solution of low concentration showed a value of 97.1% as opposed to middle and upper solutions that showed the values 99 and 98.5%, respectively.
In this work, the RS/PLS was evaluated by parameters as linearity, accuracy, and precision in the same way of Orkoula et al. and Rossignoli et al. [28, 29]. It has been shown that RS/PLS showed excellent linearity with 𝑟 = 0 . 9 9 , high precision, with DPR as low as 1.7%. Corroborating the method validation Orkoula et al. [28], commercial samples of the reference drug containing 200 mg/mL of acetaminophen were also used to validate the model. These samples showed very low errors, not exceeding 2.5%, suggesting that the dispersive Raman spectroscopy using optical probe can detect the amount of API presented in the commercial samples.
As the study of Szostak and Mazurek (2002) [27], the Raman model developed here using PLS was used to quantify the commercial samples containing acetaminophen, by applying the PLS model to the six samples of commercial drugs. Table 4 shows the results of the RS/PLS predicted concentrations and NIRS concentrations and the errors for the commercial medicine. The predicted concentrations were highly correlated to the concentration of the medicine presented in the label; all drugs fell below the ±2.5% error limit, indicating excellent results in determining the API concentration. The same was done with the NIRS values, and the largest error found was 3.7%. The present study used a fiber optic probe made up with low-OH fused-silica fibers, with the intention of getting reproducibility in the excitation/collection geometry of the volume sampled [20].
In order to achieve optimum yield and product quality, the global pharmaceutical industry is searching for better methods for tracking drugs throughout the production process, since most of the high-performance techniques currently used need sample preparation using solvents, are destructive and time consuming, and mostly need expensive inputs [7, 30, 31]. Recent literature indicates that the dispersive Raman spectroscopy, combined with chemometric method PLS (RS/PLS), would meet the expectations for a reliable quantitative determination of active ingredients in pharmaceutical drugs [7, 32, 33].
In complex systems, such as API in medicines and formulations, the classical univariate regression method may become unreliable because of possible band interactions and superimposition. The PLS method has been suitable for quantitative measurement in chemical analysis when applied to spectroscopic data [13], calculating the amount of API in prepared solutions or directly in the commercial medicines [14], quantification of nimesulide in excipient and commercial samples [26], and quantitative determination of acetylsalicylic acid and acetaminophen in tablets [27]. The PLS model showed a lower prediction error when using lv = 1, indicating that only one spectral variation was found in the samples and is related to the changes in the acetaminophen concentration, since most of pharmaceuticals consist of a few APIs and a lot of additives [27]. Since the intensity of a Raman band is associated with changes in the polarizability of the molecules presented in the sample, this has a relevant consequence in investigating pharmaceutical materials because APIs usually contain aromatic conjugated systems which show strong Raman transitions and, in contrast, excipients are generally weaker Raman scatterers due to a lack of C bonds in most formulations [27].
According to the PAT guide, the gains in processes quality, safety, and/or efficiency can change depending on the product and these changes come from factors such as prevention of rejections, disposal, and reprocessing of products; possibility of carrying out analysis in real time with reduced or no sample preparation [17, 34], leading to a reduction of human error and increasing the safety of the operator and no risk to the environment [5, 18, 19]. The RS/PLS model could be able to address most of these factors, because each sample was irradiated rapidly and examined directly with the use of an optical probe, without prior preparation, in real time at line, which leads to minimization of human error.

4. Conclusions

This study demonstrated the effectiveness of the near-infrared dispersive Raman Spectroscopy with fiber optic probe as a quantitative method, with high linearity, selectivity < 2.5%, precision < 5%, and accuracy < 3%, with recovery values of <5% and SEP of 1.19 mg/mL (0.6%) for quantifying the acetaminophen diluted in excipient. For the quantification of API in the commercial formulation, the model showed concentration values as closer as 200 mg/mL, with less than 2.5% error compared to the labeled concentration, while, by NIRS method, the error found was 3.7%. These results suggest that RS/PLS with fiber optic probe could replace the NIRS method in quality control of the product under study, while maintaining the safety and efficacy.

Acknowledgments

V. G. Borio thanks CAPES/PROSUP for the Master Fellowship. L. Silveira Jr. and R. A. Nicolau thank CNPq for the Productivity Fellowship (process no. 305610/2008-2 and 314455/2009-4, respectively). H. P. M. de Oliveira thanks FAPESP (2006/56701-3) and CNPq (479655/2008-1) for the financial support.