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.

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.

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

A generic whole body physiologically based pharmacokinetic model for therapeutic proteins in PK-Sim.

Proteins are an increasingly important class of drugs used as therapeutic as well as diagnostic agents. A generic physiologically based pharmacokinetic (PBPK) model was developed in order to represent at whole body level the fundamental mechanisms driving the distribution and clearance of large molecules like therapeutic proteins. The model was built as an extension of the PK-Sim model for small molecules incorporating (i) the two-pore formalism for drug extravasation from blood plasma to interstitial space, (ii) lymph flow, (iii) endosomal clearance and (iv) protection from endosomal clearance by neonatal Fc receptor (FcRn) mediated recycling as especially relevant for antibodies. For model development and evaluation, PK data was used for compounds with a wide range of solute radii. The model supports the integration of knowledge gained during all development phases of therapeutic proteins, enables translation from pre-clinical species to human and allows predictions of tissue concentration profiles which are of relevance for the analysis of on-target pharmacodynamic effects as well as off-target toxicity. The current implementation of the model replaces the generic protein PBPK model available in PK-Sim since version 4.2 and becomes part of the Open Systems Pharmacology Suite.

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

In vivo [6,6-2H2]-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) 2H2-glucose studies and very-long (9-week) 2H2O 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.

Human neutrophil kinetics: modeling of stable isotope labeling data supports short blood neutrophil half-lives.

Human neutrophils have traditionally been thought to have a short half-life in blood; estimates vary from 4 to 18 hours. This dogma was recently challenged by stable isotope labeling studies with heavy water, which yielded estimates in excess of 3 days. To investigate this disparity, we generated new stable isotope labeling data in healthy adult subjects using both heavy water (n = 4) and deuterium-labeled glucose (n = 9), a compound with more rapid labeling kinetics. To interpret results, we developed a novel mechanistic model and applied it to previously published (n = 5) and newly generated data. We initially constrained the ratio of the blood neutrophil pool to the marrow precursor pool (ratio = 0.26; from published values). Analysis of heavy water data sets yielded turnover rates consistent with a short blood half-life, but parameters, particularly marrow transit time, were poorly defined. Analysis of glucose-labeling data yielded more precise estimates of half-life (0.79 ± 0.25 days; 19 hours) and marrow transit time (5.80 ± 0.42 days). Substitution of this marrow transit time in the heavy water analysis gave a better-defined blood half-life of 0.77 ± 0.14 days (18.5 hours), close to glucose-derived values. Allowing the ratio of blood neutrophils to mitotic neutrophil precursors (R) to vary yielded a best-fit value of 0.19. Reanalysis of the previously published model and data also revealed the origin of their long estimates for neutrophil half-life: an implicit assumption that R is very large, which is physiologically untenable. We conclude that stable isotope labeling in healthy humans is consistent with a blood neutrophil half-life of less than 1 day.

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 Km and Vmax 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.