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.

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.

Quantitative Systems Pharmacology Approaches for Immuno-Oncology: Adding Virtual Patients to the Development Paradigm

Drug development in oncology commonly exploits the tools of molecular biology to gain therapeutic benefit through reprograming of cellular responses. In immuno‐oncology (IO) the aim is to direct the patient’s own immune system to fight cancer. After remarkable successes of antibodies targeting PD1/PD‐L1 and CTLA4 receptors in targeted patient populations, the focus of further development has shifted toward combination therapies. However, the current drug‐development approach of exploiting a vast number of possible combination targets and dosing regimens has proven to be challenging and is arguably inefficient. In particular, the unprecedented number of clinical trials testing different combinations may no longer be sustainable by the population of available patients. Further development in IO requires a step change in selection and validation of candidate therapies to decrease development attrition rate and limit the number of clinical trials. Quantitative systems pharmacology (QSP) proposes to tackle this challenge through mechanistic modeling and simulation. Compounds’ pharmacokinetics, target binding, and mechanisms of action as well as existing knowledge on the underlying tumor and immune system biology are described by quantitative, dynamic models aiming to predict clinical results for novel combinations. Here, we review the current QSP approaches, the legacy of mathematical models available to quantitative clinical pharmacologists describing interaction between tumor and immune system, and the recent development of IO QSP platform models. We argue that QSP and virtual patients can be integrated as a new tool in existing IO drug development approaches to increase the efficiency and effectiveness of the search for novel combination therapies.

Physiologically Based Pharmacokinetic Modeling of Monoclonal Antibodies in Pediatric Populations Using PK-Sim

Physiologically based pharmacokinetic (PBPK) models are increasingly used to support pediatric dose selection for small molecule drugs. In contrast, only a few pediatric PBPK models for therapeutic antibodies have been published recently, and the knowledge on the maturation of the processes relevant for antibody pharmacokinetics (PK) is limited compared to small molecules. The aim of this study was, thus, to evaluate predictions from antibody PBPK models for children which were scaled from PBPK models for adults in order to identify respective knowledge gaps. For this, we used the generic PBPK model implemented in PK-Sim without further modifications. Focusing on general clearance and distribution mechanisms, we selected palivizumab and bevacizumab as examples for this evaluation since they show simple, linear PK which is not governed by drug-specific target mediated disposition at usual therapeutic dosages, and their PK has been studied in pediatric populations after intravenous application. The evaluation showed that the PK of palivizumab was overall reasonably well predicted, while the clearance for bevacizumab seems to be underestimated. Without implementing additional ontogeny for antibody PK-specific processes into the PBPK model, bodyweight normalized clearance increases only moderately in young children compared to adults. If growth during aging at the time of the simulation was considered, the apparent clearance is approximately 20% higher compared to simulations for which growth was not considered for newborns due to the long half-life of antibodies. To fully understand the differences and similarities in the PK of antibodies between adults and children, further research is needed. By integrating available information and data, PBPK modeling can contribute to reveal the relevance of involved processes as well as to generate and test hypothesis.

Quantitative systems pharmacology of interferon alpha administration: A multi-scale approach

The therapeutic effect of a drug is governed by its pharmacokinetics which determine the downstream pharmacodynamic response within the cellular network. A complete understanding of the drug-effect relationship therefore requires multi-scale models which integrate the properties of the different physiological scales. Computational modelling of these individual scales has been successfully established in the past. However, coupling of the scales remains challenging, although it will provide a unique possibility of mechanistic and holistic analyses of therapeutic outcomes for varied treatment scenarios. We present a methodology to combine whole-body physiologically-based pharmacokinetic (PBPK) models with mechanistic intracellular models of signal transduction in the liver for therapeutic proteins. To this end, we developed a whole-body distribution model of IFN-α in human and a detailed intracellular model of the JAK/STAT signalling cascade in hepatocytes and coupled them at the liver of the whole-body human model. This integrated model infers the time-resolved concentration of IFN-α arriving at the liver after intravenous injection while simultaneously estimates the effect of this dose on the intracellular signalling behaviour in the liver. In our multi-scale physiologically-based pharmacokinetic/pharmacodynamic (PBPK/PD) model, receptor saturation is seen at low doses, thus giving mechanistic insights into the pharmacodynamic (PD) response. This model suggests a fourfold lower intracellular response after administration of a typical IFN-α dose to an individual as compared to the experimentally observed responses in in vitro setups. In conclusion, this work highlights clear differences between the observed in vitro and in vivo drug effects and provides important suggestions for future model-based study design.

Abstract 2791: Physiologically based pharmacokinetic modeling and simulations to estimate the efficacious dose of the CEACAM6 function-blocking antibody BAY 1834942

BAY 1834942 is an immunostimulatory function-blocking (fb) antibody (Ab) against the target carcinoembryonic antigen related cell adhesion molecule 6 (CEACAM6) expressed on tumor cells in multiple cancer indications. The suggested mode of action of BAY 1834942 is the blockade of the immunosuppressive effect of CEACAM6 on activated T cells which restores the immune response against cancer cells.

Available preclinical pharmacokinetic (PK) and in vitro pharmacodynamic (PD) data, target receptor density information and tumor (patho-)physiology were used to create a model framework taking into account the three most essential ‘pillars' (target exposure, target binding and drug activity) to estimate the human efficacious dosing of BAY 1834942. For this purpose, a physiologically-based pharmacokinetic (PBPK) model considering target binding of BAY 1834942 in tumor tissue and on blood granulocytes has been developed. The aim of the PBPK model was to estimate the dose range and regimen for humans leading to exposure level at the tumor site that allows sufficient target binding. The PBPK simulations were based on 1) an analysis of PD in vitro data in order to estimate the degree of saturation needed for maximum drug activity, 2) the assessment of CEACAM6 receptor numbers on tumor cells and blood granulocytes and 3) in vivo PK data in order to develop and evaluate the PBPK model. For the latter, plasma PK data of BAY 1834942 in monkeys were used as well as known tumor concentration-time profiles of the antibodies MOPC21 (non-targeting Ab) and ZCE025 (anti-CEA Ab) in mice and humans. Uncertainty of parameters which are relevant for CEACAM6 target saturation was considered by stochastic in silico simulations to estimate the CEACAM6 saturation at the tumor vs. dosing. This analysis revealed that the predicted human efficacious dose strongly depends on CEACAM6 density. Thus, a low CEACAM6 density scenario (25,000 CEACAM6/tumor cell) and a high CEACAM6 density scenario (250,000 CEACAM6/tumor cell) were simulated and used to support dose selection for the first-in-man (FIM) study of BAY 1834942. The FIM study is currently under preparation.

Citation Format: Sabine Wittemer-Rump, Christoph Niederalt, Joerg Willuda, Mark Trautwein, Merlin Luetke-Eversloh, Wolf-Dietrich Doecke, Clemens Guenther, Christian Scheerans. Physiologically based pharmacokinetic modeling and simulations to estimate the efficacious dose of the CEACAM6 function-blocking antibody BAY 1834942 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2791.

Physiologically Based Simulations of Deuterated Glucose for Quantifying Cell Turnover in Humans

In vivo [6,6-²H₂]-glucose labeling is a state-of-the-art technique for quantifying cell proliferation and cell disappearance in humans. However, there are discrepancies between estimates of T cell proliferation reported in short (1-day) versus long (7-day) ²H₂-glucose studies and very-long (9-week) ²H₂O studies. It has been suggested that these discrepancies arise from underestimation of true glucose exposure from intermittent blood sampling in the 1-day study. Label availability in glucose studies is normally approximated by a “square pulse” (Sq pulse). Since the body glucose pool is small and turns over rapidly, the availability of labeled glucose can be subject to large fluctuations and the Sq pulse approximation may be very inaccurate. Here, we model the pharmacokinetics of exogenous labeled glucose using a physiologically based pharmacokinetic (PBPK) model to assess the impact of a more complete description of label availability as a function of time on estimates of CD4+ and CD8+ T cell proliferation and disappearance. The model enabled us to predict the exposure to labeled glucose during the fasting and de-labeling phases, to capture the fluctuations of labeled glucose availability caused by the intake of food or high-glucose beverages, and to recalculate the proliferation and death rates of immune cells. The PBPK model was used to reanalyze experimental data from three previously published studies using different labeling protocols. Although using the PBPK enrichment profile decreased the 1-day proliferation estimates by about 4 and 7% for CD4 and CD8+ T cells, respectively, differences with the 7-day and 9-week studies remained significant. We conclude that the approximations underlying the “square pulse” approach—recently suggested as the most plausible hypothesis—only explain a component of the discrepancy in published T cell proliferation rate estimates.

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.

Computational investigation of potential dosing schedules for a switch of medication from warfarin to rivaroxaban-an oral, direct Factor Xa inhibitor

The long-lasting anticoagulant effect of vitamin K antagonists can be problematic in cases of adverse drug reactions or when patients are switched to another anticoagulant therapy. The objective of this study was to examine in silico the anticoagulant effect of rivaroxaban, an oral, direct Factor Xa inhibitor, combined with the residual effect of discontinued warfarin. Our simulations were based on the recommended anticoagulant dosing regimen for stroke prevention in patients with atrial fibrillation. The effects of the combination of discontinued warfarin plus rivaroxaban were simulated using an extended version of a previously validated blood coagulation computer model. A strong synergistic effect of the two distinct mechanisms of action was observed in the first 2–3 days after warfarin discontinuation; thereafter, the effect was close to additive. Nomograms for the introduction of rivaroxaban therapy after warfarin discontinuation were derived for Caucasian and Japanese patients using safety and efficacy criteria described previously, together with the coagulation model. The findings of our study provide a mechanistic pharmacologic rationale for dosing schedules during the therapy switch from warfarin to rivaroxaban and support the switching strategies as outlined in the Summary of Product Characteristics and Prescribing Information for rivaroxaban.

Spatio-temporal simulation of first pass drug perfusion in the liver

The liver is the central organ for detoxification of xenobiotics in the body. In pharmacokinetic modeling, hepatic metabolization capacity is typically quantified as hepatic clearance computed as degradation in well-stirred compartments. This is an accurate mechanistic description once a quasi-equilibrium between blood and surrounding tissue is established. However, this model structure cannot be used to simulate spatio-temporal distribution during the first instants after drug injection. In this paper, we introduce a new spatially resolved model to simulate first pass perfusion of compounds within the naive liver. The model is based on vascular structures obtained from computed tomography as well as physiologically based mass transfer descriptions obtained from pharmacokinetic modeling. The physiological architecture of hepatic tissue in our model is governed by both vascular geometry and the composition of the connecting hepatic tissue. In particular, we here consider locally distributed mass flow in liver tissue instead of considering well-stirred compartments. Experimentally, the model structure corresponds to an isolated perfused liver and provides an ideal platform to address first pass effects and questions of hepatic heterogeneity. The model was evaluated for three exemplary compounds covering key aspects of perfusion, distribution and metabolization within the liver. As pathophysiological states we considered the influence of steatosis and carbon tetrachloride-induced liver necrosis on total hepatic distribution and metabolic capacity. Notably, we found that our computational predictions are in qualitative agreement with previously published experimental data. The simulation results provide an unprecedented level of detail in compound concentration profiles during first pass perfusion, both spatio-temporally in liver tissue itself and temporally in the outflowing blood. We expect our model to be the foundation of further spatially resolved models of the liver in the future.

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.

A computational systems biology software platform for multiscale modeling and simulation: integrating whole-body physiology, disease biology, and molecular reaction networks

Today, in silico studies and trial simulations already complement experimental approaches in pharmaceutical R&D and have become indispensable tools for decision making and communication with regulatory agencies. While biology is multiscale by nature, project work, and software tools usually focus on isolated aspects of drug action, such as pharmacokinetics at the organism scale or pharmacodynamic interaction on the molecular level. We present a modeling and simulation software platform consisting of PK-Sim^® and MoBi^® capable of building and simulating models that integrate across biological scales. A prototypical multiscale model for the progression of a pancreatic tumor and its response to pharmacotherapy is constructed and virtual patients are treated with a prodrug activated by hepatic metabolization. Tumor growth is driven by signal transduction leading to cell cycle transition and proliferation. Free tumor concentrations of the active metabolite inhibit Raf kinase in the signaling cascade and thereby cell cycle progression. In a virtual clinical study, the individual therapeutic outcome of the chemotherapeutic intervention is simulated for a large population with heterogeneous genomic background. Thereby, the platform allows efficient model building and integration of biological knowledge and prior data from all biological scales. Experimental in vitro model systems can be linked with observations in animal experiments and clinical trials. The interplay between patients, diseases, and drugs and topics with high clinical relevance such as the role of pharmacogenomics, drug–drug, or drug–metabolite interactions can be addressed using this mechanistic, insight driven multiscale modeling approach.

Dynamically simulating the interaction of midazolam and the CYP3A4 inhibitor itraconazole using individual coupled whole-body physiologically-based pharmacokinetic (WB-PBPK) models

Background

Drug-drug interactions resulting from the inhibition of an enzymatic process can have serious implications for clinical drug therapy. Quantification of the drugs internal exposure increase upon administration with an inhibitor requires understanding to avoid the drug reaching toxic thresholds. In this study, we aim to predict the effect of the CYP3A4 inhibitors, itraconazole (ITZ) and its primary metabolite, hydroxyitraconazole (OH-ITZ) on the pharmacokinetics of the anesthetic, midazolam (MDZ) and its metabolites, 1' hydroxymidazolam (1OH-MDZ) and 1' hydroxymidazolam glucuronide (1OH-MDZ-Glu) using mechanistic whole body physiologically-based pharmacokinetic simulation models. The model is build on MDZ, 1OH-MDZ and 1OH-MDZ-Glu plasma concentration time data experimentally determined in 19 CYP3A5 genotyped adult male individuals, who received MDZ intravenously in a basal state. The model is then used to predict MDZ, 1OH-MDZ and 1OH-MDZ-Glu concentrations in an CYP3A-inhibited state following ITZ administration.

Results

For the basal state model, three linked WB-PBPK models (MDZ, 1OH-MDZ, 1OH-MDZ-Glu) for each individual were elimination optimized that resulted in MDZ and metabolite plasma concentration time curves that matched individual observed clinical data. In vivo K_m and V_max optimized values for MDZ hydroxylation were similar to literature based in vitro measures. With the addition of the ITZ/OH-ITZ model to each individual coupled MDZ + metabolite model, the plasma concentration time curves were predicted to greatly increase the exposure of MDZ as well as to both increase exposure and significantly alter the plasma concentration time curves of the MDZ metabolites in comparison to the basal state curves. As compared to the observed clinical data, the inhibited state curves were generally well described although the simulated concentrations tended to exceed the experimental data between approximately 6 to 12 hours following MDZ administration. This deviations appeared to be greater in the CYP3A5 *1/*1 and CYP3A5 *1/*3 group than in the CYP3A5 *3/*3 group and was potentially the result of assuming that ITZ/OH-ITZ inhibits both CYP3A4 and CYP3A5, whereas in vitro inhibition is due to CYP3A4.

Conclusion

This study represents the first attempt to dynamically simulate metabolic enzymatic drug-drug interactions via coupled WB-PBPK models. The workflow described herein, basal state optimization followed by inhibition prediction, is novel and will provide a basis for the development of other inhibitor models that can be used to guide, interpret, and potentially replace clinical drug-drug interaction trials.