Best Practice in the Use of Physiologically Based Pharmacokinetic Modeling and Simulation to Address Clinical Pharmacology Regulatory Questions

Use of physiologically-based pharmacokinetic modeling to understand the effect of omeprazole administration on the pharmacokinetics of oral extended-release nifedipine

Proton pump inhibitors (PPIs) can affect the release of drugs from their dosage forms in vivo by elevating the gastric pH. Our recent clinical study has demonstrated that drug-drug interactions (DDIs) exist between a PPI, omeprazole, and nifedipine extended-release formulations, where systemic exposure of nifedipine was increased in subjects after multiple-dose pretreatment of omeprazole. However, the mechanism of the observed DDIs between omeprazole and nifedipine has not been well-understood, as the DDI may also be mediated through CYP3A4 enzyme inhibition in addition to the elevated gastric pH caused by omeprazole. This study used physiologically-based pharmacokinetic (PBPK) modeling and simulations to investigate the underlying mechanism of these complex DDIs. A formulation exhibiting differences in in vitro dissolution across physiological pH range and another formulation where pH does not impact dissolution appreciably (e.g., an osmotic pump) were chosen to characterize the potential impact of pH. The PBPK models incorporated two-stage in vitro release profiles via US Pharmacopeia 2 apparatus. PBPK simulations suggest that the elevated gastric pH following multiple-dose administration of omeprazole has a minimal effect on nifedipine pharmacokinetics (PKs), whereas CYP3A4-mediated DDI is likely the main driver to the observed change of nifedipine PKs in the presence of omeprazole. Compared to the osmotic formulation, the slightly increased exposure of nifedipine can be accounted for by the enhanced drug release in the pH-dependent formulation. The reported model-based approach may be useful in DDI risk assessments, product formulation designs, and bioequivalence evaluations.

A tutorial for model-based evaluation and translation of cardiovascular safety in preclinical trials

Assessment of drug-induced effects on the cardiovascular (CV) system remains a critical component of the drug discovery process enabling refinement of the therapeutic index. Predicting potential drug-related unintended CV effects in the preclinical stage is necessary for first-in-human dose selection and preclusion of adverse CV effects in the clinical stage. According to the current guidelines for small molecules, nonclinical CV safety assessment conducted via telemetry analyses should be included in the safety pharmacology core battery studies. However, the manual for quantitative evaluation of the CV safety signals in animals is available only for electrocardiogram parameters (i.e., QT interval assessment), not for hemodynamic parameters (i.e., heart rate, blood pressure, etc.). Various model-based approaches, including empirical pharmacokinetic-toxicodynamic analyses and systems pharmacology modeling could be used in the framework of telemetry data evaluation. In this tutorial, we provide a comprehensive workflow for the analysis of nonclinical CV safety on hemodynamic parameters with a sequential approach, highlight the challenges associated with the data, and propose respective solutions, complemented with a reproducible example. The work is aimed at helping researchers conduct model-based analyses of the CV safety in animals with subsequent translation of the effect to humans seamlessly and efficiently.

Computational Amendment of Parenteral In Situ Forming Particulates' Characteristics: Design of Experiment and PBPK Physiological Modeling

Lipid and/or polymer-based drug conjugates can potentially minimize side effects by increasing drug accumulation at target sites and thus augment patient compliance. Formulation factors can present a potent influence on the characteristics of the obtained systems. The selection of an appropriate solvent with satisfactory rheological properties, miscibility, and biocompatibility is essential to optimize drug release. This work presents a computational study of the effect of the basic formulation factors on the characteristics of the obtained in situ-forming particulates (IFPs) encapsulating a model drug using a 21.31 full factorial experimental design. The emulsion method was employed for the preparation of lipid and/or polymer-based IFPs. The IFP release profiles and parameters were computed. Additionally, a desirability study was carried out to choose the optimum formulation for further morphological examination, rheological study, and PBPK physiological modeling. Results revealed that the type of particulate forming agent (lipid/polymer) and the incorporation of structure additives like Brij 52 and Eudragit RL can effectively augment the release profile as well as the burst of the drug. The optimized formulation exhibited a pseudoplastic rheological behavior and yielded uniformly spherical-shaped dense particulates with a PS of 573.92 ± 23.5 nm upon injection. Physiological modeling simulation revealed the pioneer pharmacokinetic properties of the optimized formulation compared to the observed data. These results assure the importance of controlling the formulation factors during drug development, the potentiality of the optimized IFPs for the intramuscular delivery of piroxicam, and the reliability of PBPK physiological modeling in predicting the biological performance of new formulations with effective cost management.

Physiologically based pharmacokinetic modelling of semaglutide in children and adolescents with healthy and obese body weights

AIMS

To develop paediatric physiologically based pharmacokinetic modelling (PBPK) models of semaglutide to estimate the pharmacokinetic profile for subcutaneous injections in children and adolescents with healthy and obese body weights.

METHODS

Pharmacokinetic modelling and simulations of semaglutide subcutaneous injections were performed using the Transdermal Compartmental Absorption & Transit model implemented in GastroPlus v.9.5 modules. A PBPK model of semaglutide was developed and verified in the adult population, by comparing the simulated plasma exposure with the observed data, and further scaled to the paediatric populations with normal and obese body weight.

RESULTS

The semaglutide PBPK model was successfully developed in adults and scaled to the paediatric population. Our paediatric PBPK simulations indicated a significant increase in maximum plasma concentrations for the 10-14 years' paediatric population with healthy body weights, which was higher than the observed values in adults at the reference dose. Since gastrointestinal adverse events are related to increased semaglutide concentrations, peak concentrations outside the target range may represent a safety risk for this paediatric age group. Besides, paediatric PBPK models indicated that body weight was inversely related to semaglutide maximum plasma concentration, corroborating the consensus on the influence of body weight on semaglutide PK in adults.

CONCLUSION

Paediatric PBPK was successfully achieved using a top-down approach and drug-related parameters. The development of unprecedented PBPK models will support paediatric clinical therapy for applying aid-safe dosing regimens for the paediatric population in diabetes treatment.

Development and application of a PBPK modeling strategy to support antimalarial drug development

As part of a collaboration between Medicines for Malaria Venture (MMV), Certara UK and Monash University, physiologically-based pharmacokinetic (PBPK) models were developed for 20 antimalarials, using data obtained from standardized in vitro assays and clinical studies within the literature. The models have been applied within antimalarial drug development at MMV for more than 5 years. During this time, a strategy for their impactful use has evolved. All models are described in the supplementary material and are available to researchers. Case studies are also presented, demonstrating real-world development and clinical applications, including the assessment of the drug-drug interaction liability between combination partners or with co-administered drugs. This work emphasizes the benefit of PBPK modeling for antimalarial drug development and decision making, and presents a strategy to integrate it into the research and development process. It also provides a repository of shared information to benefit the global health research community.

Clinical Pharmacology of Brigatinib: A Next-Generation Anaplastic Lymphoma Kinase Inhibitor

Brigatinib, a next-generation anaplastic lymphoma kinase (ALK) inhibitor designed to overcome mechanisms of resistance associated with crizotinib, is approved for the treatment of ALK-positive advanced or metastatic non-small cell lung cancer. After oral administration of single doses of brigatinib 30-240 mg, the median time to reach maximum plasma concentration ranged from 1 to 4 h. In patients with advanced malignancies, brigatinib showed dose linearity over the dose range of 60-240 mg once daily. A high-fat meal had no clinically meaningful effect on systemic exposures of brigatinib (area under the plasma concentration-time curve); thus, brigatinib can be administered with or without food. In a population pharmacokinetic analysis, a three-compartment pharmacokinetic model with transit absorption compartments was found to adequately describe brigatinib pharmacokinetics. In addition, the population pharmacokinetic analyses showed that no dose adjustment is required based on body weight, age, race, sex, total bilirubin (< 1.5× upper limit of normal), and mild-to-moderate renal impairment. Data from dedicated phase I trials have indicated that no dose adjustment is required for patients with mild or moderate hepatic impairment, while a dose reduction of approximately 40% (e.g., from 180 to 120 mg) is recommended for patients with severe hepatic impairment, and a reduction of approximately 50% (e.g., from 180 to 90 mg) is recommended when administering brigatinib to patients with severe renal impairment. Brigatinib is primarily metabolized by cytochrome P450 (CYP) 3A, and results of clinical drug-drug interaction studies and physiologically based pharmacokinetic analyses have demonstrated that coadministration of strong or moderate CYP3A inhibitors or inducers with brigatinib should be avoided. If coadministration with a strong or moderate CYP3A inhibitor cannot be avoided, the dose of brigatinib should be reduced by approximately 50% (strong CYP3A inhibitor) or approximately 40% (moderate CYP3A inhibitor), respectively. Brigatinib is a weak inducer of CYP3A in vivo; data from a phase I drug-drug interaction study showed that coadministration of brigatinib 180 mg once daily reduced the oral midazolam area under the plasma concentration-time curve from time zero to infinity by approximately 26%. Brigatinib did not inhibit CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, or CYP2D6 at clinically relevant concentrations in vitro. Exposure-response analyses based on data from the ALTA (ALK in Lung Cancer Trial of AP26113) and ALTA-1L pivotal trials of brigatinib confirm the favorable benefit versus risk profile of the approved titration dosing regimen of 180 mg once daily (after a 7-day lead-in at 90 mg once daily).

Current State and Challenges of Physiologically Based Biopharmaceutics Modeling (PBBM) in Oral Drug Product Development

Physiologically based biopharmaceutics modeling (PBBM) emphasizes the integration of physicochemical properties of drug substance and formulation characteristics with system physiological parameters to predict the absorption and pharmacokinetics (PK) of a drug product. PBBM has been successfully utilized in drug development from discovery to postapproval stages and covers a variety of applications. The use of PBBM facilitates drug development and can reduce the number of preclinical and clinical studies. In this review, we summarized the major applications of PBBM, which are classified into six categories: formulation selection and development, biopredictive dissolution method development, biopharmaceutics risk assessment, clinically relevant specification settings, food effect evaluation and pH-dependent drug-drug-interaction risk assessment. The current state of PBBM applications is illustrated with examples from published studies for each category of application. Despite the variety of PBBM applications, there are still many hurdles limiting the use of PBBM in drug development, that are associated with the complexity of gastrointestinal and human physiology, the knowledge gap between the in vitro and the in vivo behavior of drug products, the limitations of model interfaces, and the lack of agreed model validation criteria, among other issues. The challenges and essential considerations related to the use of PBBM are discussed in a question-based format along with the scientific thinking on future research directions. We hope this review can foster open discussions between the pharmaceutical industry and regulatory agencies and encourage collaborative research to fill the gaps, with the ultimate goal to maximize the applications of PBBM in oral drug product development.

Regulatory utility of mechanistic modeling to support alternative bioequivalence approaches: A workshop overview

On September 30 and October 1, 2021, the US Food and Drug Administration (FDA) and the Center for Research on Complex Generics cosponsored a live virtual workshop titled "Regulatory Utility of Mechanistic Modeling to Support Alternative Bioequivalence Approaches." The overall aims of the workshop included (i) engaging the generic drug industry and other involved stakeholders regarding how mechanistic modeling and simulation can support their product development and regulatory submissions; (ii) sharing the current state of mechanistic modeling for bioequivalence (BE) assessment through case studies; (iii) establishing a consensus on best practices for using mechanistic modeling approaches, such as physiologically based pharmacokinetic modeling and computational fluid dynamics modeling, for BE assessment; and (iv) introducing the concept of a Model Master File to improve model sharing between model developers, industry, and the FDA. More than 1500 people registered for the workshop. Based on a postworkshop survey, the majority of participants reported that their fundamental scientific understanding of mechanistic models was enhanced, there was greater consensus on model validation and verification, and regulatory expectations for mechanistic modeling submitted in abbreviated new drug applications were clarified by the workshop.

Design and conduct considerations for studies in patients with hepatic impairment

Despite the liver being the primary site for clearance of xenobiotics utilizing a myriad of mechanisms ranging from cytochrome P450 enzyme pathways, glucuronidation, and biliary excretion, there is a dearth of information available as to how the severity of hepatic impairment (HI) can alter drug absorption and disposition (i.e., pharmacokinetics [PK]) as well as their efficacy and safety or pharmacodynamics (PD). In general, regulatory agencies recommend conducting PK studies in subjects with HI when hepatic metabolism/excretion accounts for more than 20% of drug elimination or if the drug has a narrow therapeutic range. In this tutorial, we provide an overview of the global regulatory landscape, clinical measures for hepatic function assessment, methods to stage HI severity, and consequently the impact on labeling. In addition, we provide an in-depth practical guidance for designing and conducting clinical trials for patients with HI and on the application of modeling and simulation strategies in lieu of dedicated trials for dosing recommendations in patients with HI.

Investigation of CYP3A induction by PF-05251749 in early clinical development: comparison of linear slope physiologically based pharmacokinetic prediction and biomarker response

PF-05251749 is a dual inhibitor of casein kinase 1 δ/ε under clinical development to treat disruption of circadian rhythm in Alzheimer's and Parkinson's diseases. In vitro, PF-05251749 (0.3-100 μM) induced CYP3A in cryopreserved human hepatocytes, demonstrating non-saturable, dose-dependent CYP3A mRNA increases, with induction slopes in the range 0.036-0.39 μM^-1 . In a multiple-dose study (B8001002) in healthy participants, CYP3A activity was explored by measuring changes in 4β-hydroxycholesterol/cholesterol ratio. Following repeated oral administration of PF-05251749, up to 400 mg q.d., no significant changes were observed in 4β-hydroxycholesterol/cholesterol ratio; this ratio increased significantly (~1.5-fold) following administration of PF-05251749 at 750 mg q.d., suggesting potential CYP3A induction at this dose. Physiologically based pharmacokinetic (PBPK) models were developed to characterize the observed clinical pharmacokinetics (PK) of PF-05251749 at 400 and 750 mg q.d.; the PBPK induction model was calibrated using the in vitro linear fit induction slope, with rifampin as reference compound (Ind_max  = 8, EC₅₀  = 0.32 μM). Clinical trial simulation following co-administration of PF-05251749, 400 mg q.d. with oral midazolam 2 mg, predicted no significant drug interaction risk. PBPK model predicted weak drug interaction following co-administration of PF-05251749, 750 mg q.d. with midazolam 2 mg. In conclusion, good agreement was obtained between CYP3A drug interaction risk predicted using linear-slope PBPK model and exploratory biomarker trends. This agreement between two orthogonal approaches enabled assessment of drug interaction risks of PF-05251749 in early clinical development, in the absence of a clinical drug-drug interaction study.

Development and Evaluation of a Physiologically Based Pharmacokinetic Model for Predicting Haloperidol Exposure in Healthy and Disease Populations

The physiologically based pharmacokinetic (PBPK) approach can be used to develop mathematical models for predicting the absorption, distribution, metabolism, and elimination (ADME) of administered drugs in virtual human populations. Haloperidol is a typical antipsychotic drug with a narrow therapeutic index and is commonly used in the management of several medical conditions, including psychotic disorders. Due to the large interindividual variability among patients taking haloperidol, it is very likely for them to experience either toxic or subtherapeutic effects. We intend to develop a haloperidol PBPK model for identifying the potential sources of pharmacokinetic (PK) variability after intravenous and oral administration by using the population-based simulator, PK-Sim. The model was initially developed and evaluated to predict the PK of haloperidol and its reduced metabolite in adult healthy population after intravenous and oral administration. After evaluating the developed PBPK model in healthy adults, it was used to predict haloperidol–rifampicin drug–drug interaction and was extended to tuberculosis patients. The model evaluation was performed using visual assessments, prediction error, and mean fold error of the ratio of the observed-to-predicted values of the PK parameters. The predicted PK values were in good agreement with the corresponding reported values. The effects of the pathophysiological changes and enzyme induction associated with tuberculosis and its treatment, respectively, on haloperidol PK, have been predicted precisely. For all clinical scenarios that were evaluated, the predicted values were within the acceptable two-fold error range.

Insights into the mapping of green synthesis conditions for ZnO nanoparticles and their toxicokinetics

Research on ZnO nanoparticles (NPs) has broad medical applications. However, the green synthesis of ZnO NPs involves a wide range of properties requiring optimization. ZnO NPs show toxicity at lower doses. This toxicity is a function of NP properties and pharmacokinetics. Moreover, NP toxicity and pharmacokinetics are affected by the species type and age of the animals tested. Physiologically based pharmacokinetic (PBPK) modeling offers a mechanistic platform to scrutinize the colligative effect of the interplay between these factors, which reduces the need for in vivo studies. This review provides a guide to choosing green synthesis conditions that result in minimal toxicity using a mechanistic tool, namely PBPK modeling.

Guide to development of compound files for PBPK modeling in the Simcyp population-based simulator

The Simcyp Simulator is a software platform for population physiologically‐based pharmacokinetic (PBPK) modeling and simulation. It links in vitro data to in vivo absorption, distribution, metabolism, excretion and pharmacokinetic/pharmacodynamic outcomes to explore clinical scenarios and support drug development decisions, including regulatory submissions and drug labels. This tutorial describes the different input parameters required, as well as the considerations needed when developing a PBPK model within the Simulator, for a small molecule intended for oral administration. A case study showing the development and application of a PBPK model for ondansetron is herein used to aid the understanding of different PBPK model development concepts.

Prediction of CYP‐mediated DDIs involving inhibition: Approaches to address the requirements for system qualification of the Simcyp Simulator

Physiologically‐based pharmacokinetic (PBPK) modeling is being increasingly used in drug development to avoid unnecessary clinical drug–drug interaction (DDI) studies and inform drug labels. Thus, regulatory agencies are recommending, or indeed requesting, more rigorous demonstration of the prediction accuracy of PBPK platforms in the area of their intended use. We describe a framework for qualification of the Simcyp Simulator with respect to competitive and mechanism‐based inhibition (MBI) of CYP1A2, CYP2D6, CYP2C8, CYP2C9, CYP2C19, and CYP3A4/5. Initially, a DDI matrix, consisting of a range of weak, moderate, and strong inhibitors and substrates with varying fraction metabolized by specific CYP enzymes that were susceptible to different degrees of inhibition, were identified. Simulations were run with 123 clinical DDI studies involving competitive inhibition and 78 clinical DDI studies involving MBI. For competitive inhibition, the overall prediction accuracy was good with an average fold error (AFE) of 0.91 and 0.92 for changes in the maximum plasma concentration (C_max) and area under the plasma concentration (AUC) time profile, respectively, as a consequence of the DDI. For MBI, an AFE of 1.03 was determined for both C_max and AUC. The prediction accuracy was generally comparable across all CYP enzymes, irrespective of the isozyme and mechanism of inhibition. These findings provide confidence in application of the Simcyp Simulator (V19 R1) for assessment of the DDI potential of drugs in development either as inhibitors or victim drugs of CYP‐mediated interactions. The approach described herein and the identified DDI matrix can be used to qualify subsequent versions of the platform.

Recent advances in the translation of drug metabolism and pharmacokinetics science for drug discovery and development

Drug metabolism and pharmacokinetics (DMPK) is an important branch of pharmaceutical sciences. The nature of ADME (absorption, distribution, metabolism, excretion) and PK (pharmacokinetics) inquiries during drug discovery and development has evolved in recent years from being largely descriptive to seeking a more quantitative and mechanistic understanding of the fate of drug candidates in biological systems. Tremendous progress has been made in the past decade, not only in the characterization of physiochemical properties of drugs that influence their ADME, target organ exposure, and toxicity, but also in the identification of design principles that can minimize drug-drug interaction (DDI) potentials and reduce the attritions. The importance of membrane transporters in drug disposition, efficacy, and safety, as well as the interplay with metabolic processes, has been increasingly recognized. Dramatic increases in investments on new modalities beyond traditional small and large molecule drugs, such as peptides, oligonucleotides, and antibody-drug conjugates, necessitated further innovations in bioanalytical and experimental tools for the characterization of their ADME properties. In this review, we highlight some of the most notable advances in the last decade, and provide future perspectives on potential major breakthroughs and innovations in the translation of DMPK science in various stages of drug discovery and development.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

PBPK Modeling to Simulate the Fate of Compounds in Living Organisms

Pharmacokinetics study the fate of xenobiotics in a living organism. Physiologically based pharmacokinetic (PBPK) models provide realistic descriptions of xenobiotics' absorption, distribution, metabolism, and excretion processes. They model the body as a set of homogeneous compartments representing organs, and their parameters refer to anatomical, physiological, biochemical, and physicochemical entities. They offer a quantitative mechanistic framework to understand and simulate the time-course of the concentration of a substance in various organs and body fluids. These models are well suited for performing extrapolations inherent to toxicology and pharmacology (e.g., between species or doses) and for integrating data obtained from various sources (e.g., in vitro or in vivo experiments, structure-activity models). In this chapter, we describe the practical development and basic use of a PBPK model from model building to model simulations, through implementation with an easily accessible free software.

Integration of advanced methods and models to study drug absorption and related processes: An UNGAP perspective

This collection of contributions from the European Network on Understanding Gastrointestinal Absorption-related Processes (UNGAP) community assembly aims to provide information on some of the current and newer methods employed to study the behaviour of medicines. It is the product of interactions in the immediate pre-Covid period when UNGAP members were able to meet and set up workshops and to discuss progress across the disciplines. UNGAP activities are divided into work packages that cover special treatment populations, absorption processes in different regions of the gut, the development of advanced formulations and the integration of food and pharmaceutical scientists in the food-drug interface. This involves both new and established technical approaches in which we have attempted to define best practice and highlight areas where further research is needed. Over the last months we have been able to reflect on some of the key innovative approaches which we were tasked with mapping, including theoretical, in silico, in vitro, in vivo and ex vivo, preclinical and clinical approaches. This is the product of some of us in a snapshot of where UNGAP has travelled and what aspects of innovative technologies are important. It is not a comprehensive review of all methods used in research to study drug dissolution and absorption, but provides an ample panorama of current and advanced methods generally and potentially useful in this area. This collection starts from a consideration of advances in a priori approaches: an understanding of the molecular properties of the compound to predict biological characteristics relevant to absorption. The next four sections discuss a major activity in the UNGAP initiative, the pursuit of more representative conditions to study lumenal dissolution of drug formulations developed independently by academic teams. They are important because they illustrate examples of in vitro simulation systems that have begun to provide a useful understanding of formulation behaviour in the upper GI tract for industry. The Leuven team highlights the importance of the physiology of the digestive tract, as they describe the relevance of gastric and intestinal fluids on the behaviour of drugs along the tract. This provides the introduction to microdosing as an early tool to study drug disposition. Microdosing in oncology is starting to use gamma-emitting tracers, which provides a link through SPECT to the next section on nuclear medicine. The last two papers link the modelling approaches used by the pharmaceutical industry, in silico to Pop-PK linking to Darwich and Aarons, who provide discussion on pharmacometric modelling, completing the loop of molecule to man.

Physiologically based pharmacokinetic (PBPK) modelling of tamsulosin related to CYP2D6*10 allele

Tamsulosin, a selective [Formula: see text]-adrenoceptor blocker, is commonly used for alleviation of lower urinary tract symptoms related to benign prostatic hyperplasia. Tamsulosin is predominantly metabolized by CYP3A4 and CYP2D6 enzymes, and several studies reported the effects of CYP2D6 genetic polymorphism on the pharmacokinetics of tamsulosin. This study aims to develop and validate the physiologically based pharmacokinetic (PBPK) model of tamsulosin in CYP2D6*wt/*wt, CYP2D6*wt/*10, and CYP2D6*10/*10 genotypes, using Simcyp® simulator. Physicochemical, and formulation properties and data for absorption, distribution, metabolism and excretion were collected from previous publications, predicted in the simulator, or optimized in different CYP2D6 genotypes. The tamsulosin PBPK model in CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes were developed based on the clinical pharmacokinetic study where a single oral dose of 0.2 mg tamsulosin was administered to 25 healthy Korean male volunteers with CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes. A previous pharmacokinetic study was used to develop the model in CYP2D6*10/*10 genotype. The developed model was validated using other clinical pharmacokinetic studies not used in development. The predicted exposures via the PBPK model in CYP2D6*wt/*10 and CYP2D6*10/*10 genotype was 1.23- and 1.76-fold higher than CYP2D6*wt/*wt genotype, respectively. The simulation profiles were visually similar to the observed profiles, and fold errors of all development and validation datasets were included within the criteria. Therefore, the tamsulosin PBPK model in different CYP2D6 genotypes with regards to CYP2D6*10 alleles was appropriately established. Our model can contribute to the implementation of personalized pharmacotherapy of patients, appropriately predicting the pharmacokinetics of tamsulosin reflecting their demographic and CYP2D6 genotype characteristics without unnecessary drug exposure.

Expression dynamics of pregnane X receptor-controlled genes in 3D primary human hepatocyte spheroids

The pregnane X receptor (PXR) is a ligand-activated nuclear receptor controlling hepatocyte expression of numerous genes. Although expression changes in xenobiotic-metabolizing, lipogenic, gluconeogenic and bile acid synthetic genes have been described after PXR activation, the temporal dynamics of their expression is largely unknown. Recently, 3D spheroids of primary human hepatocytes (PHHs) have been characterized as the most phenotypically relevant hepatocyte model. We used 3D PHHs to assess time-dependent expression profiles of 12 prototypic PXR-controlled genes in the time course of 168 h of rifampicin treatment (1 or 10 µM). We observed a similar bell-shaped time-induction pattern for xenobiotic-handling genes (CYP3A4, CYP2C9, CYP2B6, and MDR1). However, we observed either biphasic profiles for genes involved in endogenous metabolism (FASN, GLUT2, G6PC, PCK1, and CYP7A1), a decrease for SHP or oscillation for PDK4 and PXR. The rifampicin concentration determined the expression profiles for some genes. Moreover, we calculated half-lives of CYP3A4 and CYP2C9 mRNA under induced or basal conditions and we used a mathematical model to describe PXR-mediated regulation of CYP3A4 expression employing 3D PHHs. The study shows the importance of long-term time-expression profiling of PXR target genes in phenotypically stable 3D PHHs and provides insight into PXR function in liver beyond our knowledge from conventional 2D in vitro models.

A Bayesian population physiologically based pharmacokinetic absorption modeling approach to support generic drug development: application to bupropion hydrochloride oral dosage forms

We propose a Bayesian population modeling and virtual bioequivalence assessment approach to establishing dissolution specifications for oral dosage forms. A generalizable semi-physiologically based pharmacokinetic absorption model with six gut segments and liver, connected to a two-compartment model of systemic disposition for bupropion hydrochloride oral dosage forms was developed. Prior information on model parameters for gut physiology, bupropion physicochemical properties, and drug product properties were obtained from the literature. The release of bupropion hydrochloride from immediate-, sustained- and extended-release oral dosage forms was described by a Weibull function. In vitro dissolution data were used to assign priors to the in vivo release properties of the three bupropion formulations. We applied global sensitivity analysis to identify the influential parameters for plasma bupropion concentrations and calibrated them. To quantify inter- and intra-individual variability, plasma concentration profiles in healthy volunteers that received the three dosage forms, each at two doses, were used. The calibrated model was in good agreement with both in vitro dissolution and in vivo exposure data. Markov Chain Monte Carlo samples from the joint posterior parameter distribution were used to simulate virtual crossover clinical trials for each formulation with distinct drug dissolution profiles. For each trial, an allowable range of dissolution parameters ("safe space") in which bioequivalence can be anticipated was established. These findings can be used to assure consistent product performance throughout the drug product life-cycle and to support manufacturing changes. Our framework provides a comprehensive approach to support decision-making in drug product development.

Physiologically-based pharmacokinetic model predictions of inter-ethnic differences in imatinib pharmacokinetics and dosing regimens

AIMS

This study implements a physiologically-based pharmacokinetic (PBPK) modelling approach to investigate inter-ethnic differences in imatinib pharmacokinetics and dosing regimens.

METHODS

A PBPK model of imatinib was built in the Simcyp Simulator (version 17) integrating in vitro drug metabolism and clinical pharmacokinetic data. The model accounts for ethnic differences in body size and abundance of drug-metabolising enzymes and proteins involved in imatinib disposition. Utility of this model for prediction of imatinib pharmacokinetics was evaluated across different dosing regimens and ethnic groups. The impact of ethnicity on imatinib dosing was then assessed based on the established range of trough concentrations (C_ss,min ).

RESULTS

The PBPK model of imatinib demonstrated excellent predictive performance in describing pharmacokinetics and the attained C_ss,min in patients from different ethnic groups, shown by prediction differences that were within 1.25-fold of the clinically-reported values in published studies. PBPK simulation suggested a similar dose of imatinib (400-600 mg/d) to achieve the desirable range of C_ss,min (1000-3200 ng/mL) in populations of European, Japanese and Chinese ancestry. The simulation indicated that patients of African ancestry may benefit from a higher initial dose (600-800 mg/d) to achieve imatinib target concentrations, due to a higher apparent clearance (CL/F) of imatinib compared to other ethnic groups; however, the clinical data to support this are currently limited.

CONCLUSION

PBPK simulations highlighted a potential ethnic difference in the recommended initial dose of imatinib between populations of European and African ancestry, but not populations of Chinese and Japanese ancestry.

Physiologically-Based Pharmacokinetic Modelling of Entrectinib Parent and Active Metabolite to Support Regulatory Decision-Making

BACKGROUND AND OBJECTIVE

Entrectinib is a selective inhibitor of ROS1/TRK/ALK kinases, recently approved for oncology indications. Entrectinib is predominantly cleared by cytochrome P450 (CYP) 3A4, and modulation of CYP3A enzyme activity profoundly alters the pharmacokinetics of both entrectinib and its active metabolite M5. We describe development of a combined physiologically based pharmacokinetic (PBPK) model for entrectinib and M5 to support dosing recommendations when entrectinib is co-administered with CYP3A4 inhibitors or inducers.

METHODS

A PBPK model was established in Simcyp^® Simulator. The initial model based on in vitro-in vivo extrapolation was refined using sensitivity analysis and non-linear mixed effects modeling to optimize parameter estimates and to improve model fit to data from a clinical drug-drug interaction study with the strong CYP3A4 inhibitor, itraconazole. The model was subsequently qualified against clinical data, and the final qualified model used to simulate the effects of moderate to strong CYP3A4 inhibitors and inducers on entrectinib and M5 pharmacokinetics.

RESULTS

The final model showed good predictive performance for entrectinib and M5, meeting commonly used predictive performance acceptance criteria in each case. The model predicted that co-administration of various moderate CYP3A4 inhibitors (verapamil, erythromycin, clarithromycin, fluconazole, and diltiazem) would result in an average increase in entrectinib exposure between 2.2- and 3.1-fold, with corresponding average increases for M5 of approximately 2-fold. Co-administration of moderate CYP3A4 inducers (efavirenz, carbamazepine, phenytoin) was predicted to result in an average decrease in entrectinib exposure between 45 and 79%, with corresponding average decreases for M5 of approximately 50%.

CONCLUSIONS

The model simulations were used to derive dosing recommendations for co-administering entrectinib with CYP3A4 inhibitors or inducers. PBPK modeling has been used in lieu of clinical studies to enable regulatory decision-making.

A Mechanistic Site-Of-Action Model: A Tool for Informing Right Target, Right Compound, And Right Dose for Therapeutic Antagonistic Antibody Programs

Quantitative modeling is increasingly utilized in the drug discovery and development process, from the initial stages of target selection, through clinical studies. The modeling can provide guidance on three major questions–is this the right target, what are the right compound properties, and what is the right dose for moving the best possible candidate forward. In this manuscript, we present a site-of-action modeling framework which we apply to monoclonal antibodies against soluble targets. We give a comprehensive overview of how we construct the model and how we parametrize it and include several examples of how to apply this framework for answering the questions postulated above. The utilities and limitations of this approach are discussed.

Comparison between physiologically based pharmacokinetic and population pharmacokinetic modelling to select paediatric doses of gepotidacin in plague

Aims

To develop physiologically based pharmacokinetic (PBPK) and population pharmacokinetic (PopPK) models to predict effective doses of gepotidacin in paediatrics for the treatment of pneumonic plague (Yersinia pestis).

Methods

A gepotidacin PBPK model was constructed using a population‐based absorption, distribution, metabolism and excretion simulator, Simcyp®, with physicochemical and in vitro data, optimized with clinical data from a dose‐escalation intravenous (IV) study and a human mass balance study. A PopPK model was developed with pooled PK data from phase 1 studies with IV gepotidacin in healthy adults.

Results

For both the PopPK and PBPK models, body weight was found to be a key covariate affecting gepotidacin clearance. With PBPK, ~90% of the predicted PK for paediatrics fell between the 5th and 95th percentiles of adult values except for subjects weighing ≤5 kg. PopPK‐simulated paediatric means for C_max and AUC_(0‐τ) were similar to adult exposures across various weight brackets. The proposed dosing regimens were weight‐based for subjects ≤40 kg and fixed‐dose for subjects >40 kg. Comparison of observed and predicted exposures in adults indicated that both PBPK and PopPK models achieved similar AUC and C_max for a given dose, but the C_max predictions with PopPK were slightly higher than with PBPK. The two models differed on dose predictions in children <3 months old. The PopPK model may be suboptimal for low age groups due to the absence of maturation characterization of drug‐metabolizing enzymes involved with clearance in adults.

Conclusions

Both PBPK and PopPK approaches can reasonably predict gepotidacin exposures in children.

Scientific and regulatory evaluation of mechanistic in silico drug and disease models in drug development: Building model credibility

The value of in silico methods in drug development and evaluation has been demonstrated repeatedly and convincingly. While their benefits are now unanimously recognized, international standards for their evaluation, accepted by all stakeholders involved, are still to be established. In this white paper, we propose a risk‐informed evaluation framework for mechanistic model credibility evaluation. To properly frame the proposed verification and validation activities, concepts such as context of use, regulatory impact and risk‐based analysis are discussed. To ensure common understanding between all stakeholders, an overview is provided of relevant in silico terminology used throughout this paper. To illustrate the feasibility of the proposed approach, we have applied it to three real case examples in the context of drug development, using a credibility matrix currently being tested as a quick‐start tool by regulators. Altogether, this white paper provides a practical approach to model evaluation, applicable in both scientific and regulatory evaluation contexts.

Physiologically-Based Pharmacokinetic Modeling to Support Determination of Bioequivalence for Dermatological Drug Products: Scientific and Regulatory Considerations

Physiologically-based pharmacokinetic (PBPK) modeling and simulation provides mechanism-based predictions of the pharmacokinetics of an active ingredient following its administration in humans. Dermal PBPK models describe the skin permeation and disposition of the active ingredient following the application of a dermatological product on the skin of virtual healthy and diseased human subjects. These models take into account information on product quality attributes, physicochemical properties of the active ingredient and skin (patho)physiology, and their interplay with each other. Regulatory and product development decision makers can leverage these quantitative tools to identify factors impacting local and systemic exposure. In the realm of generic drug products, the number of US Food and Drug Administratioin (FDA) interactions that use dermal PBPK modeling to support alternative bioequivalence (BE) approaches is increasing. In this report, we share scientific considerations on the development, verification and validation (V&V), and application of PBPK models within the context of a virtual BE assessment for dermatological drug products. We discuss the challenges associated with model V&V for these drug products stemming from the fact that target-site active ingredient concentrations are typically not measurable. Additionally, there are no established relationships between local and systemic PK profiles, when the latter are quantifiable. To that end, we detail a multilevel model V&V approach involving validation for the model of the drug product of interest coupled with the overall assessment of the modeling platform in use while leveraging in vitro and in vivo data related to local and systemic bioavailability.

In vitro and In silico biopharmaceutic regulatory guidelines for generic bioequivalence for oral products: Comparison among various regulatory agencies

Generic drug development is a complex process that involves development of formulation similar to reference product. Because of the complexity associated with generic drug development, many regulatory agencies have come up with various guidelines. Out of many guidelines, the biopharmaceutics classification system that was introduced in 1995 based on aqueous solubility and permeability helped many pharmaceutical scientists across the globe to utilize the tool for formulation development, waiver of in vivo studies. Later on in vitro guidelines based on dissolution and in vitro in vivo correlation were introduced by many regulatory agencies with an intent to reduce number of in vivo human testing thereby facilitating shorter development time and faster approvals and launch. Most recently, understanding the importance in silico approaches such as physiologically based pharmacokinetic modelling, regulatory agencies such as United States Food and Drug Administration (USFDA) and European Middle East and Africa (EMA) came up with modelling guidance documents. Even though consensus exists between guidance documents from various regulatory agencies, still there are many minor to major differences exists between these guidance documents that needs to be considered while submitting a generic drug application. This review aims to compare all the in vitro and in silico guidance documents from major regulatory agencies with emphasis on latest trends and technologies combined with regulatory acceptability with an intention to harmonize regulations. Guidance documents from major regulatory agencies such as USFDA, EMA, World Health Organization, International Council for Harmonization and other emerging markets were compared. Similarities &differences among these guidance documents are critically reviewed to provide the reader a detailed overview of these guidance documents at one place.

Characteristics of mathematical modeling languages that facilitate model reuse in systems biology: a software engineering perspective

Reuse of mathematical models becomes increasingly important in systems biology as research moves toward large, multi-scale models composed of heterogeneous subcomponents. Currently, many models are not easily reusable due to inflexible or confusing code, inappropriate languages, or insufficient documentation. Best practice suggestions rarely cover such low-level design aspects. This gap could be filled by software engineering, which addresses those same issues for software reuse. We show that languages can facilitate reusability by being modular, human-readable, hybrid (i.e., supporting multiple formalisms), open, declarative, and by supporting the graphical representation of models. Modelers should not only use such a language, but be aware of the features that make it desirable and know how to apply them effectively. For this reason, we compare existing suitable languages in detail and demonstrate their benefits for a modular model of the human cardiac conduction system written in Modelica.

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.

Design and conduct considerations for studies in patients with impaired renal function

An impaired renal function, including acute and chronic kidney disease and end‐stage renal disease, can be the result of aging, certain disease conditions, the use of some medications, or as a result of smoking. In patients with renal impairment (RI), the pharmacokinetics (PKs) of drugs or drug metabolites may change and result in increased safety risks or decreased efficacy. In order to make specific dose recommendations in the label of drugs for patients with RI, a clinical trial may have to be conducted or, when not feasible, modeling and simulations approaches, such as population PK modeling or physiologically‐based PK modelling may be applied. This tutorial aims to provide an overview of the global regulatory landscape and a practical guidance for successfully designing and conducting clinical RI trials or, alternatively, on applying modeling and simulation tools to come to a dose recommendation for patients with RI in the most efficient manner.

Current Evidence, Challenges, and Opportunities of Physiologically Based Pharmacokinetic Models of Atorvastatin for Decision Making

Atorvastatin (ATS) is the gold-standard treatment worldwide for the management of hypercholesterolemia and prevention of cardiovascular diseases associated with dyslipidemia. Physiologically based pharmacokinetic (PBPK) models have been positioned as a valuable tool for the characterization of complex pharmacokinetic (PK) processes and its extrapolation in special sub-groups of the population, leading to regulatory recognition. Several PBPK models of ATS have been published in the recent years, addressing different aspects of the PK properties of ATS. Therefore, the aims of this review are (i) to summarize the physicochemical and pharmacokinetic characteristics involved in the time-course of ATS, and (ii) to evaluate the major highlights and limitations of the PBPK models of ATS published so far. The PBPK models incorporate common elements related to the physicochemical aspects of ATS. However, there are important differences in relation to the analyte evaluated, the type and effect of transporters and metabolic enzymes, and the permeability value used. Additionally, this review identifies major processes (lactonization, P-gp contribution, ATS-Ca solubility, simultaneous management of multiple analytes, and experimental evidence in the target population), which would enhance the PBPK model prediction to serve as a valid tool for ATS dose optimization.

Food constituent- and herb-drug interactions in oncology: Influence of quantitative modelling on Drug labelling

AIMS

Herbal products, spices and/or fruits are perceived as inherently healthy; for instance, St. John's wort (SJW) is marketed as a natural antidepressant and patients often self-administer it concomitantly with oncology medications. However, food constituents/herbs can interfere with drug pharmacokinetics, with risk of altering pharmacodynamics and efficacy. The objective of this work was to develop a strategy to prioritize herb- or food constituent-drug interactions (FC-DIs) to better assess oncology drug clinical risk.

METHODS

Physiologically based pharmacokinetic (PBPK) models were developed by integrating in vitro parameters with the clinical pharmacokinetics of food constituents in grapefruit juice (bergamottin), turmeric (curcumin) or SJW (hyperforin). Perpetrator files were linked to verified victim PBPK models through appropriate interaction mechanisms (cytochrome P450 3A, breast cancer resistance protein, P-glycoprotein) and applied in prospective PBPK simulations to inform the likelihood and magnitude of changes in exposure to osimertinib, olaparib or acalabrutinib.

RESULTS

Reported FC-DIs with oncology drugs were well recovered, with absolute average fold error values of 1.10 (bergamottin), 1.05 (curcumin) and 1.01 (hyperforin). Prospective simulations with grapefruit juice and turmeric showed clinically minor to insignificant changes in exposure (<1.50-fold) to acalabrutinib, osimertinib and olaparib, but predicted 1.57-fold FC-DI risk between acalabrutinib and curcumin. Moderate DDI risk was expected when acalabrutinib, osimertinib or olaparib were dosed with SJW.

CONCLUSIONS

A model-informed decision tree based on mechanistic understanding of transporter and/or enzyme-mediated FC-DI is proposed based on bergamottin, curcumin and hyperforin FC-DI clinical data. Adopting this quantitative modelling approach should streamline herbal product safety assessments, assist in FC-DI management, and ultimately promote safe clinical use of oncology drugs.

Physiological-based pharmacokinetic modeling trends in pharmaceutical drug development over the last 20-years; in-depth analysis of applications, organizations, and platforms

We assess the advancement of physiologically based pharmacokinetic (PBPK) modeling and simulation (M&S) over the last 20 years (start of 2000 to end of 2019) focusing on the trends in each decade with the relative contributions from different organizations, areas of applications, and software tools used. Unlike many of the previous publications which focused on regulatory applications, our analysis is based on PBPK publications in peer-reviewed journals based on a large sample (>700 original articles). We estimated a rate of growth for PBPK (>40 fold/20 years) that was much steeper than the general pharmacokinetic modeling (<3 fold/20 years) or overall scientific publications (∼3 fold/20 years). The analyses demonstrated that contrary to commonly held belief, commercial specialized PBPK platforms with graphical-user interface were a much more popular choice than open-source alternatives even within academic organizations. These platforms constituted 81% of the whole set of the sample we assessed. The major PBPK applications (top 3) were associated with the study design, predicting formulation effects, and metabolic drug-drug interactions, while studying the fate of drugs in special populations, predicting kinetics in early drug development, and investigating transporter drug interactions have increased proportionally over the last decade. The proportions of application areas based on published research were distinctively different from those shown previously for the regulatory submissions and impact on labels. This may demonstrate the lag time between the research applications versus verified usage within the regulatory framework. The report showed the trend of overall PBPK publications in pharmacology drug development from the past 2 decades stratified by the organizations involved, software used, and area of applications. The analysis showed a more rapid increase in PBPK than that of the pharmacokinetic space itself with an equal contribution from academia and industry. By establishing and recording the journey of PBPK modeling in the past and looking at its current status, the analysis can be used for devising plans based on the anticipated trajectory of future regulatory applications.

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

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

Regulatory Considerations for the Mother, Fetus and Neonate in Fetal Pharmacology Modeling

The regulatory framework for considering the fetal effects of new drugs is limited. This is partially due to the fact that pediatric regulations (21 CFR subpart D) do not apply to the fetus, and only US Health and Human Service (HHS) regulations apply to the fetus. The HHS regulation 45 CFR Part 46 Subpart B limits research approvable by an institutional review board to research where the risk to the fetus is minimal unless the research holds out the prospect of a direct benefit to the fetus or the pregnant woman (45 CFR 46.204). Research that does not meet these requirements, but presents an opportunity to understand, prevent, or alleviate a serious problem affecting the health of pregnant women, fetuses, or neonates, may be permitted by the Secretary of the HHS after expert panel consultation and opportunity for public review and comment (45 CFR 46.407). If the product is regulated by the US Food and Drug Administration (FDA), FDA may get involved in the review process. The FDA does however have a Reviewer Guidance on Evaluating the Risks of Drug Exposure in Human Pregnancies from 2005 and this guidance does discuss the intensity of drug exposure. Estimation of that exposure using physiologically based pharmacokinetic (PBPK) modeling has been suggested by some investigators. Given that drug exposure during pregnancy will impact the fetus, a number of new guidances in the last 2 years also address inclusion of pregnant women in clinical drug trials. Therefore, the drug-specific information on fetal pharmacology will increase dramatically in the next decade due to interest in drugs administered in pregnancy and with the assistance of model-informed drug development.

A Physiologically Based Pharmacokinetic Model to Predict Potential Drug-Drug Interactions and Inform Dosing of Acumapimod, an Oral p38 MAPK Inhibitor

Acumapimod, an investigational oral p38 mitogen‐activated protein kinase inhibitor for treatment during severe acute exacerbations of chronic obstructive pulmonary disease, is metabolized primarily by cytochrome P450 3A4 (CYP3A4) and is a P‐glycoprotein (P‐gp) substrate. Concerns about drug–drug interactions (DDIs) have meant patients receiving drugs that inhibit CYP3A4 were ineligible for acumapimod trials. We report on how 2 acumapimod clinical DDI studies and a physiologically‐based pharmacokinetic (PBPK) model assessing how co‐administration of a weak (azithromycin) and strong (itraconazole) CYP3A4 inhibitor affected acumapimod systemic exposure, informed decision making and supported concomitant use of CYP3A4 and P‐gp inhibitors. Studies MBCT102 and MBCT103, respectively, demonstrated that co‐administration of azithromycin or itraconazole had no clinically meaningful impact on acumapimod pharmacokinetics. Findings were consistent with PBPK model results. Safety profiles were similar when acumapimod was co‐administered with azithromycin or itraconazole. These studies highlight the value of PBPK modeling in drug development, and its potential to inform DDI investigations.

Physiological parameter values for physiologically based pharmacokinetic models in food-producing animals. Part III: Sheep and goat

This report is the third in a series of studies that aimed to compile physiological parameters related to develop physiologically based pharmacokinetic (PBPK) models for drugs and environmental chemicals in food‐producing animals including swine and cattle (Part I), chickens and turkeys (Part II), and finally sheep and goats (the focus of this manuscript). Literature searches were conducted in multiple databases (PubMed, Google Scholar, ScienceDirect, and ProQuest), with data on relevant parameters including body weight, relative organ weight (% of body weight), cardiac output, relative organ blood flow (% of cardiac output), residual blood volume (% of organ weight), and hematocrit reviewed and statistically summarized. The mean and standard deviation of each parameter are presented in tables. Equations describing the growth curves of sheep and goats are presented in figures. When data are sufficient, parameter values are reported for different ages or production classes of sheep, including fetal sheep, lambs, and market‐age sheep (mature sheep). These data provide a reference database for developing standardized PBPK models to predict drug withdrawal intervals in sheep and goats, and also provide a basis for extrapolating PBPK models from major species such as cattle to minor species such as sheep and goats.

A Physiologically-Based Pharmacokinetic (PBPK) Model Network for the Prediction of CYP1A2 and CYP2C19 Drug-Drug-Gene Interactions with Fluvoxamine, Omeprazole, S-mephenytoin, Moclobemide, Tizanidine, Mexiletine, Ethinylestradiol, and Caffeine

Physiologically-based pharmacokinetic (PBPK) modeling is a well-recognized method for quantitatively predicting the effect of intrinsic/extrinsic factors on drug exposure. However, there are only few verified, freely accessible, modifiable, and comprehensive drug–drug interaction (DDI) PBPK models. We developed a qualified whole-body PBPK DDI network for cytochrome P450 (CYP) CYP2C19 and CYP1A2 interactions. Template PBPK models were developed for interactions between fluvoxamine, S-mephenytoin, moclobemide, omeprazole, mexiletine, tizanidine, and ethinylestradiol as the perpetrators or victims. Predicted concentration–time profiles accurately described a validation dataset, including data from patients with genetic polymorphisms, demonstrating that the models characterized the CYP2C19 and CYP1A2 network over the whole range of DDI studies investigated. The models are provided on GitHub (GitHub Inc., San Francisco, CA, USA), expanding the library of publicly available qualified whole-body PBPK models for DDI predictions, and they are thereby available to support potential recommendations for dose adaptations, support labeling, inform the design of clinical DDI trials, and potentially waive those.

Dose Optimization in Kidney Disease: Opportunities for PBPK Modeling and Simulation

Kidney disease affects pharmacokinetic (PK) profiles of not only renally cleared drugs but also nonrenally cleared drugs. The impact of kidney disease on drug disposition has not been fully elucidated, but describing the extent of such impact is essential for conducting dose optimization in kidney disease. Accurate evaluation of kidney function has been a clinical interest for dose optimization, and more scientists pay attention and conduct research for clarifying the role of drug transporters, metabolic enzymes, and their interplay in drug disposition as kidney disease progresses. Physiologically based pharmacokinetic (PBPK) modeling and simulation can provide valuable insights for dose optimization in kidney disease. It is a powerful tool to integrate discrete knowledge from preclinical and clinical research and mechanistically investigate system- and drug-dependent factors that may contribute to the changes in PK profiles. PBPK-based prediction of drug exposures may be used a priori to adjust dosing regimens and thereby minimize the likelihood of drug-related toxicity. With real-time clinical studies, parameter estimation may be performed with PBPK approaches that can facilitate identification of sources of interindividual variability. PBPK modeling may also facilitate biomarker research that aids dose optimization in kidney disease. U.S. Food and Drug Administration guidances related to conduction of PK studies in kidney impairment and PBPK documentation provide the foundation for facilitating model-based dose-finding research in kidney disease.

Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling to Predict the Impact of CYP2C9 Genetic Polymorphisms, Co-Medication and Formulation on the Pharmacokinetics and Pharmacodynamics of Flurbiprofen

Physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models can serve as a powerful framework for predicting the influence as well as the interaction of formulation, genetic polymorphism and co-medication on the pharmacokinetics and pharmacodynamics of drug substances. In this study, flurbiprofen, a potent non-steroid anti-inflammatory drug, was chosen as a model drug. Flurbiprofen has absolute bioavailability of ~95% and linear pharmacokinetics in the dose range of 50–300 mg. Its absorption is considered variable and complex, often associated with double peak phenomena, and its pharmacokinetics are characterized by high inter-subject variability, mainly due to its metabolism by the polymorphic CYP2C9 (fmCYP2C9 ≥ 0.71). In this study, by leveraging in vitro, in silico and in vivo data, an integrated PBPK/PD model with mechanistic absorption was developed and evaluated against clinical data from PK, PD, drug-drug and gene-drug interaction studies. The PBPK model successfully predicted (within 2-fold) 36 out of 38 observed concentration-time profiles of flurbiprofen as well as the CYP2C9 genetic effects after administration of different intravenous and oral dosage forms over a dose range of 40–300 mg in both Caucasian and Chinese healthy volunteers. All model predictions for C_max, AUC_inf and CL/F were within two-fold of their respective mean or geometric mean values, while 90% of the predictions of C_max, 81% of the predictions of AUC_inf and 74% of the predictions of Cl/F were within 1.25 fold. In addition, the drug-drug and drug-gene interactions were predicted within 1.5-fold of the observed interaction ratios (AUC, C_max ratios). The validated PBPK model was further expanded by linking it to an inhibitory E_max model describing the analgesic efficacy of flurbiprofen and applying it to explore the effect of formulation and genetic polymorphisms on the onset and duration of pain relief. This comprehensive PBPK/PD analysis, along with a detailed translational biopharmaceutic framework including appropriately designed biorelevant in vitro experiments and in vitro-in vivo extrapolation, provided mechanistic insight on the impact of formulation and genetic variations, two major determinants of the population variability, on the PK/PD of flurbiprofen. Clinically relevant specifications and potential dose adjustments were also proposed. Overall, the present work highlights the value of a translational PBPK/PD approach, tailored to target populations and genotypes, as an approach towards achieving personalized medicine.

Use of a physiologically based pharmacokinetic-pharmacodynamic model for initial dose prediction and escalation during a paediatric clinical trial

AIMS

To build and verify a physiologically based pharmacokinetic (PBPK) model for radiprodil in adults and link this to a pharmacodynamic (PD) receptor occupancy (RO) model derived from in vitro data. Adapt this model to the paediatric population and predict starting and escalating doses in infants based on RO. Use the model to guide individualized dosing in a clinical trial in 2- to 14-month-old children with infantile spasms.

METHODS

A PBPK model for radiprodil was developed to investigate the systemic exposure of the drug after oral administration in fasted and fed adults; this was then linked to RO via a PD model. The model was then expanded to include developmental physiology and ontogeny to predict escalating doses in infants that would result in a specific RO of 20, 40 and 60% based on average unbound concentration following a twice daily (b.i.d.) dosing regimen. Dose progression in the clinical trial was based on observed concentration-time data against PBPK predictions.

RESULTS

For paediatric predictions, the elimination of radiprodil, based on experimental evidence, had no ontogeny. Predicted b.i.d. doses ranged from 0.04 mg/kg for 20% RO, 0.1 mg/kg for 40% RO to 0.21 mg/kg for 60% RO. For all infants recruited in the study, observed concentration-time data following the 0.04 mg/kg and subsequent doses were within the PBPK model predicted 5^th and 95^th percentiles.

CONCLUSION

To our knowledge, this is the first time a PBPK model linked to RO has been used to guide dose selection and escalation in the live phase of a paediatric clinical trial.

Current Approaches and Techniques in Physiologically Based Pharmacokinetic (PBPK) Modelling of Nanomaterials

There have been efforts to develop physiologically based pharmacokinetic (PBPK) models for nanomaterials (NMs). Since NMs have quite different kinetic behaviors, the applicability of the approaches and techniques that are utilized in current PBPK models for NMs is warranted. Most PBPK models simulate a size-independent endocytosis from tissues or blood. In the lungs, dosimetry and the air-liquid interface (ALI) models have sometimes been used to estimate NM deposition and translocation into the circulatory system. In the gastrointestinal (GI) tract, kinetics data are needed for mechanistic understanding of NM behavior as well as their absorption through GI mucus and their subsequent hepatobiliary excretion into feces. Following absorption, permeability (P_t) and partition coefficients (PCs) are needed to simulate partitioning from the circulatory system into various organs. Furthermore, mechanistic modelling of organ- and species-specific NM corona formation is in its infancy. More recently, some PBPK models have included the mononuclear phagocyte system (MPS). Most notably, dissolution, a key elimination process for NMs, is only empirically added in some PBPK models. Nevertheless, despite the many challenges still present, there have been great advances in the development and application of PBPK models for hazard assessment and risk assessment of NMs.

Pharmacokinetics and population pharmacokinetics in pediatric oncology

Pharmacokinetic research has become increasingly important in pediatric oncology as it can have direct clinical implications and is a crucial component in individualized medicine. Population pharmacokinetics has become a popular method especially in children, due to the potential for sparse sampling, flexible sampling times, computing of heterogeneous data, and identification of variability sources. However, population pharmacokinetic reports can be complex and difficult to interpret. The aim of this article is to provide a basic explanation of population pharmacokinetics, using clinical examples from the field of pediatric oncology, to facilitate the translation of pharmacokinetic research into the daily clinic.

Predicting Clinical Effects of CYP3A4 Modulators on Abemaciclib and Active Metabolites Exposure Using Physiologically Based Pharmacokinetic Modeling

Abemaciclib, a selective inhibitor of cyclin‐dependent kinases 4 and 6, is metabolized mainly by cytochrome P450 (CYP)3A4. Clinical studies were performed to assess the impact of strong inhibitor (clarithromycin) and inducer (rifampin) on the exposure of abemaciclib and active metabolites. A physiologically based pharmacokinetic (PBPK) model incorporating the metabolites was developed to predict the effect of other strong and moderate CYP3A4 inhibitors and inducers. Clarithromycin increased the area under the plasma concentration‐time curve (AUC) of abemaciclib and potency‐adjusted unbound active species 3.4‐fold and 2.5‐fold, respectively. Rifampin decreased corresponding exposures 95% and 77%, respectively. These changes influenced the fraction metabolized via CYP3A4 in the model. An absolute bioavailability study informed the hepatic and gastric availability. In vitro data and a human radiolabel study determined the fraction and rate of formation of the active metabolites as well as absorption‐related parameters. The predicted AUC ratios of potency‐adjusted unbound active species with rifampin and clarithromycin were within 0.7‐ and 1.25‐fold of those observed. The PBPK model predicted 3.78‐ and 7.15‐fold increases in the AUC of the potency‐adjusted unbound active species with strong CYP3A4 inhibitors itraconazole and ketoconazole, respectively; and 1.62‐ and 2.37‐fold increases with the concomitant use of moderate CYP3A4 inhibitors verapamil and diltiazem, respectively. The model predicted modafinil, bosentan, and efavirenz would decrease the AUC of the potency‐adjusted unbound active species by 29%, 42%, and 52%, respectively. The current PBPK model, which considers changes in unbound potency‐adjusted active species, can be used to inform dosing recommendations when abemaciclib is coadministered with CYP3A4 perpetrators.

Recent advances in physiologically based pharmacokinetic and pharmacodynamic models for anticancer nanomedicines

Nanoparticles (NPs) have distinct pharmacokinetic (PK) properties and can potentially improve the absorption, distribution, metabolism, and elimination (ADME) of small-molecule drugs loaded therein. Owing to the unwanted toxicities of anticancer agents in healthy organs and tissues, their precise delivery to the tumor is an essential requirement. There have been numerous advancements in the development of nanomedicines for cancer therapy. Physiologically based PK (PBPK) models serve as excellent tools for describing and predicting the ADME properties and the efficacy and toxicity of drugs, in combination with pharmacodynamic (PD) models. The recent preliminary application of these modeling approaches to NPs demonstrated their potential benefits in research and development processes relevant to the ADME and pharmacodynamics of NPs and nanomedicines. Here, we comprehensively review the pharmacokinetics of NPs, the developed PBPK models for anticancer NPs, and the developed PD model for anticancer agents.