Abstract
Keywords
Introduction
Oral contraceptives such as levonorgestrel – ethinyl estradiol (Levo-EE) are deemed essential medications by the World Health Organization (WHO). 1 Oral contraceptives are a common form of family planning and supplied through global health initiatives. 2 There are approximately 1.9 billion women of reproductive age globally, and an estimated 842 million are using contraceptive methods. 3 Clinicians commonly prescribe two low-dose active pharmaceutical indigents (API)s with one being a progesterone and the other an estrogen, which may offer other advantages in preventing venous events such as blood clots. 4 By the year 2030, the number of age bearing women will increase, especially in sub-Saharan Africa. 3
Poor-quality and falsified pharmaceuticals are an ongoing problem in supplying essential medications to those in resource limited areas including medications that are substandard, falsified, or degraded.5,6 These substandard medications risk prolonging illness, failed treatment, and possibly mortality. A previous meta-analysis found the prevalence of substandard medications reported in the literature continues to rise, however, estimating an economic impact was difficult and could range between $10 and $200 billion USD annually. 7 Rapid diagnostic tools for substandard medications are necessary to the response and key to quality compliance. 8
Diffuse reflectance spectroscopy (DRS) methods have been used for a number of years in pharmaceutical quality compliance of incoming raw materials (Igne and Ciurczak, 2021a), 9 and occasionally to monitor finished pharmaceutical products compliance.10,11 DRS methods use light interactions with tablets and suspensions to capture information pertaining to the physical and chemical nature of the samples. 12 In terms of a solid dosage tablet, this can relay information relating to the API and excipients. 13
In a global pharmaceutical supply chain, careful attention to detail should be paid to the influence of adverse sample holding and storage conditions. 14 For example, if a shipment of solid dosage tablets is stored in adverse conditions, there is a question if the product still can offer the same intended health benefit and if the product still is regarded as safe for consumption. The influence of heat and humidity are key factors in product stability. 15 Given the sensitivity of the near infrared region to moisture, there is potential for detecting tablets stored under adverse conditions.
In recent years, technological advances have positioned handheld spectrometers in the market that are compact and come with a much lower cost-point than traditional benchtop spectrometers. 16 However, these spectrometers are limited in the spectral range of data collection, the optical bandwidth, and typically come with the unknown influence of changes to the operating environment, such as increased temperature and humidity. Several detailed recent reviews can be found that describe the current benefits, conflicts, and unknowns concerning applied usage of these instruments in various fields such as agriculture, textiles, and pharmaceuticals.17–19 Additionally, the use of qualitative analysis for falsified pharmaceutical tablets was undertaken with handheld Raman and near infrared (NIR) spectrometers. 20
In resource limited areas where it may be restrictive to screen essential medications for quality compliance with a benchtop DRS that requires a higher initial investment in funding, training, and software, it may be advantageous to implement handheld spectrometers for compliance screenings. Handheld and portable spectroscopy tools are not intended to replace traditional laboratory-based chemical testing and their confirmatory power, but are intended to compliment traditional methods when resources are limited. 21 In this instance there is concern regards brand differentiation and determining whether environmentally stressed samples can be detected and flagged with DRS. A handheld spectrometer is used (in comparison to a benchtop spectrometer) to sample Levo-EE tablets containing two low-dose active ingredients, to determine if (a) brand discrimination of solid-dosage tablets is possible and (b) to determine if environmentally stressed tablets are rejected from a single-class screening method for one of the Levo-EE brands.
Materials and methods
Oral dosage tablets
Three brands of oral dosage tablets containing levonorgestrel and ethinyl estradiol were obtained from manufacturing sources for the initial brand qualitative discrimination. Oral contraceptives were obtained through routine screening from a donation based global health supply chain. Manufacturing sources are blinded for confidentiality purposes and are referred to as brands “A”, “B”, and “C” from here onwards. Products are approximately 75 mg per tablet, containing 30 µg of ethinyl estradiol and 150 µg of levonorgestrel. For the brand discrimination portion, 233 tablets of brand A were obtained from 8 lots, 120 tablets of brand B were obtained from 4 lots, additionally, 10 samples from one lot of brand A under a different label were collected, 108 tablets from 4 lots of brand C. Tablets come packaged with placebos that are typically taken after 21 days of taking the tablets containing the API, and are commonly either ferrous fumarate or lactose. Here, nine ferrous fumarate placebo tablets from brand A, B, C, and 3 lactose placebo tablets from brand A were also sampled and included as negative controls.
Levonorgestrel – ethinyl estradiol tablet samples scanned by a handheld spectrometer (900 – 1700 nm) and a benchtop spectrometer (350 – 2500 nm). Samples used in construction of the brand specific qualitative model for brand A were from tablets stored at ambient conditions.
Levo = levonorgestrel, EE = ethinyl estradiol, and RH = relative humidity.
Brand A tablets held under adverse environmental conditions were stored for 7 years and 1 month, Brand B tablets held under adverse environmental conditions were stored for 5 years and 6 months.
Spectrometers and data collection
The two NIR instruments used for spectral acquisition are shown in Figure 1 with the handheld spectrometer, a Tellspec Food Sensor (Tellspec Inc., Toronto, ON, Canada), shown in Figure 1(a) and the ASD LabSpec 5000 with Muglight® assembly (Malvern Panalytical, Malvern, UK) shown in Figure 1(b). The portable handheld spectrometer operates in the short-wave infrared range (900 – 1700 nm) collecting data at a variable 2 – 4 nm for 256 points was used for scanning tablets. The Tellspec spectrometer consists of the NIR-S-G1 spectrometer (Innospectra Corp., Hsinchu, Taiwan), and operates through the free and downloadable version of the DC&M 2.0 software (Tellspec, Inc., Montreal, Canada), which was installed on to a mobile device for data collection, after being paired through Bluetooth® with data stored in the app. For sample collection, a 12 × 35 mm borosilicate glass vial (FisherbrandTM 51, Fisher Scientific, Hampton, NH, USA) was cut in half using a motorized saw with a 19 mm or 38 mm rotary diamond blades (Dremel 4000, Racine, MI, USA). After cleaning the cut vial, it was affixed to the handheld spectrometer. Inside the vial, a 1.6 mm thick piece of Teflon sheeting (McMaster Carr, Elmhurst, IL, USA) was cut to the inside diameter of the vial and inserted. A #1 cork borer was used to cut a small hole (4 mm) in the center of the Teflon before inserting and was positioned directly over the scanning window of the handheld spectrometer. Background scans were collected with the empty vial and Teflon insert, and a second Teflon insert inserted into the vial to cover the small opening over the scanning window. For sample scans, the top piece of Teflon was removed, and tablets were positioned on top of the opening and scanning window. Tablets covered the opening in the Teflon completely. Triplicate scans were collected, and tablets were rotated approximately 45 degrees between scans. After scanning was complete, CSV files were exported as raw reflectance spectra. The background scans were used to adjust the raw reflectance and convert to absorbance Log10(1/R) for further analysis, where R = I/Io with I for the sample data and Io for the reference. NIR spectra of standards for levonorgestrel, ethinyl estradiol, and lactose (purity 97.6% or higher) (Sigma Aldrich, St Louis, MO, USA) were acquired by placing the powders in individual 12 mm × 35 mm borosilicate glass vials (Fisher Scientific, Waltham, MA, USA). A background was taken by another of the previously mentioned vials and cutting a piece of 1.6 mm thick Teflon to the inside diameter and placing in the bottom of the vial, then scanning in triplicate against the sample window. (a) Tablet scanning adapter and setup for use with the handheld spectrometer (900 – 1700 nm). (b) Tablet scanning setup for the benchtop spectrometer using the Muglight® attachment, sample tray adapter, and Spectralon cover for tablet data collection.
The Labspec 5000 was used for scanning the tablets from 350 to 2500 nm collecting data every 1 nm for 2150 points. The Muglight® accessory was used to collect data with a small sample tray adapter. First, a background was collected as the baseline using a 25.4 mm Spectralon puck inserted into the sample tray. Tablets were centered on the sample tray adapter and covered with a 89 mm Spectralon disk (Malvern Panalytical, Malvern, UK). Scans were collected in triplicate with tablets rotated approximately 45 degrees between scans. Each spectra was an average of 50 scans, resulting in a final spectra averaging 150 scans. Spectra were saved as Log10(1/R) and exported as an ASCII file. The three secondary reference standards for levonorgestrel, ethinyl estradiol, and lactose were scanned in a similar fashion as described with the handheld, but with the previously mentioned Muglight® attachment for the benchtop instrument. Here, another piece of 1.6 mm thick Teflon was cut to the inside diameter of the sample tray adapter with a hole cut in the center, the same diameter as the sample vial. References vials (previously described for the handheld) were also scanned in this manner for the benchtop spectrometer.
Multivariate data analysis
After converting raw reflectance into Log10(1/R), data were imported into R-studio (2022.07.01) (Boston, MA, USA), where all analyses were conducted. A mean-centered singular value decomposition (SVD) based principal component analysis (PCA) with a leave-one-out cross validation was conducted on samples from both the benchtop spectrometer dataset and separately with the handheld spectrometer’s dataset using an R-Studio Shiny App. 22 A single product specific multivariate classification model was developed for “brand A”. A classification model built around Mahalanobis distances (M-dist.) based on principal component (PC) scores (mean-centered) and residuals was established for brand A as the reference dataset. In brief, reference samples were used to build a spectral library, and samples of other brands, placebos, and tablets stored under environmental stress were projected against the model to determine if the benchtop and portable spectrometers could differentiate negative controls from the reference, and verify the model with positive controls of new tablets from brand A. Reference and positive/negative controls were split amongst lots to approximately 70 % for references and 30 % for positive controls. More information on the M-dist. classification approach and establishing the thresholding value for classification at a 95% confidence interval can be found in Gemperline and Boyer and Kaale et al.23,24
Results and discussion
Brand comparison
The mean spectra by brand for tablets stored at ambient conditions for the benchtop is seen in Figure 2(a) and handheld spectrometer is shown in Figure 2(b). The spectra for brands A, B, and C are similar in pattern and intensity, whereas the ferrous fumarate and lactose placebos are noticeably different with spectra obtained from the benchtop spectrometer. Spectral differences are more pronounced with the handheld spectrometer’s mean spectra showing that brands B and C are different around the 1000 nm region from brand A. The ferrous fumarate and lactose placebos differ from the spectra of brands A, B, and C by a substantial amount, and are easily differentiated, visibly. The spectra for reference standards are shown for the benchtop spectrometer in Figure 2(c) and the handheld spectrometer in Figure 2(d). The variable 2 – 4 nm bandwidth of the handheld spectrometer is noticeably different from the benchtop’s 1 nm bandwidth, as the spectral peaks are broader for the handheld than the benchtop. Each of the three reference standards have definitive spectral peaks that are also noticed in the benchtop spectra. Lactose peaks at 1543 nm and 1593 nm are noticeable, while levonorgestrel peaks found at 1185 nm, 1386 nm, and 1558 nm. While ethinyl estradiol peaks can be found at 1420 nm, 1468 nm, and 1543 nm are observed with both spectrometers. However, there are multiple key spectral features for each of the three standards found beyond the 1700 nm limit of the handheld spectrometer that are captured by the benchtop spectrometer. Mean spectra of tablets stored at ambient conditions by brand/placebo for both the benchtop spectrometer (a) and the handheld spectrometer (b). Reference spectra for lactose (a common excipient), levenorgesterol (levo), and ethinyl-estradiol (EE) are shown for the benchtop spectrometer (c), and the handheld spectrometer (d).
Principal component analysis was applied to the dataset containing the three brands of Levo-EE stored at ambient conditions. The PCA scores plot can be seen in Figure 3 with the benchtop spectrometer shown in the top row, Figure 3(a) and (b), and handheld spectrometer shown in the bottom row, Figure 3(c) and (d). The PCA scores plots represent similarities and dissimilarities between samples of Levo-EE. Brand A was established as the reference samples with positive controls also from brand A. Negative controls 1 (NC1) (brand B) and 2 (NC2) (brand C) are two differently labeled brands but are understood to be produced by the same manufacturer and packaged with two different labels. Considering this, it explains why the NC1 and NC2 samples overlap each other on the PCA scores plots for both spectrometers. Also, the PC samples overlap the reference (Ref) samples showing consistency and potential for single class model sensitivity. Score values from the benchtop spectrometer show an approximately normal distribution of all samples and cluster overlapping when comparing PC1 (72.23%) and PC2 (16.75%). Plotting lower PCs shows a clear separation between Ref/PC and the NC1/NC2 samples, which can be seen by the plot in Figure 3(b). The handheld spectrometer resulted in PCA score values being separated with PC1 (71.5%) and PC2 (24.2%). Whereas the benchtop spectrometer relates more influence from the visible light range, with tablet color likely having a greater impact on PC1. Scores plots from principal component analyses of data collected with the benchtop spectrometer (a and b) and the handheld spectrometer (c and d) for brands A, B, and C Benchtop spectrometer results shown in the top row and handheld spectrometer shown in the bottom row. Ref = reference, PC = positive control, NC1/NC2 = negative controls 1 and 2.
Comparing PC loading vectors gives an indication to the variable-to-variable implications on the PCA scores. The benchtop spectrometer’s PCA loading vectors for PCs 1 – 4 (explaining 98.7% of the PCA model variance) are shown in Figure 4(a) with PC1 shown as an insert, due to scaler differences in the Loading vectors for principal components 1 – 4 for principal component analysis of benchtop spectrometer with the PC1 loading vector inserted (a) and the handheld spectrometer (b) for PCs 1 - 3. Samples included in PCA are for the brand discrimination section consisting of in-date and ambient stored brands A, B, and C only. Proposed spectral peaks from levonorgestrel (Levo), ethinyl-estradiol (EE), and lactose secondary reference standards collected across the Labspec range (350 – 2500 nm).
Environmental stress and brand classification
Once tablets were removed from the environmental chambers, they were visually inspected while still inside of the blister packs. Samples for both brands stored at 30°C and 40°C at 0% RH were visually similar, while brand B stored at 50 °C–54°C and 75% RH appeared dark red, and once the blister packs opened, tablets turned to a paste upon removing tablets from the package, having lost physical stability. These tablets were unable to be scanned with either spectrometer. Brand A stored under the same conditions had a slight yellowish hue, but tablets maintained their physical integrity and were able to be scanned by benchtop and handheld spectrometers. The mean spectra from tablets that were stored under adverse environmental conditions can be seen in Figure 5 (under desiccation) and Figure 6 (75 % RH). Spectra shown here are not preprocessed and we can see that there is a scaler difference for brand A and B with samples stored at 50 °C–54°C resulting in a lower Absorbance Log10(1/R) spectra for both the benchtop and handheld spectra. Samples stored at 30°C and 40°C showed more similarities between mean spectra. The mean spectra for brand B at 30 and 40°C are shown, but spectra for 50 °C–54°C are missing from brand B, because the tablets absorbed enough moisture to deteriorate the tablets in a way that rendered them unusable. Brand B shows that there is a greater change in spectra at and above the peak shown for moisture peak at approximately 1940 nm, however, this change is less pronounced for brand A. Brand A samples stored at 75% RH also shows a greater difference between tablets stored at 50 °C–54°C than 30 or 40°C, similar to the temperature correlations for desiccated samples. Samples stored at 50 °C–54°C also resulted in lower spectral peak resolutions with broader peaks noticed. This would also suggest that tablets stored in the range of 30 °C–40°C do exhibit noticeable changes in their spectra, but tablets stored at 50 °C–54°C have a noticeably larger change in spectral pattern. With these tablets stored in environmental chambers, we are observing spectral changes under long-term storage with adverse heat and humidity conditions, representing a possible worst-case scenario. Mean spectra for benchtop spectrometer (a and b) and handheld spectrometer (c and d) for tablets stored under desiccation at 30, 40, or 50°C–54°C. Brand A shown in a and c, while brand B is shown in (b and d). Mean spectra for benchtop spectrometer (a and b) and handheld spectrometer (c and d) for tablets stored at 75 % relative humidity and at temperatures of 30, 40, or 50 °C–54°C. Brand A shown in a and c, while brand B is shown in (b and d).

Moisture influence
Solid dosage tablets can have a shelf life of several years, and therefore preventing moisture pickup during storage is a primary concern during manufacturing and packaging. Previously, our group obtained blister packs of swollen and cracked oral contraceptive tablets provided through USAID’s program. Tablets were scanned with the benchtop spectrometer and results are shown in Figure S1. These results were previously unpublished but offer insight into the moisture retention in the tablets and changes to the spectra. Figure S1(a) shows that the visibly noticeable cracked tablets inside the blister packs, while Figure S1(b) shows the number of tablets cracked or damaged for each of the four lots, where the vast majority of cracked tablets occurred within cells that showed no evidence of indentation. The mean raw spectra of damaged and non-damaged tablets can be seen in Figure S1(c). The water overtone peak at approximately 1940 nm, which is emphasized by performing a 1st derivative Savitzky-Golay 1st preprocessing step shown in Figure S1(d).
After removing the tablets for brands, A and B from storage in the environmental chambers and scanning, spectra were also derivatized (1st derivative Savitzky-Golay of Absorbance Log10(1/R)) and plotted to observe the spectral changes in each brand by changes in storage conditions, with brand A shown in Figure 7 and brand B shown in Figure 8. For each figure, the benchtop spectrometer’s resulting spectra is shown on the left (a) (1000 – 2500 nm), while the handheld spectrometer’s results are shown on the right (b) (900 – 1700 nm). Similar to previous work (Figure S1), derivatized spectra show evidence of moisture at different temperatures/humidity conditions (relative to ambient conditions) in the regions of 1930-1940 nm for both brands and spectrometers. Recently, Patel, et al., (2023) showed that moisture in tablets could be quantified using a handheld near-infrared spectrometer (1596 – 2396 nm).
25
Mean spectra of samples stored under 75 % relative humidity for brand A tablets. Boxes highlight the water absorbance ranges. (a) The benchtop spectra with a 1st derivative Savitzky-Golay transform (1st order polynomial, and a window size of 5). (b) The handheld spectra with a 1st derivative Savitzky-Golay transform (1st order polynomial, and a window size of 3). Mean spectra of samples stored under 75 % relative humidity for brand B tablets. Boxes highlight the water absorbance ranges. (a) The benchtop spectra with a 1st derivative Savitzky-Golay transform (1st order polynomial, and a window size of 5). (b) The handheld spectra with 1st derivative Savitzky-Golay transform (1st order polynomial, and a window size of 3).

Mahalanobis distance results for a qualitative discriminatory method for product A using the benchtop and handheld spectrometers.
A or B refer to the brand of levonorgestrel and ethinyl estradiol.
All the benchtop reference samples, and 159 of the 160 handheld samples resulted in a passing classification for the brand A specific method, while all of the PCs (
Conclusions
In low- and middle-income countries, sub-standard and counterfeit pharmaceuticals are an ongoing concern. Portable and low-cost screening methods for essential medicines have the potential for impactful quality compliance systems. Here, a handheld spectrometer collecting data in the short-wave infra-red spectral range was able to qualitatively differentiate between brands of Levo-EE. When additional Levo-EE samples stored under environmental stress for an extended period of time were compared against the model for brand “A”, the stressed tablets were easily identified as different from the reference model for tablets stored under ambient conditions. The spectra collected from the environmentally stressed tablets had noticeable increases at known near infrared water absorbance peaks at approximately 1450 nm and 1940 nm. Not only were spectral changes noted on the benchtop spectrometer, but was also picked up by the handheld spectrometer, suggesting that the handheld spectrometer is not only sensitive enough to detect brand differences in tablets containing two low-dose active ingredients, but also sensitive enough to detect tablets that have been under environmental stress, potentially reducing the efficaciousness of the medicine.
Supplemental Material
Supplemental Material - Qualitative discrimination of low-dose oral contraceptives and identification of environmentally stressed tablets: Comparing a handheld near-infrared spectrometer (900 – 1700 nm) to a benchtop spectrometer (350 – 2500 nm)
Supplemental Material for Qualitative discrimination of low-dose oral contraceptives and identification of environmentally stressed tablets: Comparing a handheld near-infrared spectrometer (900 – 1700 nm) to a benchtop spectrometer (350 – 2500 nm) by Matthew Eady, Noah Peters and David Jenkins in Journal of Near Infrared Spectroscopy.
Footnotes
Acknowledgements
Declaration of conflicting interests
Funding
Supplemental Material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
