Pharmacokinetics, pharmacodynamics, and safety of fesomersen, a novel antisense inhibitor of factor XI, in healthy Chinese, Japanese, and Caucasian volunteers

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 (Cmax) in plasma observed within a few hours. After reaching Cmax, 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.

A Phase II randomized controlled trial evaluated antithrombotic treatment with fesomersen in patients with kidney failure on hemodialysis

Patients with kidney failure on hemodialysis (KF-HD) are at high risk for both atherothrombotic 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), and 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 = 0.002) and AV-access thrombosis (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. REGISTRATION: URL: https://www.clinicaltrials.gov; unique identifier: NCT04534114.

The Importance of Assessing Drug Pharmacokinetics and Pharmacodynamics in the Obese Population During Drug Development

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.

Longitudinal observations of TFPI levels in paediatric Haemophilia A patients

INTRODUCTION

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

AIM

Knowledge about potential developmental changes of TFPI levels during childhood are a prerequisite to translate adult doses of TFPI inhibitors to doses in paediatric patients.

METHODS

In this study we present longitudinal data for total TFPI concentrations (TFPI-T) and TFPI activity (TFPI-A) from 48 paediatric Haemophilia A patients in the age range from 3 to 18 years (2-12 observations per patient).

RESULTS

TFPI-T and TFPI-A tend to decrease over age during childhood. Lowest values were observed between 12 and <18 years. On average, TFPI-T and TFPI-A were lower in adolescent haemophilia patients than in adult haemophilia patients.

CONCLUSION

In summary, the presented information on TFPI levels in children adds to the current knowledge of developmental haemostasis and it can be helpful in evaluating how children respond to haemophilia treatment including the new class of anti-TFPI compounds.

Population pharmacokinetics of riociguat in a pediatric population (aged ≥ 6 years) with pulmonary arterial hypertension

BACKGROUND

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.

METHODS

Using the data set from PATENT-CHILD and building on existing population pharmacokinetic (PK) models for riociguat and its major metabolite (M1) in adults with PAH, a coupled riociguat-M1 PK model was developed. The final model developed incorporated a one-compartment model for riociguat, coupled to a one-compartment model for M1, allowing for presystemic formation of M1. It included allometric scaling exponents for bodyweight.

RESULTS

Apparent clearance of riociguat was similar in children and adult patients with PAH (median [interquartile range] 2.20 [1.75-3.44] and 2.08 L/h [1.55-2.97]). Factors contributing to lower PK exposure were lower riociguat maintenance dose in PATENT-CHILD, and a higher riociguat clearance in some adolescent patients, compared with adult patients. No effects of formulation, sex, or age on riociguat PK were observed. An exploratory PK/pharmacodynamics analysis found the increase in 6-min walking distance in pediatric patients treated with riociguat was not related to riociguat PK.

CONCLUSIONS

Body size is the main determinant of PK in growing children, and the model supports clinical data that, for children weighing < 50 kg, a bodyweight-adjusted dose of riociguat should be used to achieve a similar exposure to that observed in adults with PAH.

Riociguat in children with pulmonary arterial hypertension:the PATENT‐CHILD study

Riociguat, 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.

Nifurtimox for Treatment of Chagas Disease in Pediatric Patients: the Challenges of Applying Pharmacokinetic-Pharmacodynamic Principles to Dose Finding

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.

Graphical abstract

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.

Population Pharmacokinetics of Nifurtimox in Adult and Pediatric Patients with Chagas' Disease

Nifurtimox (LAMPIT®) has been used for decades for the treatment of Chagas' Disease, a chronic and potentially life-threatening disease caused by the parasite Trypanonosma Cruzi. The pharmacokinetics (PK) information on nifurtimox in humans derived from controlled clinical studies is very limited. The objective was to investigate and compare the Population PK (PopPK) 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 (WT). 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 PopPK model for pediatric patients, significant covariates were WT and age. As compared to adults, children older than 2 years were estimated to have 50.6 % higher apparent clearance (CL/F). 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 to 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. This article is protected by copyright. All rights reserved.

Leveraging translational approaches for accelerated clinical development of vericiguat

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

Dosing Regimen Prediction and Confirmation with Rivaroxaban for Thromboprophylaxis in Children after the Fontan Procedure: Insights from the Phase III UNIVERSE Study

Thrombosis remains an important complication for children with single ventricle physiology post-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 utilizing a model-informed drug development approach. A physiologically based pharmacokinetic (PBPK) 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 3 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 UNIVERSE. The safety and efficacy in the rivaroxaban group were compared to those in the acetylsalicylic acid group for 12 months. Pharmacodynamic endpoints were assessed in both parts of the study. Rivaroxaban exposures achieved in Part 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 PBPK 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. This article is protected by copyright. All rights reserved.

Riociguat for the treatment of Phe508del homozygous adults with cystic fibrosis

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, placebo-controlled 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₁), 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₁) 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.

PK/PD modeling of FXI antisense oligonucleotides to bridge the dose-FXI activity relation from healthy volunteers to end-stage renal disease patients

IONIS-FXIRX (BAY2306001) is an antisense oligonucleotide that inhibits the synthesis of coagulation factor XI (FXI) and has been investigated in healthy volunteers and patients with end-stage renal disease (ESRD). FXI-LICA (BAY2976217) shares the same RNA sequence as IONIS-FXIRX but contains a GalNAc-conjugation that facilitates asialoglycoprotein receptor (ASGPR)-mediated uptake into hepatocytes. FXI-LICA has been studied in healthy volunteers and is currently investigated in patients with ESRD on hemodialysis. We present a model-informed bridging approach that facilitates the extrapolation of the dose-exposure-FXI relationship from IONIS-FXIRX to FXI-LICA in patients with ESRD and, thus, supports the selection of FX-LICA doses being investigated in patients with ESRD. A two-compartment pharmacokinetic (PK) model, with mixed first- and zero-order subcutaneous absorption and first-order elimination, was combined with an indirect response model for the inhibitory effect on the FXI synthesis rate via an effect compartment. This PK/pharmacodynamic model adequately described the median trends, as well as the interindividual variabilities for plasma drug concentration and FXI activity in healthy volunteers of IONIS-FXIRX and FXI-LICA, and in patients with ESRD of IONIS-FXIRX . The model was then used to predict dose-dependent steady-state FXI activity following repeat once-monthly doses of FXI-LICA in a virtual ESRD patient population. Under the assumption of similar ASGPR expression in patients with ESRD and healthy volunteers, doses of 40 mg, 80 mg, and 120 mg FXI-LICA are expected to cover the target range of clinical interest for steady-state FXI activity in the phase IIb study of FXI-LICA in patients with ESRD undergoing hemodialysis.

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

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.

Clinical investigation of the biopharmaceutical characteristics of nifurtimox tablets - Implications for quality control and application

Nifurtimox is approved in Chagas disease and has been used in endemic countries since the 1960s. Nifurtimox, available as a 120 mg tablet, is administered with food typically three times daily, and dose is adjusted for age and bodyweight. Accurately or reproducibly fragmenting the 120 mg tablet for dose adjustment in young children and those with low bodyweight is problematic. Based on the existing tablet formulation, new nifurtimox 30 mg and 120 mg tablets have been developed in a format that can be divided accurately into 15 mg and 60 mg fragments. In adults with chronic Chagas disease, we investigated whether nifurtimox bioavailability is affected by tablet dissolution rate, and whether different diets affect nifurtimox bioavailability. In an open-label, three-period cross-over study (n=36; ClinicalTrials.gov, NCT03350295), patients randomly received three 30 mg tablet formulations (slow, medium, or fast dissolution; a 4 × 30 mg dose of one formulation per period). In an open-label, four-period cross-over study (n=24; ClinicalTrials.gov, NCT03334838) patients randomly fasted or received one of three meal types (high-fat/high-calorie, low-fat, dairy-based) before ingesting nifurtimox (a 4 × 30 mg dose per period). Acceptance criteria for no difference between groups were 90% confidence intervals (CIs) of exposure ratios in the range 0.8-1.25. Nifurtimox bioavailability was unaffected by tablet dissolution kinetics. Ratios of area under the curve at final assessment (AUC_(0-tlast) [90% CI]) were: fast/medium dissolution, 1.061 (0.990-1.137); slow/medium dissolution, 0.964 (0.900-1.033); fast/slow dissolution, 1.100 (1.027-1.179). Compared with a fasting state, nifurtimox bioavailability increased by 73% after a high-fat/high-calorie meal (AUC_(0-tlast) ratio [90% CI], 1.732 [1.581-1.898]); smaller increases were seen with the other meal types (low-fat: 1.602 [1.462-1.755]; dairy-based: 1.340 [1.222-1.468]). Although type of diet can affect bioavailability, taking nifurtimox with food is most important.

Predictive Performance of Physiology-Based Pharmacokinetic Dose Estimates for Pediatric Trials: Evaluation With 10 Bayer Small-Molecule Compounds in Children

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.

Comparing Predictions of a PBPK Model for Cyclosporine With Drug Levels From Therapeutic Drug Monitoring

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.

Associations between model-predicted rivaroxaban exposure and patient characteristics and efficacy and safety outcomes in the prevention of venous thromboembolism

Anticoagulant plasma concentrations and patient characteristics might affect the benefit–risk balance of therapy. The study objective was to assess the impact of model-predicted rivaroxaban exposure and patient characteristics on outcomes in patients receiving rivaroxaban for venous thromboembolism (VTE) prophylaxis (VTE-P) after hip/knee replacement surgery. Post hoc exposure–response analyses were conducted using data from the phase 3 RECORD1–4 studies, in which 12,729 patients were randomized to rivaroxaban 10 mg once daily or enoxaparin for ≤ 39 days. Multivariate regression approaches were used to correlate model-predicted individual rivaroxaban exposures and patient characteristics with outcomes. In the absence of measured rivaroxaban exposure, exposure estimates were predicted based on individual increases in prothrombin time (PT) and by making use of the known correlation between rivaroxaban plasma concentration and dynamics of PT. No significant associations between rivaroxaban exposure and total VTE or major bleeding were identified. A significant association between exposure and a composite of major or non-major clinically relevant (NMCR) bleeding from day 4 after surgery was observed. The relationship was shallow, with an approximate predicted absolute increase in a composite of major or NMCR bleeding from 1.08 [95% confidence interval (CI) 0.76–1.54] to 2.18% (95% CI 1.51–3.17) at the 5th and 95th percentiles of trough plasma concentration, respectively. In conclusion, based on the underlying data and analysis, no reliable target window for exposure with improved benefit–risk could be identified within the investigated exposure range. Hence, monitoring rivaroxaban levels is unlikely to be beneficial in VTE-P.

Associations between model-predicted rivaroxaban exposure and patient characteristics and efficacy and safety outcomes in patients with non-valvular atrial fibrillation

Rivaroxaban exposure and patient characteristics may affect the rivaroxaban benefit–risk balance. This study aimed to quantify associations between model-predicted rivaroxaban exposure and patient characteristics and efficacy and safety outcomes in patients with non-valvular atrial fibrillation (NVAF), using data from the phase 3 ROCKET AF trial (NCT00403767). In ROCKET AF, 14,264 patients with NVAF were randomized to rivaroxaban (20 mg once daily [OD], or 15 mg OD if creatinine clearance was 30–49 mL/min) or dose-adjusted warfarin (median follow-up: 707 days); rivaroxaban plasma concentration was measured in a subset of 161 patients. In this post hoc exposure–response analysis, a multivariate Cox model was used to correlate individual predicted rivaroxaban exposures and patient characteristics with time-to-event efficacy and safety outcomes in 7061 and 7111 patients, respectively. There was no significant association between model-predicted rivaroxaban trough plasma concentration (C_trough) and efficacy outcomes. Creatinine clearance and history of stroke were significantly associated with efficacy outcomes. C_trough was significantly associated with the composite of major or non-major clinically relevant (NMCR) bleeding (hazard ratio [95th percentile vs. median]: 1.26 [95% confidence interval 1.13–1.40]) but not with major bleeding alone. The exposure–response relationship for major or NMCR bleeding was shallow with no clear threshold for an acceleration in risk. History of gastrointestinal bleeding had a greater influence on safety outcomes than C_trough. These results support fixed rivaroxaban 15 mg and 20 mg OD dosages in NVAF. Therapeutic drug monitoring is unlikely to offer clinical benefits in this indication beyond evaluation of patient characteristics.

Applications of Physiologically Based Pharmacokinetic Modeling of Rivaroxaban-Renal and Hepatic Impairment and Drug-Drug Interaction Potential

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

Predictive Pediatric Modeling and Simulation Using Ontogeny Information

Food and Drug Administration submissions of physiologically based pharmacokinetic (PBPK) modeling and simulation of small-molecule drugs document the relevance of pediatric drug development and, in particular, information on dosing strategies in children. The most relevant prerequisite for reliable PBPK-based translation of adult pharmacokinetics of a small molecule to children is knowledge of the drug-specific absorption, distribution, metabolism, and elimination (ADME) processes in adults together with existing information about ontogeny of ADME processes relevant for the drug. All mechanisms driving a drug's clearance are of specific importance. For other drug modalities, our knowledge of ADME processes and ontogeny is still limited. More research is required, for example, to understand why some therapeutic proteins show complex differences in pharmacokinetics between adults and children, whereas other proteins seem to follow simple allometric scaling rules. Ontogeny information originates from various sources, such as (semi)quantitative mRNA expression, in vitro activity data, and deconvolution of in vivo pharmacokinetic data. The workflow for pediatric predictions is well described in several articles documenting successful translation from adults to children. The technical hurdles for PBPK modeling are low. State-of-the-art PBPK modeling software tools provide integrated pediatric translation workflows. For example, PK-Sim and MoBi are freely available as fully transparent open-source software via Open Systems Pharmacology (OSP). With the latest 2019 software release, version 8.0, OSP even provides a fully integrated technical framework for the qualification (and requalification) of any specific intended PBPK use in line with Food and Drug Administration and European Medicines Agency PBPK guidance. Qualification packages for pediatric translation are available on the OSP platform.

Rivaroxaban for treatment of pediatric venous thromboembolism. An Einstein-Jr phase 3 dose-exposure-response evaluation

BACKGROUND

Recently, the randomized EINSTEIN-Jr study showed similar efficacy and safety for rivaroxaban and standard anticoagulation for treatment of pediatric venous thromboembolism (VTE). The rivaroxaban dosing strategy was established based on phase 1 and 2 data in children and through pharmacokinetic (PK) modeling.

METHODS

Rivaroxaban treatment with tablets or the newly developed granules-for-oral suspension formulation was bodyweight-adjusted and administered once-daily, twice-daily, or thrice-daily for children with bodyweights of ≥30, ≥12 to <30, and <12 kg, respectively. Previously, these regimens were confirmed for children weighing ≥20 kg but only predicted in those <20 kg. Based on sparse blood sampling, the daily area under the plasma concentration-time curve [AUC(0-24)ss ] and trough [Ctrough,ss ] and maximum [Cmax,ss ] steady-state plasma concentrations were derived using population PK modeling. Exposure-response graphs were generated to evaluate the potential relationship of individual PK parameters with recurrent VTE, repeat imaging outcomes, and bleeding or adverse events. A taste-and-texture questionnaire was collected for suspension-recipients.

RESULTS

Of the 335 children (aged 0-17 years) allocated to rivaroxaban, 316 (94.3%) were evaluable for PK analyses. Rivaroxaban exposures were within the adult exposure range. No clustering was observed for any of the PK parameters with efficacy, bleeding, or adverse event outcomes. Results were similar for the tablet and suspension formulation. Acceptability and palatability of the suspension were favorable.

DISCUSSION

Based on this analysis and the recently documented similar efficacy and safety of rivaroxaban compared with standard anticoagulation, we conclude that bodyweight-adjusted pediatric rivaroxaban regimens with either tablets or suspension are validated and provide for appropriate treatment of children with VTE.

A Physiologically-Based Pharmacokinetic Model to Describe Ciprofloxacin Pharmacokinetics Over the Entire Span of Life

Background

Physiologically-based pharmacokinetic (PBPK) modeling has received growing interest as a useful tool for the assessment of drug pharmacokinetics by continuous knowledge integration.

Objective

The objective of this study was to build a ciprofloxacin PBPK model for intravenous and oral dosing based on a comprehensive literature review, and evaluate the predictive performance towards pediatric and geriatric patients.

Methods

The aim of this report was to establish confidence in simulations of the ciprofloxacin PBPK model along the development process to facilitate reliable predictions outside of the tested adult age range towards the extremes of ages. Therefore, mean data of 69 published clinical trials were identified and integrated into the model building, simulation and verification process. The predictive performance on both ends of the age scale was assessed using individual data of 258 subjects observed in own clinical trials.

Results

Ciprofloxacin model verification demonstrated no concentration-related bias and accurate simulations for the adult age range, with only 4.8% of the mean observed data points for intravenous administration and 12.1% for oral administration being outside the simulated twofold range. Predictions towards the extremes of ages for the area under the plasma concentration–time curve (AUC) and the maximum plasma concentration (C_max) over the entire span of life revealed a reliable estimation, with only two pediatric AUC observations outside the 90% prediction interval.

Conclusion

Overall, this ciprofloxacin PBPK modeling approach demonstrated the predictive power of a thoroughly informed middle-out approach towards age groups of interest to potentially support the decision-making process.

Electronic supplementary material

The online version of this article (10.1007/s40262-018-0661-6) contains supplementary material, which is available to authorized users.

Bodyweight-adjusted rivaroxaban for children with venous thromboembolism (EINSTEIN-Jr): results from three multicentre, single-arm, phase 2 studies

BACKGROUND

Rivaroxaban has been shown to be efficacious for treatment of venous thromboembolism in adults, and has a reduced risk of bleeding compared with standard anticoagulants. We aimed to develop paediatric rivaroxaban regimens for the treatment of venous thromboembolism in children and adolescents.

METHODS

In this phase 2 programme, we did three studies to evaluate rivaroxaban treatment in children younger than 6 months, aged 6 months to 5 years, and aged 6-17 years. Our studies used a multicentre, single-arm design at 54 sites in Australia, Europe, Israel, Japan, and north America. We included children with objectively confirmed venous thromboembolism previously treated with low-molecular weight heparin, fondaparinux, or a vitamin K antagonist for at least 2 months or, in children who had catheter-related venous thromboembolism for at least 6 weeks. We administered rivaroxaban orally in a bodyweight-adjusted 20 mg-equivalent dose, based on physiologically-based pharmacokinetic modelling predictions and EINSTEIN-Jr phase 1 data in young adults, in either a once-daily (tablets; for those aged 6-17 years), twice-daily (in suspension; for those aged 6 months to 11 years), or three times-daily (in suspension; for those younger than 6 months) dosing regimen for 30 days (or 7 days for those younger than 6 months). The primary aim was to define rivaroxaban treatment regimens that match the target adult exposure range. The principal safety outcome was major bleeding and clinically relevant non-major bleeding. Analyses were per-protocol. The predefined efficacy outcomes were symptomatic recurrent venous thromboembolism, asymptomatic deterioration on repeat imaging at the end of the study treatment period. These trials are registered at ClinicalTrials.gov, numbers NCT02564718, NCT02309411, and NCT02234843.

FINDINGS

Between Feb 11, 2013, and Dec 20, 2017, we enrolled 93 children (ten children younger than 6 months; 15 children aged 6 months to 1 year; 25 children aged 2-5 years; 32 children aged 6-11 years; and 11 children aged 12-17 years) into our study. 89 (96%) children completed study treatment (30 days of treatment, or 7 days in those younger than 6 months), and 93 (100%) children received at least one dose of study treatment and were evaluable for the primary endpoints. None of the children had a major bleed, and four (4%, 95% CI 1·2-10·6) of these children had a clinically relevant non-major bleed (three children aged 12-17 years with menorrhagia and one child aged 6-11 years with gingival bleeding). We found no symptomatic recurrent venous thromboembolism in any patients (0%, 0·0-3·9). 24 (32%) of 75 patients with repeat imaging had their thrombotic burden resolved, 43 (57%) patients improved, and eight (11%) patients were unchanged. No patient deteriorated. We confirmed therapeutic rivaroxaban exposures with once-daily dosing in children with bodyweights of at least 30 kg and with twice-daily dosing in children with bodyweights of at least 20 kg and less than 30 kg. Children with low bodyweights (<20 kg, particularly <12 kg) showed low exposures so, for future studies, rivaroxaban dosages were revised for these weight categories, to match the target adult exposure range. 61 (66%) of 93 children had adverse events during the study. Pyrexia was the most common adverse event (ten [11%] events), and anaemia and neutropenia or febrile neutropenia were the most frequent grade 3 or worse events (four [4%] events each). No children died or were discontinued from rivaroxaban because of adverse events.

INTERPRETATION

Treatment with bodyweight-adjusted rivaroxaban appears to be safe in children. The treatment regimens that we confirmed in children with bodyweights of at least 20 kg and the revised treatment regimens that we predicted in those with bodyweights less than 20 kg will be evaluated in the EINSTEIN-Jr phase 3 trial in children with acute venous thromboembolism.

FUNDING

Bayer AG, Janssen Research and Development.

Influence of model-predicted rivaroxaban exposure and patient characteristics on efficacy and safety outcomes in patients with acute coronary syndrome

Background:

This analysis aimed to evaluate the impact of rivaroxaban exposure and patient characteristics on efficacy and safety outcomes in patients with acute coronary syndrome (ACS) and to determine whether therapeutic drug monitoring might provide additional information regarding rivaroxaban dose, beyond what patient characteristics provide.

Methods:

A post hoc exposure–response analysis was conducted using data from the phase III ATLAS ACS 2 Thrombolysis in Myocardial Infarction (TIMI) 51 study, in which 15,526 randomized ACS patients received rivaroxaban (2.5 mg or 5 mg twice daily) or placebo for a mean of 13 months (maximum follow up: 31 months). A multivariate Cox model was used to correlate individual predicted rivaroxaban exposures and patient characteristics with time-to-event clinical outcomes.

Results:

For the incidence of myocardial infarction (MI), ischemic stroke, or nonhemorrhagic cardiovascular death, hazard ratios (HRs) for steady-state maximum plasma concentration (C_max) in the 5th and 95th percentiles versus the median were statistically significant but close to 1 for both rivaroxaban doses. For TIMI major bleeding events, a statistically significant association was observed with C_max [HR, 1.08; 95% CI, 1.06–1.11 (95th percentile versus median, 2.5 mg twice daily)], sex [HR, 0.56; 95% CI, 0.38–0.84 (female versus male)], and previous revascularization [HR, 0.62; 95% CI, 0.44–0.87 (no versus yes)].

Conclusions:

The shallow slopes of the exposure–response relationships and the lack of a clear therapeutic window render it unlikely that therapeutic drug monitoring in patients with ACS would provide additional information regarding rivaroxaban dose beyond that provided by patient characteristics.

Application of Physiologically-Based and Population Pharmacokinetic Modeling for Dose Finding and Confirmation During the Pediatric Development of Moxifloxacin

Moxifloxacin is a widely used fluoroquinolone for the treatment of complicated intra‐abdominal infections. We applied physiologically‐based pharmacokinetic (PBPK) and population pharmacokinetic (popPK) modeling to support dose selection in pediatric patients. We scaled an existing adult PBPK model to children based on prior physiological knowledge. The resulting model proposed an age‐dependent dosing regimen that was tested in a phase I study. Refined doses were then tested in a phase III study. A popPK analysis of all clinical pediatric data confirmed the PBPK predictions, including the proposed dosing schedule in children, and supported pharmacokinetics‐related safety/efficacy questions. The pediatric PBPK model adequately predicted the doses necessary to achieve antimicrobial efficacy while maintaining safety in the phase I and III pediatric studies. Altogether, this study retroactively demonstrated the robustness and utility of modeling to support dose finding and confirmation in pediatric drug development for moxifloxacin.

Enhancing the Quality of Rivaroxaban Exposure Estimates Using Prothrombin Time in the Absence of Pharmacokinetic Sampling

Prothrombin time (PT) is a measure of coagulation status and was assessed in the majority of patients in the rivaroxaban phase II and III clinical trials as a pharmacodynamic marker. In the absence of sufficient phase III pharmacokinetic (PK) data to provide individual exposure measures for input into rivaroxaban exposure–response analyses, the aim of the present study was to investigate the use of PT‐adjustment approaches (i.e., the use of observed individual PT measurements) to enhance the prediction of individual rivaroxaban exposure metrics (derived using a previously developed integrated population PK model) based on the observed linear relationship between PT and rivaroxaban plasma concentrations. The PT‐adjustment approaches were established using time‐matched PK and PT measurements, which were available from 1,779 patients across four phase II trials and one phase III trial of rivaroxaban. PT‐adjusted exposure estimates improved the identification of statistically significant effects when compared with covariate‐only exposure estimates.

Exploratory evaluation of pharmacodynamics, pharmacokinetics and safety of rivaroxaban in children and adolescents: an EINSTEIN-Jr phase I study

BACKGROUND

The EINSTEIN-Jr program will evaluate rivaroxaban for the treatment of venous thromboembolism (VTE) in children, targeting exposures similar to the 20 mg once-daily dose for adults.

METHODS

This was a multinational, single-dose, open-label, phase I study to describe the pharmacodynamics (PD), pharmacokinetics (PK) and safety of a single bodyweight-adjusted rivaroxaban dose in children aged 0.5-18 years. Children who had completed treatment for a venous thromboembolic event were enrolled into four age groups (0.5-2 years, 2-6 years, 6-12 years and 12-18 years) receiving rivaroxaban doses equivalent to 10 mg or 20 mg (either as a tablet or oral suspension). Blood samples for PK and PD analyses were collected within specified time windows.

RESULTS

Fifty-nine children were evaluated. In all age groups, PD parameters (prothrombin time, activated partial thromboplastin time and anti-Factor Xa activity) showed a linear relationship versus rivaroxaban plasma concentrations and were in line with previously acquired adult data, as well as in vitro spiking experiments. The rivaroxaban pediatric physiologically based pharmacokinetic model, used to predict the doses for the individual body weight groups, was confirmed. No episodes of bleeding were reported, and treatment-emergent adverse events occurred in four children and all resolved during the study.

CONCLUSIONS

Bodyweight-adjusted, single-dose rivaroxaban had predictable PK/PD profiles in children across all age groups from 0.5 to 18 years. The PD assessments based on prothrombin time and activated partial thromboplastin time demonstrated that the anticoagulant effect of rivaroxaban was not affected by developmental hemostasis in children.

TRIAL REGISTRATION

ClinicalTrials.gov number, NCT01145859.

Pharmacokinetics of rivaroxaban in children using physiologically based and population pharmacokinetic modelling: an EINSTEIN-Jr phase I study

BACKGROUND

The EINSTEIN-Jr program will evaluate rivaroxaban for the treatment of venous thromboembolism (VTE) in children, targeting exposures similar to the 20 mg once-daily dose for adults. A physiologically based pharmacokinetic (PBPK) model for pediatric rivaroxaban dosing has been constructed.

METHODS

We quantitatively assessed the pharmacokinetics (PK) of a single rivaroxaban dose in children using population pharmacokinetic (PopPK) modelling and assessed the applicability of the PBPK model. Plasma concentration-time data from the EINSTEIN-Jr phase I study were analysed by non-compartmental and PopPK analyses and compared with the predictions of the PBPK model. Two rivaroxaban dose levels, equivalent to adult doses of rivaroxaban 10 mg and 20 mg, and two different formulations (tablet and oral suspension) were tested in children aged 0.5-18 years who had completed treatment for VTE.

RESULTS

PK data from 59 children were obtained. The observed plasma concentration-time profiles in all subjects were mostly within the 90% prediction interval, irrespective of dose or formulation. The PopPK estimates and non-compartmental analysis-derived PK parameters (in children aged ≥6 years) were in good agreement with the PBPK model predictions.

CONCLUSIONS

These results confirmed the applicability of the rivaroxaban pediatric PBPK model in the pediatric population aged 0.5-18 years, which in combination with the PopPK model, will be further used to guide dose selection for the treatment of VTE with rivaroxaban in EINSTEIN-Jr phase II and III studies.

TRIAL REGISTRATION

ClinicalTrials.gov number, NCT01145859; registration date: 17 June 2010.

Physiologically Based Pharmacokinetic Modeling of Renally Cleared Drugs in Pregnant Women

BACKGROUND

Since pregnant women are considerably underrepresented in clinical trials, information on optimal dosing in pregnancy is widely lacking. Physiologically based pharmacokinetic (PBPK) modeling may provide a method for predicting pharmacokinetic changes in pregnancy to guide subsequent in vivo pharmacokinetic trials in pregnant women, minimizing associated risks.

OBJECTIVES

The goal of this study was to build and verify a population PBPK model that predicts the maternal pharmacokinetics of three predominantly renally cleared drugs (namely cefazolin, cefuroxime, and cefradine) at different stages of pregnancy. It was further evaluated whether the fraction unbound (f _u) could be estimated in pregnant women using a proposed scaling approach.

METHODS

Based on a recent literature review on anatomical and physiological changes during pregnancy, a pregnancy population PBPK model was built using the software PK-Sim^®/MoBi^®. This model comprised 27 compartments, including nine pregnancy-specific compartments. The PBPK model was verified by comparing the predicted maternal pharmacokinetics of cefazolin, cefuroxime, and cefradine with observed in vivo data taken from the literature. The proposed scaling approach for estimating the f _u in pregnancy was evaluated by comparing the predicted f _u with experimentally observed f _u values of 32 drugs taken from the literature.

RESULTS

The pregnancy population PBPK model successfully predicted the pharmacokinetics of cefazolin, cefuroxime, and cefradine at all tested stages of pregnancy. All predicted plasma concentrations fell within a 2-fold error range and 85% of the predicted concentrations within a 1.25-fold error range. The f _u in pregnancy could be adequately predicted using the proposed scaling approach, although a slight underestimation was evident in case of drugs bound to α₁-acidic glycoprotein.

CONCLUSION

Pregnancy population PBPK models can provide a valuable tool to predict a priori the pharmacokinetics of predominantly renally cleared drugs in pregnant women. These models can ultimately support informed decision making regarding optimal dosing regimens in this vulnerable special population.

Gestation-Specific Changes in the Anatomy and Physiology of Healthy Pregnant Women: An Extended Repository of Model Parameters for Physiologically Based Pharmacokinetic Modeling in Pregnancy

BACKGROUND

In the past years, several repositories for anatomical and physiological parameters required for physiologically based pharmacokinetic modeling in pregnant women have been published. While providing a good basis, some important aspects can be further detailed. For example, they did not account for the variability associated with parameters or were lacking key parameters necessary for developing more detailed mechanistic pregnancy physiologically based pharmacokinetic models, such as the composition of pregnancy-specific tissues.

OBJECTIVES

The aim of this meta-analysis was to provide an updated and extended database of anatomical and physiological parameters in healthy pregnant women that also accounts for changes in the variability of a parameter throughout gestation and for the composition of pregnancy-specific tissues.

METHODS

A systematic literature search was carried out to collect study data on pregnancy-related changes of anatomical and physiological parameters. For each parameter, a set of mathematical functions was fitted to the data and to the standard deviation observed among the data. The best performing functions were selected based on numerical and visual diagnostics as well as based on physiological plausibility.

RESULTS

The literature search yielded 473 studies, 302 of which met the criteria to be further analyzed and compiled in a database. In total, the database encompassed 7729 data. Although the availability of quantitative data for some parameters remained limited, mathematical functions could be generated for many important parameters. Gaps were filled based on qualitative knowledge and based on physiologically plausible assumptions.

CONCLUSION

The presented results facilitate the integration of pregnancy-dependent changes in anatomy and physiology into mechanistic population physiologically based pharmacokinetic models. Such models can ultimately provide a valuable tool to investigate the pharmacokinetics during pregnancy in silico and support informed decision making regarding optimal dosing regimens in this vulnerable special population.

Integrated Population Pharmacokinetic Analysis of Rivaroxaban Across Multiple Patient Populations

The population pharmacokinetics (PK) of rivaroxaban have been evaluated in several population‐specific models. We developed an integrated population PK model using pooled data from 4,918 patients in 7 clinical trials across all approved indications. Effects of gender, age, and weight on apparent clearance (CL/F) and apparent volume of distribution (V/F), renal function, and comedication on CL/F, and relative bioavailability as a function of dose (F) were analyzed. Virtual subpopulations for exposure simulations were defined by age, creatinine clearance (CrCL) and body mass index (BMI). Rivaroxaban PK were adequately described by a one‐compartment disposition model with a first‐order absorption rate constant. Significant effects of CrCL, use of comedications, and study population on CL/F, age, weight, and gender on V/F, and dose on F were identified. CrCL had a modest influence on exposure, whereas age and BMI had a minor influence. The model was suitable to predict rivaroxaban exposure in patient subgroups of special interest.

Development of a Whole-Body Physiologically Based Pharmacokinetic Approach to Assess the Pharmacokinetics of Drugs in Elderly Individuals

Background

Because of the vulnerability and frailty of elderly adults, clinical drug development has traditionally been biased towards young and middle-aged adults. Recent efforts have begun to incorporate data from paediatric investigations. Nevertheless, the elderly often remain underrepresented in clinical trials, even though persons aged 65 years and older receive the majority of drug prescriptions. Consequently, a knowledge gap exists with regard to pharmacokinetic (PK) and pharmacodynamic (PD) responses in elderly subjects, leaving the safety and efficacy of medicines for this population unclear.

Objectives

The goal of this study was to extend a physiologically based pharmacokinetic (PBPK) model for adults to encompass the full course of healthy aging through to the age of 100 years, to support dose selection and improve pharmacotherapy for the elderly age group.

Methods

For parameterization of the PBPK model for healthy aging individuals, the literature was scanned for anthropometric and physiological data, which were consolidated and incorporated into the PBPK software PK-Sim^®. Age-related changes that occur from 65 to 100 years of age were the main focus of this work. For a sound and continuous description of an aging human, data on anatomical and physiological changes ranging from early adulthood to old age were included. The capability of the PBPK approach to predict distribution and elimination of drugs was verified using the test compounds morphine and furosemide, administered intravenously. Both are cleared by a single elimination pathway. PK parameters for the two compounds in younger adults and elderly individuals were obtained from the literature. Matching virtual populations—with regard to age, sex, anthropometric measures and dosage—were generated. Profiles of plasma drug concentrations over time, volume of distribution at steady state (V_ss) values and elimination half-life (t_½) values from the literature were compared with those predicted by PBPK simulations for both younger adults and the elderly.

Results

For most organs, the age-dependent information gathered in the extensive literature analysis was dense. In contrast, with respect to blood flow, the literature study produced only sparse data for several tissues, and in these cases, linear regression was required to capture the entire elderly age range. On the basis of age-informed physiology, the predicted PK profiles described age-associated trends well. The root mean squared prediction error for the prediction of plasma concentrations of furosemide and morphine in the elderly were improved by 32 and 49 %, respectively, by use of age-informed physiology. The majority of the individual V_ss and t_½ values for the two model compounds, furosemide and morphine, were well predicted in the elderly population, except for long furosemide half-lifes.

Conclusion

The results of this study support the feasibility of using a knowledge-driven PBPK aging model that includes the elderly to predict PK alterations throughout the entire course of aging, and thus to optimize drug therapy in elderly individuals. These results indicate that pharmacotherapy and safety-related control of geriatric drug therapy regimens may be greatly facilitated by the information gained from PBPK predictions.

Electronic supplementary material

The online version of this article (doi:10.1007/s40262-016-0422-3) contains supplementary material, which is available to authorized users.

Addressing Adherence Using Genotype-Specific PBPK Modeling-Impact of Drug Holidays on Tamoxifen and Endoxifen Plasma Levels

Introduction: Tamoxifen is one of the most common treatment opportunities for hormonal positive breast cancer. Despite its good tolerability, patients demonstrate decreasing adherence over years impacting on therapeutic success. PBPK modeling was applied to demonstrate the impact of drug holidays on plasma levels of tamoxifen and its active metabolite endoxifen for different CYP2D6 genotypes.

Materials and Methods: A virtual study with 24,000 patients was conducted in order to investigate the development of tamoxifen steady-state kinetics in patient groups of different CYP2D6 genotypes. The impact of drug holidays on steady-state kinetics was investigated assuming changing drug holiday scenarios.

Results: Drug holidays in CYP2D6 extensive and intermediate metabolizers (EMs, IMs) exceeding 1 month lead to a decrease of endoxifen steady-state trough levels below the 5th percentile of the control group. Assuming drug holidays of 1, 2, or 3 months and administering a fixed-dose combination of 20 mg tamoxifen and 3 mg endoxifen EMs demonstrated re-established endoxifen steady-state trough levels after 5, 8, and 9 days. IMs receiving the same fixed-dose combination demonstrated re-established endoxifen steady-state trough levels after 7, 10, and 11 days.

Discussion: The PBPK model impressively demonstrates the impact of drug holidays in different CYP2D6 genotypes on PK. Population simulation results indicate that drug holidays of more than 2 weeks cause a tremendous decrease of plasma levels despite the long half-life of tamoxifen. To improve therapeutic success, PBPK modeling allows identifying genotype-specific differences in PK following drug holidays and adequate treatment with loading doses.

Applied Concepts in PBPK Modeling: How to Build a PBPK/PD Model

The aim of this tutorial is to introduce the fundamental concepts of physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling with a special focus on their practical implementation in a typical PBPK model building workflow. To illustrate basic steps in PBPK model building, a PBPK model for ciprofloxacin will be constructed and coupled to a pharmacodynamic model to simulate the antibacterial activity of ciprofloxacin treatment.

Development of a Physiologically-Based Pharmacokinetic Model for Preterm Neonates: Evaluation with In Vivo Data

Among pediatric patients, preterm neonates and newborns are the most vulnerable subpopulation. Rapid developmental changes of physiological factors affecting the pharmacokinetics of drug substances in newborns require extreme care in dose and dose regimen decisions. These decisions could be supported by in silico methods such as physiologically-based pharmacokinetic (PBPK) modeling. In a comprehensive literature search, the physiological information of preterm neonates that is required to establish a PBPK model has been summarized and implemented into the database of a generic PBPK software. Physiological parameters include the organ weights and blood flow rates, tissue composition, as well as ontogeny information about metabolic and elimination processes in the liver and kidney. The aim of this work is to evaluate the model's accuracy in predicting the pharmacokinetics following intravenous administration of two model drugs with distinct physicochemical properties and elimination pathways based on earlier reported in vivo data. To this end, PBPK models of amikacin and paracetamol have been set up to predict their plasma levels in preterm neonates. Predicted plasma concentration-time profiles were compared to experimentally obtained in vivo data. For both drugs, plasma concentration time profiles following single and multiple dosing were appropriately predicted for a large range gestational and postnatal ages. In summary, PBPK simulations in preterm neonates appear feasible and might become a useful tool in the future to support dosing decisions in this special patient population.

Evaluation of changes in oral drug absorption in preterm and term neonates for Biopharmaceutics Classification System (BCS) class I and II compounds

AIMS

Evidence suggests that the rate of oral drug absorption changes during early childhood. Yet, respective clinical implications are currently unclear, particularly for preterm neonates. The objective of this study was to evaluate changes in oral drug absorption after birth for different Biopharmaceutics Classification System (BCS) class I and II compounds to better understand respective implications for paediatric pharmacotherapy.

METHODS

Two paradigm compounds were selected for BCS class I (paracetamol (acetaminophen) and theophylline) and II (indomethacin and ibuprofen), respectively, based on the availability of clinical literature data following intravenous and oral dosing. A comparative population pharmacokinetic analysis was performed in a step-wise manner in NONMEM® 7.2 to characterize and predict changes in oral drug absorption after birth for paracetamol, theophylline and indomethacin.

RESULTS

A one compartment model with an age-dependent maturation function for oral drug absorption was found appropriate to characterize the pharmacokinetics of paracetamol. Our findings indicate that the rate at which a drug is absorbed from the GI tract reaches adult levels within about 1 week after birth. The maturation function for paracetamol was found applicable to theophylline and indomethacin once solubility limitations were overcome via drug formulation. The influence of excipients on solubility and, hence, oral bioavailability was confirmed for ibuprofen, a second BCS class II compound.

CONCLUSIONS

The findings of our study suggest that the processes underlying changes in oral drug absorption after birth are drug-independent and that the maturation function identified for paracetamol may be generally applicable to other BCS class I and II compounds for characterizing drug absorption in preterm as well as term neonates.

Physiology-based pharmacokinetic modeling: ready to be used

Physiology-based pharmacokinetic (PBPK) modeling is well recognized as a technology for mechanistically simulating and predicting the fate of substances in a mammalian body. Today, the demand for this methodology is higher than ever. The pharma industry and regulatory agencies are looking for new methods, which help to speed up and increase the efficiency of the development process for new drugs. Implementing PBPK modeling in the drug research and development workflow contributes significantly to reach this goal.:

Concomitant use of tamoxifen and endoxifen in postmenopausal early breast cancer: prediction of plasma levels by physiologically-based pharmacokinetic modeling

Purpose

To overcome cytochrome P450 2D6 (CYP2D6) mediated tamoxifen resistance in postmenopausal early breast cancer, CYP2D6 phenotype-adjusted tamoxifen dosing in patients with impaired CYP2D6 metabolism and/or the application of endoxifen, the most potent tamoxifen metabolite, are alternative treatment options. To elucidate both strategies comprehensively we used a physiologically-based pharmacokinetic (PBPK) modeling approach.

Methods

Firstly simulation of increasing tamoxifen dosages was performed by a virtual clinical trial including populations of CYP2D6 poor (PM), intermediate (IM) and extensive metabolizers (EM) (N = 8,000). Secondly we performed PBPK-simulations under consideration of tamoxifen use plus concomitant increasing dosages of endoxifen (N = 7,000).

Results

Our virtual study demonstrates that dose escalation of tamoxifen in IMs resulted in endoxifen steady-state plasma concentrations similar to CYP2D6 EMs whereas PMs did not reach EM endoxifen levels. Steady-state plasma concentrations of tamoxifen, N-desmethyl-tamoxifen, 4-hydroxy-tamoxifen and endoxifen were similar in CYP2D6 IMs and PMs versus EMs using once daily dosing of 20 mg tamoxifen and concomitant CYP2D6 phenotype-adjusted endoxifen dosing in IMs and PMs (1 mg/d and 3 mg/d, respectively).

Conclusion

In conclusion, we suggest that co-administration of endoxifen in tamoxifen treated early breast cancer women with impaired CYP2D6 metabolism is a promising alternative to reach plasma concentrations comparable to CYP2D6 EM patients.

A Generic Integrated Physiologically based Whole-body Model of the Glucose-Insulin-Glucagon Regulatory System

Models of glucose metabolism are a valuable tool for fundamental and applied medical research in diabetes. Use cases range from pharmaceutical target selection to automatic blood glucose control. Standard compartmental models represent little biological detail, which hampers the integration of multiscale data and confines predictive capabilities. We developed a detailed, generic physiologically based whole-body model of the glucose-insulin-glucagon regulatory system, reflecting detailed physiological properties of healthy populations and type 1 diabetes individuals expressed in the respective parameterizations. The model features a detailed representation of absorption models for oral glucose, subcutaneous insulin and glucagon, and an insulin receptor model relating pharmacokinetic properties to pharmacodynamic effects. Model development and validation is based on literature data. The quality of predictions is high and captures relevant observed inter- and intra-individual variability. In the generic form, the model can be applied to the development and validation of novel diabetes treatment strategies.

Utilizing in vitro and PBPK tools to link ADME characteristics to plasma profiles: case example nifedipine immediate release formulation

One of the most prominent food-drug interactions is the inhibition of intestinal cytochrome P450 (CYP) 3A enzymes by grapefruit juice ingredients, and, as many drugs are metabolized via CYP 3A, this interaction can be of clinical importance. Calcium channel-blocking agents of the dihydropyridine type, such as felodipine and nifedipine, are subject to extensive intestinal first pass metabolism via CYP 3A, thus resulting in significantly enhanced in vivo exposure of the drug when administered together with grapefruit juice. Physiologically based pharmacokinetic (PBPK) modeling was used to simulate pharmacokinetics of a nifedipine immediate release formulation following concomitant grapefruit juice ingestion, that is, after inhibition of small intestinal CYP 3A enzymes. For this purpose, detailed data about CYP 3A levels were collected from the literature and implemented into commercial PBPK software. As literature reports show that grapefruit juice (i) leads to a marked delay in gastric emptying, and (ii) rapidly lowers the levels of intestinal CYP 3A enzymes, inhibition of intestinal first pass metabolism following ingestion of grapefruit juice was simulated by altering the intestinal CYP 3A enzyme levels and simultaneously decelerating the gastric emptying rate. To estimate the in vivo dispersion and dissolution behavior of the formulation, dissolution tests in several media simulating both the fasted and fed state stomach and small intestine were conducted, and the results from the in vitro dissolution tests were used as input function to describe the in vivo dissolution of the drug. Plasma concentration-time profiles of the nifedipine immediate release formulation both with and without simultaneous CYP 3A inhibition were simulated, and the results were compared with data gathered from the literature. Using this approach, nifedipine plasma profiles could be simulated well both with and without enzyme inhibition. A reduction in small intestinal CYP 3A levels by 60% was found to yield the best results, with simulated nifedipine concentration-time profiles within 20% of the in vivo observed results. By additionally varying the dissolution input of the PBPK model, a link between the dissolution characteristics of the formulation and its in vivo performance could be established.

Using Bayesian-PBPK modeling for assessment of inter-individual variability and subgroup stratification

Purpose

Inter-individual variability in clinical endpoints and occurrence of potentially severe adverse effects represent an enormous challenge in drug development at all phases of (pre-)clinical research. To ensure patient safety it is important to identify adverse events or critical subgroups within the population as early as possible. Hence, a comprehensive understanding of the processes governing pharmacokinetics and pharmacodynamics is of utmost importance. In this paper we combine Bayesian statistics with detailed mechanistic physiologically-based pharmacokinetic (PBPK) models. On the example of pravastatin we demonstrate that this combination provides a powerful tool to investigate inter-individual variability in groups of patients and to identify clinically relevant homogenous subgroups in an unsupervised approach. Since PBPK models allow the identification of physiological, drug-specific and genotype-specific knowledge separately, our approach supports knowledge-based extrapolation to other drugs or populations.

Methods

PBPK models are based on generic distribution models and extensive collections of physiological parameters and allow a mechanistic investigation of drug distribution and drug action. To systematically account for parameter variability within patient populations, a Bayesian-PBPK approach is developed rigorously quantifying the probability of a parameter given the amount of information contained in the measured data. Since these parameter distributions are high-dimensional, a Markov chain Monte Carlo algorithm is used, where the physiological and drug-specific parameters are considered in separate blocks.

Results

Considering pravastatin pharmacokinetics as an application example, Bayesian-PBPK is used to investigate inter-individual variability in a cohort of 10 patients. Correlation analyses infer structural information about the PBPK model. Moreover, homogeneous subpopulations are identified a posteriori by examining the parameter distributions, which can even be assigned to a polymorphism in the hepatic organ anion transporter OATP1B1.

Conclusions

The presented Bayesian-PBPK approach systematically characterizes inter-individual variability within a population by updating prior knowledge about physiological parameters with new experimental data. Moreover, clinically relevant homogeneous subpopulations can be mechanistically identified. The large scale PBPK model separates physiological and drug-specific knowledge which allows, in combination with Bayesian approaches, the iterative assessment of specific populations by integrating information from several drugs.

Electronic supplementary material

The online version of this article (doi:10.1186/2193-9616-1-6) contains supplementary material, which is available to authorized users.

A detailed physiologically based model to simulate the pharmacokinetics and hormonal pharmacodynamics of enalapril on the circulating endocrine Renin-Angiotensin-aldosterone system

The renin-angiotensin-aldosterone system (RAAS) plays a key role in the pathogenesis of cardiovascular disorders including hypertension and is one of the most important targets for drugs. A whole body physiologically based pharmacokinetic (wb PBPK) model integrating this hormone circulation system and its inhibition can be used to explore the influence of drugs that interfere with this system, and thus to improve the understanding of interactions between drugs and the target system. In this study, we describe the development of a mechanistic RAAS model and exemplify drug action by a simulation of enalapril administration. Enalapril and its metabolite enalaprilat are potent inhibitors of the angiotensin-converting-enzyme (ACE). To this end, a coupled dynamic parent-metabolite PBPK model was developed and linked with the RAAS model that consists of seven coupled PBPK models for aldosterone, ACE, angiotensin 1, angiotensin 2, angiotensin 2 receptor type 1, renin, and prorenin. The results indicate that the model represents the interactions in the RAAS in response to the pharmacokinetics (PK) and pharmacodynamics (PD) of enalapril and enalaprilat in an accurate manner. The full set of RAAS-hormone profiles and interactions are consistently described at pre- and post-administration steady state as well as during their dynamic transition and show a good agreement with literature data. The model allows a simultaneous representation of the parent-metabolite conversion to the active form as well as the effect of the drug on the hormone levels, offering a detailed mechanistic insight into the hormone cascade and its inhibition. This model constitutes a first major step to establish a PBPK-PD-model including the PK and the mode of action (MoA) of a drug acting on a dynamic RAAS that can be further used to link to clinical endpoints such as blood pressure.

Development of a physiologically based computational kidney model to describe the renal excretion of hydrophilic agents in rats

A physiologically based kidney model was developed to analyze the renal excretion and kidney exposure of hydrophilic agents, in particular contrast media, in rats. In order to study the influence of osmolality and viscosity changes, the model mechanistically represents urine concentration by water reabsorption in different segments of kidney tubules and viscosity dependent tubular fluid flow. The model was established using experimental data on the physiological steady state without administration of any contrast media or drugs. These data included the sodium and urea concentration gradient along the cortico-medullary axis, water reabsorption, urine flow, and sodium as well as urea urine concentrations for a normal hydration state. The model was evaluated by predicting the effects of mannitol and contrast media administration and comparing to experimental data on cortico-medullary concentration gradients, urine flow, urine viscosity, hydrostatic tubular pressures and single nephron glomerular filtration rate. Finally the model was used to analyze and compare typical examples of ionic and non-ionic monomeric as well as non-ionic dimeric contrast media with respect to their osmolality and viscosity. With the computational kidney model, urine flow depended mainly on osmolality, while osmolality and viscosity were important determinants for tubular hydrostatic pressure and kidney exposure. The low diuretic effect of dimeric contrast media in combination with their high intrinsic viscosity resulted in a high viscosity within the tubular fluid. In comparison to monomeric contrast media, this led to a higher increase in tubular pressure, to a reduction in glomerular filtration rate and tubular flow and to an increase in kidney exposure. The presented kidney model can be implemented into whole body physiologically based pharmacokinetic models and extended in order to simulate the renal excretion of lipophilic drugs which may also undergo active secretion and reabsorption.

Pharmacogenomics of codeine, morphine, and morphine-6-glucuronide: model-based analysis of the influence of CYP2D6 activity, UGT2B7 activity, renal impairment, and CYP3A4 inhibition

BACKGROUND AND OBJECTIVE

The analgesic effect of codeine depends on the formation of the opioid metabolites morphine and morphine-6-glucuronide. Different factors have been shown or suspected to affect the safety and efficacy of codeine treatment. The objective of the current study is to assess and quantify the impact of important pharmacokinetic factors, using a mechanistic modeling approach.

METHODS

By means of a generic modeling approach integrating prior physiologic knowledge, we systematically investigated the complex dependence of opioid exposure on cytochrome P450 2D6 and 3A4 (CYP2D6 and CYP3A4), and uridine diphosphate glucuronosyltransferase 2B7 (UGT2B7) activity, as well as renal function, by means of a virtual clinical trial.

RESULTS

First, the known dominant role of CYP2D6 activity for morphine exposure was reproduced. Second, the model demonstrated that mild and moderate renal impairment and co-administration of CYP3A4 inhibitors have only minor influences on opioid exposure. Third, the model showed - in contrast to current opinion - that increased UGT2B7 activity is associated with a decrease in active opioid exposure.

CONCLUSION

Overall, the model-based analysis predicts a wide range of morphine levels after codeine administration and supports recent doubts about safe and efficacious use of codeine for analgesia in non-genotyped individuals.

Physiologically Based Pharmacokinetic Modeling of Tamoxifen and its Metabolites in Women of Different CYP2D6 Phenotypes Provides New Insight into the Tamoxifen Mass Balance

Tamoxifen is a first-line endocrine agent in the mechanism-based treatment of estrogen receptor positive (ER⁺) mammary carcinoma and applied to breast cancer patients all over the world. Endoxifen is a secondary and highly active metabolite of tamoxifen that is formed among others by the polymorphic cytochrome P450 2D6 (CYP2D6). It is widely accepted that CYP2D6 poor metabolizers exert a pronounced decrease in endoxifen steady-state plasma concentrations compared to CYP2D6 extensive metabolizers. Nevertheless, an in-depth understanding of the chain of cause and effect between CYP2D6 genotype, endoxifen steady-state plasma concentration, and subsequent tamoxifen treatment benefit still remains to be evolved. In this study, physiologically based pharmacokinetic (PBPK)-modeling was applied to mechanistically investigate the impact of CYP2D6 phenotype on endoxifen formation in female breast cancer patients undergoing tamoxifen therapy. A PBPK-model of tamoxifen and its pharmacologically important metabolites N-desmethyltamoxifen (NDM-TAM), 4-hydroxytamoxifen (4-OH-TAM), and endoxifen was developed and validated. This model is able to simulate the pharmacokinetics (PK) after single and repeated oral tamoxifen doses in female breast cancer patients in dependence of the CYP2D6 phenotype. A detailed model-based analysis of the mass balance offered support for a recent hypothesis stating a more prominent role for endoxifen formation from 4-OH-TAM. In the future this model provides a good basis to further investigate the linkage of PK, mode of action, and treatment outcome in dependence of factors such as phenotype, ethnicity, or co-treatment with CYP2D6 inhibitors.