From physicochemistry to absorption and distribution: predictive mechanistic modelling and computational tools

A physiologically-based quantitative systems pharmacology model for mechanistic understanding of the response to alogliptin and its application in patients with renal impairment

Alogliptin is a highly selective inhibitor of dipeptidyl peptidase-4 and primarily excreted as unchanged drug in the urine, and differences in clinical outcomes in renal impairment patients increase the risk of serious adverse reactions. In this study, we developed a comprehensive physiologically-based quantitative systematic pharmacology model of the alogliptin-glucose control system to predict plasma exposure and use glucose as a clinical endpoint to prospectively understand its therapeutic outcomes with varying renal function. Our model incorporates a PBPK model for alogliptin, DPP-4 activity described by receptor occupancy theory, and the crosstalk and feedback loops for GLP-1-GIP-glucagon, insulin, and glucose. Based on the optimization of renal function-dependent parameters, the model was extrapolated to different stages renal impairment patients. Ultimately our model adequately describes the pharmacokinetics of alogliptin, the progression of DPP-4 inhibition over time and the dynamics of the glucose control system components. The extrapolation results endorse the dose adjustment regimen of 12.5 mg once daily for moderate patients and 6.25 mg once daily for severe and ESRD patients, while providing additional reflections and insights. In clinical practice, our model could provide additional information on the in vivo fate of DPP4 inhibitors and key regulators of the glucose control system.

Polymer-interface-tissue model to estimate leachable release from medical devices

The ability to predict clinically relevant exposure to potentially hazardous compounds that can leach from polymeric components can help reduce testing needed to evaluate the biocompatibility of medical devices. In this manuscript, we compare two physics-based exposure models: 1) a simple, one-component model that assumes the only barrier to leaching is the migration of the compound through the polymer matrix and 2) a more clinically relevant, two-component model that also considers partitioning across the polymer-tissue interface and migration in the tissue away from the interface. Using data from the literature, the variation of the model parameters with key material properties were established, enabling the models to be applied to a wide range of combinations of leachable compound, polymer matrix and tissue type. Exposure predictions based on the models suggest that the models are indistinguishable over much of the range of clinically relevant scenarios. However, for systems with low partitioning and/or slow tissue diffusion, the two-component model predicted up to three orders of magnitude less mass release over the same time period. Thus, despite the added complexity, in some scenarios it can be beneficial to use the two-component model to provide more clinically relevant estimates of exposure to leachable substances from implanted devices.

Development and validation of PBPK models for genistein and daidzein for use in a next-generation risk assessment

Introduction

All cosmetic ingredients must be evaluated for their safety to consumers. In the absence of in vivo data, systemic concentrations of ingredients can be predicted using Physiologically based Pharmacokinetic (PBPK) models. However, more examples are needed to demonstrate how they can be validated and applied in Next-Generation Risk Assessments (NGRA) of cosmetic ingredients. We used a bottom-up approach to develop human PBPK models for genistein and daidzein for a read-across NGRA, whereby genistein was the source chemical for the target chemical, daidzein.

Methods

An oral rat PBPK model for genistein was built using PK-Sim® and in vitro ADME input data. This formed the basis of the daidzein oral rat PBPK model, for which chemical-specific input parameters were used. Rat PBPK models were then converted to human models using human-specific physiological parameters and human in vitro ADME data. In vitro skin metabolism and penetration data were used to build the dermal module to represent the major route of exposure to cosmetics.

Results

The initial oral rat model for genistein was qualified since it predicted values within 2-fold of measured in vivo PK values. This was used to predict plasma concentrations from the in vivo NOAEL for genistein to set test concentrations in bioassays. Intrinsic hepatic clearance and unbound fractions in plasma were identified as sensitive parameters impacting the predicted Cmax values. Sensitivity and uncertainty analyses indicated the developed PBPK models had a moderate level of confidence. An important aspect of the development of the dermal module was the implementation of first-pass metabolism, which was extensive for both chemicals. The final human PBPK model for daidzein was used to convert the in vitro PoD of 33 nM (from an estrogen receptor transactivation assay) to an external dose of 0.2% in a body lotion formulation.

Conclusion

PBPK models for genistein and daidzein were developed as a central component of an NGRA read-across case study. This will help to gain regulatory confidence in the use of PBPK models, especially for cosmetic ingredients.

PBPK-based translation from preclinical species to humans for the full-size IgG therapeutic efalizumab

Introduction

Animal models play a vital role in pharmaceutical research and development by supporting the planning and design of later clinical studies. To improve confidence and reliability of first in human dose estimates it is essential to assess the comparability of animal studies with the human situation. In the context of large molecules, it is particularly important to evaluate the cross-species-translatability of parameters related to neonatal fragment crystallizable receptor (FcRn) binding and target mediated drug disposition (TMDD), as they greatly influence distribution and disposition of proteins in the body of an organism.

Methods

Plasma pharmacokinetic data of the therapeutic protein efalizumab were obtained from literature. Physiologically based pharmacokinetic (PBPK) models were built for three different species (rabbit, non-human primate (NHP), human). Target binding was included in the NHP and human models. The assumption of similar target turnover and target-binding in NHP and human was explored, to gain insights into how these parameters might be translated between species.

Results

Efalizumab PBPK models were successfully developed for three species and concentration-time-profiles could be described appropriately across different intravenously administered doses. The final NHP and human models feature a common set of parameters for target turnover and drug-target-complex internalization, as well as comparable target-binding parameters. Our analyses show that different parameter values for FcRn affinity are crucial to accurately describe the concentration-time profiles.

Discussion

Based on the available data in rabbits, NHP and humans, parameters for FcRn affinity cannot be translated between species, but parameters related to target mediated drug disposition can be translated from NHP to human. The inclusion of additional pharmacokinetic (PK) data including different efalizumab doses would further support and confirm our findings on identifying TMDD and, thus, binding kinetics of efalizumab in NHPs. Furthermore, we suggest that information on target expression and internalization rates could make it possible to develop comprehensive human PBPK models with minimal animal testing. In this project, we compared the pharmacokinetics of a therapeutic protein in rabbit, NHP and human using an open PBPK modeling platform (Open Systems Pharmacology Suite, http://www.open-systems-pharmacology.org). Our findings could support similar translatory studies for first in human dose predictions in the future.

Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans

Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration-time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration-time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.

Application of partition coefficient methods to predict tissue:plasma affinities in common farm animals: Influence of ionisation state

Tissue affinities are conventionally determined from in vivo steady-state tissue and plasma or plasma-water chemical concentration data. In silico approaches were initially developed for preclinical species but standardly applied and tested in human physiologically-based kinetic (PBK) models. Recently, generic PBK models for farm animals have been made available and require partition coefficients as input parameters. In the current investigation, data for species-specific tissue compositions have been collected, and prediction of chemical distribution in various tissues of livestock species for cattle, chicken, sheep and swine have been performed. Overall, tissue composition was very similar across the four farm animal species. However, small differences were observed in moisture, fat and protein content in the various organs within each species. Such differences could be attributed to factors such as variations in age, breed, and weight of the animals and general conditions of the animal itself. With regards to the predictions of tissue:plasma partition coefficients, 80 %, 71 %, 77 % of the model predictions were within a factor 10 using the methods of Berezhkovskiy (2004), Rodgers and Rowland (2006) and Schmitt (2008). The method of Berezhkovskiy (2004) was often providing the most reliable predictions except for swine, where the method of Schmitt (2008) performed best. In addition, investigation of the impact of chemical classes on prediction performance, all methods had very similar reliability. Notwithstanding, no clear pattern regarding specific chemicals or tissues could be detected for the values predicted outside a 10-fold change in certain chemicals or specific tissues. This manuscript concludes with the need for future research, particularly focusing on lipophilicity and species differences in protein binding.

Physiologically based pharmacokinetics modeling and transporter proteomics to predict systemic and local liver and muscle disposition of statins

Statins are used to reduce liver cholesterol levels but also carry a dose-related risk of skeletal muscle toxicity. Concentrations of statins in plasma are often used to assess efficacy and safety, but because statins are substrates of membrane transporters that are present in diverse tissues, local differences in intracellular tissue concentrations cannot be ruled out. Thus, plasma concentration may not be an adequate indicator of efficacy and toxicity. To bridge this gap, we used physiologically based pharmacokinetic (PBPK) modeling to predict intracellular concentrations of statins. Quantitative data on transporter clearance were scaled from in vitro to in vivo conditions by integrating targeted proteomics and transporter kinetics data. The developed PBPK models, informed by proteomics, suggested that organic anion-transporting polypeptide 2B1 (OATP2B1) and multidrug resistance-associated protein 1 (MRP1) play a pivotal role in the distribution of statins in muscle. Using these PBPK models, we were able to predict the impact of alterations in transporter function due to genotype or drug-drug interactions on statin systemic concentrations and exposure in liver and muscle. These results underscore the potential of proteomics-guided PBPK modeling to scale transporter clearance from in vitro data to real-world implications. It is important to evaluate the role of drug transporters when predicting tissue exposure associated with on- and off-target effects.

Machine learning framework to predict pharmacokinetic profile of small molecule drugs based on chemical structure

Accurate prediction of a new compound's pharmacokinetic (PK) profile is pivotal for the success of drug discovery programs. An initial assessment of PK in preclinical species and humans is typically performed through allometric scaling and mathematical modeling. These methods use parameters estimated from in vitro or in vivo experiments, which although helpful for an initial estimation, require extensive animal experiments. Furthermore, mathematical models are limited by the mechanistic underpinning of the drugs' absorption, distribution, metabolism, and elimination (ADME) which are largely unknown in the early stages of drug discovery. In this work, we propose a novel methodology in which concentration versus time profile of small molecules in rats is directly predicted by machine learning (ML) using structure-driven molecular properties as input and thus mitigating the need for animal experimentation. The proposed framework initially predicts ADME properties based on molecular structure and then uses them as input to a ML model to predict the PK profile. For the compounds tested, our results demonstrate that PK profiles can be adequately predicted using the proposed algorithm, especially for compounds with Tanimoto score greater than 0.5, the average mean absolute percentage error between predicted PK profile and observed PK profile data was found to be less than 150%. The suggested framework aims to facilitate PK predictions and thus support molecular screening and design earlier in the drug discovery process.

Physiologically Based Pharmacokinetic Modeling to Unravel the Drug-gene Interactions of Venlafaxine: Based on Activity Score-dependent Metabolism by CYP2D6 and CYP2C19 Polymorphisms

BACKGROUND

Venlafaxine (VEN) is a commonly utilized medication for alleviating depression and anxiety disorders. The presence of genetic polymorphisms gives rise to considerable variations in plasma concentrations across different phenotypes. This divergence in phenotypic responses leads to notable differences in both the efficacy and tolerance of the drug.

PURPOSE

A physiologically based pharmacokinetic (PBPK) model for VEN and its metabolite O-desmethylvenlafaxine (ODV) to predict the impact of CYP2D6 and CYP2C19 gene polymorphisms on VEN pharmacokinetics (PK).

METHODS

The parent-metabolite PBPK models for VEN and ODV were developed using PK-Sim® and MoBi®. Leveraging prior research, derived and implemented CYP2D6 and CYP2C19 activity score (AS)-dependent metabolism to simulate exposure in the drug-gene interactions (DGIs) scenarios. The model's performance was evaluated by comparing predicted and observed values of plasma concentration-time (PCT) curves and PK parameters values.

RESULTS

In the base models, 91.1%, 94.8%, and 94.6% of the predicted plasma concentrations for VEN, ODV, and VEN + ODV, respectively, fell within a twofold error range of the corresponding observed concentrations. For DGI scenarios, these values were 81.4% and 85% for VEN and ODV, respectively. Comparing CYP2D6 AS = 2 (normal metabolizers, NM) populations to AS = 0 (poor metabolizers, PM), 0.25, 0.5, 0.75, 1.0 (intermediate metabolizers, IM), 1.25, 1.5 (NM), and 3.0 (ultrarapid metabolizers, UM) populations in CYP2C19 AS = 2.0 group, the predicted DGI AUC0-96 h ratios for VEN were 3.65, 3.09, 2.60, 2.18, 1.84, 1.56, 1.34, 0.61, and for ODV, they were 0.17, 0.35, 0.51, 0.64, 0.75, 0.83, 0.90, 1.11, and the results were similar in other CYP2C19 groups. It should be noted that PK differences in CYP2C19 phenotypes were not similar across different CYP2D6 groups.

CONCLUSIONS

In clinical practice, the impact of genotyping on the in vivo disposition process of VEN should be considered to ensure the safety and efficacy of treatment.

Physiologically based pharmacokinetic modeling of levetiracetam to predict the exposure in hepatic and renal impairment and elderly populations

Levetiracetam (LEV) is an anti-epileptic drug approved for use in various populations. The pharmacokinetic (PK) behavior of LEV may be altered in the elderly and patients with renal and hepatic impairment. Thus, dosage adjustment is required. This study was conducted to investigate how the physiologically-based PK (PBPK) model describes the PKs of LEV in adult and elderly populations, as well as to predict the PKs of LEV in patients with renal and hepatic impairment in both populations. The whole-body PBPK models were developed using the reported physicochemical properties of LEV and clinical data. The models were validated using data from clinical studies with different dose ranges and different routes and intervals of administration. The fit performance of the models was assessed by comparing predicted and observed blood concentration data and PK parameters. It is recommended that the doses be reduced to ~70%, 60%, and 45% of the adult dose for the mild, moderate, and severe renal impairment populations and ~95%, 80%, and 57% of the adult dose for the Child Pugh-A (CP-A), Child Pugh-B (CP-B), and Child Pugh-C (CP-C) hepatic impairment populations, respectively. No dose adjustment is required for the healthy elderly population, but dose reduction is required for the elderly with organ dysfunction accordingly, on a scale similar to that of adults. A PBPK model of LEV was successfully developed to optimize dosing regimens for special populations.

Model-informed approach to support pediatric dosing for the pan-PI3K inhibitor copanlisib in children and adolescents with relapsed/refractory solid tumors

Copanlisib is an intravenously administered phosphatidylinositol 3-kinase (PI3K) inhibitor which was investigated in pediatric patients with relapsed/refractory solid tumors. A model-informed approach was undertaken to support and confirm an empirically selected starting dose of 28 mg/m2 for pediatric patients ≥1 year old, corresponding to 80% of the adult recommended dose adjusted for body surface area. An adult physiologically based pharmacokinetic (PBPK) model was initially established using copanlisib physicochemical and disposition properties and clinical pharmacokinetics (PK) data and was shown to adequately capture clinical PK across a range of copanlisib doses in adult cancer patients. The adult PBPK model was then extended to the pediatric population through incorporation of age-dependent anatomical and physiological changes and used to simulate copanlisib exposures in pediatric cancer patient age groups. The pediatric PBPK model predicted that the copanlisib 28 mg/m2 dose would achieve similar copanlisib exposures across pediatric ages when compared with historical adult exposures following the approved copanlisib 60 mg dose administered on Days 1, 8, and 15 of a 28-day cycle. Clinical PK were collected from a phase I study in pediatric patients with relapsed/refractory solid tumors (aged ≥4 years). An established adult population PK model was extended to incorporate an allometrically-scaled effect of body surface area and confirmed that the copanlisib maximum tolerated dose of 28 mg/m2 was appropriate to achieve uniform copanlisib exposures across the investigated pediatric age range and consistent exposures to historical data in adult cancer patients. The model-informed approach successfully supported and confirmed the copanlisib pediatric dose recommendation.

In-depth analysis of patterns in selection of different physiologically based pharmacokinetic modeling tools: PartI - Applications and rationale behind the use of open source-code software

PBPK applications published in the literature support a greater adoption of non-open source-code (NOSC) software as opposed to open source-code (OSC) alternatives. However, a significant number of PBPK modelers are still using OSC software, understanding the rationale for the use of this modality is important and may help those embarking on PBPK modeling. No previous analysis of PBPK modeling trends has included the rationale of the modeler. An in-depth analysis of PBPK applications of OSC software is warranted to determine the true impact of OSC software on the rise of PBPK. Publications focussing on PBPK modeling applications, which used OSC software, were identified by systematically searching the scientific literature for original articles. A total of 171 articles were extracted from the narrowed subset. The rise in the use of OSC software for PBPK applications was greater than the general discipline of pharmacokinetics (9 vs. 4), but less than the overall growth of the PBPK area (9 vs. 43). Our report demonstrates conclusively that the surge in PBPK usage is primarily attributable to the availability and implementations of NOSC software. Modelers preferred not to share the reasons for their selection of certain modeling software and no 'explicit' rationale was given to support the use of OSC analysed by this study. As the preference for NOSC versus OSC software tools in the PBPK area continues to be divided, initiatives to add the rationale in using one form over another to every future PBPK modeling report will be a welcomed and informative addition.

Develop adult extrapolation to pediatrics and pediatric dose optimization based on the physiological pharmacokinetic model of azithromycin

Physiologically-based pharmacokinetic (PBPK) models are more frequently used for supporting pediatric dose selection in small-molecule drugs. Through literature research, drug parameters of azithromycin and clinical data from different studies were obtained. Through parameter optimization of the absorption and dissolution process, the adult intravenous model was extended to the adult oral model. The adult intravenous and oral PBPK models are precise to meet the AAFE<2 standard, and the pharmacokinetic parameters of the predicted values of the model are all within the mean standard deviation of the clinical observations. The values of plasma protein unbound fraction, renal clearance, and gastric juice pH between adults and pediatrics were changed by using the age-dependent pediatric organ maturity formula, and the adult model was extrapolated to the pediatric model. The final developed pediatric PBPK model was used to evaluate optimal dosing for children of different developmental ages. The relationship between the frist dose and age was as follows: 8.8 mg/kg/day from 0.5 to 2 years old, 9.2 mg/kg/day from 3 to 6 years old, 9.4 mg/kg/day from 7 to 12 years old, and 8.2 mg/kg/day from 13 to 18 years old, taken in half for 2-5 days. Simultaneously, the simulated exposures achieved with the dosing regimen proposed were comparable to adult plasma exposures for treatment of community-acquired pneumonia. A reasonable azithromycin pharmacokinetic-pharmacodynamic model for adults and pediatrics has been established, which can be demonstrated by the use of literature pediatric data to develop pediatric PBPK models, expanding the scope of this powerful modeling tool.

The investigation of the complex population-drug-drug interaction between ritonavir-boosted lopinavir and chloroquine or ivermectin using physiologically-based pharmacokinetic modeling

OBJECTIVES

Therapy failure caused by complex population-drug-drug (PDDI) interactions including CYP3A4 can be predicted using mechanistic physiologically-based pharmacokinetic (PBPK) modeling. A synergy between ritonavir-boosted lopinavir (LPVr), ivermectin, and chloroquine was suggested to improve COVID-19 treatment. This work aimed to study the PDDI of the two CYP3A4 substrates (ivermectin and chloroquine) with LPVr in mild-to-moderate COVID-19 adults, geriatrics, and pregnancy populations.

METHODS

The PDDI of LPVr with ivermectin or chloroquine was investigated. Pearson's correlations between plasma, saliva, and lung interstitial fluid (ISF) levels were evaluated. Target site (lung epithelial lining fluid [ELF]) levels of ivermectin and chloroquine were estimated.

RESULTS

Upon LPVr coadministration, while the chloroquine plasma levels were reduced by 30, 40, and 20%, the ivermectin plasma levels were increased by a minimum of 425, 234, and 453% in adults, geriatrics, and pregnancy populations, respectively. The established correlation equations can be useful in therapeutic drug monitoring (TDM) and dosing regimen optimization.

CONCLUSIONS

Neither chloroquine nor ivermectin reached therapeutic ELF levels in the presence of LPVr despite reaching toxic ivermectin plasma levels. PBPK modeling, guided with TDM in saliva, can be advantageous to evaluate the probability of reaching therapeutic ELF levels in the presence of PDDI, especially in home-treated patients.

Predictive performance and verification of physiologically based pharmacokinetic model of propylthiouracil

Background: Propylthiouracil (PTU) treats hyperthyroidism and thyroid crisis in all age groups. A variety of serious adverse effects can occur during clinical use and require attention to its pharmacokinetic and pharmacodynamic characteristics in various populations.

Objective: To provide information for individualized dosing and clinical evaluation of PTU in the clinical setting by developing a physiologically based pharmacokinetic (PBPK) model, predicting ADME characteristics, and extrapolating to elderly and pediatric populations.

Methods: Relevant databases and literature were retrieved to collect PTU’s pharmacochemical properties and ADME parameters, etc. A PBPK model for adults was developed using PK-Sim® software to predict tissue distribution and extrapolated to elderly and pediatric populations. The mean fold error (MFE) method was used to compare the differences between predicted and observed values to assess the accuracy of the PBPK model. The model was validated using PTU pharmacokinetic data in healthy adult populations.

Result: The MFE ratios of predicted to observed values of AUC0-t, Cmax, and Tmax were mainly within 0.5 and 2. PTU concentrations in various tissues are lower than venous plasma concentrations. Compared to healthy adults, the pediatric population requires quantitative adjustment to the appropriate dose to achieve the same plasma exposure levels, while the elderly do not require dose adjustments.

Conclusion: The PBPK model of PTU was successfully developed, externally validated, and applied to tissue distribution prediction and special population extrapolation, which provides a reference for clinical individualized drug administration and evaluation.

Does the choice of applied physiologically‐based pharmacokinetics platform matter? A case study on simvastatin disposition and drug–drug interaction

Physiologically‐based pharmacokinetic (PBPK) models have an important role in drug discovery/development and decision making in regulatory submissions. This is facilitated by predefined PBPK platforms with user‐friendly graphical interface, such as Simcyp and PK‐Sim. However, evaluations of platform differences and the potential implications for disposition‐related applications are still lacking. The aim of this study was to assess how PBPK model development, input parameters, and model output are affected by the selection of PBPK platform. This is exemplified via the establishment of simvastatin PBPK models (workflow, final models, and output) in PK‐Sim and Simcyp as representatives of established whole‐body PBPK platforms. The major finding was that the choice of PBPK platform influenced the model development strategy and the final model input parameters, however, the predictive performance of the simvastatin models was still comparable between the platforms. The main differences between the structure and implementation of Simcyp and PK‐Sim were found in the absorption and distribution models. Both platforms predicted equally well the observed simvastatin (lactone and acid) pharmacokinetics (20–80 mg), BCRP and OATP1B1 drug–gene interactions (DGIs), and drug–drug interactions (DDIs) when co‐administered with CYP3A4 and OATP1B1 inhibitors/inducers. This study illustrates that in‐depth knowledge of established PBPK platforms is needed to enable an assessment of the consequences of PBPK platform selection. Specifically, this work provides insights on software differences and potential implications when bridging PBPK knowledge between Simcyp and PK‐Sim users. Finally, it provides a simvastatin model implemented in both platforms for risk assessment of metabolism‐ and transporter‐mediated DGIs and DDIs.

Quality Assurance of PBPK Modeling Platforms and Guidance on Building, Evaluating, Verifying and Applying PBPK Models Prudently under the Umbrella of Qualification: Why, When, What, How and By Whom?

Modeling and simulation emerges as a fundamental asset of drug development. Mechanistic modeling builds upon its strength to integrate various data to represent a detailed structural knowledge of a physiological and biological system and is capable of informing numerous drug development and regulatory decisions via extrapolations outside clinically studied scenarios. Herein, physiologically based pharmacokinetic (PBPK) modeling is the fastest growing branch, and its use for particular applications is already expected or explicitly recommended by regulatory agencies. Therefore, appropriate applications of PBPK necessitates trust in the predictive capability of the tool, the underlying software platform, and related models. That has triggered a discussion on concepts of ensuring credibility of model-based derived conclusions. Questions like 'why', 'when', 'what', 'how' and 'by whom' remain open. We seek for harmonization of recent ideas, perceptions, and related terminology. First, we provide an overview on quality assurance of PBPK platforms with the two following concepts. Platform validation: ensuring software integrity, security, traceability, correctness of mathematical models and accuracy of algorithms. Platform qualification: demonstrating the predictive capability of a PBPK platform within a particular context of use. Second, we provide guidance on executing dedicated PBPK studies. A step-by-step framework focuses on the definition of the question of interest, the context of use, the assessment of impact and risk, the definition of the modeling strategy, the evaluation of the platform, performing model development including model building, evaluation and verification, the evaluation of applicability to address the question, and the model application under the umbrella of a qualified platform.

Use of Real-World Data and Physiologically-Based Pharmacokinetic Modeling to Characterize Enoxaparin Disposition in Children With Obesity

Dosing guidance for children with obesity is often unknown despite the fact that nearly 20% of US children are classified as obese. Enoxaparin, a commonly prescribed low-molecular-weight heparin, is dosed based on body weight irrespective of obesity status to achieve maximum concentration within a narrow therapeutic or prophylactic target range. However, whether children with and without obesity experience equivalent enoxaparin exposure remains unclear. To address this clinical question, 2,825 anti-activated factor X (anti-Xa) surrogate concentrations were collected from the electronic health records of 596 children, including those with obesity. Using linear mixed-effects regression models, we observed that 4-hour anti-Xa concentrations were statistically significantly different in children with and without obesity, even for children with the same absolute dose (P = 0.004). To further mechanistically explore obesity-associated differences in anti-Xa concentration, a pediatric physiologically-based pharmacokinetic (PBPK) model was developed in adults, and then scaled to children with and without obesity. This PBPK model incorporated binding of enoxaparin to antithrombin to form anti-Xa and elimination via heparinase-mediated metabolism and glomerular filtration. Following scaling, the PBPK model predicted real-world pediatric concentrations well, with an average fold error (standard deviation of the fold error) of 0.82 (0.23) and 0.87 (0.26) in children with and without obesity, respectively. PBPK model simulations revealed that children with obesity have at most 20% higher 4-hour anti-Xa concentrations under recommended, total body weight-based dosing compared to children without obesity owing to reduced weight-normalized clearance. Enoxaparin exposure was better matched across age groups and obesity status using fat-free mass weight-based dosing.

Physiologically based pharmacokinetic modeling of brivaracetam and its interactions with rifampin based on CYP2C19 phenotypes

Brivaracetam (BRV), a third-generation antiepileptic drug (AED), is primarily metabolized through amidase hydrolysis and CYP2C19-mediated hydroxylation in vivo. This study utilized physiologically based pharmacokinetic (PBPK) modeling to explore the pharmacokinetics of BRV and drug interactions between BRV and rifampin (RIF), a CYP2C19 inducer, based on CYP2C19 genetic polymorphisms. A PBPK model of BRV was developed in the general population and in individuals with different CYP2C19 phenotypes by adjusting catalytic rate constants (k_cat), and the model was validated with observed clinical data. The model was then extrapolated to predict BRV steady-state plasma concentration in individuals with different CYP2C19 phenotypes, with or without coadministration of RIF. The developed model adequately described BRV exposure in the abovementioned populations. The predicted steady-state area under the curve (AUC_τ-ss) increases by 20% in heterozygous extensive metabolizers (hEMs) and 55% in poor metabolizers (PMs), compared to homozygous extensive metabolizer (EMs). When coadministered with RIF, the model predicted the most significant magnitude of drug-drug interaction (DDI) in EMs, while the exposure change of BRV was minimal in PMs. Referencing the recommended concentration for therapeutic drug monitoring (TDM), we concluded that the current clinical maintenance dose of BRV is acceptable regardless of CYP2C19 polymorphisms and coadministration with RIF.

Determining the Effects of Chronic Kidney Disease on Organic Anion Transporter1/3 Activity Through Physiologically Based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The renal excretion of drugs via organic anion transporters 1 and 3 (OAT1/3) is significantly decreased in patients with renal impairment. This study uses physiologically based pharmacokinetic models to quantify the reduction in OAT1/3-mediated secretion of drugs throughout varying stages of chronic kidney disease.

METHODS

Physiologically based pharmacokinetic models were constructed for four OAT1/3 substrates in healthy individuals: acyclovir, meropenem, furosemide, and ciprofloxacin. Observed data from drug-drug interaction studies with probenecid, a potent OAT1/3 inhibitor, were used to parameterize the contribution of OAT1/3 to the renal elimination of each drug. The models were then translated to patients with chronic kidney disease by accounting for changes in glomerular filtration rate, kidney volume, renal blood flow, plasma protein binding, and hematocrit. Additionally, a relationship was derived between the estimated glomerular filtration rate and the reduction in OAT1/3-mediated secretion of drugs based on the renal extraction ratios of ƿ-aminohippuric acid in patients with varying degrees of renal impairment. The relationship was evaluated in silico by evaluating the predictive performance of each final model in describing the pharmacokinetics (PK) of drugs across stages of chronic kidney disease.

RESULTS

OAT1/3-mediated renal excretion of drugs was found to be decreased by 27-49%, 50-68%, and 70-96% in stage 3, stage 4, and stage 5 of chronic kidney disease, respectively. In support of the parameterization, physiologically based pharmacokinetic models of four OAT1/3 substrates were able to adequately characterize the PK in patients with different degrees of renal impairment. Total exposure after intravenous administration was predicted within a 1.5-fold error and 85% of the observed data points fell within a 1.5-fold prediction error. The models modestly under-predicted plasma concentrations in patients with end-stage renal disease undergoing intermittent hemodialysis. However, results should be interpreted with caution because of the limited number of molecules analyzed and the sparse sampling in observed chronic kidney disease pharmacokinetic studies.

CONCLUSIONS

A quantitative understanding of the reduction in OAT1/3-mediated excretion of drugs in differing stages of renal impairment will contribute to better predictive accuracy for physiologically based pharmacokinetic models in drug development, assisting with clinical trial planning and potentially sparing this population from unnecessary toxic exposures.

The Development of a PBPK Model for Atomoxetine Using Levels in Plasma, Saliva and Brain Extracellular Fluid in Patients with Normal and Deteriorated Kidney Function

BACKGROUND

Atomoxetine is a treatment for attention-deficit hyperactivity disorder. It inhibits Norepinephrine Transporters (NET) in the brain. Renal impairment can reduce hepatic CYP2D6 activity and atomoxetine elimination which may increase its body exposure. Atomoxetine can be secreted in saliva.

OBJECTIVE

The objective of this work was to test the hypothesis that atomoxetine saliva levels (sATX) can be used to predict ATX brain Extracellular Fluid (bECF) levels and their pharmacological effects in healthy subjects and those with End-Stage Renal Disease (ESRD).

METHODS

The pharmacokinetics of atomoxetine after intravenous administration to rats with chemically induced acute and chronic renal impairments were investigated. A physiologically-based pharmacokinetic (PBPK) model was built and verified in rats using previously published measured atomoxetine levels in plasma and brain tissue. The rat PBPK model was then scaled to humans and verified using published measured atomoxetine levels in plasma, saliva, and bECF.

RESULTS

The rat PBPK model predicted the observed reduced atomoxetine clearance due to renal impairment in rats. The PBPK model predicted atomoxetine exposure in human plasma, sATX and bECF. Additionally, it predicted that ATX bECF levels needed to inhibit NET are achieved at 80 mg dose. In ESRD patients, the developed PBPK model predicted that the previously reported 65% increase in plasma exposure in these patients can be associated with a 63% increase in bECF. The PBPK simulations showed that there is a significant correlation between sATX and bECF in human.

CONCLUSION

Saliva levels can be used to predict atomoxetine pharmacological response.

Opportunities and challenges of physiologically based pharmacokinetic modeling in drug delivery

Physiologically based pharmacokinetic (PBPK) modeling is an important in silico tool to bridge drug properties and in vivo PK behaviors during drug development. Over the recent decade, the PBPK method has been largely applied to drug delivery systems (DDS), including oral, inhaled, transdermal, ophthalmic, and complex injectable products. The related therapeutic agents have included small-molecule drugs, therapeutic proteins, nucleic acids, and even cells. Simulation results have provided important insights into PK behaviors of new dosage forms, which strongly support drug regulation. In this review, we comprehensively summarize recent progress in PBPK applications in drug delivery, which shows large opportunities for facilitating drug development. In addition, we discuss the challenges of applying this methodology from a practical viewpoint.

Use of Physiologically-Based Kinetics Modelling to Reliably Predict Internal Concentrations of the UV Filter, Homosalate, After Repeated Oral and Topical Application

Ethical and legal considerations have led to increased use of non-animal methods to evaluate the safety of chemicals for human use. We describe the development and qualification of a physiologically-based kinetics (PBK) model for the cosmetic UV filter ingredient, homosalate, to support its safety without the need of generating further animal data. The intravenous (IV) rat PBK model, using PK-Sim®, was developed and validated using legacy in vivo data generated prior to the 2013 EU animal-testing ban. Input data included literature or predicted physicochemical and pharmacokinetic properties. The refined IV rat PBK model was subject to sensitivity analysis to identify homosalate-specific sensitive parameters impacting the prediction of Cmax (more sensitive than AUC(0-∞)). These were then considered, together with population modeling, to calculate the confidence interval (CI) 95% Cmax and AUC(0-∞). Final model parameters were established by visual inspection of the simulations and biological plausibility. The IV rat model was extrapolated to oral administration, and used to estimate internal exposures to doses tested in an oral repeated dose toxicity study. Next, a human PBK dermal model was developed using measured human in vitro ADME data and a module to represent the dermal route. Model performance was confirmed by comparing predicted and measured values from a US-FDA clinical trial (Identifier: NCT03582215, https://clinicaltrials.gov/). Final exposure estimations were obtained in a virtual population and considering the in vitro and input parameter uncertainty. This model was then used to estimate the Cmax and AUC(0-24 h) of homosalate according to consumer use in a sunscreen. The developed rat and human PBK models had a good biological basis and reproduced in vivo legacy rat and human clinical kinetics data. They also complied with the most recent WHO and OECD recommendations for assessing the confidence level. In conclusion, we have developed a PBK model which predicted reasonably well the internal exposure of homosalate according to different exposure scenarios with a medium to high level of confidence. In the absence of in vivo data, such human PBK models will be the heart of future completely non-animal risk assessments; therefore, valid approaches will be key in gaining their regulatory acceptance. Clinical Trial Registration: https://clinicaltrials.gov/, identifier, NCT03582215.

A Physiological Approach to Pharmacokinetics in Chronic Kidney Disease

The conventional approach to approximating the pharmacokinetics of drugs in patients with chronic kidney disease (CKD) only accounts for changes in the estimated glomerular filtration rate. However, CKD is a systemic and multifaceted disease that alters many body systems. Therefore, the objective of this exercise was to develop and evaluate a whole-body mechanistic approach to predicting pharmacokinetics in patients with CKD. Physiologically based pharmacokinetic models were developed in PK-Sim v8.0 (www.open-systems-pharmacology.org) to mechanistically represent the disposition of 7 compounds in healthy human adults. The 7 compounds selected were eliminated by glomerular filtration and active tubular secretion by the organic cation transport system to varying degrees. After a literature search, the healthy adult models were adapted to patients with CKD by numerically accounting for changes in glomerular filtration rate, kidney volume, renal perfusion, hematocrit, plasma protein concentrations, and gastrointestinal transit. Literature-informed interindividual variability was applied to the physiological parameters to facilitate a population approach. Model performance in CKD was evaluated against pharmacokinetic data from 8 clinical trials in the literature. Overall, integration of the CKD parameterization enabled exposure predictions that were within 1.5-fold error across all compounds and patients with varying stages of renal impairment. Notable improvement was observed over the conventional approach to scaling exposure, which failed in all but 1 scenario in patients with advanced CKD. Further research is required to qualify its use for first-in-CKD dose selection and clinical trial planning for a wider selection of renally eliminated compounds, including those subject to anion transport.

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.

Quantification of the Impact of Partition Coefficient Prediction Methods on Physiologically Based Pharmacokinetic Model Output Using a Standardized Tissue Composition

Tissue:plasma partition coefficients are key parameters in physiologically based pharmacokinetic (PBPK) models, yet the coefficients are challenging to measure in vivo. Several mechanistic-based equations have been developed to predict partition coefficients using tissue composition information and the compound's physicochemical properties, but it is not clear which, if any, of the methods is most appropriate under given circumstances. Complicating the evaluation, each prediction method was developed, and is typically employed, using a different set of tissue composition information, thereby making a controlled comparison impossible. This study proposed a standardized tissue composition for humans that can be used as a common input for each of the five frequently used prediction methods. These methods were implemented in R and were used to predict partition coefficients for 11 drugs, classified as strong bases, weak bases, acids, neutrals, and zwitterions. PBPK models developed in R (mrgsolve) for each drug and each set of partition coefficient predictions were compared with respective observed plasma concentration data. Percent root mean square error and half-life percent error were used to evaluate the accuracy of the PBPK model predictions using each partition coefficient method as summarized by strong bases, weak bases, acids, neutrals, and zwitterions characterization. The analysis indicated that no partition coefficient method consistently yielded the most accurate PBPK model predictions. As such, PBPK model predictions using all partition coefficient methods should be considered during drug development. SIGNIFICANCE STATEMENT: Several mechanistic-based methods exist to predict tissue:plasma partition coefficients critical to PBPK modeling. Controlled comparisons are confounded by the use of different tissue composition values for each method; a standardized tissue composition was proposed. Resulting assessments indicated that no method was consistently superior; therefore, sensitivity of PBPK predictions to each method may be warranted prior to model optimization.

Saliva versus Plasma Therapeutic Drug Monitoring of Gentamicin in Jordanian Preterm Infants. Development of a Physiologically-Based Pharmacokinetic (PBPK) Model and Validation of Class II Drugs of Salivary Excretion Classification System

Gentamicin has proven to be a very successful treatment for bacterial infection, but it also can cause adverse effects, especially ototoxicity, which is irreversible. Therapeutic drug monitoring (TDM) in saliva is a more convenient non-invasive alternative compared to plasma. A physiologically-based pharmacokinetic (PBPK) model of gentamicin was built and validated using previously-published plasma and saliva data. The validated model was then used to predict experimentally-observed plasma and saliva gentamicin TDM data in Jordanian pediatric preterm infant patients measured using sensitive LCMS/MS method. A correlation was established between plasma and saliva exposures. The developed PBPK model predicted previously reported gentamicin levels in plasma, saliva and those observed in the current study. A good correlation was found between plasma and saliva exposures. The PBPK model predicted that gentamicin in saliva is 5-7 times that in plasma, which is in agreement with observed results. Saliva can be used as an alternative for TDM of gentamicin in preterm infant patients. Exposure to gentamicin in plasma and saliva can reliably be predicted using the developed PBPK model in patients.

Physiologically based pharmacokinetic-pharmacodynamic evaluation of meropenem plus fosfomycin in paediatrics

AIMS

The objective of the current study was to evaluate paediatric dosing regimens for meropenem plus fosfomycin that generate sufficient coverage against multidrug-resistant bacteria.

METHODS

The physiologically based pharmacokinetic (PBPK) models of meropenem and fosfomycin were developed from previously published pharmacokinetic studies in five populations: healthy subjects of Japanese origin, and healthy adults, geriatric, paediatric and renally impaired of primarily Caucasian origins. Pharmacodynamic (PD) analyses were carried out by evaluating dosing regimens that achieved a ≥90% joint probability of target attainment (PTA), which was defined as the minimum of the marginal probabilities to achieve the target PD index of each antibiotic. For meropenem, the percentage of time over a 24-hour period wherein the free drug concentration was above the minimum inhibitory concentration (fT > MIC) of at least 40% was its PD target. The fosfomycin PD index was described by fAUC/MIC of at least 40.8.

RESULTS

For coadministration consisting of 20 mg/kg meropenem q8h as a 3-hour infusion and 35 mg/kg fosfomycin q8h also as a 3-hour infusion in a virtual paediatric population between 1 month and 12 years of age with normal renal function and a corresponding body weight between 3 and 50 kg, a joint PTA ≥ 90% is achieved at MICs of 16 and 64 mg/L for meropenem and fosfomycin coadministration, respectively, against Klebsiella pneumoniae and Pseudomonas aeruginosa.

CONCLUSION

The current study identified potentially effective paediatric dosing regimens for meropenem plus fosfomycin coadministration against multidrug-resistant bacteria.

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.

Development of a physiologically based pharmacokinetic model for intravenous lenalidomide in mice

PURPOSE

Lenalidomide is used widely in B-cell malignancies for its immunomodulatory activity. It is primarily eliminated via the kidneys, with a significant proportion of renal elimination attributed to active processes. Lenalidomide is a weak substrate of P-glycoprotein (P-gp), though it is unclear whether P-gp is solely responsible for lenalidomide transport. This study aimed to determine whether the current knowledge of lenalidomide was sufficient to describe the pharmacokinetics of lenalidomide in multiple tissues.

METHODS

A physiologically based pharmacokinetic model was developed using the Open Systems Pharmacology Suite to explore the pharmacokinetics of lenalidomide in a variety of tissues. Data were available for mice dosed intravenously at 0.5, 1.5, 5, and 10 mg/kg, with concentrations measured in plasma, brain, heart, kidney, liver, lung, muscle, and spleen. P-gp expression and activity were sourced from the literature.

RESULTS

The model predictions in plasma, liver, and lung were representative of the observed data (median prediction error 13%, - 10%, and 30%, respectively, with 90% confidence intervals including zero), while other tissue predictions showed sufficient similarity to the observed data. Contrary to the data, model predictions for the brain showed no drug reaching brain tissue when P-gp was expressed at the blood-brain barrier. The data were better described by basolateral transporters at the intracellular wall. Local sensitivity analysis showed that transporter activity was the most sensitive parameter in these models for exposure.

CONCLUSION

As P-gp transport at the blood-brain barrier did not explain the observed brain concentrations alone, there may be other transporters involved in lenalidomide disposition.

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.

In silico approach in reveal traditional medicine plants pharmacological material basis

In recent years, studies of traditional medicinal plants have gradually increased worldwide because the natural sources and variety of such plants allow them to complement modern pharmacological approaches. As computer technology has developed, in silico approaches such as virtual screening and network analysis have been widely utilized in efforts to elucidate the pharmacological basis of the functions of traditional medicinal plants. In the process of new drug discovery, the application of virtual screening and network pharmacology can enrich active compounds among the candidates and adequately indicate the mechanism of action of medicinal plants, reducing the cost and increasing the efficiency of the whole procedure. In this review, we first provide a detailed research routine for examining traditional medicinal plants by in silico techniques and elaborate on their theoretical principles. We also survey common databases, software programs and website tools that can be used for virtual screening and pharmacological network construction. Furthermore, we conclude with a simple example that illustrates the whole methodology, and we present perspectives on the development and application of this in silico methodology to reveal the pharmacological basis of the effects of traditional medicinal plants.

Development and validation of a physiology-based model for the prediction of pharmacokinetics/toxicokinetics in rabbits

The environmental fates of pharmaceuticals and the effects of crop protection products on non-target species are subjects that are undergoing intense review. Since measuring the concentrations and effects of xenobiotics on all affected species under all conceivable scenarios is not feasible, standard laboratory animals such as rabbits are tested, and the observed adverse effects are translated to focal species for environmental risk assessments. In that respect, mathematical modelling is becoming increasingly important for evaluating the consequences of pesticides in untested scenarios. In particular, physiologically based pharmacokinetic/toxicokinetic (PBPK/TK) modelling is a well-established methodology used to predict tissue concentrations based on the absorption, distribution, metabolism and excretion of drugs and toxicants. In the present work, a rabbit PBPK/TK model is developed and evaluated with data available from the literature. The model predictions include scenarios of both intravenous (i.v.) and oral (p.o.) administration of small and large compounds. The presented rabbit PBPK/TK model predicts the pharmacokinetics (Cmax, AUC) of the tested compounds with an average 1.7-fold error. This result indicates a good predictive capacity of the model, which enables its use for risk assessment modelling and simulations.

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.

A physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) model of the histone deacetylase (HDAC) inhibitor vorinostat for pediatric and adult patients and its application for dose specification

PURPOSE

This study aimed at recommending pediatric dosages of the histone deacetylase (HDAC) inhibitor vorinostat and potentially more effective adult dosing regimens than the approved standard dosing regimen of 400 mg/day, using a comprehensive physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling approach.

METHODS

A PBPK/PD model for vorinostat was developed for predictions in adults and children. It includes the maturation of relevant metabolizing enzymes. The PBPK model was expanded by (1) effect compartments to describe vorinostat concentration-time profiles in peripheral blood mononuclear cells (PBMCs), (2) an indirect response model to predict the HDAC inhibition, and (3) a thrombocyte model to predict the dose-limiting thrombocytopenia. Parameterization of drug and system-specific processes was based on published and unpublished in silico, in vivo, and in vitro data. The PBPK modeling software used was PK-Sim and MoBi.

RESULTS

The PBPK/PD model suggests dosages of 80 and 230 mg/m² for children of 0-1 and 1-17 years of age, respectively. In comparison with the approved standard treatment, in silico trials reveal 11 dosing regimens (9 oral, and 2 intravenous infusion rates) increasing the HDAC inhibition by an average of 31%, prolonging the HDAC inhibition by 181%, while only decreasing the circulating thrombocytes to a tolerable 53%. The most promising dosing regimen prolongs the HDAC inhibition by 509%.

CONCLUSIONS

Thoroughly developed PBPK models enable dosage recommendations in pediatric patients and integrated PBPK/PD models, considering PD biomarkers (e.g., HDAC activity and platelet count), are well suited to guide future efficacy trials by identifying dosing regimens potentially superior to standard dosing regimens.

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.

Development of Physiologically Based Organ Models to Evaluate the Pharmacokinetics of Drugs in the Testes and the Thyroid Gland

We extended a generic whole-body physiologically based pharmacokinetic (PBPK) model for rats and humans for organs of the reproductive and endocrine systems (i.e., the testes and the thyroid gland). An extensive literature search was performed, first, to determine the most generic organ model structures for testes and thyroid across species, and, second, to identify the corresponding anatomic and physiological parameters in rats and humans. The testes and thyroid organ models were implemented in the PBPK modeling software PK-Sim and MoBi. The capability of the PBPK approach to simulate the testes and thyroid tissue concentration data was demonstrated using a series of test compounds. The presented organ model structures and parameterization yielded a close agreement between observed and simulated tissue concentrations over time. The organ models are ready to be used to predict the pharmacokinetics of passively entering drugs in the testes and thyroid tissue in a generic PBPK modeling framework.

Translational learning from clinical studies predicts drug pharmacokinetics across patient populations

Early indication of late-stage failure of novel candidate drugs could be facilitated by continuous integration, assessment, and transfer of knowledge acquired along pharmaceutical development programs. We here present a translational systems pharmacology workflow that combines drug cocktail probing in a specifically designed clinical study, physiologically based pharmacokinetic modeling, and Bayesian statistics to identify and transfer (patho-)physiological and drug-specific knowledge across distinct patient populations. Our work builds on two clinical investigations, one with 103 healthy volunteers and one with 79 diseased patients from which we systematically derived physiological information from pharmacokinetic data for a reference probe drug (midazolam) at the single-patient level. Taking into account the acquired knowledge describing (patho-)physiological alterations in the patient cohort allowed the successful prediction of the population pharmacokinetics of a second, candidate probe drug (torsemide) in the patient population. In addition, we identified significant relations of the acquired physiological processes to patient metadata from liver biopsies. The presented prototypical systems pharmacology approach is a proof of concept for model-based translation across different stages of pharmaceutical development programs. Applied consistently, it has the potential to systematically improve predictivity of pharmacokinetic simulations by incorporating the results of clinical trials and translating them to subsequent studies.

Physiologically based modeling together with Bayesian statistics allows the prediction of drug pharmacokinetics in specific patient populations. An interdisciplinary group of clinicians and computational scientists led by Dr. Lars Kuepfer from Bayer developed a generic workflow consisting of several consecutive learning steps where knowledge about both individual physiology as well as drug physicochemistry can be efficiently derived from plasma concentration profiles. The acquired information is then be used for the prediction of the pharmacokinetic behavior of a new drug candidate in a diseased population. This allows to simulate the variability in drug exposure virtually before starting clinical investigation in real patients in order to evaluate drug safety or efficacy through the simulation of virtual populations. Further development of this workflow could improve the safety of clinical development programs to assess the risk-benefit ratio of novel drug candidates in silico.

A Comparative Analysis of Drug-Induced Hepatotoxicity in Clinically Relevant Situations

Drug-induced toxicity is a significant problem in clinical care. A key problem here is a general understanding of the molecular mechanisms accompanying the transition from desired drug effects to adverse events following administration of either therapeutic or toxic doses, in particular within a patient context. Here, a comparative toxicity analysis was performed for fifteen hepatotoxic drugs by evaluating toxic changes reflecting the transition from therapeutic drug responses to toxic reactions at the cellular level. By use of physiologically-based pharmacokinetic modeling, in vitro toxicity data were first contextualized to quantitatively describe time-resolved drug responses within a patient context. Comparatively studying toxic changes across the considered hepatotoxicants allowed the identification of subsets of drugs sharing similar perturbations on key cellular processes, functional classes of genes, and individual genes. The identified subsets of drugs were next analyzed with regard to drug-related characteristics and their physicochemical properties. Toxic changes were finally evaluated to predict both molecular biomarkers and potential drug-drug interactions. The results may facilitate the early diagnosis of adverse drug events in clinical application.

Liver toxicity may occur at drug levels above the therapeutic range and is thus a crucial problem in clinical care. However, the cellular changes induced by drug administration of therapeutic and toxic doses in humans are still not well understood. We here coupled patient-specific drug concentration-time profiles following oral administration of therapeutic and toxic doses with in vitro drug response data to predict toxic changes that quantitatively reflect the transition from desired drug effects to undesired toxic reactions. These toxic changes were comparatively evaluated for fifteen hepatotoxic drugs to identify subsets of drugs, which show similar drug effects on key cellular processes, functional classes of genes, and individual genes, respectively. In addition, analyzing toxic changes for individual genes allowed the prediction of molecular biomarkers and potential drug-drug interactions. Our results may hence support the early diagnosis of liver toxicity in clinical care in the future and may, moreover, help to assess potential risks of drug combination therapies.

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.

Building a 3D Virtual Liver: Methods for Simulating Blood Flow and Hepatic Clearance on 3D Structures

In this paper, we develop a spatio-temporal modeling approach to describe blood and drug flow, as well as drug uptake and elimination, on an approximation of the liver. Extending on previously developed computational approaches, we generate an approximation of a liver, which consists of a portal and hepatic vein vasculature structure, embedded in the surrounding liver tissue. The vasculature is generated via constrained constructive optimization, and then converted to a spatial grid of a selected grid size. Estimates for surrounding upscaled lobule tissue properties are then presented appropriate to the same grid size. Simulation of fluid flow and drug metabolism (hepatic clearance) are completed using discretized forms of the relevant convective-diffusive-reactive partial differential equations for these processes. This results in a single stage, uniformly consistent method to simulate equations for blood and drug flow, as well as drug metabolism, on a 3D structure representative of a liver.

Modeling and Simulation of In Vivo Drug Effects

The concept of a pharmacokinetics-pharmacodynamics (PK/PD) assessment of drug development candidates is well established in pharmaceutical research and development, and PK/PD modeling is common practice in all pharmaceutical companies. A recent analysis (Morgan et al., Drug Discov Today 17(9-10):419-424, 2012) revealed however that insufficient certainty in the integrity of the causal chain of fundamental pharmacological steps from drug dosing through systemic exposure, target tissue exposure, and engagement of molecular target to pharmacological response is still the major driver of failure in phase II of clinical drug development. Despite the rise of molecular biomarkers, ethical, scientific, and practical constraints very often still prevent a direct assessment of each necessary step ultimately leading to an intended drug effect or an unintended adverse reaction. Yet, incomplete investigation of the causality of drug responses is a major risk for translational assessments and the prediction of drug responses in different species or other populations. Mechanism-based modeling and simulation (M&S) offers a means to investigate complex physiological and pharmacological processes and to complement experimental data for non-accessible steps in the pharmacological causal chain. With the help of two examples, it is illustrated, what level of physiological detail, state-of-the-art models can represent, how predictive these models are and how mechanism-based approaches can be combined with empirical correlation-based concepts.

Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations

Interindividual variability in anatomical and physiological properties results in significant differences in drug pharmacokinetics. The consideration of such pharmacokinetic variability supports optimal drug efficacy and safety for each single individual, e.g. by identification of individual-specific dosings. One clear objective in clinical drug development is therefore a thorough characterization of the physiological sources of interindividual variability. In this work, we present a Bayesian population physiologically-based pharmacokinetic (PBPK) approach for the mechanistically and physiologically realistic identification of interindividual variability. The consideration of a generic and highly detailed mechanistic PBPK model structure enables the integration of large amounts of prior physiological knowledge, which is then updated with new experimental data in a Bayesian framework. A covariate model integrates known relationships of physiological parameters to age, gender and body height. We further provide a framework for estimation of the a posteriori parameter dependency structure at the population level. The approach is demonstrated considering a cohort of healthy individuals and theophylline as an application example. The variability and co-variability of physiological parameters are specified within the population; respectively. Significant correlations are identified between population parameters and are applied for individual- and population-specific visual predictive checks of the pharmacokinetic behavior, which leads to improved results compared to present population approaches. In the future, the integration of a generic PBPK model into an hierarchical approach allows for extrapolations to other populations or drugs, while the Bayesian paradigm allows for an iterative application of the approach and thereby a continuous updating of physiological knowledge with new data. This will facilitate decision making e.g. from preclinical to clinical development or extrapolation of PK behavior from healthy to clinically significant populations.

Network pharmacodynamic models for customized cancer therapy

Pharmacokinetics (PKs) and pharmacodynamics (PDs) have always been integral to the design of rational drug dosing regimens. Early on PK-driven approaches came under the auspices of therapeutic drug monitoring that progressed into population-based PK and PK/PD modeling analyses. As the availability of tissue samples for measurement of drug concentrations is limited in patients, the bulk of such model-based methods relied on plasma drug concentrations to both build models and monitor therapy. The continued advances in systems biology and the spawning of systems pharmacology propelled the creation of enhanced PD (ePD) models. One of the main characteristic of ePD models is that they are derived from mechanistically grounded biochemical reaction networks. These models are commonly represented as systems of coupled ordinary differential equations with the ability to tailor each reaction and protein concentration to an individual's genomic/proteomic profile. As patient genomic analyses become more common, many genetic and protein abnormalities can be represented in the ePD models, and thus offer a path toward personalized anticancer therapies. By linking PK models to ePD models, a full spectrum of pharmacological simulation tools is available to design sophisticated multidrug regimens. However, ePD models are not a panacea and face challenges in model identifiability, scaling and parameter estimation. Nonetheless, as new technologies evolve and are coupled with fresh ideas on model implementation, it is likely that ePD and PK/ePD models will be considered a viable enterprise to customize anticancer drug therapy.

Development of a paediatric population-based model of the pharmacokinetics of rivaroxaban

Background

Venous thromboembolism has been increasingly recognised as a clinical problem in the paediatric population. Guideline recommendations for antithrombotic therapy in paediatric patients are based mainly on extrapolation from adult clinical trial data, owing to the limited number of clinical trials in paediatric populations. The oral, direct Factor Xa inhibitor rivaroxaban has been approved in adult patients for several thromboembolic disorders, and its well-defined pharmacokinetic and pharmacodynamic characteristics and efficacy and safety profiles in adults warrant further investigation of this agent in the paediatric population.

Objective

The objective of this study was to develop and qualify a physiologically based pharmacokinetic (PBPK) model for rivaroxaban doses of 10 and 20 mg in adults and to scale this model to the paediatric population (0–18 years) to inform the dosing regimen for a clinical study of rivaroxaban in paediatric patients.

Methods

Experimental data sets from phase I studies supported the development and qualification of an adult PBPK model. This adult PBPK model was then scaled to the paediatric population by including anthropometric and physiological information, age-dependent clearance and age-dependent protein binding. The pharmacokinetic properties of rivaroxaban in virtual populations of children were simulated for two body weight-related dosing regimens equivalent to 10 and 20 mg once daily in adults. The quality of the model was judged by means of a visual predictive check. Subsequently, paediatric simulations of the area under the plasma concentration–time curve (AUC), maximum (peak) plasma drug concentration (C_max) and concentration in plasma after 24 h (C_24h) were compared with the adult reference simulations.

Results

Simulations for AUC, C_max and C_24h throughout the investigated age range largely overlapped with values obtained for the corresponding dose in the adult reference simulation for both body weight-related dosing regimens. However, pharmacokinetic values in infants and preschool children (body weight <40 kg) were lower than the 90 % confidence interval threshold of the adult reference model and, therefore, indicated that doses in these groups may need to be increased to achieve the same plasma levels as in adults. For children with body weight between 40 and 70 kg, simulated plasma pharmacokinetic parameters (C_max, C_24h and AUC) overlapped with the values obtained in the corresponding adult reference simulation, indicating that body weight-related exposure was similar between these children and adults. In adolescents of >70 kg body weight, the simulated 90 % prediction interval values of AUC and C_24h were much higher than the 90 % confidence interval of the adult reference population, owing to the weight-based simulation approach, but for these patients rivaroxaban would be administered at adult fixed doses of 10 and 20 mg.

Conclusion

The paediatric PBPK model developed here allowed an exploratory analysis of the pharmacokinetics of rivaroxaban in children to inform the dosing regimen for a clinical study in paediatric patients.

Electronic supplementary material

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

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.