Next Article in Journal
A Review on the Application of Machine Learning in Gamma Spectroscopy: Challenges and Opportunities
Next Article in Special Issue
Spectroscopic Analyses of Blue Pigments in the Manoscritto Parmense 3285 from the 14th Century
Previous Article in Journal / Special Issue
Interpreting the Microwave Spectra of Diatomic Molecules—Part II: Nuclear Quadrupole Coupling of One Nucleus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development and Validation of a Near Infra-Red (NIR) Hand-held Spectrophotometric Method Using PCA Approaches and Chemometric Tools: Application for Qualitative and Quantitative Determination of Tadalafil Marketed in Kinshasa—D.R. Congo

by
Jocelyn Mankulu Kakumba
1,2,3,*,
Patient Ciza Hamuli
1,3,
Merdie Mpemba Luyinama
1,
Freddy Mugisho Kasago
2,3,
Malachie Tembo Monyele
1,3,
Dadit Ive Kitenge
1,3,
Trésor Bayebila Menanzambi
1,3,
Trésor Kimbeni Malongo
1,2,3,
Didi Mana Kialengila
2,3 and
Jérémie Mbinze Kindenge
1,2,3
1
Laboratory of Drug Analysis and Health Products, Faculty of Pharmaceutical Sciences, University of Kinshasa, Kinshasa P.O. Box 212, Congo
2
Laboratory of Pharmaceutical Organic and Analytical Chemistry, Faculty of Pharmaceutical Sciences, University of Kinshasa, Kinshasa P.O. Box 212, Congo
3
Laboratory of Drug Analysis and Diets Control, Faculty of Pharmaceutical Sciences, University of Kinshasa, Kinshasa P.O. Box 212, Congo
*
Author to whom correspondence should be addressed.
Spectrosc. J. 2024, 2(3), 105-122; https://doi.org/10.3390/spectroscj2030007
Submission received: 3 March 2024 / Revised: 29 May 2024 / Accepted: 17 June 2024 / Published: 3 July 2024
(This article belongs to the Special Issue Feature Papers in Spectroscopy Journal)

Abstract

:
A hand-held NIR spectrophotometric method was developed, validated, and applied for the determination of tadalafil in tablets. The aim of our work was to develop analytical methods based on vibrational techniques using low-cost portable equipment. Based on different chemometric modeling, we attempted to validate the method, which gave encouraging results from the principal component analysis (PCA), DD-SIMCA, and PLS modeling. Following this, we optimized the method using an appropriate experiment plan. For validation, we used the total error approach with acceptance limits set at ±10% with a risk level of 5%. The method showed that it was possible to perform both qualitative and quantitative analysis of pharmaceutical products using low-cost portable NIR systems with chemometric tools. The developed approach enabled the completion of the first step in implementing an NIR method for quality control of tadalafil-based drugs in the DRC. Validation difficulties of the PLS method resulted from the lack of information about inter-day serial variations of spectral responses. It would be interesting to extend the study to a larger calibration interval in order to correct uncertainties that may result from the variability observed under different conditions and to verify robustness. These are the limitations of this work, but the results are nevertheless very encouraging.

1. Introduction

Globally, sub-standard and falsified medicines pose a serious risk to public health. Pharmaceutical falsification affects wealthy countries such as the United States, European countries, Australia, Japan, and others. This includes developing countries where falsified medicines are frequently distributed. According to the World Health Organization (WHO), any medication that has been fraudulently and purposefully mislabeled regarding its identity or source is considered falsified. Both branded and generic products are subject to falsification; examples of falsified products include those with the proper or incorrect active ingredients, no active ingredients, insufficient active ingredients, or phony packaging [1]. Phosphodiesterase (PDE) 5 inhibitors for erectile dysfunction, such as tadalafil (Cialis), sildenafil (Viagra), and vardenafil (Levitra), are among the most popular categories of personal-import medications. In Japan, reports from 2011 and 2012 indicated that tadalafil that had been fabricated posed health risks. Compared to 2005, the number of counterfeit medications that failed customs inspection surged one hundredfold in 2015, with PDE 5 inhibitors accounting for a large portion of such medications. A collaborative survey conducted in 2016 by four producers and distributors of erectile dysfunction medications revealed that 40% of the test samples were fabricated. There have also been reports of tadalafil and other phosphodiesterase type 5 inhibitors that were not real from Europe, the US, Canada, and Australia [2,3,4]. These medications are frequently made in unregulated street labs that do not adhere to appropriate manufacturing practices. The effectiveness and safety of falsified medicines are extremely unreliable since their source is unknown, and no quality control is carried out. Ingredients can range from harmless components to toxic and dangerous elements. The number of medications that require testing has increased due to growing concerns about falsified products, which emphasizes the need for quick and low-cost preliminary screening techniques. The problem of fake medications makes it abundantly evident that analytical methods are required in order to identify these false products and differentiate them from real medications. Several analytical methods have already been documented in the literature and can be broadly categorized into two categories: spectroscopic and chromatographic methods. Large quantities of these products arrive on the European market, which requires the development of simple and rapid methods to help customs officers detect them; this is the case with the near-infrared technique [1,2,5].
In their Roadmap for Supply Chain Security, the Asia Pacific Economic Cooperation stated that it is critical to create strategies and procedures for identifying and stopping the production of falsified products. However, the advancement of detection techniques has been sluggish, and no trustworthy techniques have surfaced. Pharmaceutical analysis has traditionally relied on high-performance liquid chromatography (HPLC) and thin-layer chromatography (TLC), which are capable of precisely identifying and measuring active pharmaceutical ingredients (API). Chromatographic techniques, particularly liquid chromatography coupled with mass spectrometry (LC-MS, HPLC-MS, UPLC-MS), are widely used because they produce consistent results. Nevertheless, there are a number of disadvantages to their use, such as the time-consuming nature of sample preparation and analysis, the high cost of the equipment, and the requirement for a skilled analyst. Furthermore, chromatography-based techniques may be harmful to the environment, primarily because solvents are used in the sample preparation and chromatographic separation processes [6,7,8,9]. Thus, several alternative methods are used to verify the identity of batches of starting materials in accordance with the Pharmaceutical Inspection Convention and Pharmaceutical Inspection Co-operation Scheme, such as portable spectroscopic analyzers. Spectroscopic techniques, including near-infrared (NIR), Fourier transform IR, nuclear magnetic resonance (NMR), and Raman spectroscopy, have been developed as quick, non-destructive methods of analysis. The main non-destructive technique used in the pharmaceutical industry that does not require sample preparation is NIR spectroscopy [10,11]. One such technique is gas chromatography-mass spectrometry (GC–MS). They do not cause destruction [12]. Since longer wavelengths of fluorescence are generally acknowledged to result in lower levels of external interference, attempts to alter the fluorescence structure (such as by creating probes with near-infrared emission) may be successful. Furthermore, a ratio metric pattern, which simultaneously monitors the fluorescence intensity of two wavelengths, can effectively eliminate external interference and improve detection sensitivity. However, vibrational spectroscopy is a preferred detection and quantification technology in the crop industry and can address some of the limitations of the aforementioned methods when combined with multivariate statistical analysis techniques (chemometrics). In particular, near-infrared (NIR) spectroscopy showed that it could accurately quantify a wide range of elements and compounds [2,13,14]. Infrared spectroscopy has shown potential as a tool for control in different environments [15,16,17,18], and it also demonstrated benefits for non-invasive and immediate detection, which allowed the technique to gain acceptance and a wide range of use in industries, including biomedicine and more. It also performs well when applied to other tasks such as adulteration detection, drug forgery identification, crop quality testing, etc. Miniaturized NIR spectrometers are becoming more and more common due to advancements in electronic technologies and microelectromechanical systems; portable models have been commercialized up to this point. Portable NIR spectrometers are very popular because of their flexibility, and their range of uses has grown from traditional laboratory settings to spot inspection [19,20,21,22]. Combined with chemometric tools, this technique can be more performant in evaluating products from any matrix, with possible online analytical applications. One of its advantages is the greater penetration of radiation in the sample; this makes it possible to analyze the thicker solid samples by reflectance to obtain more representative information [23,24,25,26].
Several measurement modes can be used depending on the sample shape, reflection, transmission, or transflection [13]. It is important to emphasize that the development of NIR methods for the quantification of solid dosage forms (tablets) presents major advantages in that they allow measurements to be carried out with little or no prior treatment of the samples. However, these methods present several challenges due to the influence of physical properties on the spectra. Due to the presence of a variety of excipients and dosage forms, these methods generally require the analyst to develop specific predictive models for each formulation. In order to solve this problem, this is a continuation of the work of P.H. Ciza et al. [13,14], except that here, for the first time, we tested a non-soluble compound and used HCl 0.1 N as solvent.
Therefore, this research aimed to develop and validate an analytical method using near-infrared for qualitative and quantitative determination of tadalafil marketed in Kinshasa and for its market surveillance in D.R. Congo with the use of chemometric tools combining PCA application and data pretreating. After development, the PLS method was validated and applied for quantitative application.

2. Materials and Methods

2.1. Chemicals

2.1.1. Pharmaceutical Substances and Brands

Tadalafil (Standard, 99.3%) was from ZM Technologies Limited, Liverpool, England. Magnesium stearate, Monohydrate lactose, cellulose microcrystalline, cross carmelose, and tadalafil samples were from Kim Pharma, Kinshasa, D.R. Congo. For reasons of confidentiality, the Tadalafil brand is coded as MCC, MMP, MXX, and MTT. The Table 1 groups all the analyzed brands giving in extension for each of them the applicant name, brand name, manufacturer, origin country, the authorization number, the registration date, validity period and expiration date.
  • Brand MCC
Dark yellow, biconvex, diamond-shaped, film-coated tablets in rectangular blister packs, each containing two tablets. Each blister is packaged in a yellow, green, and white cardboard box, with two blister packs in the latter. Each box also includes a leaflet containing composition, pharmacological properties, indication, contraindication, warnings, dosage, side effects, drug interactions, overdose, presentation, and storage.
Dosage: 20 mg
Manufacturer: Eli Lilly
Preemption date: June 2025
Batch number: 09968
b.
Brand MTT
Dark pink, rounded, lozenge-shaped, film-coated tablet contained in a rectangular blister pack. Each blister pack contains four tablets. Each blister is packaged in a blue-gold cardboard box containing a leaflet with composition, pharmacological properties, indication, contraindication, warning, dosage, side effects, drug interactions, overdose, presentation, and storage.
Dosage: 20 mg
Manufacturer: Lincoln Pharmaceuticals Ltd
Preemption date: January 2024
Batch number: NX 1001
c.
Brand MXX
Dosage: 20 mg
Manufacturer: Kim Pharma
Preemption date: July 2025
Batch number: XT-11
d.
Brand MMP
Pink, biconvex, diamond-shaped, film-coated tablet contained in a rectangular blister pack. Each blister pack contains two tablets. Each blister is packaged in a black, yellow cardboard box containing a leaflet with composition, pharmacological properties, indication, contraindication, warning, dosage, side effects, drug interactions, overdose, presentation, and storage.
Dosage: 10 mg
Manufacturer: Paloma
Preemption date: September 2025
Batch number: M-03

2.1.2. Solvents and Reagents

Analytical grade methanol used for sample preparation was obtained from Kim Pharma-Kinshasa, D.R. Congo. 37% analytical grade hydrochloric acid was used to dissolve components provided by VWR in Leuven, Belgium. The demineralized water was obtained from the faculty MilliQ Plus 185 (Waters, UK) system and purchased from Molshein in France and Acetone from Merck in Germany.

2.2. Near Infrared Instrumentations

Two NIR spectrophotometric devices were used in our research, both dispersive portable NIR spectrophotometer systems, namely NIR-S-G1, Innospectra Corp branded (Hsinchu, Taiwan), and NIR-M-T1, Innospectra Corp branded (Hsinchu, Taiwan). Measurements were conducted in two modes (reflection and transmission modes). The NIR devices were checked for their wavenumbers accuracy, and the X-axis was recalibrated as necessary as possible to ensure a lower difference (±1 cm−1).

2.3. Glassware and Small Equipment

An Eppendorf (5 and 50 mL), a volumetric pipette (1 mL), an analytical balance SV gram precision branded beaker, a bowl, gloves, and paper towels were used to facilitate experiments.

2.4. Preparation of Solutions

2.4.1. NIR Calibration and Validation Standards

Calibration and validation samples were prepared by dissolving reference tadalafil in methanol. Tadalafil standards were prepared in the presence of a sufficient quantity of mixing excipients to mimic tablets using gravimetric data as a reference. Three independent sample runs were performed for calibration and validation purposes with five concentration levels (5, 7.5, 10, 12.5, and 15 mg/mL tadalafil) in which the target concentration for the samples was 10 mg/mL.
Before analysis, tadalafil standard solutions containing excipients were shaken by hand for a maximum of 5 min, then left to settle for 10 min. The supernatant was finally filtered through a 4.5 μm filter. All calibration and validation solutions were prepared in triplicate for each concentration level and measured three times. For validation calculations, the three predictions were averaged per sample.
For qualitative model validation, solutions containing different brands of Tadalafil tablets were prepared at the target concentration of 10 mg/mL. In addition, a placebo solution with the excipient mixture was prepared at the target level.

2.4.2. Sample Preparation

Test solutions were prepared by simple dissolution or dilution of the unit dosage form with methanol to obtain a final solution of 1% w/v of the declared strength. The suspensions obtained from the tablet test sample were then shaken by hand for a maximum of 5 min, allowed to stand for 10 min, and passed through a 4.5 μm filter as for Tadalafil standards.

2.4.3. Sample Essays and Spectral Measurements

The target concentration of 10 mg/mL was obtained by dilution in HCl 0.1 N for each collected sample of Tadalafil. This was performed for every single sample, and each sample solution was filtered to obtain three measurements. For all the tests, predictions were averaged for n equals 3.
Tadalafil spectral measurements were performed using a NIR-M-T1-A device in transmission (sample in solution). Chemometric techniques and the spectral range were helpful for optimization and allowed the building of PLS models. Table 2 describes how the NIR spectra were processed, including all applied models to smooth, filter, and remove the noise.

3. Results and Discussion

3.1. Data Acquisition and Chemometric Pretreating

Chemometric calculations and analyses were performed in MATLAB (R2018a) (The Mathworks, Inc., Natick, MA, USA). For PLS models, the PLS Toolbox v.4.11 (EigenvectorResearch, Inc., Wenatchee, WA, USA) was used.
To reduce variability and improve chemical spectral characteristics, the raw NIR spectra were pre-processed.
Tadalafil spectra recorded with the NIR-S-G1 device in reflection (solid-state) and with NIR-M-T1 in transmission (dissolved in methanol) are shown in Figure 1 and Figure 2, respectively. Specific spectral features of Tadalafil are present in the spectral range where the absorption of O-H bonds in methanol is low (between 1500 and 1640 nm). For the spectrum of tadalafil in methanol, no difference is observable to the naked eye. However, chemometric techniques can detect very small differences in absorption. Figure 3 shows comparative spectra recorded using NIR-S-T1 in transmission during determinations of tadalafil brand samples, tadalafil standard, matrix (made of excipients mixture), and the blank (methanol).

3.2. Chemometric Data Analysis

3.2.1. PCA Models (Principal Component Analysis)

The PLS models were obtained using tadalafil-branded sample spectra and were compared once obtained from the two NIR Spectrometers. As can be seen in Figure 4, there is considerable variability between tadalafil tablet samples (green) and tadalafil standard samples (blue). When the results are compared between the different samples, there is a significant difference in terms of spectral signature.

3.2.2. DD-SIMCA Models

Due to the nature of NIR spectra (broad, highly overlapping bands), it is recommended to use a qualitative model before the PLS regression model.
In this study, Data-Driven Soft Independent Modelling Class Analogy (DD-SIMCA) was used to build qualitative identification models. DD-SIMCA is based on the construction of a principal component analysis (PCA) model of the target class, which, in our case, corresponds to the calibration spectra. In this way, the score distance (SD) and orthogonal distance (OD) can be calculated for each future spectrum, making it possible to determine, at a given confidence level, the acceptance zone for authentication of a specific brand. In our case, DD-SIMCA models were also used to distinguish X-1 Cialis from related (tadalafil-based) formulations but with different matrices.
A placebo solution (excipient and blank) was also projected onto the DD-SIMCA model to ensure that X-1’s excipients did not interfere with the qualitative models.
Sensitivity was assessed using the new set (tadalafil validation set) on the DD-SIMCA models, using the following formula:
S e n s i t i v i t y   ( % ) = T r u e   P o s i t i v e P o s i t i v e × 100
S p e c i f i c i t y = T r u e   N e g a t i v e N e g a t i v e × 100
As a first validation criterion, the specificity of the method must be demonstrated. It is generally good practice to build a qualitative model before the quantitative one. This ensures that only spectra similar to those used for calibration are projected onto the regression model.
As can be seen, visual spectral correlation cannot be a reliable method of selectivity. Due to the very small differences between spectra, pre-processing methods using Savitzky-Golay derivatives improve spectral characteristics and facilitate discrimination of different APIs in SIMCA models. In addition, chemometric methods were used to help detect spectral dissimilarities between these samples. A single-class classification model was developed using the DD-SIMCA approach. The model parameters selected are listed in Table 1, together with the associated sensitivity and specificity.
DD-SIMCA acceptance plots for calibration and validation data are shown in Figure 5 and Figure 6, while Figure 7 and Figure 8 show DD-SIMCA acceptance plots for X-1 and Cialis tablet data. The DD-SIMCA (Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9) models applied to NIR data enabled perfect recognition of X-1 samples and perfect discrimination of placebo and, in part, Cialis samples. This confirms its applicability to systematically reject non-X-1 Tadalafil samples prior to quantitative analysis.
Nevertheless, performance showed a sensitivity of 97.7% for tadalafil validation samples, a specificity of 100% for matrix, and a specificity of 33.3% for Cialis.

3.3. PLS Analysis

Several PLS models have been built using different preprocessing methods or combinations of them, taking into account different numbers of latent variables (see Figure 10 and Figure 11).
Selecting an appropriate number of latent variables avoids under- or over-fitting the model.
The PLS Toolbox model optimizer was used to preselect spectral ranges, preprocess, and a number of latent variables, using RMSEP as a quality criterion. The last few models selected were compared based on accuracy profiles reflecting the method’s current use.
The wide dispersion of the relative error of different concentration levels can be explained by the fact that the matrix used contained excipients that could pass into the methanolic solution (due to the higher excipient/API ratio). This has a random impact on the amount of tadalafil present in the solution after the filtration step. After being developed, the PLS model was tested for its predictive ability on spectral data taken under multi-source environmental conditions (change in temperature and relative humidity). Unfortunately, a bias was observed when the developed model was used to predict validation samples measured under different conditions (see Figure 6).
When looking at the score graphs, the validation data are considered to be outliers in relation to the model and, unfortunately, cannot be analyzed directly without correction (Figure 12).
The slope-bias correction SBC performed a linear-univariate model study between the predicted data and the standard data; then the linear model (taking into account the slope and y-intercept) was used to correct the predicted data using the following equation:
Y ( v a l ) C o r r = B i a s + S l o p e   Y ( v a l )
In some cases, this may be one of the alternative solutions for compensating for NIR spectral variations that may result from instrumental or environmental variability, for example, when samples and/or spectral measurements are significantly affected by changes in temperature, relative humidity, or instrumental variability. NIR-M-T1 portable instruments are still in the development phase and do not have waterproof shells to protect the system from temperature and/or humidity variations. As a result, particular attention needs to be paid to their use in tropical zones, where there are strong seasonal variations in temperature and relative humidity.
For SBC, the values of R2 (coefficient of determination), (the slope of the linear fit), and y-intercept, using predicted versus calibration data, are presented in the linear equations in Figure 8. The SBC corrected the bias observed when predicting standard tadalafil samples measured under different environmental conditions. The corrected predicted data enabled the model to be successfully validated in the tadalafil concentration range.
However, it would be interesting to test under different environmental conditions and check when convergence is reached, indicating the optimal condition to include in the overall modeling. Another approach could be global modeling. This approach consists of adding the new variability resulting from the new data taken under different environmental conditions to the calibration set.
Validation
The NIR predictive model was validated using the total error approach with acceptance limits of ±10% and a risk level of 5%.
All validation calculations were performed with E-noval 4.0b (PharmalexBelgium, Mont-saint-Guibert, Belgium).
Experiment plan
Validation standards are samples reconstituted in the matrix containing a known concentration and whose value is considered true by consensus.
Table 3 shows the number of validation standards per concentration level, the concentration levels considered, and the different series performed.

3.4. Validation Criteria Study

3.4.1. Trueness

Trueness expresses the closeness of agreement between the mean value obtained from a large series of test results and an accepted reference value. Trueness gives an indication of systematic errors.
As shown in Table 4, trueness is expressed in terms of absolute bias (mg/mL), relative bias (%), or recovery rate (%) for each concentration level of the validation standards.

3.4.2. Precision

Precision expresses the closeness of agreement between a series of measurements taken from multiple replicates of the same homogeneous sample under prescribed conditions. It provides information on random error and is assessed at two levels: repeatability and intermediate precision.
As shown in Table 5, precision (repeatability and intermediate precision) can be expressed in terms of standard deviation (SD) and coefficient of variation (CV).
CVs in % of Repeatability and Intermediate precision were obtained by dividing the standard deviation (SD) obtained by the corresponding mean of the introduced concentrations.
The precision of this method was assessed at two levels: repeatability and intermediate precision. The coefficient of variation (CV) was used as an expression of this precision.
For good method fidelity, the percentage CV must not exceed 2.000 at all concentration levels studied. This explains the good precision only at levels 4.0 and 5.0; repeatability and inter-day intermediate precision are very good considering only these two levels.

3.4.3. Accuracy

Accuracy expresses the closeness of agreement between the test result and the reference value accepted as such, also known as the “conventionally true value”. Accuracy takes into account the total error, i.e., the systematic error and the random error associated with the result. Consequently, accuracy is the sum of trueness and precision. It is estimated from the accuracy profile shown in Figure 13.
Acceptance limits have been set at ±10%, in line with the objective of the analytical procedure (USP).
The solid red line represents the relative bias, the dashed blue lines define the limits of the tolerance interval expected at the beta level, and the dashed black lines are the acceptance limits. The dots represent the relative error of the results and are plotted against their target concentrations.

3.4.4. Linearity

The linearity of an analytical method is its ability, within a certain assay interval, to obtain results directly proportional to the analyte concentration in the sample. A linear regression model (see Figure 14) was fitted to the results calculated as a function of the concentrations introduced in order to obtain the following equation:
Y = 0.3650 + 0.9655 X
where Y = results (mg/mL)
And X = concentrations introduced (mg/mL)
The coefficient of determination (r2) is 0.9890.
The correlation coefficient (r) is 0.9945.
Figure 14. Relationship between introduced concentrations and results (The solid black line depicts the linea-rity of the relationship between different concentrations as they progress from one level to the next, as well as the number of analysis repetitions per level, whereas the blue dots represent the various concentration levels, with the number of dots corresponding to the number of repetitions).
Figure 14. Relationship between introduced concentrations and results (The solid black line depicts the linea-rity of the relationship between different concentrations as they progress from one level to the next, as well as the number of analysis repetitions per level, whereas the blue dots represent the various concentration levels, with the number of dots corresponding to the number of repetitions).
Spectroscj 02 00007 g014
The solid black line is the identity line (Y = X). The limits represented by the dashed blue lines on the graph correspond to the linearity profile, the blue dots represent the various concentration levels, with the number of dots corresponding to the number of repetitions and dotted black line represents the dispersion limit from the lower to the upper concentration values (Figure 15), i.e., the “beta-expectation” tolerance limits expressed in absolute values.
In order to demonstrate the linearity of the method, the approach based on the expected tolerance interval at the beta level, expressed as an absolute value, can be used and is illustrated in the previous figure. As the figure shows, the linearity of the method is not exactly valid. This is due to the fact that R2 (0.989) is not close to 1. In this case, we say that there is no good agreement between the concentrations introduced and the concentrations calculated for certain levels.

3.5. Method Application to the Drug Determination

It is important to remember that routine application is the final phase in the life cycle of an analytical procedure. As the developed method is not validated, it is important to cover all environmental variability and revalidate the method before considering its routine use. Nevertheless, a preliminary test was carried out on samples of X-1 tablets (Tadalafil 20 mg). Test solutions were prepared by simply dissolving the unit dosage form with aqueous methanol solution to obtain a final solution of 1% w/v of the declared content. The resulting test tablet sample suspensions were then shaken by hand for a maximum of 5 min and then allowed to stand for 10 min before passing through a 4.5 μm filter as per the tadalafil standards. The analytical results are presented in the Table 6 showed below and including for each sample the batches (for those where they were available, predicted concentration, obtained concentration and the content in percentage).

4. Conclusions

Hand-held spectrophotometers have now become essential tools for quality control of pharmaceutical products, particularly for the field detection of substandard and/or falsified drugs. The aim of our work was to contribute to the development of analytical methods based on vibrational techniques and using low-cost portable equipment. Based on the results of our research, we can conclude that the developed method has shown that it is possible to perform both qualitative and quantitative analysis of pharmaceutical products using low-cost portable PIR systems combined with chemometric tools.
Difficulties in validating the PLS method were encountered due to the lack of information on variations in spectral responses between the different series prepared on different days. However, it would be interesting to extend this study to a larger number of calibration data in order to correct measurement uncertainties that may result from the variability observed under different environmental conditions, and to verify its robustness. These are the limitations of this work, but the results are nevertheless very encouraging.

Author Contributions

Conceptualization, J.M.K. (Jocelyn Mankulu Kakumba) and P.C.H.; methodology, J.M.K. (Jeremie Mbinze Kindenge), J.M.K. (Jocelyn Mankulu Kakumba) and M.M.L.; software, P.C.H. and M.M.L.; validation, M.M.L., D.I.K., M.T.M. and J.M.K. (Jeremie Mbinze Kindenge); formal analysis, J.M.K. (Jocelyn Mankulu Kakumba), D.M.K. and F.M.K.; investigation, M.M.L., D.I.K. and T.B.M.; resources, D.M.K., T.K.M. and J.M.K. (Jeremie Mbinze Kindenge); data curation, P.C.H.; writing—original draft preparation, J.M.K. (Jocelyn Mankulu Kakumba); writing—review and editing, J.M.K. (Jocelyn Mankulu Kakumba), D.M.K. and J.M.K. (Jeremie Mbinze Kindenge); visualization, P.C.H.; supervision, J.M.K. (Jeremie Mbinze Kindenge); project administration, J.M.K. (Jeremie Mbinze Kindenge). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the research data produced during this research are available to all the readers and can be found within this research article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Custers, D.; Krakowska, B.; De Beer, J.O.; Courselle, P.; Daszykowski, M.; Apers, S.; Deconinck, E. Chromatographic impurity fingerprinting of genuine and counterfeit Cialis® as a means to compare the discriminating ability of PDA and MS detection. Talanta 2016, 146, 540–548. [Google Scholar] [CrossRef] [PubMed]
  2. Sanada, T.; Yoshida, N.; Matsushita, R.; Kimura, K.; Tsuboi, H. Falsified tadalafil tablets distributed in Japan via the internet. Forensic Sci. Int. 2020, 307, 110143. [Google Scholar] [CrossRef]
  3. Tajanko, E.; Sein Anand, J.; Roszkowska-Kranc, K.; Kozłowska-Boszko, B.; Chodorowski, Z.; Korolkiewicz, R.P. Drug counterfeiting--The risk to the public health. Przegl. Lek. 2007, 64, 357–359. [Google Scholar] [PubMed]
  4. Fleshner, N.; Harvey, M.; Adomat, H.; Wood, C.; Eberding, A.; Hersey, K.; Guns, E. Evidence for contamination of herbal erectile dysfunction products with phosphodiesterase type 5 inhibitors. J. Urol. 2005, 174, 636–641; discussion 641; quiz 801. [Google Scholar] [CrossRef] [PubMed]
  5. Custers, D.; Vandemoortele, S.; Bothy, J.-L.; De Beer, J.O.; Courselle, P.; Apers, S.; Deconinck, E. Physical profiling and IR spectroscopy: Simple and effective methods to discriminate between genuine and counterfeit samples of Viagra® and Cialis®. Drug Test. Anal. 2016, 8, 378–387. [Google Scholar] [CrossRef] [PubMed]
  6. Caramês, E.T.D.S.; Piacentini, K.C.; Aparecida Almeida, N.; Lopes Pereira, V.; Azevedo Lima Pallone, J.; de Oliveira Rocha, L. Rapid assessment of enniatins in barley grains using near infrared spectroscopy and chemometric tools. Food Res. Int. 2022, 161, 111759. [Google Scholar] [CrossRef] [PubMed]
  7. Zhao, Z.; Wang, L.; Chen, J.; Zhang, N.; Zhou, W.; Song, Y. Altitudinal variation of dragon fruit metabolite profiles as revealed by UPLC-MS/MS-based widely targeted metabolomics analysis. BMC Plant Biol. 2024, 24, 344. [Google Scholar] [CrossRef] [PubMed]
  8. Sharma, D.; Wangoo, N.; Sharma, R.K. Ultrasensitive NIR fluorometric assay for inorganic pyrophosphatase detection via Cu2+-PPi interaction using bimetallic Au-Ag nanoclusters. Anal. Chim. Acta 2024, 1305, 342584. [Google Scholar] [CrossRef] [PubMed]
  9. Patyra, E.; Osiński, Z.; Kwiatek, K. The Identification and Quantification of 21 Antibacterial Substances by LC-MS/MS in Natural and Organic Liquid Fertilizer Samples. Molecules 2024, 29, 1644. [Google Scholar] [CrossRef]
  10. Sanada, T.; Yoshida, N.; Kimura, K.; Tsuboi, H. Detection Method of Falsified Medicines by Using a Low-Cost Raman Scattering Spectrometer Combined with Soft Independent Modeling of Class Analogy and Partial Least Squares Discriminant Analysis. Biol. Pharm. Bull. 2021, 44, 691–700. [Google Scholar] [CrossRef]
  11. Sanada, T.; Yoshida, N.; Kimura, K.; Tsuboi, H. Discrimination of Falsified Erectile Dysfunction Medicines by Use of an Ultra-Compact Raman Scattering Spectrometer. Pharmacy 2020, 9, 3. [Google Scholar] [CrossRef]
  12. Marín-San Román, S.; Fernández-Novales, J.; Cebrián-Tarancón, C.; Sánchez-Gómez, R.; Diago, M.P.; Garde-Cerdán, T. Application of near-infrared spectroscopy for the estimation of volatile compounds in Tempranillo Blanco grape berries during ripening. J. Sci. Food Agric. 2023, 103, 6317–6329. [Google Scholar] [CrossRef]
  13. Yao, S.; Ball, C.; Miyagusuku-Cruzado, G.; Giusti, M.M.; Aykas, D.P.; Rodriguez-Saona, L.E. A novel handheld FT-NIR spectroscopic approach for real-time screening of major cannabinoids content in hemp. Talanta 2022, 247, 123559. [Google Scholar] [CrossRef] [PubMed]
  14. Yuan, W.; Wan, C.; Zhang, J.; Li, Q.; Zhang, P.; Zheng, K.; Zhang, Q.; Ding, C. Near-infrared ratiometric fluorescent strategy for butyrylcholinesterase activity and its application in the detection of pesticide residue in food samples and biological imaging. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2023, 297, 122719. [Google Scholar] [CrossRef]
  15. Silva, F.S.R.; da Silva, Y.J.A.B.; Maia, A.J.; Biondi, C.M.; Araújo, P.R.M.; Barbosa, R.S.; Silva, C.M.C.A.C.; Luiz, T.C.S.; Araújo, A.F.V. Prediction of heavy metals in polluted mangrove soils in Brazil with the highest reported levels of mercury using near-infrared spectroscopy. Environ. Geochem. Health 2023, 45, 8337–8352. [Google Scholar] [CrossRef] [PubMed]
  16. Kachalkin, M.N.; Ryazanova, T.K.; Sokolova, I.V. Quantitative determination of ademetionine in tablets utilizing ATR-FTIR and partial least squares methods approaches. J. Pharm. Biomed. Anal. 2024, 241, 115991. [Google Scholar] [CrossRef] [PubMed]
  17. Maulidya, V.; Hasanah, A.N.; Rijai, L.; Muchtaridi, M. Quality Control and Authentication of Black Betel Leaf Extract (Piper acre Blume) from East Kalimantan as an Antimicrobial Agent Using a Combination of High-Performance Liquid Chromatography and Chemometric Fourier Transform Infrared. Molecules 2023, 28, 5666. [Google Scholar] [CrossRef]
  18. Attia, K.A.M.; El-Olemy, A.; Serag, A.; Abbas, A.E.F.; Eid, S.M. Environmentally sustainable DRS-FTIR probe assisted by chemometric tools for quality control analysis of cinnarizine and piracetam having diverged concentration ranges: Validation, greenness, and whiteness studies. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2023, 302, 123161. [Google Scholar] [CrossRef]
  19. Ren, S.; Jia, Y. Near-Infrared data classification at phone terminal based on the combination of PCA and CS-RBFSVC algorithms. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2023, 287, 122080. [Google Scholar] [CrossRef]
  20. Balage, J.M.; Amigo, J.M.; Antonelo, D.S.; Mazon, M.R.; da Luz e Silva, S. Shear force analysis by core location in Longissimus steaks from Nellore cattle using hyperspectral images—A feasibility study. Meat Sci. 2018, 143, 30–38. [Google Scholar] [CrossRef]
  21. Ripoll, G.; Lobón, S.; Joy, M. Use of visible and near infrared reflectance spectra to predict lipid peroxidation of light lamb meat and discriminate dam’s feeding systems. Meat Sci. 2018, 143, 24–29. [Google Scholar] [CrossRef] [PubMed]
  22. Qu, J.-H.; Liu, D.; Cheng, J.-H.; Sun, D.-W.; Ma, J.; Pu, H.; Zeng, X.-A. Applications of near-infrared spectroscopy in food safety evaluation and control: A review of recent research advances. Crit. Rev. Food Sci. Nutr. 2015, 55, 1939–1954. [Google Scholar] [CrossRef] [PubMed]
  23. Campos, M.I.; Debán, L.; Antolín, G.; Pardo, R. A quantitative on-line analysis of salt in cured ham by near-infrared spectroscopy and chemometrics. Meat Sci. 2023, 200, 109167. [Google Scholar] [CrossRef] [PubMed]
  24. Mokgalaka-Matlala, N.S.; Regnier, T.; Combrinck, S.; Kouekam, C.R.; Weiersbye, I.M. Chemometrics and vibrational spectroscopy as green tools for mine phytoremediation strategies. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2013, 100, 138–143. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, Y.; Li, M.; Li, L.; Ning, J.; Zhang, Z. Green analytical assay for the quality assessment of tea by using pocket-sized NIR spectrometer. Food Chem. 2021, 345, 128816. [Google Scholar] [CrossRef]
  26. Yu, F.-L.; Zhao, N.; Wu, Z.-S.; Huang, M.; Wang, D.; Zhang, Y.-B.; Hu, X.; Chen, X.-L.; Huang, L.-Q.; Pang, Y.-X. NIR Rapid Assessments of Blumea balsamifera (Ai-na-xiang) in China. Molecules 2017, 22, 1730. [Google Scholar] [CrossRef]
Figure 1. Tadalafil spectra were recorded with the NIR-S-G1 device in reflection (with pink curve representing the blank, capri blue curve representing excipients, blue color for standard and all the other colors as different tadalafil brands used for testing).
Figure 1. Tadalafil spectra were recorded with the NIR-S-G1 device in reflection (with pink curve representing the blank, capri blue curve representing excipients, blue color for standard and all the other colors as different tadalafil brands used for testing).
Spectroscj 02 00007 g001
Figure 2. Transmission spectra of NIR-M-T1 recorded Tadalafil (with pink curve representing the blank, capri blue curve representing excipients, blue color for standard and all the other colors as different tadalafil brands used for testing).
Figure 2. Transmission spectra of NIR-M-T1 recorded Tadalafil (with pink curve representing the blank, capri blue curve representing excipients, blue color for standard and all the other colors as different tadalafil brands used for testing).
Spectroscj 02 00007 g002
Figure 3. Spectra of tadalafil tablet samples (green curve), matrix (excipient mixture, in capri blue curve), blank (represented by the pink curve), and tadalafil SCR (blue curve) recorded NIR-M-T1 in transmission.
Figure 3. Spectra of tadalafil tablet samples (green curve), matrix (excipient mixture, in capri blue curve), blank (represented by the pink curve), and tadalafil SCR (blue curve) recorded NIR-M-T1 in transmission.
Spectroscj 02 00007 g003
Figure 4. Principal component analysis of spectral data obtained. Score diagram representing the variability between different samples.
Figure 4. Principal component analysis of spectral data obtained. Score diagram representing the variability between different samples.
Spectroscj 02 00007 g004
Figure 5. DD-SIMCA acceptance diagram for calibration samples. Samples below the threshold (green col-ored dots) are considered as standards of the target class and were used to build the model (The green dots reflect samples that are considered regular in reference to the calibration model and so authentic for the analysis performed. The green line defines the acceptance zone, within which the result can be deemed good and authentic, while the red line indicates the impassable boundary beyond which the product would be inauthentic, and the gap between the two lines is considered the warning zone).
Figure 5. DD-SIMCA acceptance diagram for calibration samples. Samples below the threshold (green col-ored dots) are considered as standards of the target class and were used to build the model (The green dots reflect samples that are considered regular in reference to the calibration model and so authentic for the analysis performed. The green line defines the acceptance zone, within which the result can be deemed good and authentic, while the red line indicates the impassable boundary beyond which the product would be inauthentic, and the gap between the two lines is considered the warning zone).
Spectroscj 02 00007 g005
Figure 6. DD-SIMCA acceptance diagram for validation samples. Samples below the threshold (green) are considered as belonging to the target class. Samples above the threshold (red) are considered out-liers.
Figure 6. DD-SIMCA acceptance diagram for validation samples. Samples below the threshold (green) are considered as belonging to the target class. Samples above the threshold (red) are considered out-liers.
Spectroscj 02 00007 g006
Figure 7. DD-SIMCA acceptance diagram for X-1 samples. All samples below the threshold (green color) are considered authentic.
Figure 7. DD-SIMCA acceptance diagram for X-1 samples. All samples below the threshold (green color) are considered authentic.
Spectroscj 02 00007 g007
Figure 8. DD-SIMCA acceptance chart for Cialis samples. Only a third of the samples below the threshold (green colored dots) are considered authentic.
Figure 8. DD-SIMCA acceptance chart for Cialis samples. Only a third of the samples below the threshold (green colored dots) are considered authentic.
Spectroscj 02 00007 g008
Figure 9. DD-SIMCA acceptance diagram for matrix and blank. All samples above the threshold (red color dots) are considered outliers as the green line represents the acceptability region limitation.
Figure 9. DD-SIMCA acceptance diagram for matrix and blank. All samples above the threshold (red color dots) are considered outliers as the green line represents the acceptability region limitation.
Spectroscj 02 00007 g009
Figure 10. Graph showing the characteristics of the PLS calibration model for NIR data (The blue dots repre-sent the number of analysis repetitions on each level of concentration as they progress from one to the next, as well as the dispersion of the responses obtained each time; the green line represents what would be a better characteristic of the calibration model and what would be desirable to obtain with the samples for acceptance of the analysis results; and the red line corresponds exactly to the inverse of what the green line provides as information according to this calibration model).
Figure 10. Graph showing the characteristics of the PLS calibration model for NIR data (The blue dots repre-sent the number of analysis repetitions on each level of concentration as they progress from one to the next, as well as the dispersion of the responses obtained each time; the green line represents what would be a better characteristic of the calibration model and what would be desirable to obtain with the samples for acceptance of the analysis results; and the red line corresponds exactly to the inverse of what the green line provides as information according to this calibration model).
Spectroscj 02 00007 g010
Figure 11. Graph showing the characteristics of the PLS calibration model developed for the prediction of validation data.
Figure 11. Graph showing the characteristics of the PLS calibration model developed for the prediction of validation data.
Spectroscj 02 00007 g011
Figure 12. Characteristics of the linear Y model.
Figure 12. Characteristics of the linear Y model.
Spectroscj 02 00007 g012
Figure 13. Accuracy profile.
Figure 13. Accuracy profile.
Spectroscj 02 00007 g013
Figure 15. Linearity profile graph of acceptability.
Figure 15. Linearity profile graph of acceptability.
Spectroscj 02 00007 g015
Table 1. Different tadalafil brand information.
Table 1. Different tadalafil brand information.
Applicant Product NameManufacturer Origin CountryAuthorization N°Registration DateValidity (Year)Preemption
KIM PHARMAX-1 20mg, capsule, bxe of 2.ETS KIM PHARMADR.CongoMS.1253/10/05/DGM/0108/2016May 20165May 2021
ELI LILLY NEDERLANDCialis20mgEli lilly NederlandNetherlandEU/1/02/237/006Nov 201210Nov 2022
Mon plaisir------
Tadalanique------
Table 2. Selected model parameters.
Table 2. Selected model parameters.
MetricNIR-M-T1
Spectral range1530–1642 nm
PreprocessingSG(1,2,15) + SNV + MC
PC3
A0.001
Sn (VAL)%97.7
Sp (Cialis)%33.3
Sp (placebo)%100.0
SG: Savitzky-Golay (derivative, polynomial order, window size). MC: Mean centering. SNV: Standard Normal Variable. Sn: Sensitivity. Sp: Specificity.
Table 3. Experiment plan.
Table 3. Experiment plan.
SeriesConcentration Level (mg/mL)No of Independent Replicates
11.03
12.03
13.03
14.03
15.03
21.03
22.03
23.03
24.03
25.03
31.03
32.03
33.03
34.03
35.03
Total number of observations: 45.
Table 4. The trueness.
Table 4. The trueness.
Concentration Level (mg/mL)Mean of Entered Concentrations (mg/mL)Mean of Results (mg/mL)Absolute Bias (mg/mL)Relative Bias (%)Recovery Rate (%)95% Confidence Interval Recoveries (%)
1.005.004.94−0.06−1.2498.76[92.48, 105.00]
2.007.507.690.192.51102.50[99.64, 105.40]
3.0010.0010.350.353.49103.50[101.70, 105.20]
4.0012.5012.550.050.42100.40[99.32, 101.50]
5.0015.0014.58−0.42−2.8397.17[95.57, 98.77]
Table 5. Repeatability and relative intermediate precision.
Table 5. Repeatability and relative intermediate precision.
Concentration Level (mg/mL)Mean of Entered Concentrations (mg/mL)RepeatabilityIntermediate Precision
SDCV%SDCV%
1.005.000.203.980.469.15
2.007.500.202.670.304.02
3.0010.000.232.290.232.29
4.0012.500.151.200.181.49
5.0015.000.120.800.352.36
Table 6. Analytical results.
Table 6. Analytical results.
SamplesBatchesPredicted Conc. (mg/mL)Corrected Conc. (mg/mL)Content%
1XT-1117.2410.57105.75
2XT-1117.4210.75107.54
3XT-1116.7510.10100.77
4 17.9811.31113.13
5 17.2410.57105.72
6 17.6510.98109.83
7 17.8411.17111.66
8 17.6711.00110.01
9 17.5410.87108.72
Target concentration: 10 mg/mL.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mankulu Kakumba, J.; Hamuli, P.C.; Luyinama, M.M.; Kasago, F.M.; Tembo Monyele, M.; Ive Kitenge, D.; Bayebila Menanzambi, T.; Kimbeni Malongo, T.; Mana Kialengila, D.; Mbinze Kindenge, J. Development and Validation of a Near Infra-Red (NIR) Hand-held Spectrophotometric Method Using PCA Approaches and Chemometric Tools: Application for Qualitative and Quantitative Determination of Tadalafil Marketed in Kinshasa—D.R. Congo. Spectrosc. J. 2024, 2, 105-122. https://doi.org/10.3390/spectroscj2030007

AMA Style

Mankulu Kakumba J, Hamuli PC, Luyinama MM, Kasago FM, Tembo Monyele M, Ive Kitenge D, Bayebila Menanzambi T, Kimbeni Malongo T, Mana Kialengila D, Mbinze Kindenge J. Development and Validation of a Near Infra-Red (NIR) Hand-held Spectrophotometric Method Using PCA Approaches and Chemometric Tools: Application for Qualitative and Quantitative Determination of Tadalafil Marketed in Kinshasa—D.R. Congo. Spectroscopy Journal. 2024; 2(3):105-122. https://doi.org/10.3390/spectroscj2030007

Chicago/Turabian Style

Mankulu Kakumba, Jocelyn, Patient Ciza Hamuli, Merdie Mpemba Luyinama, Freddy Mugisho Kasago, Malachie Tembo Monyele, Dadit Ive Kitenge, Trésor Bayebila Menanzambi, Trésor Kimbeni Malongo, Didi Mana Kialengila, and Jérémie Mbinze Kindenge. 2024. "Development and Validation of a Near Infra-Red (NIR) Hand-held Spectrophotometric Method Using PCA Approaches and Chemometric Tools: Application for Qualitative and Quantitative Determination of Tadalafil Marketed in Kinshasa—D.R. Congo" Spectroscopy Journal 2, no. 3: 105-122. https://doi.org/10.3390/spectroscj2030007

APA Style

Mankulu Kakumba, J., Hamuli, P. C., Luyinama, M. M., Kasago, F. M., Tembo Monyele, M., Ive Kitenge, D., Bayebila Menanzambi, T., Kimbeni Malongo, T., Mana Kialengila, D., & Mbinze Kindenge, J. (2024). Development and Validation of a Near Infra-Red (NIR) Hand-held Spectrophotometric Method Using PCA Approaches and Chemometric Tools: Application for Qualitative and Quantitative Determination of Tadalafil Marketed in Kinshasa—D.R. Congo. Spectroscopy Journal, 2(3), 105-122. https://doi.org/10.3390/spectroscj2030007

Article Metrics

Back to TopTop