Physiologically-Based Pharmacokinetic Modeling-Guided Dose Management of Oral Anticoagulants when Initiating Nirmatrelvir/Ritonavir (Paxlovid) for COVID-19 Treatment.

Patients with coronavirus disease 2019 (COVID-19) with cardiovascular diseases who are at higher risk of progressing to critical illness should be treated with nirmatrelvir/ritonavir (Paxlovid). Ritonavir, the booster in nirmatrelvir/ritonavir, modulates multiple drug metabolizing enzymes and transporters, complicating its use in real-world clinics. We aimed to apply physiologically-based pharmacokinetic (PBPK) modeling to simulate the complex drug-drug interactions (DDIs) of ritonavir with two anticoagulants, rivaroxaban and racemic warfarin, to address this important clinical conundrum. Simulations were implemented within Simcyp Simulator. Compound and population models were adopted from Simcyp and our previous studies. Upon verification and validation of the PBPK model of ritonavir, prospective DDI simulations with the anticoagulants were performed in both the general population (20-65 years) and geriatric subjects (65-85 years) with or without moderate renal impairment. Elevated rivaroxaban concentrations were simulated with nirmatrelvir/ritonavir treatment, where the impact was more profound among geriatric subjects with renal impairment. The overexposure of rivaroxaban was restored to normal range on day 4 post-discontinuation of nirmatrelvir/ritonavir, corroborating with the recovery of enzyme activity. A lower 10 mg daily dose of rivaroxaban could effectively maintain acceptable systemic exposure of rivaroxaban during nirmatrelvir/ritonavir treatment. Treatment of ritonavir marginally declined simulated S-warfarin concentrations, but substantially elevated that of R-warfarin, resulting in a decrease in the international normalized ratio (INR). As INR only recovered 2 weeks post-nirmatrelvir/ritonavir treatment, a longer surveillance INR for warfarin becomes important. Our PBPK-guided simulations evaluated clinically important yet untested DDIs and supports clinical studies to ensure proper anticoagulation management of patients with COVID-19 with chronic coagulative abnormalities when initiating nirmatrelvir/ritonavir therapy.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

Patients with coronavirus disease 2019 (COVID‐19) with cardiovascular diseases who are at higher risk of progressing to critical illness should be treated with nirmatrelvir/ritonavir (Paxlovid). Ritonavir, the booster in nirmatrelvir/ritonavir, modulates multiple drug metabolizing enzymes and transporters, complicating its concomitant use with two widely prescribed oral anticoagulants, rivaroxaban and warfarin.

WHAT QUESTION DID THIS STUDY ADDRESS?

For patients with COVID‐19 who undergo chronic anticoagulant therapy, what are the recommendations for appropriate oral anticoagulation management when initiating nirmatrelvir/ritonavir therapy?

WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

With the application of physiologically‐based pharmacokinetic modeling and simulations of complex disease‐drug–drug interactions, we provide scientific evidence of dose adjustment of rivaroxaban to 10 mg daily to mitigate the potential overexposure and up to 2‐weeks surveillance of international normalized ratio for warfarin treatment after nirmatrelvir/ritonavir discontinuation to avoid warfarin‐related thrombosis.

HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE?

Our novel findings are pertinent in supporting further clinical studies and/or systematic analyses of real‐world data to ensure efficacious and safe anticoagulant therapy while preventing progression to severe illness from COVID‐19 based on nirmatrelvir/ritonavir therapy.

Nirmatrelvir/ritonavir (Paxlovid) has received emergency use authorizations (EUAs) from the US Food and Drug Administration (FDA) for the treatment of patients with mild to moderate coronavirus disease (COVID‐19) at risk of progression to severe infection. COVID‐19 Treatment Guidelines also strongly recommended the use of nirmatrelvir/ritonavir for treating nonhospitalized patients. Treatment with nirmatrelvir/ritonavir within 5 days of the appearance of symptoms could effectively reduce the risk of hospitalization or death for high‐risk patients by 88%. Nirmatrelvir/ritonavir comprises nirmatrelvir and ritonavir, where ritonavir, a potent synthetic HIV protease inhibitor, was given at a sub‐antiviral‐dose of 100 mg twice daily to boost the systemic exposure of the novel severe acute respiratory syndrome‐coronavirus 2 (SARS‐CoV‐2) Mpro inhibitor, nirmatrelvir. ,

Ritonavir yields potent mechanism‐based inactivation (MBI) of cytochrome P450 3A4 enzyme (CYP3A4) and circumvents the rapid and extensive elimination of nirmatrelvir by CYP3A4. Additionally, in vitro and in vivo studies have shown that ritonavir could remarkedly inhibit or induce several other CYP450 enzymes (e.g., CYP2C9 and CYP2J2) as well as multiple transporters (e.g., P‐glycoprotein (P‐gp) and organic anion transporters (OATs)). , , Due to the unique co‐formulation of nirmatrelvir and ritonavir, the concomitant use of nirmatrelvir/ritonavir and certain other drugs may result in potentially significant drug–drug interactions (DDIs) and subsequent serious and/or life‐threatening reactions associated with elevated drug exposures. Such DDIs complicate the prescription of nirmatrelvir/ritonavir in real‐world clinics. ,

Patients with COVID‐19 have shown high incidence of thromboembolism and cardiac arrhythmias. , , Several large cohort clinical studies evaluating anticoagulation therapy in hospitalized patients with COVID‐19 recommended the use of heparin analogues as effective interventions, whereas oral anticoagulants should be used with caution. , , Hence, a clinical conundrum arises where patients with cardiovascular diseases who undergo chronic anticoagulant therapy are the most vulnerable cohort when diagnosed with COVID‐19, necessitating nirmatrelvir/ritonavir treatment to prevent the progression of COVID‐19 to critical illness, whereas it is recommended that these patients should continue taking the oral anticoagulants after they receive a diagnosis of COVID‐19 infection.

Rivaroxaban and warfarin are two widely prescribed oral anticoagulants. S‐warfarin is predominantly metabolized by CYP2C9, whereas CYP3A4 was responsible for the metabolism of R‐warfarin and rivaroxaban. Additionally, CYP2J2 and OAT3 are pivotally involved in hepatic metabolism and renal secretion of rivaroxaban, respectively. EUA recommended close monitoring of warfarin if its co‐administration with nirmatrelvir/ritonavir is necessary, whereas concomitant use of nirmatrelvir/ritonavir with rivaroxaban should be avoided. Clinical studies have reported that higher therapeutic doses of ritonavir (600 mg twice daily) increased systemic exposure of rivaroxaban by 2.53‐fold in healthy subjects, and 100 mg once or twice daily dose of ritonavir could modestly decrease S‐warfarin concentration. Furthermore, individual cases reported significant decline of international normalized ratio (INR) values by ritonavir during warfarin treatment. Nevertheless, there remains a lack of direct clinical evidence in guiding appropriate oral anticoagulation management for patients with COVID‐19 during and post‐nirmatrelvir/ritonavir therapy.

To address this important clinical conundrum in the pharmacotherapy of patients with COVID‐19 with nirmatrelvir/ritonavir who are already prescribed with rivaroxaban or warfarin, we applied physiologically‐based pharmacokinetic (PBPK) and PBPK/pharmacodynamic (PD) modeling and simulations to investigate complex disease‐drug–drug interactions (DDDIs) and recommend dosage adjustment and surveillance in clinical scenarios where the DDDIs are potentially unfavorable. Our findings are pertinent in supporting clinicians in treating patients with COVID‐19 efficaciously and safely using both nirmatrelvir/ritonavir and oral anticoagulants.

METHODS

The models and simulations were implemented within the population‐based Simcyp Simulator (version 19, Sheffield, UK). Compound profiles for midazolam and digoxin were directly adopted from Simcyp. The ritonavir model provided by Simcyp was optimized based on previous studies. , , Rivaroxaban PBPK model and racemic warfarin PBPK/PD model were judiciously constructed and directly utilized for our simulations. Details of model development are summarized in and Table . Nirmatrelvir was not incorporated into the simulations due to its low DDI potential. Healthy population (Healthy Volunteers), general White population (NEurCaucasian) and geriatric population (Geriatric NEC) provided in Simcyp were adopted as population models. The moderate renal impairment (creatinine clearance (CrCL) 30 to 49 mL/min) population model modified previously , were utilized for simulations. Each simulation was designed as a total size of 100 subjects (10 trials with 10 subjects for each).

The pharmacokinetic (PK) profiles of ritonavir following multiple oral dose administrations were first simulated to verify the performance of the modified ritonavir model. The model was further validated with clinically available DDI studies to assess the model liability for DDI simulations against different enzymes and transporters. The dosage regimens of drugs were matched to the design of the corresponding clinical studies. Details of model verification and validation are provided in .

Age is a critical factor that contributed to the COVID‐19‐related severe illness and deaths, and could substantially affect the PK of compounds as revealed by our previous studies. , Prospective DDDI simulations were performed using both general population (20–65 years) and geriatric subjects (65–85 years). Because nirmatrelvir/ritonavir is prescribed for the treatment of mild to moderate COVID‐19, we postulated that the impact of COVID‐19 related cytokine storm on enzymes and transporters (e.g., CYP3A4) is marginal. The effect was not incorporated into our modeling. Ritonavir dose was set as 100 mg twice daily for 5 days. In patients with moderate renal impairment, the dosage of nirmatrelvir/ritonavir is adjusted as 150 mg nirmatrelvir and 100 mg ritonavir twice daily, and a dose reduction of rivaroxaban from 20 to 15 mg daily is also required for patients with moderate renal impairment during the treatment of non‐valvular atrial fibrillation. Rivaroxaban was simulated with these two doses with or without moderate renal dysfunction, respectively. A lower 10 mg daily dose of rivaroxaban, which was approved for prophylaxis of deep vein thrombosis and venous thromboembolism, was prospectively simulated as a dose adjustment strategy, starting from day 1 of nirmatrelvir/ritonavir treatment to day 3 post‐nirmatrelvir/ritonavir treatment. The daily dose of racemic warfarin was fixed as 5 mg.

The DDI potential with rivaroxaban was assessed based on simulated plasma concentrations of rivaroxaban and the fold‐change of area under curve (AUC) on the fifth day of nirmatrelvir/ritonavir treatment (AUCD5,0–24 hours), which was further calculated for risk of major bleeding ( ). INR was utilized for the safety evaluation of warfarin therapy.

RESULTS

Ritonavir PBPK model verification

As illustrated in Table and Figure , the ritonavir PBPK model successfully reproduced the nonlinear PKs of ritonavir with daily doses ranging from 100 to 1200 mg, where 23 out of 27 simulated PK parameters fell within fold‐error criterion of 0.67–1.50. DDI predictive performance of ritonavir was further validated using available PK data sets with multiple substrates, where digoxin, midazolam, S‐warfarin, and rivaroxaban validated the effects of ritonavir on P‐gp, CYP3A4, CYP2C9, and OAT3, respectively (Table ).

Prospective PBPK DDDI simulation for rivaroxaban

DDDIs between rivaroxaban and ritonavir were simulated among virtual subjects with normal or impaired renal functions with clinically recommended dosage regimens. Elevated rivaroxaban concentrations were simulated for the four virtual populations along with nirmatrelvir/ritonavir treatment (Figure ), resulting in AUCD5,0‐24h fold‐change exceeding the 1.43 threshold (Table ). The impact was more profound among geriatric subjects with renal impairment, where simulated risk of major bleeding was as high as 24.03% (Table ). As illustrated in Figure , the overexposure of rivaroxaban normalized on day 4 post‐discontinuation of nirmatrelvir/ritonavir treatment, which corroborated with the recovery of simulated enzyme activities to 62.44, 82.12, and 77.58% for CYP3A4 (liver), CYP3A4 (gut), and CYP2J2, respectively (Figure ). Simulation results of the lower 10 mg dose of rivaroxaban showed that the dose adjustment could effectively mitigate the bleeding tendency by half for all four cohorts while maintaining acceptable systemic exposure of rivaroxaban (Figure , Table ).

Figure 1

Plasma concentration profile over time of rivaroxaban in the general population (20–65 years) (a, c) and geriatric subjects (65–85 years) (b, d) with normal renal function a and b or moderate renal impairment (creatinine clearance (CrCL) 30 to 49 mL/min) c and d. Normal rivaroxaban daily dose of 20 mg and 15 mg was simulated for subjects with normal renal function or moderate renal impairment, respectively. Gray background represents 5‐days treatment of ritonavir of 100 mg twice daily. Data are presented as mean value. Plasma concentrations of 249 and 44 ng/mL were obtained from previous rivaroxaban Phase II studies as maximum and minimum plasma concentration at steady‐state, respectively.

Prospective PBPK/PD simulation for warfarin

Expectedly, 5‐days’ treatment of ritonavir marginally declined S‐warfarin concentrations while substantially elevated that of R‐warfarin (Figure ). The polarized PKs of warfarin enantiomers resulted in a decrease in INR, which was more notable among older adults (Figure ). The simulated steady‐state INRs before nirmatrelvir/ritonavir treatment were 1.9 and 2.4, whereas the lowest INRs during treatment were 1.3 and 1.6 for the general population and older adults, respectively. In view of the indirect response of INR toward the changes of warfarin concentrations, INR values was simulated to return to pretreatment levels by week 2 post‐nirmatrelvir/ritonavir treatment (Figure ).

Figure 2

Plasma concentration profile of warfarin enantiomers (a, b) and INR profile (c, d) over time in the general population (20–65 years) a and c and geriatric subjects (65–85 years) b and d. Racemic warfarin daily dose of 5 mg was simulated. Gray background represents 5‐days treatment of ritonavir of 100 mg twice daily. Data are presented as mean value. INR, international normalized ratio.

DISCUSSION

COVID‐19 has been associated with inflammation and a prothrombotic state, resulting in high prevalence of cardiac arrhythmias and venous thromboembolism among patients with COVID‐19. , , Meanwhile, patients with heart conditions and stroke disease are of high risk with COVID‐19 infection, as well as for progression to severe COVID‐19 condition. Nirmatrelvir/ritonavir treatment should be initiated with these nonhospitalized patients as soon as possible while appropriately maintaining their anticoagulant therapy at the same time. , However, direct clinical data between nirmatrelvir/ritonavir and oral anticoagulants are still limited. , Notably, application of PBPK models in diseased and specific populations, in addition to DDIs, has demonstrated promising performance in the evaluation of clinically relevant yet untested scenarios. PBPK‐guided dose recommendations, as an alternative to real‐world studies, have been approved in several drug labels in recent years, , and has been applied for nirmatrelvir/ritonavir‐related study. In view of this urgent conundrum where clinical studies cannot be conducted in a timely manner, our study aimed to utilize a ritonavir PBPK model, which was validated against multiple DDI data and elimination pathways, to investigate the DDI potential of ritonavir with rivaroxaban and warfarin, and offer dosing and risk management strategy, during and after 5‐days of nirmatrelvir/ritonavir treatment.

Ritonavir inhibits OAT3‐ and P‐gp‐mediated transport and causes MBI of CYP3A4 and CYP2J2, , , where these enzymes and transporters are implicated in the metabolism and excretion of rivaroxaban. Considering the complex DDDIs, it remains critical to establish the currently arcane risk between rivaroxaban and nirmatrelvir/ritonavir comprising lower subtherapeutic dose of ritonavir, despite clinical DDI that had already been reported for rivaroxaban with 1,200 mg daily dose of ritonavir. Our quantitative and mechanistical simulations verified the potent inhibition by both dosage regimens of ritonavir against enzymes and transporters associated with elimination of rivaroxaban, especially the comparable inhibition potency against CYP3A4 (Figure ). Corroborating our previous studies, , prospective simulations revealed that geriatric patients with moderate renal dysfunction were more susceptible to the DDDIs between rivaroxaban and ritonavir, who could potentially benefit from a further dose adjustment to 10 mg daily rivaroxaban while receiving nirmatrelvir/ritonavir treatment. Considering the enzyme half‐lives of 20–30 hours for CYP3A4 and CYP2J2 and simulated recoveries of CYP3A4 (liver), CYP3A4 (gut), and CYP2J2 activities post‐MBI by ritonavir were 62.4%, 82.1%, and 77.6%, respectively, the reduced dose of rivaroxaban could be maintained for 3 days post‐nirmatrelvir/ritonavir treatment, which was in accordance with recommendations by Marzolini et al. and Hong et al.

As a racemate, S‐warfarin is metabolized by CYP2C9 whereas R‐warfarin is primarily metabolized by CYP3A4. Unlike rivaroxaban, warfarin achieves its PD effect by indirectly interfering with vitamin K dependent clotting factors. Several case reports found significantly decreased INR for patients taking warfarin after the initiation of antiretroviral therapy that included ritonavir, whereas 15 days of oral treatment of ritonavir mildly decreased 24% of S‐warfarin systemic exposure in a clinical study, attributing to the induction of CYP2C9 by ritonavir. Our simulations utilizing racemic warfarin additionally demonstrated up to 3.4‐fold increase of R‐warfarin concentrations during nirmatrelvir/ritonavir treatment, resulting from the inhibition of CYP3A4 by ritonavir. Because R‐warfarin was incorporated as a competitive antagonist of S‐warfarin against vitamin K epoxide reductase in our mechanistic PD model, the elevated R‐warfarin concentration possibly contributed to the remarkable INR decline. As the increase in R‐warfarin concentrations was more profound than the decrease in S‐warfarin concentrations, it reasonably explained the clinically observed discrepancy between PKs and PDs of warfarin when interacting with ritonavir. , Because the appropriate dose of warfarin is determined on an individual basis and adjusted via INR monitoring, only a fixed 5 mg daily dose was simulated in current study without proposing further dose adjustment. Our prospective simulations highlighted the increased tendency of thrombosis during nirmatrelvir/ritonavir treatment, particularly for older adults. As 2 weeks were anticipated for the recovery of INR, our findings underscored the importance for a longer surveillance of INR of warfarin post‐nirmatrelvir/ritonavir therapy.

Taken together, for patients with chronic coagulative abnormalities and diagnosis of COVID‐19, our prospective simulations provided scientific evidence regarding rational and prudent oral anticoagulation management when initiating nirmatrelvir/ritonavir treatment. Clinical studies and/or systematic analyses of real‐world data are warranted to further confirm our recommendations and ensure the safety of anticoagulant therapy while preventing severe illness from COVID‐19 based on nirmatrelvir/ritonavir therapy.

FUNDING

Z.W. received postdoctoral fellowship from the National University of Singapore. E.C.Y.C. received research grant funding from the Joseph Lim Boon Tiong Urology Cancer Research Initiative (Grant: A‐0002678‐01‐00).

CONFLICT OF INTEREST

The authors declared no competing interests for this work.

AUTHOR CONTRIBUTIONS

Z.W. and E.C.Y.C. designed the research and wrote the manuscript. Z.W. performed the research and analyzed the data.

Appendix S1

Click here for additional data file.

Figure S1

Click here for additional data file.

Figure S2

Click here for additional data file.

Table S1

Click here for additional data file.

Table S2

Click here for additional data file.

Table S3

Click here for additional data file.

Table S4

Click here for additional data file.