- Open Access
Green adherent degradation kinetics study of Nirmatrelvir, an oral anti-COVID-19: characterization of degradation products using LC–MS with insilico toxicity profile
BMC Chemistry volume 17, Article number: 23 (2023)
The SARS-CoV-2 virus sets up a global catastrophe, and countries all around the world made significant efforts to halt the spread. Nirmatrelvir (NMV) was lately approved by the FDA as a safe and well-tolerated oral direct-acting antiviral medication for SARS-CoV-2 virus infection. Therefore, a fast completely validated stability indicating method was established-for the first time- for NMV determination. The study used NaOH, HCl, neutral, H2O2, and sunlight to test NMV stability under various stress conditions followed by kinetics degradation investigation and derivation of Arrhenius plot. The analysis was performed using Agilent Zorbax Eclipse-C18 column (5 µm, 4.6 × 250 mm) with a mobile phase consisting of acetonitrile: 50 mM ammonium acetate, pH = 5 (50:50, v/v, respectively) at a flow rate of 1.0 mL/min with 5 min run time. Diode array detector (DAD) was set at 225 nm to quantify NMV at the concentration range of 5–500 µg/mL with LOD and LOQ of 0.6 and 2 µg/mL, respectively. Method’s greenness was assessed using different metrics including Analytical Eco-Scale, Greenness Assessment Procedure Index, GAPI, and Analytical Greenness, AGREE. A thorough study of stress stability revealed that NMV was more susceptible to alkaline hydrolysis compared with acid hydrolysis. In contrast, it was found that NMV remained stable when subjected to oxidative, neutral, and sun-induced degradation conditions. Moreover, acid and alkali-induced hydrolysis were found to follow pseudo first order kinetics. Consequently, the half lifetime of the studied degradation conditions at room temperature were calculated using the Arrhenius plot. The mechanism of the degradation pathways under stress circumstances was proposed using LC–MS-UV. Toxicities of the proposed degradation products were assessed using ProTox-II, along with the parent medication NMV, and were shown to be hardly hazardous.
Coronavirus disease 2019 (Covid-19) with its severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued to be a concern to world health. Elderly people, smokers, or those with underlying clinical diseases including cardiovascular diseases, diabetes, obesity, or cancer, are at a higher risk of severe Covid-19 and its consequences compared with normal people. It was reported that people with co-existing diseases are almost two times more probable to develop severe Covid-19 infections with a mortality rate five times more than those without such underlying diseases [1,2,3,4,5,6].
There is a need for safe and effective oral Covid-19 therapies that can accelerate the recovery and prevent the progression of the infection to more serious symptoms including hospitalisation, and probably death. Such therapies would aid the health care system to manipulate the problems of overcrowded hospitals with inadequate facilities. Monoclonal antibodies were widely used for non-hospitalized individuals with mild-to-moderate Covid-19 infections. However, they have several drawbacks, such as the requirement for clinical investigation and the possibility of diminished efficacy against developing SARS-CoV-2 variants .
Nirmatrelvir (NMV) is an antiviral drug that targets the SARS-CoV-2 enzyme, 3-chymotrypsin–like cysteine protease. NMV effectively inhibits the enzymatic activity with subsequent inhibition of virus multiplication across a wide range of coronaviruses. Being a CYP3A4 substrate, it was found that coadministration of 300 mg NMV with a low dose (100 mg) of ritonavir, a CYP3A4 inhibitor, twice daily enhances NMV pharmacokinetics [8, 9]. In December 2021, Pfizer’s Paxlovid® (nirmatrelvir tablets and ritonavir tablets, co-packaged for oral use) has received an emergency use authorization (EUA) from the US Food and Drug Administration for the treatment of mild-to-moderate COVID-19 in adults and paediatric patients with positive results of direct SARS-CoV-2 testing and who are at a high risk of progressing to severe COVID-19, including hospitalisation or death .
Drug stability affects product safety and efficacy since degradation products (DPs) may cause a loss of potency and result in potential toxic effect. Accordingly, chemical and physical stability are important to ensure the quality and safety of active ingredients .
To the best of our knowledge, only one fully validated method was published for analysing NMV in human plasma along with ritonavir using LC–MS/MS . Another article has been recently published for the characterization of NMV’s DPs using LC-UV-HRMSn . Although this published work  provided detailed postulation of degradation pathways, it did not cover any quantitative aspects concerning method validation, determination of the intact drug in presence of its DPs, or studying the degradation kinetics. Considering the high expenses of LC–MS/MS instrumentation besides its need to high expertise in its operation, and its unavailability in many analytical labs, HPLC–DAD is considered a simple alternative for quantitative applications. To the best of our knowledge, neither a validated HPLC–DAD method for the determination of NMV nor information about the degradation kinetics was found in the literature so far. This work aims at developing and validating a simple stability-indicating HPLC–DAD method-for the first time-for the determination of NMV in presence of its DPs generated under different stress conditions. The application of this method was extended to studying the degradation kinetics of NMV which is extremely important in determining the order of the reaction along with the half lifetime using Arrhenius plot. Moreover, it is a crucial target to study the safety profile of drugs and their DPs, accordingly insilico toxicity prediction software was used in this work with NMV and its DPs which adds to the novelty of the proposed study. In addition, the eco-friendlessness of the validated HPLC–DAD method was assessed using different greenness tools.
Materials and reagents
Nirmatrelvir (NMV) was kindly supplied from KU Leuven Research & Development, Belgium with approximate purity of 94%. El-Nasr Chemical Industry Company in Egypt provided Analytical quality grade hydrochloric acid (HCl) and sodium hydroxide (NaOH). Hydrogen peroxide (H2O2), 50%, was also obtained from El-Nasr Chemical Industry Company.
Both HPLC- grade acetonitrile (ACN) and methanol (MeOH) were purchased from Baker's, Ireland. Water was deionized and filtered using Ultrapure PTFE syringe filters, 0.2 µm pore size.
The HPLC–DAD system consisted of Agilent 1200 series (auto-injector, quaternary pump, vacuum degasser and diode array with multiple wavelength detector G1315 C/D and G1365 C/D connected to a computer loaded with Agilent ChemStation Software (Agilent Technologies, Santa Clara, CA, USA).
Shimadzu® UFLC series (Shimadzu Corporation, Kyoto, Japan) was used with its LC–MS-UV system hyphenated with two detectors: UV detector SPD-20A and MS detector LCMS-2020 with mass range of m/z 10 to 2000. The detector was supplied with single quadrupole system and two ion sources ESI and APCI (Shimadzu Corporation, Kyoto, Japan). The scan rate was up to 15,000 amu/sec. A Shim-pack XR-ODS II® C18 analytical column (100 × 3 mm, 2 µm) was used as the stationary phase (Shimadzu, Japan).
The separation was performed on an Agilent Zorbax Eclipse-C18 analytical column (250 × 4.6 mm, 5 µm). The mobile phase system consisted of solvent A (50 mM ammonium acetate, pH = 5) and solvent B (ACN) in an isocratic elution made of 50:50, v/v, respectively. Detection wavelength was set to 225 nm. A flow rate of 1 mL/min and 20-µL injection volume were used. The entire run time was 5 min.
A Shim-pack XR-ODS II® C18 analytical column (100 × 3 mm, 2 m), Shimadzu, Japan, was employed for the separation. The mobile phase used was the same as that for HPLC–DAD. Detection was performed using an electrospray ionisation (ESI) interface set to the selected ion monitoring (SIM) mode with positive and negative polarity. The flow rate was kept at 0.4 mL/min. Injection volumes of 10 µL were used all over the study. The run time was 7 min.
Preparation of stock standard solutions and construction of calibration graphs
Stock solutions of NMV 1000 µg/mL were prepared in HPLC-grade methanol. The solutions were kept refrigerated at 0 °C. Dilution of the stock solutions using deionized water was performed to get working standards within the concentration range 5–500 µg/mL. Each concentration was injected in triplicate into the HPLC–DAD system under the optimized conditions using 20 µL injection volume. The peak area of NMV recorded at its Rt of 4.01 ± 0.01 min was related to each concentration to get the regression equation and the calibration graph.
Analytical method validation of NMV was appraised for accuracy, precision, linearity, specificity, limit of quantitation and robustness in agreement with ICH guidelines .
Forced degradation and stability-indicating study
Forced degradation studies were done using separate volumes of 1 mL of NMV stock solution (1000 µg/mL). Following degradation, the solutions were then completed to volume in a series of 10-mL volumetric flask with deionized water to reach a final concentration of 100 µg/mL.
Acidic and alkaline hydrolysis
NMV solutions were treated with 1 mL volumes of 1 M HCl or 0.1 M NaOH, for acidic and alkaline hydrolysis, respectively. Solutions were then heated in a thermostated water bath kept at 90 ℃ for 1 h for acidic hydrolysis or at 60 °C for 30 min for alkaline hydrolysis. After the specified time intervals, the solutions were cooled then neutralized to pH 7. Degraded solutions were then diluted with deionized water in 10-mL volumetric flasks to reach a specified final concentration of 100 µg/mL.
Volume of 1 mL of deionized water was added to NMV solutions. The mixture was heated in a thermostated water bath at 80 ℃ for 1 h, then cooled and completed to volume as mentioned above.
NMV solution was treated with 1 mL of hydrogen peroxide (30% v/v), heated in a thermostated water bath at 60 ℃ for ½ h. The solutions were then cooled and diluted as above.
Photo-stability study was performed by exposing NMV stock solution (1 mL of 1000 µg/mL) to sunlight during daytime for 4 h. After the specified time, the volume was completed as shown above.
Kinetics investigation and Arrhenius plot
Degradation of NMV in an alkaline medium using 0.1 M NaOH and in an acidic medium with 2 M HCl was studied as a function of temperature.
For alkaline hydrolysis, volumes of 1 mL 0.1 M NaOH were separately added to NMV stock solution (1 mL of 1000 µg/mL) and mixed well. The solutions were kept in a thermostated water bath at different temperatures (50, 60 and 70 °C). After specified time intervals (0.5, 1, 1.5 and 2 h), the solutions were taken, left to cool, then neutralized to pH 7 before being diluted with the mobile phase to reach a final concentration of 100 μg/mL of NMV.
For acidic hydrolysis, separate volumes of 1 mL 2 M HCl were added to 1 mL of NMV stock solution and mixed well. The solutions were kept in a thermostated water bath at different temperatures (60, 70 and 80 °C). After the specified time intervals, 0.5, 1, 1.5 and 2 h, the solutions were then neutralized and diluted with the mobile phase to reach final concentrations of 100 μg/mL of NMV.
The regression equation for NMV determination was then used to determine the recovered concentration of the intact drug in the treated solutions. By monitoring the decline in NMV concentration over time, the degradation kinetics was determined.
LC–MS-UV study of NMV and its degradation products
LC–MS-UV was used to elucidate the structure of NMV DPs’ under alkaline and acidic stress conditions. Alkaline-induced DPs were generated using 0.1 M NaOH, 25 °C, 4 h or 0.1 M NaOH, 100 °C, 1 h, for partial and complete hydrolysis, respectively. Acidic-induced DPs were generated through partial and complete degradation of NMV using 1 M HCl at 90 °C for 1 h and 2 M HCl at 100 °C for 4 h, respectively. Degraded solutions were neutralised before being diluted to 100 μg/mL of NMV. NMV degraded solutions were initially injected into the LC–MS-UV system for the purpose of identification using the settings listed under section “Chromatographic conditions”. Using ESI with a mass range of 100 to 1000 m/z, the peaks of NMV and its probable DPs were monitored in both positive and negative scan mode.
In-silico toxicity studies
Using the Pro Tox-II web server tool, the toxicity of NMV and its main DPs was predicted in silico. Leave-one-out cross-validation was used to validate the ProTox-II acute toxicity model. The parameters of ProTox-prediction II had been adjusted to enhance toxicity class and LD50 prediction hit rates. Generally, after being exposed to a chemical, the median lethal dose (LD50) was calculated .
Greenness of the method
Different methodologies were used to assess the degree of greenness of the HPLC–DAD involved in the study. These methodologies include the Analytical Eco-Scale tool, the green analytical procedure index (GAPI), and the Analytical greenness metric (AGREE).
Results and discussion
A green stability-indicating assay method that can separate NMV and its DPs from each other is required to conduct a detailed kinetic degradation investigation of NMV. DPs should be structurally characterized to define their in-silico toxicity. As a result, Pro Tox- II was used to determine their potential toxicity. Finally, the proposed HPLC–DAD method’s degree of greenness was evaluated using different tools.
Optimization of the HPLC–DAD chromatographic conditions
Two reversed phase columns were tried, namely Agilent Zorbax eclipse-C18 (150 × 4.6 mm, 5-µm) and Agilent Zorbax eclipse-C18 (250 × 4.6 mm, 5-µm). The separation efficiency of the two columns was evaluated using mobile phases of ACN and ammonium acetate, pH = 5 in the ratio of 50:50, v/v, respectively. Experimental trials revealed that the longer column achieved better separation between NMV and its DPs. Thus, the longer column, of 250 mm length, was used for the study.
Different pH values of the aqueous phase were tried (3, 5 and 7). It was found that pH 3 and 5 gave better peak shape since the pka value of NMV is 7.1 . While pH 3 and 5 gave similar results concerning peak area, resolution and tailing factor, pH 5 was chosen as it is more suitable with respect to the buffer capacity of ammonium acetate . Both MeOH and ACN were tried as organic modifiers. ACN produced more symmetrical and sharper peaks compared with MeOH, thus ACN was used in subsequent trials. Isocratic elution was suitable for the development of a stability indicating method of NMV with adequate resolution between NMV and its nearest DP. Pure peaks were confirmed by calculating purity factors using DAD and they were within the appropriate automatically calculated noise acceptable limit, Additional file 1: Fig. S1. The run time was 5 min, with a 4.01-min retention time for NMV determination and the detection wavelength was 225 nm which corresponded to NMV maximum wavelength as shown in Fig. 1.
Forced degradation and stability-indicating study
NMV was subjected to different stress conditions, including acid, base, neutral, oxidative, and photo-degradation, as per ICH guidelines . The proposed HPLC–DAD method was used to track the degradation of NMV. It is noteworthy to mention that the choice of degradation conditions is a matter of “trial and error” to get reasonable % degradation of the intact drug, 30–70%. The stability-indicating potential of the developed methods is then assessed by its ability to separate the intact drug from its DPs. Accordingly, studying the degradation kinetics requires the use of stress conditions which result in stepwise degradation of the intact drug. On the other hand, structure elucidation of DPs requires the use of violent degradation conditions resulting optimally in 100% degradation of the intact drug.
Table 1 shows the different testing stress conditions that were tried during method development. Preliminary experiments showed that NMV remained stable with nearly 100% recovery under the following conditions: oxidative hydrolysis with 10% H2O2 at room temperature for 24 h, neutral hydrolysis with H2O, 60 °C, ½ h, and acidic hydrolysis using 0.1 M HCl at room temperature for 24 h. However, under alkaline degradation conditions, 0.1 M NaOH at 60 °C for ½ h, excessive degradation of NMV was noticed, with 39.56% recovery for NMV, Fig. 2. Further experimentation showed that NMV remained stable following the application of stronger oxidative/neutral degradation conditions, Table 1. On the other hand, stronger acidic conditions, 1 M HCl, 90 °C for 1 h, resulted in NMV degradation with 66.26% recovery. These findings revealed higher susceptibility of NMV to alkaline compared with acidic hydrolysis and that the drug was stable to oxidative, neutral, and sun light degradation conditions, as shown in Table 1, Fig. 2.
In conclusion, the proposed stability-indicating HPLC–DAD method enabled the separation of NMV from its DP under different degradation conditions. Peak purity was assessed by DAD, Additional file 1: Fig. S1.
Validation of the proposed HPLC–DAD method
As a crucial step in developing analytical methods, system suitability parameters were determined to assess the chromatographic system's performance. The obtained results, shown in Additional file 1: Table S1, included resolution, tailing factor, column efficiency, and theoretical plates, and they all fell within FDA acceptance standards .
Different validation parameters were evaluated as per FDA guidelines  including linearity and range, detection and quantitation limits, accuracy and precision, specificity, robustness, and stability of solutions.
Linearity and concentrations ranges
The linearity of the proposed HPLC–DAD method was evaluated in the concentration range of 5 to 500 µg/mL NMV and the linear regression equation Y = 4.14 + 1.10X was derived. Good linearity was shown by the high correlation coefficient (r = 0.9998). Additionally, the slope's RSD percent (Sb%) was determined to be 0.01 (less than 2%) indicating a small deviation around the slope. In addition, the small significance F (6.3 × 10–8) values indicated little scattering of the experimental data sets around the regression line [19, 20].
Detection and quantitation limits
Limits of detection (LOD) and quantitation (LOQ) were evaluated using a signal-to-noise ratio of 3:1 for LOD and 10:1 for LOQ. Low LOD and LOQ values (0.6 and 2 µg/mL, respectively) of the proposed method provided evidence of its good sensitivity.
Accuracy and precision
Three selected concentrations of NMV were analysed three times on the same day (intra-day) and on three successive days (inter-day). To assess the method's accuracy and precision, relative error (Er) and relative standard deviations (RSD) were calculated, respectively, Additional file 1: Table S2. The low values of Er and RSD illustrated the method's high accuracy and precision .
According to FDA guidelines , the specificity was evaluated by the method's capacity to successfully quantify NMV while being separated from its forced degradants with a reasonable resolution (more than 2). Additionally, NMV peak purity, assessed by DAD, was an index of NMV purity during its analysis with the potential degradants (stress degradation study), Additional file 1: Fig. S1.
The working wavelength (± 1 nm), mobile phase ratio (± 2%) and buffer pH (± 0.2 pH units) of the working chromatographic conditions were slightly altered to test the method's robustness. It was found that no significant changes were recorded on the measured responses or retention times, RSD% less than 2%, Additional file 1: Table S3.
Stability of solutions
No chromatographic changes were observed for NMV standards when kept at room temperature for 6 h. Additionally, it was found that NMV stock solution, prepared in HPLC-grade MeOH, was stable when kept refrigerated at 0 °C for a week. This was supported by the absence of any foreign peaks besides unchanged peak area and retention times of NMV.
Investigation of degradation kinetics and Arrhenius plot
The susceptibility of NMV to degradation in alkaline and acidic conditions paid our attention to extend the applicability of the proposed HPLC–DAD method to study the degradation kinetics, a novel aspect in our study. This was done by calculating the concentration of NMV remaining over a period (t) at different temperatures (50, 60, and 70 °C) and (60, 70 and 80 °C) for alkaline (0.1 M NaOH) and acidic (2 M HCl) degradation, respectively. The half-life (t1/2) and rate constant (K) of NMV degradation at the specified temperatures were calculated. In addition, the Arrhenius plot enabled the prediction of the on-shelf stability of NMV at ambient temperature. Figures 3 and 4 showed a linear decline in the log of the percent NMV remaining over time under the alkaline and acidic conditions, respectively. This indicated pseudo first-order degradation kinetics. The following equations were used .
where, Ct is NMV remaining concentration, after time (t), Co is the initial concentration of NMV,
K, with negative sign, is the apparent first order rate constant.
For NMV alkaline hydrolysis, the (t1/2) values calculated at 70, 60 and 50 ℃ were found to be 0.28, 0.35 and 0.52 h, with corresponding observed reaction rate constants (kobs) of 2.46, 1.95 and 1.33 h−1, respectively. On the other hand, NMV acid hydrolysis resulted in (t1/2) values of 0.78, 1.11 and 1.72 h calculated at 80, 70 and 60 ℃, respectively with corresponding (kobs) values of 0.89, 0.62 and 0.40 h−1, respectively, Table 2.
The Arrhenius plot was then derived by plotting the logarithm of degradation rate constant of NMV (kobs) against the reciprocals of the absolute temperature  as illustrated in Figs. 5 and 6 for the alkaline (0.1 M NaOH) and acidic (2 M HCl) degradation of NMV, respectively. The stability of NMV to the applied degradation conditions was predicted at room temperature by extrapolation to get logK25, Figs. 5 and 6. The t1/2 value of NMV degradation was then calculated, 0.56 h and 5.39 h for alkaline and acidic degradation, respectively.
LC–MS-UV study and postulation of degradation mechanisms
Different acidic and alkaline degraded solutions were injected into the LC–MS-UV instrument, and the m/z of the eluted peaks was recorded to investigate the degradation process of NMV. DPs produced under acidic and alkaline conditions with their equivalent m/z values using the scan mode were shown in Additional file 1: Figs. S2–S4.
Partial alkaline degradation of NMV was accomplished using 0.1 M NaOH at RT for 4 h, mild conditions, with the formation of different intermediates (INs) along with NMV peak, m/z 501. The corresponding m/z of the formed INs were 519, 537 and 555, Additional file 1: Fig. S3a. INs peaks were probably formed by the hydrolysis of amide group in the oxopyrrolidine ring followed by opening of cyclopropane ring and hydrolysis of the nitrile group. The mechanism of base-induced amide hydrolysis was previously studied [23, 24]. In addition, ring opening of cyclopropane using NaOH was studied later by Wong et al. . Due to steric hindrance in NMV structure, it was postulated that the first hydrolysis step was the oxopyrrolidine ring through the addition of one molecule of water to NMV resulting in the formation of IN-1 with m/z = 501 + 18 = 519, then opening of the cyclopropane ring by the addition of the second water molecule resulting in IN-2 peak with m/z = 519 + 18 = 537. Afterwards, the nitrile group would undergo hydrolysis to the corresponding amide and an IN-3 of m/z = 555 (537 + 18) was formed. On the other hand, complete degradation of NMV was practically possible using 0.1 M NaOH at 100 °C for 1 h, hard conditions. In this case, the INs peaks disappeared, and other DPs peaks, DP-1, DP-2 and DP-3 were formed at m/z 423, 441 and 459, Additional file 1: Fig. S3b, as a result of the hydrolysis of the trifluoroacetamide group from IN-1, IN-2 and IN-3, respectively. Moreover, hydrolysis of the amide of the azabicyclo ring resulted in the appearance of DP-4 with m/z 291. In addition, Hydrolysis of the two aliphatic amide groups yielded DP-5 with m/z 269. Figure 7 summarizes the proposed pathways of the alkaline degradation of NMV. It is noteworthy to mention that the proposed degradation pathway was different from that postulated by Secreatan et al. . The latter attributed the alkaline degradation to the hydrolysis of the nitrile group and/or trifluoroacetamide group only. However, in the current study, it was proposed that the degradation occurred due to the hydrolysis of the amide group in the oxopyrrolidine ring followed by opening of cyclopropane ring and hydrolysis of the nitrile group. Moreover, hydrolysis of the trifluoroacetamide group and the amide of the azabicyclo ring took place under hard conditions. The difference between the two pathways, our proposed work and ref. , could be attributed to the variable conditions used in the degradation, the latter  used 1 M NaOH for one week at room temperature (long term stability study) while our current study used 0.1 M NaOH at RT for 4 h and 0.1 M NaOH at 100 °C for 1 h (short term stability study).As per alkaline degradation, mild acidic conditions, 1 M HCl at 90 °C for 1 h, resulted in the formation of INs with m/z 519 and 537, Additional file 1: Fig. S4a. The formation of IN 1 was previously explained under alkaline degradation, while the formation of IN 2 was attributed to the opening of the cyclopropane ring using HCl. Hard acidic conditions with 2 M HCl at 100 °C for 4 h resulted in the formation of DPs peaks, DP-1, DP-4, DP-5, and DP-6 corresponding to m/z 423, 291, 269 and 441, Additional file 1: Fig. S4b. Further hydrolysis of INs with removal of the trifluoroacetamide group resulted in the formation of DP-1 and DP-6 of m/z 423 and 441, respectively. Additionally, the hydrolysis of the amide of the azabicyclo ring resulted in the formation of DP-4 with m/z 291, while hydrolysis of the two aliphatic amide groups yielded DP-5 with m/z 269. The mechanism of acid-induced amide hydrolysis was previously postulated [23, 24]. Moreover, the addition of hydrogen chloride to cyclopropane was discussed by Lee et al. and Lambert et al. [26, 27]. The proposed pathways for acidic hydrolysis of NMV was shown in Fig. 8. As per alkaline hydrolysis, degradation pathways suggested in our study was different from that proposed by Secreatan et al.  which attributed the acidic degradation to the hydrolysis of the nitrile group or the trifluoroacetamide group only. Our proposed pathway assumed that the degradation was due to hydrolysis of amide group in the oxopyrrolidine ring or by opening of cyclopropane ring followed by hydrolysis of the trifluoroacetamide group under hard conditions. Moreover, the mechanism of hydrolysis of the azabicyclo ring and the two aliphatic amide groups were differently postulated which could be attributed to the different degradation conditions in the two studies; the published study  used 1 M HCl for one week at RT (long term stability study) while the current one used 1 M HCl at 90 °C for 1 h or 2 M HCl at 100 °C for 4 h (short term stability study).
In-silico toxicity studies
Numerous toxicity end points were predicted using the popular computational toxicity web server, ProTox-II. It featured details about chemical and molecular targets, including computer models, pharmacophores, fragment tolerance, and molecular similarity. ProTox-II models were validated utilising the probability-based CLUSTER cross-validation with selective oversampling . In this study, this database was applied to predict the hazardous potential of NMV and its DPs including cytotoxicity, immunotoxicity, hepatotoxicity, mutagenicity, and carcinogenicity. Table 3 showed that NMV and the proposed DPs (1 to 6) did not exhibit hepatotoxicity, carcinogenicity, mutagenicity, or cytotoxicity with confidence levels ranging from 0.57 to 0.87 . Additionally, the ability of NMV and its DPs to kill 50% of test animals (rodents) after 24 h of exposure (LD50) was used to assess their acute toxicity . NMV, DP-1, DP-2, DP-3, and DP-4 each had an LD50 of 3000 mg/kg, while the LD50 for DP-5 and DP-6 were found to be 5000 and 826 mg/kg, respectively. Based on LD50 (mg/kg) and with reference to ‘classification and labelling of compounds” methodology, NMV and all its DPs, except for DP-6, were classified as being in toxicity class V, while DP-6, being in class IV, is harmful if swallowed .
Greenness of the method
The Analytical Eco-Scale , the Green Analytical Procedure Index (GAPI) , and the innovative Analytical Greenness metric (AGREE)  were used in our study to evaluate the degree of greenness of the proposed HPLC–DAD method. Analytical Eco-Scale tool quantifies the method’s degree of greenness  based on assigning penalty points (PPs) to each category of the method that deviates from the notion of greenness. The volume of any created waste, the amount of hazardous chemicals that were utilised, and the amount of electricity and energy spent were the categories to be examined. The total determined penalty points were then subtracted from 100. The method is then deemed to be of good greenness if the value is above 50 and outstanding greenness if it is above 75. This tool determined that the estimated value of our proposed method was greater than 75 (Eco-Scale score = 84), Table 4, indicating that it was extremely environmentally friendly with the least adverse effects.
In addition, our method’s greenness was evaluated using the GAPI and AGREE tools. The GAPI tool assesses the entire analytical procedure, from sample collection and preparation to the conclusion. A vibrant pictogram with 15 elements that refer to 15 assessed process criteria, such as sample preparation and extraction, solvent and reagent properties, energy consumption, waste generation, and occupational dangers, is used to illustrate the method’s eco-friendliness. The non-green characters are given a red tint, and each component is either green (for the greenest element) or yellow (for the least green element). By evaluating our method, it is considered green with 7 green, 6 yellow and 2 red zones, Fig. 9.
Lastly, the more recent AGREE metric was used to provide a more detailed, inclusive, and succinct evaluation of the greenness of the proposed method, as well as an interdisciplinary grade based on the fundamental ideas of GAC . This application could be downloaded, and after inputting the 12 characters for the method, it generates a vivid pictogram with 12 segments that are each coloured differently, ranging from deep green to deep red. Additionally, it displays a fractional number in the middle of the circular pictogram; the greener the technique, the closer the final score is to 1. Utilizing this novel strategy led to a more accurate and quantifiable assessment of the methodologies. The AGREE value of the proposed HPLC–DAD method was 0.81 which indicated the greenness of the method as shown in Fig. 10. (The AGREE report was displayed in the Additional file 2).
The aforementioned outcomes showed that the proposed HPLC–DAD method was successfully green.
A quick, reliable, and precise stability-indicating HPLC–DAD method was developed and validated, for the first time, followed by a kinetics degradation investigation, to measure NMV in the presence of its forced DPs. It was found that NMV was highly vulnerable to alkali hydrolysis, compared with acid hydrolysis, and that it was stable to conditions that promote oxidative, neutral, and sun-induced degradation. Studying the rates of acid and alkali hydrolysis revealed that the kinetics followed pseudo-first order. Subsequent Arrhenius plot was used to determine the half lifetime and reaction constants of degradation under the studied conditions at room temperature using the proposed method. LC–MS-UV was used to suggest the mechanism of the acid and base degradation pathways. Toxicities of NMV and its DPs were evaluated using ProTox-II and they were determined to be negligibly harmful. Finally, greenness appraisal proved that the proposed HPLC–DAD method was of excellent greenness.
Availability of data and materials
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
High performance liquid chromatography with diode array detector
Limit of detection
Limit of quantification
Ultra-fast liquid chromatography-mass spectrometry
Green analytical procedure index
Analytical Greenness metric
- t1/2 :
- kobs :
Reaction rate constants
Docherty AB, Harrison EM, Green CA, Hardwick HE, Pius R, Norman L, Holden KA, Read JM, Dondelinger F, Carson G, Merson L, Lee J, Plotkin D, Sigfrid L, Halpin S, Jackson C, Gamble C, Horby PW, Nguyen-Van-Tam JS, Ho A, Russell CD, Dunning J, Openshaw PJ, Baillie JK, Semple MG. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO clinical characterisation protocol: prospective observational cohort study. BMJ. 2020;369:1985. https://doi.org/10.1136/bmj.m1985.
https://coronavirus.jhu.edu/map.html. Accessed 25 July 2022.
Kim L, Garg S, O’Halloran A, Whitaker M, Pham H, Anderson EJ, Armistead I, Bennett NM, Billing L, Como-Sabetti K, Hill M, Kim S, Monroe ML, Muse A, Reingold AL, Schaffner W, Sutton M, Talbot HK, Torres SM, Yousey-Hindes K, Holstein R, Cummings C, Brammer L, Hall AJ, Fry AM, Langley GE. Risk factors for intensive care unit admission and In-hospital mortality among hospitalized adults identified through the US coronavirus disease 2019 (COVID-19)-associated hospitalization surveillance network (COVID-NET). Clin Infect Dis. 2021;72(9):e206–14. https://doi.org/10.1093/cid/ciaa1012.
Thakur B, Dubey P, Benitez J, Torres JP, Reddy S, Shokar N, Aung K, Mukherjee D, Dwivedi AK. A systematic review and meta-analysis of geographic differences in comorbidities and associated severity and mortality among individuals with COVID-19. Sci Rep. 2021;11(1):8562. https://doi.org/10.1038/s41598-021-88130-w.
Zheng Z, Peng F, Xu B, Zhao J, Liu H, Peng J, Li Q, Jiang C, Zhou Y, Liu S, Ye C, Zhang P, Xing Y, Guo H, Tang W. Risk factors of critical & mortal COVID-19 cases: a systematic literature review and meta-analysis. J Infect. 2020;81(2):e16–25. https://doi.org/10.1016/j.jinf.2020.04.021.
Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) Outbreak in China: summary of a report of 72 314 cases from the Chinese center for disease control and prevention. JAMA. 2020;323(13):1239–42. https://doi.org/10.1001/jama.2020.2648.
Hammond J, Leister-Tebbe H, Gardner A, Abreu P, Bao W, Wisemandle W, Baniecki M, Hendrick VM, Damle B, Simón-Campos A, Pypstra R, Rusnak JM. Oral nirmatrelvir for high-risk, nonhospitalized adults with Covid-19. N Engl J Med. 2022;386(15):1397–408. https://doi.org/10.1056/NEJMoa2118542.
Owen DR, Allerton CMN, Anderson AS, Aschenbrenner L, Avery M, Berritt S, Boras B, Cardin RD, Carlo A, Coffman KJ, Dantonio A, Di L, Eng H, Ferre R, Gajiwala KS, Gibson SA, Greasley SE, Hurst BL, Kadar EP, Kalgutkar AS, Lee JC, Lee J, Liu W, Mason SW, Noell S, Novak JJ, Obach RS, Ogilvie K, Patel NC, Pettersson M, Rai DK, Reese MR, Sammons MF, Sathish JG, Singh RSP, Steppan CM, Stewart AE, Tuttle JB, Updyke L, Verhoest PR, Wei L, Yang Q, Zhu Y. An oral SARS-CoV-2 M(pro) inhibitor clinical candidate for the treatment of COVID-19. Science. 2021;374(6575):1586–93. https://doi.org/10.1126/science.abl4784.
Abdelnabi R, Foo CS, Jochmans D, Vangeel L, De Jonghe S, Augustijns P, Mols R, Weynand B, Wattanakul T, Hoglund RM, Tarning J, Mowbray CE, Sjö P, Escudié F, Scandale I, Chatelain E, Neyts J. The oral protease inhibitor (PF-07321332) protects Syrian hamsters against infection with SARS-CoV-2 variants of concern. Nat Commun. 2022;13(1):719. https://doi.org/10.1038/s41467-022-28354-0.
https://www.fda.gov/drugs/news-events-human-drugs/fda-updates-paxlovid-health-care-providers. Accessed 25 July 2022.
Garnero C, Zoppi A, Aloisio C, Longhi MR. Chapter 7-Technological delivery systems to improve biopharmaceutical properties. In: Grumezescu AM, editor. Nanoscale fabrication, optimization scale-up and biological aspects of pharmaceutical nanotechnology. Norwich: William Andrew Publishing; 2018. p. 253–99.
Liu C, Zhu M, Cao L, Boucetta H, Song M, Hang T, Lu Y. Simultaneous determination of nirmatrelvir and ritonavir in human plasma using LC-MS/MS and its pharmacokinetic application in healthy Chinese volunteers. Biomed Chromatogr BMC. 2022. https://doi.org/10.1002/bmc.5456.
Secretan P-H, Annereau M, Kini-Matondo W, Prost B, Prudhomme J, Bournane L, Paul M, Yagoubi N, Sadou-Yayé H, Do B. Unequal behaviour between hydrolysable functions of Nirmatrelvir under stress conditions: structural and theoretical approaches in support of preformulation studies. Pharmaceutics. 2022;14(8):1720.
ICH. Stability testing of new drug substances and products Q1A(R2), Geneva (February 2003). https://www.ich.org/fileadmin/Public_Web_Site/ICH_Products/Guidelines/Quality/Q1A_R2/Step4/Q1A_R2__Guideline.pdf. Accessed 25 July 2022.
Banerjee P, Eckert AO, Schrey AK, Preissner R. ProTox-II: a webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018;46(W1):W257–63. https://doi.org/10.1093/nar/gky318%JNucleicAcidsResearch.
https://go.drugbank.com/drugs/DB16691. Accessed 25 July 2022.
Konermann L. Addressing a common misconception: ammonium acetate as neutral pH “buffer” for native electrospray mass spectrometry. J Am Soc Mass Spectrom. 2017;28(9):1827–35. https://doi.org/10.1007/s13361-017-1739-3.
C.o.D.E.R.C. FDA, Reviewer Guidance: validation of chromatographic methods, Washington, USA. 2018. https://www.fda.gov/files/drugs/published/Bioanalytical-Method-Validation-Guidance-for-Industry.pdf. Accessed 25 July 2022.
Miller JNM, Miller JC. Statistics and chemometrics for analytical chemistry. 5th ed. Harlow: Pearson Prentice Hall; 2005.
Armitage P, Berry G. Statistical methods in medical research. 3rd ed. Oxford: Blackwell; 1994.
Bhangare D, Rajput N, Jadav T, Sahu AK, Tekade RK, Sengupta P. Systematic strategies for degradation kinetic study of pharmaceuticals: an issue of utmost importance concerning current stability analysis practices. J Anal Sci Technol. 2022;13(1):7. https://doi.org/10.1186/s40543-022-00317-6.
Attwood D, Florence AT. Physical pharmacy. 1st ed. London: Pharmaceutical Press; 2008.
Duan P, Dai L, Savage PE. Kinetics and mechanism of N-substituted amide hydrolysis in high-temperature water. J Supercrit Fluids. 2010;51(3):362–8. https://doi.org/10.1016/j.supflu.2009.09.012.
Meresaar U, Bratt L, Alakuijala P, Oinonen L. Hydrolysis of amides. Alkaline and general acid catalyzed alkaline hydrolysis of some substituted acetamides and benzamides. Acta Chem Scand. 1974;28:715–22.
Wong HNC, Hon MY, Tse CW, Yip YC, Tanko J, Hudlicky T. Use of cyclopropanes and their derivatives in organic synthesis. Chem Rev. 1989;89(1):165–98. https://doi.org/10.1021/cr00091a005.
Lee CC, Hahn B-S, Wan K-M, Woodcock DJ. Reactions of cyclopropane with hydrochloric acid and of bromocyclopropane with hydrobromic acid. J Org Chem. 1969;34(10):3210–1. https://doi.org/10.1021/jo01262a097.
Lambert JB, Napoli JJ, Johnson KK, Taba KN. Scope, limitations, and mechanism of the homoconjugate electrophilic addition of hydrogen halides. J Org Chem. 1985;50(8):1291–5. https://doi.org/10.1021/jo00208a027.
https://tox-new.charite.de/protox_II/. Accessed 25 July 2022.
https://tox-new.charite.de/protox_II/index.php?site=faq. Accessed 25 July 2022.
Gałuszka A, Migaszewski ZM, Konieczka P, Namieśnik J. Analytical Eco-Scale for assessing the greenness of analytical procedures. TrAC Trends Anal Chem. 2012;37:61–72. https://doi.org/10.1016/j.trac.2012.03.013.
Płotka-Wasylka J. A new tool for the evaluation of the analytical procedure: green analytical procedure index. Talanta. 2018;181:204–9. https://doi.org/10.1016/j.talanta.2018.01.013.
Pena-Pereira F, Wojnowski W, Tobiszewski M. AGREE—analytical GREEnness metric approach and software. Anal Chem. 2020;92(14):10076–82. https://doi.org/10.1021/acs.analchem.0c01887.
Gałuszka A, Migaszewski Z, Namieśnik J. The 12 principles of green analytical chemistry and the SIGNIFICANCE mnemonic of green analytical practices. TrAC Trends Anal Chem. 2013;50:78–84. https://doi.org/10.1016/j.trac.2013.04.010.
We would like to thank Prof. Johan Neyts and Dr. Rana Abdelnabi at KU Leuven Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Laboratory of Virology and Chemotherapy, 3000, Leuven, Belgium, for their supply of samples of standard Nirmatrelvir.
Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1
: Figure S1 Purity plots of NMV peaks obtained from DAD for standard solution, a), and following exposure to different degradation conditions, alkaline, b) acidic, c), oxidative, d) neutral, e), and sunlight, f), degradations. Figure S2 LC-MS-UV chromatogram of NMV standard solution scanned at different m/z, indicating presence of a peak at m/z 501 for NMV and absence of peaks at other m/z proposed for INs and DPs. Figure S3 LC-MS-UV chromatograms of MLN using different alkaline stress degradation conditions, mild conditions: 0.1 M NaOH at RT for 4 h indicating presence of a peak at m/z 501 for NMV and m/z for IN-1, 2 & 3, a), and hard conditions: 0.1 M NaOH at 100 °C for 1 h, b) indicating absence of a peak at m/z 501 for NMV due to complete degradation and m/z for DP-1, 2, 3, 4 & 5. Figure S4 LC-MS-UV chromatograms of MLN using different acidic stress conditions degradation, mild conditions: 1 M HCl at 90 for 1 h indicating presence of a peak at m/z 501 for NMV and m/z for IN-1, 4 & 3, a), and hard conditions: 2 M HCl at 100 °C for 4 h, b) indicating absence of a peak at m/z 501 for NMV due to complete degradation and m/z for DP-1, 6, 3, 4 & 5. Table S1 System suitability parameters for the HPLC-DAD determination of NMVa. Table S2 Intra-day and inter-day accuracy and precision for the determination of NMV. Table S3 Evaluation of the robustness of the proposed HPLC method for the determination of NMV.
Additional file 2
: Analytical Greenness report sheet.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Aboras, S.I., Maher, H.M. Green adherent degradation kinetics study of Nirmatrelvir, an oral anti-COVID-19: characterization of degradation products using LC–MS with insilico toxicity profile. BMC Chemistry 17, 23 (2023). https://doi.org/10.1186/s13065-023-00928-z