Pediatric extrapolation strategies issued by health authorities have streamlined pediatric drug development and reduced the unnecessary burden of conducting pediatric clinical studies. In line with these strategies, physiologically based pharmacokinetic (PBPK) models have been utilized extensively for initial dosing regimen and sampling timepoint selection for pediatric studies, as well as dose validation throughout pediatric drug development. Here, the status and challenges of PBPK modeling in pediatric drug development have been summarized by the IQ Pediatric PBPK Working Group. Our work reviews current practices for pediatric PBPK modeling across various therapeutic areas. To enable best practice, we propose an optimized workflow for pediatric PBPK modeling recommendations. Two selected key pediatric PBPK case examples are also described, where modeling impacted the drug label extension to pediatric patients. Moreover, we analyze the current gaps and challenges in our understanding of drug absorption, distribution, metabolism, and elimination in pediatric PBPK model development. Since neonates are the least studied and the most medically fragile, the depth of our understanding of their rapidly evolving physiological processes is limited and so there exist significant modeling gaps which we summarize here. Finally, we provide recommendations, including building a public data repository, leveraging real-world data, and implementing microdose studies for addressing pediatric PBPK modeling challenges.

The inhibition of coagulation factor XI (FXI) presents an attractive approach for anticoagulation as it is not expected to increase the risk of clinically relevant bleeding and is anticipated to be at least as effective as currently available anticoagulants. Fesomersen is a conjugated antisense oligonucleotide that selectively inhibits the expression of FXI. The article describes three clinical studies that investigated the safety, pharmacokinetic (PK), and pharmacodynamic (PD) profiles of fesomersen after subcutaneous (s.c.) injection to healthy participants. The studies included participants from diverse ethnic backgrounds (Caucasian, Japanese, and Chinese). Fesomersen demonstrated good safety and tolerability in all three studies. No major bleeding events were observed. After single-dose s.c. injection, fesomersen was rapidly absorbed into the systemic circulation, with maximum fesomersen-equivalent (fesomersen-eq) concentrations (C-max) in plasma observed within a few hours. After reaching C-max, plasma fesomersen-eq concentrations declined in a biphasic fashion. The PD analyses showed that the injection of fesomersen led to dose-dependent reductions in FXI activity and increases in activated partial thromboplastin time (aPTT). The maximum observed PD effects were reached between Day 15 and 30, and FXI activity and aPTT returned to near-baseline levels by Day 90 after a single dose. The PK/PD profiles after a single injection were similar among the various ethnic groups. Collectively, the study results suggest that fesomersen has a favorable safety profile and predictable and similar PK and PD profiles across Chinese, Japanese, and Caucasian participants.

Patients with kidney failure on hemodialysis (KF-HD) are at high risk for both atherothrombo tic events and bleeding. This Phase IIb study evaluated the dose- response of fesomersen, an inhibitor of hepatic Factor XI expression, versus placebo, for bleeding and atherothrombosis in patients with KF-HD. Patients were randomized to receive fesomersen 40, 80, or 120 mg once-monthly, or matching placebo, for up to 12 months. The primary safety endpoint was a composite of major bleeding and clinically relevant non-major bleeding (MB/ CRNMB). Exploratory endpoints included post-dialysis arterio-venous (AV)-access bleeding, major atherothrombotic events (composite of fatal or non-fatal myocardial infarction, ischemic stroke, acute limb ischemia/major amputation, systemic embolism, symptomatic venous thromboembolism), AV-access thrombosis, and clotting of the hemodialysis circuit. Of 308 participants randomized, 307 received study treatment and were analyzed. Fesomersen led to a dose- dependent and sustained reduction of steady-state median FXI levels by 53.6% (40 mg group), 71.3% (80 mg group), 86.0% (120 mg group), versus 1.9% in the placebo group. MB/CRNMB events occurred in 6.5% (40 mg group), 5.1% (80 mg group), 3.9% (120 mg group), and in 4.0% of those receiving placebo (pooled fesomersen versus placebo P = 0.78). Major atherothrombotic events occurred in 1 patient (1.3%) in each treatment arm. MB/CRNMB bleeding and post- dialysis AV-access bleeding were not related to predicted FXI levels. Lower predicted FXI levels were associated with reductions in hemodialysis circuit clotting (P P = 0.002) and AV-access thrombosis (P P = 0.014). In patients with KF-HD, fesomersen produced a dose-dependent reduction in FXI levels associated with similar rates of major bleeding compared with placebo.

AbstractElinzanetant is a potent and selective dual neurokin‐1 (NK‐1) and ‐3 (NK‐3) receptor antagonist that is currently developed for the treatment of women with moderate‐to‐severe vasomotor symptoms (VMS) associated with menopause. Here, we report the development of a population pharmacokinetic (popPK) model for elinzanetant and its principal metabolites based on an integrated dataset from 366 subjects (including 197 women with VMS) collected in 10 phase I or II studies. The pharmacokinetics of elinzanetant and its metabolites could be well described by the popPK model. Within the investigated dose range of 40–160 mg, the oral bioavailability of elinzanetant was dose independent and estimated to be 36.7%. The clearance of elinzanetant was estimated to be 7.26 L/h and the central and peripheral distribution volume were 23.7 and 168 L. No intrinsic or extrinsic influencing factors have been identified in the investigated population other than the effect of a high‐fat breakfast on the oral absorption of elinzanetant. The popPK model was then coupled to a pharmacodynamic model to predict occupancies of the NK‐1 and NK‐3 receptors. After repeated once‐daily administration of the anticipated therapeutic dose of 120 mg elinzanetant, the model‐predicted median receptor occupancies are >99% for NK‐1 and >94.8% for NK‐3 during day and night‐time, indicating sustained and near‐complete inhibition of both target receptors during the dosing interval.

Inhibition of tissue factor pathway inhibitor (TFPI) is a potential new mode of action to achieve haemostasis in haemophilia A and B patients.

Obesity remains a US national health crisis and a growing concern worldwide. Concerningly, individuals who are obese are at an increased risk for comorbid diseases that include, but are not limited to, hypertension, diabetes, cardiovascular disease, and cancer. Beyond the risk for developing these conditions, obesity may also impact the pharmacological activity of the therapies being used to treat them and other disease states. The pharmacokinetics (PK), pharmacodynamics (PD), safety, and efficacy of therapies, both currently marketed and under clinical development, may be directly impacted by the physiological alterations that occur secondary to the occurrence of chronic excess body weight. The increased prevalence of this disease should not be ignored. Both private and federal institutions involved in drug research and development should consider, as appropriate, a greater inclusion of individuals who are obese in clinical trials throughout the entirety of drug development, and leverage the available PK, PD, safety, and efficacy data to make more informed dosing recommendations.

The PATENT‐CHILD study investigated riociguat in children aged ≥ 6 to <18 years with pulmonary arterial hypertension (PAH) treated with tablets or an oral pediatric suspension based on bodyweight‐adjusted dosing of up to 2.5 mg three times daily. PATENT‐CHILD demonstrated an acceptable riociguat safety profile and individual plasma concentrations in pediatric patients were consistent with those in adult patients.

Anticoagulant treatment of pediatric cancer–associated venous thromboembolism (VTE) has not been prospectively evaluated. Management of anticoagulation for cancer-associated VTE is often challenged by drug interactions and treatment interruptions. A total of 56 of the 500 children (11.2%) with VTE who participated in the recent EINSTEIN-Jr randomized study had cancer (hematologic malignancy, 64.3%, solid malignant tumor, 35.7%). Children were allocated to either therapeutic-dose bodyweight-adjusted oral rivaroxaban (n=40) or standard anticoagulation with heparins, with or without vitamin K antagonists (n=16) and received a median of 30 concomitant medications. Based on sparse blood sampling at steady-state, pharmacokinetic (PK) parameters of rivaroxaban were derived using population PK modeling. During the 3 months of treatment, no recurrent VTE or major bleeding occurred (95% confidence interval, 0.0%-6.4%), and 3-month repeat imaging showed complete or partial vein recanalization in 20 and 24 of 52 evaluable children (38.5% and 46.2%, respectively). Anticoagulant treatment was interrupted 70 times in 26 (46.4%) children because of thrombocytopenia, invasive procedures, or adverse events, for a mean individual period of 5.8 days. Anticoagulant therapy was resumed in therapeutic doses and was not associated with thrombotic or bleeding complications. Rivaroxaban exposures were within the adult exposure range and similar to those observed in children with VTE who did not have cancer-associated VTE. Rivaroxaban and standard anticoagulants appeared safe and efficacious and were associated with reduced clot burden in most children with cancer-associated VTE, including those who had anticoagulant treatment interruptions. Rivaroxaban exposures were within the adult exposure range despite significant polypharmacy use. This trial was registered at www.clinicaltrials.gov as #NCT02234843.

The antiparasitic drug nifurtimox was approved in the USA in 2020 for the treatment of patients with Chagas disease aged less than 18 years and weighing at least 2.5 kg, based on outcomes from the phase 3 CHICO study. Accordingly, pediatric patients with Chagas disease take nifurtimox thrice daily with food at one of two body weight-adjusted dose ranges. We investigated possible relationships between pharmacokinetic (PK) data, and pharmacodynamic efficacy and safety data collected in an analysis population of 111 participants in CHICO, using a published population PK model to estimate nifurtimox exposure at the patient level. Pediatric exposure to nifurtimox was benchmarked against levels of nifurtimox exposure known to be effective in adults with Chagas disease. Given the complex dosing regimen for nifurtimox, we also modeled nifurtimox exposure associated with simpler dosing strategies. We found no relationship between exposure to nifurtimox and efficacy measures (e.g., serological response to treatment), or between exposure and safety outcomes (including typical adverse events, e.g., headache, decreased appetite, nausea/vomiting). The analysis population appeared to represent the overall CHICO population based on the similarity of their baseline characteristics and the profiles of adverse events in the two groups. Modeled exposure based on the dosing regimen in CHICO was within the reference range derived from phase 1 data in adults. The relationship between nifurtimox exposure and cure is complex; a simplified pediatric dosing regimen is unlikely to be beneficial.

Thrombosis remains an important complication for children with single-ventricle physiology following the Fontan procedure, and effective thromboprophylaxis is an important unmet medical need. To obviate conventional dose-finding studies and expedite clinical development, a rivaroxaban dose regimen for this indication was determined using a model-informed drug development approach. A physiologically based pharmacokinetic rivaroxaban model was used to predict a pediatric dosing regimen that would produce drug exposures similar to that of 10 mg once daily in adults. This regimen was used in an open-label, multicenter phase III study, which investigated the use of rivaroxaban for thromboprophylaxis in post-Fontan patients 2 to 8 years of age. The pharmacokinetics (PK) of rivaroxaban was assessed in part A (n = 12) and in part B (n = 64) of the UNIVERSE study. The safety and efficacy in the rivaroxaban group were compared to those in the acetylsalicylic acid group for 12 months. Pharmacodynamic end points were assessed in both parts of the study. Rivaroxaban exposures achieved in parts A and B were similar to the adult reference exposures. Prothrombin time also showed similarity to the adult reference. Exposure-response analysis did not identify a quantitative relationship between rivaroxaban exposures and efficacy/safety outcomes within the observed exposure ranges. A body weight-based dose regimen selected by physiologically based pharmacokinetic modeling was shown in the UNIVERSE study to be appropriate for thromboprophylaxis in the post-Fontan pediatric population. Model-based dose selection can support pediatric drug development and bridge adult dose data to pediatrics, thereby obviating the need for dose-finding studies in pediatric programs.

AbstractRiociguat, a soluble guanylate cyclase stimulator, is approved for treatment of adults with pulmonary arterial hypertension (PAH). The safety, tolerability, and pharmacokinetics (PK) of oral riociguat in a pediatric population with PAH was assessed in PATENT–CHILD (NCT02562235), a multicenter, single‐arm, 24‐week, open‐label, Phase 3 study. Patients aged 6–17 years in World Health Organization functional class (WHO‐FC) I–III treated with stable endothelin receptor antagonists and/or prostacyclin analogs received riociguat equivalent to 0.5–2.5 mg three times daily in adults, as either oral pediatric suspension or tablets, based on bodyweight. Primary outcomes were safety, tolerability, and PK of riociguat. Twenty‐four patients (mean age 12.8 years), 18 of whom were in WHO‐FC II, were enrolled. Adverse events (AEs), mostly mild or moderate, were reported in 20 patients (83%). Four patients (17%) experienced a serious AE; all resolved by study end and two (8%) were considered study‐drug related. Hypotension was reported in three patients and hemoptysis in one (all mild/moderate intensity). Riociguat plasma concentrations in pediatric patients were consistent with those published in adult patients. From baseline to Week 24, mean ± standard deviation increase in 6‐minute walking distance was 23 ± 69 m (n = 19), and mean decrease in NT‐proBNP was –66 ± 585 pg/ml (n = 14). There was no change in WHO‐FC. Two patients experienced clinical worsening events of hospitalization for right heart failure. PK results confirmed a suitable riociguat dosing strategy for pediatric patients with PAH. The data suggest an acceptable safety profile with potential efficacy signals.

Abstract 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.

Nifurtimox (LAMPIT) has been used for decades for the treatment of Chagas disease, a chronic and potentially life-threatening disease caused by the parasite Trypanosoma cruzi. The pharmacokinetic (PK) information on nifurtimox in humans derived from controlled clinical studies is very limited. The objective was to investigate and compare the population PK of nifurtimox in adult and pediatric patients with Chagas disease to confirm the clinical dosing regimen in children, which was based on allometric approaches using the concept that a dose-equivalent exposure would reach equivalent antiparasitic efficacy as in adults. The resulting adult model adequately described the PK in adults. Significant predictors of the availability in PK were food intake, tablet formulation (fast- vs slow-dissolution tablet), study, and body weight. As the resulting adult model could not adequately predict the sparse sampled pediatric patient data, these data were analyzed separately to derive exposure estimates for comparison with adult exposure. In the population PK model for pediatric patients, significant covariates were body weight and age. As compared to adults, children aged >2 years were estimated to have 50.6% higher apparent clearance. No hints of dose nonlinearity were observed in a dose range of 30 to 240 mg single dose in adults and 15 to 300 mg 3 times daily (8-20 mg/kg) in children. Altogether, this study retroactively showed that the current mg/kg dosing regimen in children reached similar exposure as in adults receiving an 8 mg/kg total daily dose.

Model building is one of the fundamental knowledge generation concepts in science and medicine. We introduce the concept and its application in Translational Medicine with a series of examples. The Hodgkin & Huxley model for the neural action potential illustrates the power of models in the identification of previously unknown structural properties in biology. Compartmental and Physiology-Based Pharmacokinetics demonstrate how data can be interpreted by mathematical abstraction such that predictions for untested clinical settings are enabled. Based on a computational model of blood coagulation, the mechanistic interaction of different pharmacological principles for anticoagulation is predicted. Quantitative computer simulations are used to inform complex dosing schedules ensuring patient safety when treatment paradigms are switched. Finally, we introduce the framework of model-informed drug development, a novel paradigm in pharmaceutical industry and regulatory practice where direct experimental and clinical evidence is complemented and even substituted by model-based derivations.

The non–vitamin K antagonist oral anticoagulant rivaroxaban is used in several thromboembolic disorders. Rivaroxaban is eliminated via both metabolic degradation and renal elimination as unchanged drug. Therefore, renal and hepatic impairment may reduce rivaroxaban clearance, and medications inhibiting these clearance pathways could lead to drug‐drug interactions. This physiologically based pharmacokinetic (PBPK) study investigated the pharmacokinetic behavior of rivaroxaban in clinical situations where drug clearance is impaired. A PBPK model was developed using mass balance and bioavailability data from adults and qualified using clinically observed data. Renal and hepatic impairment were simulated by adjusting disease‐specific parameters, and concomitant drug use was simulated by varying enzyme activity in virtual populations (n = 1000) and compared with pharmacokinetic predictions in virtual healthy populations and clinical observations. Rivaroxaban doses of 10 mg or 20 mg were used. Mild to moderate renal impairment had a minor effect on area under the concentration‐time curve and maximum plasma concentration of rivaroxaban, whereas severe renal impairment caused a more pronounced increase in these parameters vs normal renal function. Area under the concentration‐time curve and maximum plasma concentration increased with severity of hepatic impairment. These effects were smaller in the simulations compared with clinical observations. AUC and Cmax increased with the strength of cytochrome P450 3A4 and P‐glycoprotein inhibitors in simulations and clinical observations. This PBPK model can be useful for estimating the effects of impaired drug clearance on rivaroxaban pharmacokinetics. Identifying other factors that affect the pharmacokinetics of rivaroxaban could facilitate the development of models that approximate real‐world pharmacokinetics more accurately.

Background: Riociguat is a first-in-class soluble guanylate cyclase stimulator for which preclinical data suggested improvements in cystic fibrosis transmembrane conductance regulator (CFTR) function. Methods: This international, multicenter, two-part, Phase II study of riociguat enrolled adults with cystic fibrosis (CF) homozygous for Phe508del CFTR . Part 1 was a 28-day, randomized, double-blind, placebocontrolled study in participants not receiving CFTR modulator therapy. Twenty-one participants were randomized 1:2 to placebo or oral riociguat (0.5 mg three times daily [tid] for 14 days, increased to 1.0 mg tid for the subsequent 14 days). The primary and secondary efficacy endpoints were change in sweat chloride concentration and percent predicted forced expiratory volume in 1 second (ppFEV 1 ), respectively, from baseline to Day 14 and Day 28 with riociguat compared with placebo. Results: Riociguat did not alter CFTR activity (change in sweat chloride) or lung function (change in ppFEV 1 ) at doses up to 1.0 mg tid after 28 days. The most common drug-related adverse event (AE) was headache occurring in three participants (21%); serious AEs occurred in one participant receiving riociguat (7%) and one participant receiving placebo (14%). This safety profile was consistent with the underlying disease and the known safety of riociguat for its approved indications. Conclusions: The Rio-CF study was terminated due to lack of efficacy and the changing landscape of CF therapeutic development. The current study , within its limits of a small sample size, did not provide evidence that riociguat could be a valid treatment option for CF. Clinical trial registration number: NCT02170025. (c) 2021 The Authors. Published by Elsevier B.V. on behalf of European Cystic Fibrosis Society. This is an open access article under the CC BY-NC-ND license ( http://creativecommons.org/licenses/by-nc-nd/4.0/ )

Development and guidance of dosing schemes in children have been supported by physiology-based pharmacokinetic (PBPK) modeling for many years. PBPK models are built on a generic basis, where compound- and system-specific parameters are separated and can be exchanged, allowing the translation of these models from adults to children by accounting for physiological differences. Owing to these features, PBPK modeling is a valuable approach to support clinical decision making for dosing in children. In this analysis, we evaluate pediatric PBPK models for 10 small-molecule compounds that were applied to support clinical decision processes at Bayer for their predictive power in different age groups. Ratios of PBPK-predicted to observed PK parameters for the evaluated drugs in different pediatric age groups were estimated. Predictive performance was analyzed on the basis of a 2-fold error range and the bioequivalence range (ie, 0.8 ≤ predicted/observed ≤ 1.25). For all 10 compounds, all predicted-to-observed PK ratios were within a 2-fold error range (n = 27), with two-thirds of the ratios within the bioequivalence range (n = 18). The findings demonstrate that the pharmacokinetics of these compounds was successfully and adequately predicted in different pediatric age groups. This illustrates the applicability of PBPK for guiding dosing schemes in the pediatric population.

Abstract Background/Introduction Vericiguat is a soluble guanylate cyclase (sGC) stimulator, like riociguat and nelociguat, and entered clinical development in 2012. Before entering Phase 2, pharmacokinetics (PK) and pharmacodynamics (PD) of vericiguat had been studied in healthy volunteers only, whereas riociguat and nelociguat had also been studied in patients with pulmonary hypertension (PH) and left ventricular dysfunction (LVD) or biventricular chronic heart failure (HF). We hypothesised that integrating all PK/PD data from these compounds into population PK/PD (popPK/PD) and physiology-based PK (PBPK) models could be used to predict optimal and safe dose ranges of vericiguat for Phase 2b studies in patients with worsening chronic HF. This novel bridging approach was applied in one of several translational stages to accelerate the development of vericiguat (Figure 1). Purpose We used prior knowledge from other sGC stimulators in a combined PK/PD and PBPK modelling approach to directly initiate Phase 2b studies of vericiguat in patients after Phase 1 studies in healthy volunteers. Methods PK, heart rate (HR) and systemic vascular resistance (SVR) data for vericiguat, nelociguat and riociguat were used to calculate PK/PD slopes of linear models, corrected with fraction unbound percentages (2.2%, 3.6% and 3.9%, respectively), to compare potency relative to riociguat based on unbound concentrations. PK estimates for nelociguat and riociguat were derived using population PK modelling (NONMEM) from patient studies with sparse PK sampling. PBPK models informed by preclinical physicochemical and PK data as well as clinical data for vericiguat were used to predict vericiguat PK in patients with HF (PK-Sim). Exposure–response data for riociguat in patients indicated the optimal range of PD responses for vericiguat (blood pressure for safety and cardiac index for efficacy). Results Vericiguat and nelociguat had lower potency than riociguat when comparing PK/PD slopes for HR and SVR (slope ratios of 0.23–0.32 for vericiguat and 0.33–0.47 for nelociguat). Plasma concentrations of vericiguat would need to be ∼3.6 times that of riociguat for equivalent responses. In patients with PH and LVD the optimal plasma concentration range for riociguat was ∼10–100 μg/l in exposure–response and safety studies, which translates to a target exposure range of ∼90–900 μg/l for vericiguat in patients with HF. PBPK modelling showed that vericiguat 2.5 mg and 10 mg would cover the target exposure range and that 1.25 mg would be a “non-effective” dose level with respect to haemodynamics. Conclusions Our novel translational approach combining popPK/PD analyses of other sGC stimulators with PBPK modelling enabled vericiguat to move directly from Phase 1 to Phase 2b, reducing development time by ∼2 years. PK and safety results from Phase 2b (SOCRATES-REDUCED) and Phase 3 (VICTORIA) trials confirmed that use of this translational approach to predict dose ranges of vericiguat was successful. Funding Acknowledgement Type of funding sources: Private company. Main funding source(s): Funding for this research was provided by Bayer AG, Berlin, Germany Figure 1

This study compared simulations of a physiologically based pharmacokinetic (PBPK) model implemented for cyclosporine with drug levels from therapeutic drug monitoring to evaluate the predictive performance of a PBPK model in a clinical population. Based on a literature search model parameters were determined. After calibrating the model using the pharmacokinetic profiles of healthy volunteers, 356 cyclosporine trough levels of 32 renal transplant outpatients were predicted based on their biometric parameters. Model performance was assessed by calculating absolute and relative deviations of predicted and observed trough levels. The median absolute deviation was 6 ng/ml (interquartile range: 30 to 31 ng/ml, minimum = −379 ng/ml, maximum = 139 ng/ml). 86% of predicted cyclosporine trough levels deviated less than twofold from observed values. The high intra-individual variability of observed cyclosporine levels was not fully covered by the PBPK model. Perspectively, consideration of clinical and additional patient-related factors may improve the model’s performance. In summary, the current study has shown that PBPK modeling may offer valuable contributions for pharmacokinetic research in clinical drug therapy.

Rivaroxaban has been investigated in the EINSTEIN-Jr program for the treatment of acute venous thromboembolism (VTE) in children aged 0 to 18 years and in the UNIVERSE program for thromboprophylaxis in children aged 2 to 8 years with congenital heart disease after Fontan-procedure. Physiologically-based pharmacokinetic (PBPK) and population pharmacokinetic (PopPK) modeling were used throughout the pediatric development of rivaroxaban according to the learn-and-confirm paradigm. The development strategy was to match pediatric drug exposures to adult exposure proven to be safe and efficacious. In this analysis, a refined pediatric PopPK model for rivaroxaban based on integrated EINSTEIN-Jr data and interim PK data from part A of the UNIVERSE phase III study was developed and the influence of potential covariates and intrinsic factors on rivaroxaban exposure was assessed. The model adequately described the observed pediatric PK data. PK parameters and exposure metrics estimated by the PopPK model were compared to the predictions from a previously published pediatric PBPK model for rivaroxaban. Ninety-one percent of the individual post hoc clearance estimates were found within the 5th to 95th percentile of the PBPK model predictions. In patients below 2 years of age, however, clearance was underpredicted by the PBPK model. The iterative and integrative use of PBPK and PopPK modeling and simulation played a major role in the establishment of the bodyweight-adjusted rivaroxaban dosing regimen that was ultimately confirmed to be a safe and efficacious dosing regimen for children aged 0 to 18 years with acute VTE in the EINSTEIN-Jr phase III study.