Population pharmacokinetic analysis of rivaroxaban in children and comparison to prospective physiologically-based pharmacokinetic predictions

Optimization of wet granulation process for manufacturing Rivaroxban generic immediate-release tablets using PBPK modeling and simulations

Granulation is the critical process for the pharmaceutical development of poorly water-soluble drug products. Poorly formulated products have challenges in dissolution and bioequivalence studies. Rivaroxaban (RXB) is a poorly soluble drug and has 66% fasting bioavailability at a high strength of 20 mg. Establishing the bioequivalence between test and reference products for high strength requires comparative dissolution profiles and bioequivalence. Improper granulation products and the rest of the batches failed in virtual bioequivalence. The present study provided insight into the optimization of the wet granulation process for manufacturing RXB generic immediate-release tablets using PBPK modeling and simulations. Furthermore, PBPK models are not only useful for formulation optimization but also for process optimization during pharmaceutical product development.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00249-6.

Potential treatment option of rivaroxaban for breastfeeding women: A case series

The use of direct oral anticoagulants (DOACs) in breastfeeding women is currently challenging due to limited safety data for breastfeeding infants, and there have been no previous studies on the drug concentration in breastfeeding infants. We treated 2 patients (one case was twin pregnancy) with venous thromboembolisms in breastfeeding women administered rivaroxaban at our institution. Blood samples from the mothers and breastmilk samples were collected at time 0 and 2 h after the rivaroxaban administration, breastfeeding was conducted 2 h after the rivaroxaban administration, and blood samples from the infants were collected 2 h after breastfeeding (4 h after maternal rivaroxaban administration). The milk-to-plasma (M:P) ratios were 0.27 in Case 1 and 0.32 in Case 2. The estimated relative infant dose (RID) was 0.82 % in Case 1 Children 1 and 2, and 1.27 % in Case 2. The rivaroxaban concentration in the infant plasma was below the lower limit of quantification in all infants. In addition, even in the high-exposure case simulation based on 5 days of breastfeeding in Case 2, the infant plasma concentration level was below the lower limit of quantification. At 3 months of follow-up, breastfeeding was continued, and all infants grew and developed without any health problems including bleeding events. The current case series showed that there were no pharmacokinetic or clinical concerns for breastfeeding women or breastfed infants, and provides support for rivaroxaban as a safe treatment option for these patients.

An extension of biorelevant fed-state dissolution tests to clinical pharmacokinetics - A study on gastrointestinal factors influencing rivaroxaban exposure and efficacy in atrial fibrillation patients

A direct oral anticoagulant rivaroxaban fails to prevent stroke and systemic embolism in one-to-several percent of patients with nonvalvular atrial fibrillation (NVAF), but the reasons are unknown. The study used semi-mechanistic in vitro-in vivo prediction (IVIVP) modeling to explore the reasons for ineffective thrombosis prevention in NVAF patients. Steady-state drug concentrations in plasma were measured at 0 h (Ctrough), 3 h (C3h), and 12 h post-dosing in thirty-four patients treated with 20 mg rivaroxaban daily. The clinical data were compared against "virtual twins" generated with a novel IVIVP model that combined drug dissolution modeling, mechanistic description of gastric drug transit, and population pharmacokinetics defining the variability of drug disposition. The nonresponders had significantly lower C3h and Ctrough than the responders (p < 0.001) and the covariates included in the population pharmacokinetic submodel did not fully explain this difference. Simulations involving varied gastrointestinal parameters in the "virtual twins" revealed that lower small intestinal effective permeability (Peff), rather than a slower stomach emptying rate, could explain low rivaroxaban exposure in the nonresponders. IVIVP modeling was effectively used for exploring pharmacotherapy failure. Low Peff, found as a major determinant of ineffective rivaroxaban treatment, encourages further research to find (pato)physiological factors influencing suboptimal absorption.

Is a Lower Dose of Rivaroxaban Required for Asians? A Systematic Review of a Population Pharmacokinetics and Pharmacodynamics Analysis of Rivaroxaban

Rivaroxaban has been widely used to prevent and treat various thromboembolic diseases for more than a decade. However, whether a lower dose of rivaroxaban is required for Asians is still debatable. This review aimed to explore the potential ethnic difference in pharmacokinetic/pharmacodynamic (PK/PD) characteristics between Asians and Caucasians. A systematic search was conducted and twenty-four studies were identified, of which 10 were conducted on Asian adults, 11 on predominantly Caucasian adults, and 3 on Caucasian pediatrics. The apparent clearance (CL/F) of rivaroxaban in Caucasian adults with non-valvular atrial fibrillation (6.45-7.64 L/h) was about 31-43% higher than that in Asians (4.46-5.98 L/h) taking 10~20 mg rivaroxaban every 24 h. Moreover, there was no obvious difference in CL/F among Japanese, Chinese, Thai, and Irani people. Regarding PK/PD relationship, prothrombin time was linked to rivaroxaban concentration in a linear or near-linear manner, and Factor Xa activity was linked with the E_max model. The exposure-response relationship was comparable between Asians and Caucasians. Renal function has a significant influence on CL/F, and no covariate was recognized for exposure-response relationship. In conclusion, a lower dose of rivaroxaban might be required for Asians, and further studies are warranted to verify this ethnic difference to facilitate optimal dosing regimens.

In-Depth Analysis of Physiologically Based Pharmacokinetic (PBPK) Modeling Utilization in Different Application Fields Using Text Mining Tools

In the past decade, only a small number of papers have elaborated on the application of physiologically based pharmacokinetic (PBPK) modeling across different areas. In this review, an in-depth analysis of the distribution of PBPK modeling in relation to its application in various research topics and model validation was conducted by text mining tools. Orange 3.32.0, an open-source data mining program was used for text mining. PubMed was used for data retrieval, and the collected articles were analyzed by several widgets. A total of 2699 articles related to PBPK modeling met the predefined criteria. The number of publications per year has been rising steadily. Regarding the application areas, the results revealed that 26% of the publications described the use of PBPK modeling in early drug development, risk assessment and toxicity assessment, followed by absorption/formulation modeling (25%), prediction of drug-disease interactions (20%), drug-drug interactions (DDIs) (17%) and pediatric drug development (12%). Furthermore, the analysis showed that only 12% of the publications mentioned model validation, of which 51% referred to literature-based validation and 26% to experimentally validated models. The obtained results present a valuable review of the state-of-the-art regarding PBPK modeling applications in drug discovery and development and related fields.

Assessment of Aging-Related Function Variations of P-gp Transporter in Old-Elderly Chinese CHF Patients Based on Modeling and Simulation

BACKGROUND AND OBJECTIVES

P-glycoprotein (P-gp) is one of the most intensely studied transporters owing to its broad tissue distribution and substrate specificity. Existing research suggests that the risk of systemic exposure to dabigatran etexilate (DABE) and digoxin, two P-gp probe substrates in vivo, has significantly increased in elderly patients. We applied a model-based quantitative pharmacological approach to assess aging-related P-gp changes in the Chinese old-elderly population.

METHODS

Population pharmacokinetic (PopPK) modeling was first performed using clinical pharmacokinetic data to explore the effect of age on the pharmacokinetic characteristics of dabigatran (DAB, the active principle of DABE) and digoxin in elderly Chinese patients. Corresponding physiologically based pharmacokinetic (PBPK) models were established to further explain the elevated systemic exposure to these two drugs. Eventually, standard dosing regimens of DABE and digoxin were assessed in Chinese old-elderly patients with chronic heart failure (CHF) with different stages of renal impairment.

RESULTS

PopPK analysis suggested that age as a covariate had an additional effect on the apparent clearance of these two drugs after correcting for creatinine clearance. PBPK simulation results suggested that disease-specific pathophysiological changes could explain DAB exposure in the young elderly. In the elderly population, 17.1% of elevated DAB exposure remained unexplained, and 25.5% of the reduced P-gp function associated with aging was ultimately obtained using sensitivity analysis. This value was further validated using digoxin data obtained by PBPK modeling. The simulation results suggest that CHF patients with advanced age and moderate-to-severe renal impairment require heightened vigilance for elevated exposure risk during the use of DABE and digoxin.

CONCLUSIONS

Aging might be a significant risk factor for elevated systemic exposure to DAB and digoxin by reducing P-gp-mediated efflux in the Chinese old elderly population.

Model‐informed bridging of rivaroxaban doses for thromboprophylaxis in pediatric patients aged 9 years and older with congenital heart disease

Rivaroxaban is approved in various regions for the treatment of acute venous thromboembolism (VTE) in children aged between 0 and 18 years and was recently investigated for thromboprophylaxis in children aged between 2 and 8 years (with body weights <30 kg) with congenital heart disease who had undergone the Fontan procedure. In the absence of clinical data, rivaroxaban doses for thromboprophylaxis in post‐Fontan children aged 9 years and older or ≥30 kg were derived by a bridging approach that used physiologically‐based pharmacokinetic (PBPK) and population pharmacokinetic (popPK) models based on pharmacokinetic (PK) data from 588 pediatric patients and from adult patients who received 10 mg once daily for thromboprophylaxis after major orthopedic surgeries as a reference. Both models showed a tendency toward underestimating rivaroxaban exposure in post‐Fontan patients aged between 2 and 5 years but accurately described rivaroxaban PK in post‐Fontan patients aged between 5 and 8 years. Under the assumption that hepatic function is not impaired in post‐Fontan patients, PBPK and popPK simulations indicated that half of the rivaroxaban doses for the same body weight given to pediatric patients treated for acute VTE would yield in pediatric post‐Fontan patients exposures similar to the exposure observed in adult patients receiving 10 mg rivaroxaban once daily for thromboprophylaxis. Simulation‐derived doses (7.5 mg rivaroxaban once daily for body weights 30–<50 kg and 10 mg once daily for body weights ≥50 kg) were therefore included in the recent US label of rivaroxaban for thromboprophylaxis in children aged 2 years and older with congenital heart disease who have undergone the Fontan procedure.

Rivaroxaban Pharmacokinetics in Obese Subjects: A Systematic Review

INTRODUCTION

Venous thromboembolism (VTE) is a leading cause of morbidity and mortality globally. The direct oral anticoagulants, including rivaroxaban, are relatively novel therapeutic options in the treatment and prevention of VTE. There is a conflicting and inconclusive evidence surrounding the pharmacokinetics (PK) of rivaroxaban in patients with VTE who are obese.

OBJECTIVES

We conducted a systematic review to provide an overview, and to synthesize the available evidence in the current literature pertaining to rivaroxaban PK in obese subjects who are healthy or diseased.

METHODS

The PubMed, Embase, ScienceDirect, Rayyan, and Cochrane Library databases were systematically searched from 1 May 2021 through 28 February 2022. Studies investigating rivaroxaban PK in adult obese subjects were included in the review. Pertinent data, including anthropometric parameters, rivaroxaban dosage regimen, PK parameters, PK model, and outcome measures were extracted. Reference values of rivaroxaban PK parameters in the general population were used for comparison purposes. The review protocol was registered in the PROSPERO database (CRD42020177770).

RESULTS

In the 11 studies included in this systematic review, over 7140 healthy or diseased subjects received rivaroxaban therapy, with varying clinical indications in the diseased population. The reported PK parameters of rivaroxaban in obese subjects compared with reference values in the general population were variable. The reported values of the volume of distribution (V_d) among obese subjects (73.4-82.8 L) fell within the range of values reported/calculated for the general population (59.4-104 L), assuming complete bioavailability. However, some of the reported values of clearance (CL) in obese subjects (7.86-16.8 L.h^-1) do not fall within the range of values reported/calculated for the general population (5.57-11.3 L.h^-1). The reported maximum plasma concentrations in obese subjects versus the general population following a 10 mg dose were 149 vs. 143-180 µg.L^-1, and following a 20 mg dose were 214-305 vs. 299-360 µg.L^-1, respectively. The area under the plasma concentration versus time curves (AUC) over different intervals in obese subjects versus the general population following a 10 mg dose were 1155 (AUC from time zero to infinity [AUC_∞]) vs. 1029 (AUC_∞) µg.h.L^-1; and 1204-2800 (AUC from time zero to 24 h [AUC₂₄]) vs. 3200 (AUC₂₄) µg.h.L^-1, respectively, following a 20 mg dose. The reported values of half-life and time to reach the maximum plasma concentration in obese subjects versus the general population were not consistent across studies.

CONCLUSION

Variable changes and inconsistencies in different rivaroxaban PK parameters were reported in obese subjects. Further well-designed studies are warranted to better characterize the PK and clinical outcomes of rivaroxaban in subjects with obesity.