Faisal K. Algethamia,
Sherif M. Eid*b,
Khadiga M. Kelanicd,
Mohamed R. Elghobashybc and
Mohamed K. Abd El-Rahmanc
aChemistry Department, Faculty of Science, Imam Mohammed ibn Saud Islamic University, Riyadh, Saudi Arabia
bAnalytical Chemistry Department, Faculty of Pharmacy, October 6 University, 6 October City, Giza, Egypt. E-mail: Sheriefmohammed@o6u.edu.eg; Sherief055@icloud.com; Tel: +201111609539
cAnalytical Chemistry Department, Faculty of Pharmacy, Cairo University, El-Kasr El-Aini Street, ET-11562 Cairo, Egypt
dAnalytical Chemistry Department, Faculty of Pharmacy, Modern University for Technology and Information, Cairo, Egypt
First published on 17th February 2020
The use of performance-enhancing drugs is prohibited in sports competitions according to the World Anti-Doping Agency (WADA) regulations. Here, ATR-FTIR spectroscopy coupled with a partial least squares regression (PLSR) chemometric tool was used for the detection of the misuse of such substances. Bambuterol and its metabolite terbutaline have been included in the list of prohibited doping agents. Therefore, we used bambuterol and terbutaline as models for the accurate and simultaneous qualitative and quantitative analysis of bambuterol and terbutaline in human urine samples. The method was straightforward and once the urine samples were collected, they could be directly applied to the surface of the ZnSe prism (ATR unit) to get the results within one minute. A calibration set with a partial factorial design was used to develop the PLSR model that could be used to predict the concentration of unknown samples containing the two drugs. The developed method was carefully validated and successfully applied to the urine sample analysis of human volunteers. The drugs were quantified at nanogram level concentrations. A side-by-side comparison of the proposed method with the routine GC-MS method was performed to demonstrate the challenges and opportunities of each method.
The WADA rules classify β2-agonists such as bambuterol HCl (BAM) and terbutaline (TER) as stimulants and anabolic agents with a restricted specific regulation.3 Athletes often use BAM for therapeutic purposes as a treatment for exercise-induced asthma.6 However, there are some adverse effects of the use of β2-agonists, such as anabolic and stimulation effects with the increase in myocardial contraction and heart rate, increase in the accretion of muscle proteins and a reduced rate of fat deposition;7 thus, this drug and its metabolite TER8 have been included in the list of prohibited doping agents.3 The detection of BAM and TER as doping drugs for routine analysis in urine samples has been accomplished mainly by chromatography techniques such as GC/MS9,10 and LC-MS/MS.11 Although these methods have high potential in the fight against doping drugs, they suffer from many disadvantages, such as the use of very expensive instruments, being non-portable, need for highly qualified operators, high operation cost, moderate throughput, limitations with respect to solvents and samples, long operation times, and too many sample preparation steps.
It has been proven that attenuated total reflectance Fourier transform infrared spectroscopy (ATR-FTIR) has ability for the quantitative and/or qualitative determination of biological factors and samples such as lung and prostate cancer biomarkers,12 protein analysis,13 and petroleum component analysis.14 Chemometrics is a well-established technique in forensic analysis,15,16 and partial least squares regression (PLSR) does not require the accurate knowledge of all the components present in urine samples. It can make concentration predictions even in the presence of interfering compounds.17 Also, the direct application of urine samples to the surface of the ATR prism as a sampling technique is sufficient for the quantitative determination of the early-stage malaria parasite,18 enzyme kinetics parameters,19 biomarkers of glomerulonephritis,20 tumor cells,21 and malignant ovarian tissues.22 A recent study compared the performances of micro-ATR-FTIR and GC-MS as identifying tools for micro-plastics and fibers isolated from different water sources.23
Therefore, this paper describes the application of ATR-FTIR in combination with a PLSR chemometric tool for the determination of the doping drugs BAM and TER. Low detection limits were obtained due to the high selectivity, sensitivity, portability, high throughput analysis, and straightforward nondestructive direct determination capabilities of PLSR/ATR-FTIR.24–26 Only a few microliters (μL) of urine samples were sufficient for the quantitative determination of the analytes. The developed anti-doping method may be considered as a preliminary step, opening the way for further investigations on other doping agents.
Coded levels | MIX no. | BAM (ng ml−1) | TER (ng ml−1) | |
---|---|---|---|---|
BAM | TER | |||
0 | 0 | 1 | 5000 | 5000 |
0 | −2 | 2 | 5000 | 50 |
−2 | −2 | 3 | 50 | 50 |
−2 | 2 | 4 | 50 | 10000 |
2 | −1 | 5 | 10000 | 2500 |
−1 | 2 | 6 | 2500 | 10000 |
2 | 0 | 7 | 10000 | 5000 |
0 | −1 | 8 | 5000 | 2500 |
−1 | −1 | 9 | 2500 | 2500 |
−1 | 1 | 10 | 2500 | 7500 |
1 | 2 | 11 | 7500 | 10000 |
2 | 1 | 12 | 10000 | 7500 |
1 | 0 | 13 | 7500 | 5000 |
0 | 2 | 14 | 5000 | 10000 |
2 | 2 | 15 | 10000 | 10000 |
2 | −2 | 16 | 10000 | 50 |
−2 | 1 | 17 | 50 | 7500 |
1 | −2 | 18 | 7500 | 50 |
−2 | 0 | 19 | 50 | 5000 |
0 | 1 | 20 | 5000 | 7500 |
1 | 1 | 21 | 7500 | 7500 |
1 | −1 | 22 | 7500 | 2500 |
−1 | −2 | 23 | 2500 | 50 |
−2 | −1 | 24 | 50 | 2500 |
−1 | 0 | 25 | 2500 | 5000 |
A representative set of 25 mixtures was prepared by taking the proper volumes from the prepared working solutions into the volumetric flasks. Then, the volumes were completed to the mark with human urine. The prepared mixtures were separated into two sets: a calibration set of 15 synthetic mixtures and a validation set of 10 mixtures (Table 1).
Fig. 1 (A) ATR-FTIR absorption spectra of the calibration set mixtures; (B) 3D plot of the 1st derivative of the calibration set mixtures. |
Complete sample dryness was ensured by the application of 10 μL of human urine over the ATR prism; then, a spectrum was obtained each minute during the air-drying process. Complete sample dryness occurred after two minutes. The cleanness of the ATR prism was tested regularly between the measurements, where an absorption spectrum was recorded without the sample to make sure that it had no peaks higher than the normal environmental noise level.
The derivative calculations were applied to the recorded spectra from the earlier steps using the following parameters: derivative order: 1, polynomial order: 2 and smoothing points: 5, as shown in Fig. 1B.
These pre-processing steps could be registered on the Unscrambler X software to be performed automatically for the prediction of unknown sample concentrations.
Finally, the PLSR chemometric algorithm29 was applied to the derivative spectra obtained from the last step (shown in Fig. 1B). The recorded absorbance spectra were related to the concentrations by Beer's law.30 The Regression calibration curves were plotted between the known and predicted concentrations.
The study involving human participants was performed in strict accordance with the institutional ethical standards of the Helsinki Declaration of 1964 and its later amendments. It was approved by the local research ethical committee of faculty of pharmacy, Cairo University, Egypt. Healthy volunteers were fully informed verbally about the objectives, nature and possible risks of this study and a written informed consent was obtained from all the volunteers involved in the study.
The analysis of BAM in the presence of its metabolite TER is a challenge as they are structurally related compounds. Fig. 2 shows the structures of BAM and TER. BAM (RS-5-(2-tert-butylamino-1-hydroxyethyl)-m-phenylene bis dimethyl carbamate) is metabolized in the liver by the cleavage of the ester bond and replacement of the branches of N,N-dialkyl carbamate by OH groups to give TER.31 This hydrolysis causes the main differences in the structures.8
Fig. 2 shows the overlay plot of the two spectra of BAM and TER and their characteristic peaks can be distinguished. For more illustration, the spectra were analyzed by the assignment function of the ACD/SpecManager software to show the peak differences between BAM and TER, as summarized in Table 1S.† The FTIR spectra of BAM and TER can be considered as their own fingerprints. The BAM spectrum has a characteristic CO group that shows a peak at 1720 cm−1, which disappears in the TER spectrum. The proposed method was convenient for the simultaneous qualitative and quantitative determination of BAM and TER.
The performance of multicomponent analysis is based on the spectral model used and the wavelength selected. PLSR with a non-linear iterative partial least squares algorithm (NIPALS) was calculated using the Unscrambler software on the Savitzky–Golay first derivative algorithm for the selected spectrum range (400–4000 cm−1), due to which all the available measured wavelengths are often used without neglecting any region that may contain important data.
The key step in building the PLSR model is the selection of optimum number of factors to achieve correct concentration prediction in PLSR calibrations to prevent data over-fitting. PLSR models with 1 to 7 factors were calculated for each wavelength region and validated using “leave-one-out cross-validation”.32,33 Briefly, the model is repeatedly recalculated leaving out a single sample individually at a time and then used to predict the concentration of the left out observation. PLSR calibration using 14-mixture spectra was performed and the regression equation was calculated. The concentration of the left-out sample was calculated using the obtained regression equation. The same procedure was employed for the 15-mixture analysis; each sample was left out once and the predicted concentration was calculated for each sample. The predicted concentrations were plotted against the known concentrations of each mixture. The obtained recoveries (between 95 and 105%) indicate the validity of the method for its application to unknown samples. The latent factors that explain the spectral matrix are sorted in decreasing order according to their contributions to the spectral features. Using the optimal pre-processing method, latent factors were selected with a cross-validation method, as mentioned above. This method depended on the selection of the model including the smallest number of factors that give the minimal root mean square error of calibration (RMSEC) and the minimum root mean square error of validation (RMSEV), by plotting the difference of RMSEV against the number of factors for the analysis of each mixture. The PLSR method proved that (2) factors were found to be optimum for BAM and TER determination. The PLSR calibration model was constructed to achieve correct quantification and prediction with minimum errors.34,35
The PLSR model's predictive ability was evaluated using two procedures: the 1st procedure was achieved by plotting the reference concentrations against the corresponding predicted concentrations, as shown in Fig. 3; then, the internal validation method was performed (estimation of BAM concentrations in its own calibration set). A satisfactory correlation coefficient (r2) (0.9988) and good recovery (101.09%) were obtained, indicating high predictive abilities of the model. The 2nd procedure for the determination of the predictive ability of the method was conducted by calculating the standard error of prediction (SEP).36
Fig. 3 Plots of reference concentrations against the predicted concentrations of the calibration sets of (A) BAM and (B) TER. |
Table 2 summarizes the characteristic parameters of PLSR chemometric urine analysis. The obtained very low BIAS (0.00048) means a very small acceptable error of prediction. Correlation is the linear relationship between the predicted and reference values in the plot, which is very close to 1, indicating that the predicted and reference concentrations are highly correlated. RMSEC measures the dispersion of the calibration samples around the regression line. BIAS is the mean of all the numbers that lie above or below the regression line and the values close to zero show a random distribution of points around the regression line. The obtained RMSEC and standard error of calibration (SEC) are very close, indicating insignificant BIAS.
Parameter | Compound | |
---|---|---|
BAM | TER | |
Elements | 15 | 15 |
Range (ng ml−1) | 50–10000 | 50–10000 |
Slope | 0.9842 | 0.9892 |
Offset | 76.98 | 53.97 |
Correlation | 0.9921 | 0.9945 |
R2 | 0.9842 | 0.9892 |
RMSEC | 456.01 | 365.67 |
SEP | 432.87 | 510.50 |
Bias | 13 × 10−4 | 13.3 × 10−4 |
External validation was achieved by a validation set of 10 mixtures and it was not implemented in model building. The concentrations of BAM and TER in the validation set solutions were selected randomly but were in the range of the calibration set. The method has high resolving power with good mean recoveries (97.3 and 99.1%) for BAM and TER, respectively.
Regression coefficient plots were obtained to summarize the relationship between all predictors and responses by analyzing the 15-mixture 1st derivative spectra to find the areas (highlighted with black marks) that the software depended on during the building of the model for the detection of either BAM or TER, as shown in Fig. 4. This plot proves and validates the ability of the developed method to differentiate between BAM and TER within the matrix of urine components.
Urine samples collected from the same volunteer before and after the oral administration of the Bambec® tablets were scanned. The urine sample that did not contain any drug was used as a blank and scanned as a background spectrum. Then, the samples collected after the oral administration of Bambec® tablets were scanned against this background to subtract any substances that may be found in urine and may interfere with the analysis results. Thus, the obtained peaks were due to the presence of BAM and TER in urine with different concentrations, as confirmed from Fig. 1S.†
Fig. 5 shows the obtained spectra of human volunteers' urine at different times (1, 2 and 6 h). It can be easily noticed that the peak intensity moderately increases with time, which may be due to the increased excretion of the drugs in urine with time. The obtained spectra were highly related to the spectra of original BAM and TER shown in Fig. 2.
Fig. 5 Typical ATR-FTIR spectra of the urine samples obtained from volunteers at various time intervals (1, 2 and 6 h) after a single oral dose of 20 mg bambuterol tablet. |
The mean cumulative excretion profiles of BAM and TER in healthy volunteers are shown in Fig. 6. It shows that the highest drug concentration in urine is detected within the first six hours; then, the concentration decreases gradually with the analysis time. Most of the drugs (BAM and TER) were eliminated within 10 hours of oral administration. The maximum elimination rates of BAM and TER were 180 and 530 ng ml−1 h−1, respectively, and this was observed within 0–5 h after their oral administration.
Fig. 6 Cumulative excretion of BAM and TER in urine samples of healthy volunteers after a single oral dose of BAMBIC® tablet. |
Parameter | PLSR chemometrics urine analysis | GC-MS reported method ref. 9 | ||
---|---|---|---|---|
BAM | TER | BAM | TER | |
a The values in parenthesis are the corresponding theoretical values of t and f at the 95% confidence level. | ||||
M.R.% | 98.92 | 99.97 | 98 | 100.2 |
S.D. | 0.982 | 1.67 | 0.4 | 0.3 |
N | 5 | 5 | 5 | 5 |
Standard error | 0.439 | 0.747 | ||
Student's t testa | 1.94 (2.31) | 0.303 (2.31) | ||
F testa | 0.166 (6.39) | 0.032 (6.39) |
Obviously, the ATR-FTIR method was straightforward; once the sample was collected (10 μL of urine), it could be directly placed over the ZnSe prism and scanned using the FTIR device to detect if the doping agent is present or not. Also, it has many advantages, which are as follows: nondestructive techniques, a wide range of samples with different masses can be analyzed, provides a spectrum of the analyte, which is considered as a fingerprint of each compound, only one minute is needed to get the results, and no chemical extraction or sample pre-treatment is required.
The disadvantages of the GC-MS method are as follows: it is a destructive technique, i.e., measurements cannot be repeated using the same sample, inorganic compounds cannot be detected, nonvolatile substances have to be chemically modified to be volatile, several sample pre-treatments have to be performed before measurements, not portable, chemical extraction steps of the analytes from the urine samples are required, analysis time is nearly 15 minutes for each sample, needs highly qualified operators, and the instruments are very expensive.
However, it gives complete data about the chemical composition of the sample and contamination that may be present in the sample.
The application of ATR-FTIR as a universal technique for the determination of doping drugs still needs much more investigation and application to all the doping drugs under all possible conditions. Our study can be considered as a preliminary study that opens the way for much more ATR-FTIR measurements for doping control.
It has many advantages in comparison with the traditional techniques by eliminating the errors during sample pre-treatment that can consume time and may lead to bad recoveries.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra10033d |
This journal is © The Royal Society of Chemistry 2020 |