Population Physiologically Based Pharmacokinetic Research Articles

Prediction of drug-drug interactions between roflumilast and CYP3A4/1A2 perpetrators using a physiologically-based pharmacokinetic (PBPK) approach.

This study aimed to develop a physiologically-based pharmacokinetic (PBPK) model to predict changes in the pharmacokinetics (PK) and pharmacodynamics (PD, PDE4 inhibition) of roflumilast (ROF) and ROF N-oxide when co-administered with eight CYP3A4/1A2 perpetrators. The population PBPK model of ROF and ROF N-oxide has been successfully developed and validated based on the four clinical PK studies and five clinical drug-drug interactions (DDIs) studies. In PK simulations, every ratio of prediction to observation for PK parameters fell within the range 0.7 to 1.5. In DDI simulations, except for tow peak concentration ratios (Cmax) of ROF with rifampicin (prediction: 0.63 vs. observation: 0.19) and with cimetidine (prediction: 1.07 vs. observation: 1.85), the remaining predicted ratios closely matched the observed ratios. Additionally, the PBPK model suggested that co-administration with the three perpetrators (cimetidine, enoxacin, and fluconazole) may use with caution, with CYP3A4 strong inhibitor (ketoconazole and itraconazole) or with dual CYP3A41A2 inhibitor (fluvoxamine) may reduce to half-dosage or use with caution, while co-administration with CYP3A4 strong or moderate inducer (rifampicin, efavirenz) should avoid. Overall, the present PBPK model can provide recommendations for adjusting dosing regimens in the presence of DDIs.

Exploring inter-ethnic and inter-patient variability and optimal dosing of osimertinib: a physiologically based pharmacokinetic modeling approach.

Purpose: This study aimed to develop and validate a physiologically based pharmacokinetic (PBPK) model for osimertinib (OSI) to predict plasma trough concentration (Ctrough) and pulmonary EGFRm+ (T790M and L858R mutants) inhibition in Caucasian, Japanese, and Chinese populations. The PBPK model was also utilized to investigate inter-ethnic and inter-patient differences in OSI pharmacokinetics (PK) and determine optimal dosing regimens. Methods: Population PBPK models of OSI for healthy and disease populations were developed using physicochemical and biochemical properties of OSI and physiological parameters of different groups. And then the PBPK models were validated using the multiple clinical PK and drug-drug interaction (DDI) study data. Results: The model demonstrated good consistency with the observed data, with most of prediction-to-observation ratios of 0.8-1.25 for AUC, Cmax, and Ctrough. The PBPK model revealed that plasma exposure of OSI was approximately 2-fold higher in patients compared to healthy individuals, and higher exposure observed in Caucasians compared to other ethnic groups. This was primarily attributed to a lower CL/F of OSI in patients and Caucasian. The PBPK model displayed that key factors influencing PK and EGFRm+ inhibition differences included genetic polymorphism of CYP3A4, CYP1A2 expression, plasma free concentration (fup), albumin level, and auto-inhibition/induction on CYP3A4. Inter-patient PK variability was most influenced by CYP3A4 variants, fup, and albumin level. The PBPK simulations indicated that the optimal dosing regimen for patients across the three populations of European, Japanese, and Chinese ancestry was OSI 80mg once daily (OD) to achieve the desired range of plasma Ctrough (328-677nmol/L), as well as 80mg and 160mg OD for desirable pulmonary EGFRm+ inhibition (>80%). Conclusion: In conclusion, this study's PBPK simulations highlighted potential ethnic and inter-patient variability in OSI PK and EGFRm+ inhibition between Caucasian, Japanese, and Chinese populations, while also providing insights into optimal dosing regimens of OSI.

Assessing Dose-Exposure–Response Relationships of Miltefosine in Adults and Children using Physiologically-Based Pharmacokinetic Modeling Approach

ObjectivesMiltefosine is the first and only oral medication to be successfully utilized as an antileishmanial agent. However, the drug is associated with differences in exposure patterns and cure rates among different population groups e.g. ethnicity and age (i.e., children v adults) in clinical trials. In this work, mechanistic population physiologically-based pharmacokinetic (PBPK) models have been developed to study the dose-exposure–response relationship of miltefosine in in silico clinical trials and evaluate the differences in population groups, particularly children and adults.MethodsThe Simcyp population pharmacokinetics platform was employed to predict miltefosine exposure in plasma and peripheral blood mononuclear cells (PBMCs) in a virtual population under different dosing regimens. The cure rate of a simulation was based on the percentage of number of the individual virtual subjects with AUCd0-28 > 535 µg⋅day/mL in the virtual population.ResultsIt is shown that both adult and paediatric PBPK models of miltefosine can be developed to predict the PK data of the clinical trials accurately. There was no significant difference in the predicted dose-exposure–response of the miltefosine treatment for different simulated ethnicities under the same dose regime and the dose-selection strategies determined the clinical outcome of the miltefosine treatment. A lower cure rate of the miltefosine treatment in paediatrics was predicted because a lower exposure of miltefosine was simulated in virtual paediatric in comparison with adult virtual populations when they received the same dose of the treatment.ConclusionsThe mechanistic PBPK model suggested that the higher fraction of unbound miltefosine in plasma was responsible for a higher probability of failure in paediatrics because of the difference in the distribution of plasma proteins between adults and paediatrics. The developed PBPK models could be used to determine an optimal miltefosine dose regime in future clinical trials.Graphical

Physiologically-based pharmacokinetic modeling for optimal dosage prediction of olaparib when co-administered with CYP3A4 modulators and in patients with hepatic/renal impairment

This study aimed to develop a physiologically-based pharmacokinetic (PBPK) model to predict the maximum plasma concentration (Cmax) and trough concentration (Ctrough) at steady-state of olaparib (OLA) in Caucasian, Japanese and Chinese. Furthermore, the PBPK model was combined with mean and 95% confidence interval to predict optimal dosing regimens of OLA when co-administered with CYP3A4 modulators and administered to patients with hepatic/renal impairment. The dosing regimens were determined based on safety and efficacy PK threshold Cmax (< 12,500 ng/mL) and Ctrough (772–2500 ng/mL). The population PBPK model for OLA was successfully developed and validated, demonstrating good consistency with clinically observed data. The ratios of predicted to observed values for Cmax and Ctrough fell within the range of 0.5 to 2.0. When OLA was co-administered with a strong or moderate CYP3A4 inhibitor, the recommended dosing regimens should be reduced to 100 mg BID and 150 mg BID, respectively. Additionally, the PBPK model also suggested that OLA could be not recommended with a strong or moderate CYP3A4 inducer. For patients with moderate hepatic and renal impairment, the dosing regimens of OLA were recommended to be reduced to 200 mg BID and 150 mg BID, respectively. In cases of severe hepatic and renal impairment, the PBPK model suggested a dosing regimen of 100 mg BID for OLA. Overall, this present PBPK model can determine the optimal dosing regimens for various clinical scenarios involving OLA.

Lamivudine and Emtricitabine Dosing Proposal for Children with HIV and Chronic Kidney Disease, Supported by Physiologically Based Pharmacokinetic Modelling.

Dose recommendations for lamivudine or emtricitabine in children with HIV and chronic kidney disease (CKD) are absent or not supported by clinical data. Physiologically based pharmacokinetic (PBPK) models have the potential to facilitate dose selection for these drugs in this population. Existing lamivudine and emtricitabine compound models in Simcyp® (v21) were verified in adult populations with and without CKD and in non-CKD paediatric populations. We developed paediatric CKD population models reflecting subjects with a reduced glomerular filtration and tubular secretion, based on extrapolation from adult CKD population models. These models were verified using ganciclovir as a surrogate compound. Then, lamivudine and emtricitabine dosing strategies were simulated in virtual paediatric CKD populations. The compound and paediatric CKD population models were verified successfully (prediction error within 0.5- to 2-fold). The mean AUC ratios in children (GFR-adjusted dose in CKD population/standard dose in population with normal kidney function) were 1.15 and 1.23 for lamivudine, and 1.20 and 1.30 for emtricitabine, with grade-3- and -4-stage CKD, respectively. With the developed paediatric CKD population PBPK models, GFR-adjusted lamivudine and emtricitabine dosages in children with CKD resulted in adequate drug exposure, supporting paediatric GFR-adjusted dosing. Clinical studies are needed to confirm these findings.

Application of population physiologically based pharmacokinetic modelling to optimize target expression and clearance mechanisms of therapeutic monoclonal antibodies.

To use population physiologically based pharmacokinetic (PopPBPK) modelling to optimize target expression, kinetics and clearance of HER1/2 directed therapeutic monoclonal antibodies (mAbs). Thus, to propose a general workflow of PopPBPK modelling and its application in clinical pharmacology. Full PBPK model of pertuzumab (PTZ) was developed in patient population using Simcyp V21R1 incorporating mechanistic targeted-mediated drug disposition process by fitting known clinical PK and sparse receptor proteomics data to optimize target expression and kinetics of HER2 receptor. Trastuzumab (TTZ) PBPK modelling was used to validate the optimized HER2 target. Additionally, the simulator was also used to develop a full PBPK model for the HER1-directed mAb cetuximab (CTX) to assess the underlying targeted-mediated drug disposition-independent elimination mechanisms. HER2 final parameterisation coming from the PBPK modelling of PTZ was successfully cross validated through PBPK modelling of TTZ with average fold error (AFE), absolute AFE and percent prediction error values for area under the concentration-time curve (AUC) and maximum plasma concentration (Cmax ) of 1.13, 1.16 and 16, and 1.01, 1.07 and 7, respectively. CTX PBPK model performance was validated after the incorporation of an additional systemic clearance of 0.033 L/h as AFE and absolute AFE showed an acceptable predictive power of AUC and Cmax with percent prediction error of 13% for AUC and 10% for Cmax . Optimisation of both system and drug related parameters were performed through PBPK modelling to improve model performance of therapeutic mAbs (PTZ, TTZ and CTX). General workflow was proposed to develop and apply PopPBPK to support clinical development of mAbs targeting same receptor.

Prediction for Plasma Trough Concentration and Optimal Dosing of Imatinib under Multiple Clinical Situations Using Physiologically Based Pharmacokinetic Modeling.

(1) Purpose: This study aimed to develop a physiologically based pharmacokinetic (PBPK) model to predict the trough concentration (C trough) of imatinib (IMA) at steady state in patients and to explore the role of free concentration (f up), α1-acid glycoprotein (AGP) level, and organic cation transporter 1 (OCT1) activity/expression in clinical efficacy. (2) Methods: The population PBPK model was built using physicochemical and biochemical properties, metabolizing and transporting kinetics, tissue distribution, and human physiological parameters. (3) Results: The PBPK model successfully predicted the C trough of IMA administered alone in chronic phase (CP) and accelerated phase (AP) patients, the C trough of IMA co-administered with six modulators, and C trough in CP patients with hepatic impairment. Most of the ratios between predicted and observed data are within 0.70-1.30. Additionally, the recommendations for dosing adjustments for IMA have been given under multiple clinical uses. The sensitivity analysis showed that exploring the f up and AGP level had a significant influence on the plasma C trough of IMA. Meanwhile, the simulations also revealed that OCT1 activity and expression had a significant impact on the intracellular C trough of IMA. (4) Conclusion: The current PBPK model can accurately predict the IMA C trough and provide appropriate dosing adjustment recommendations in a variety of clinical situations.

Physiologically-Based Pharmacokinetic Modeling of Tenofovir Disoproxil Fumarate in Pregnant Women.

Physiological changes during pregnancy can affect antiretroviral drug processes and further influence drug efficacy and safety. Physiologically-based pharmacokinetic (PBPK) modeling offers a unique modality to predict PK in pregnant women. The objective of this study was to establish a PBPK modeling of tenofovir disoproxil fumarate (TDF) in pregnant women, to provide a reference for the clinical use of TDF. A full PBPK modeling of tenofovir (TFV) and TDF following i.v. and p.o. administration was developed using the simulation software PK-Sim®. The modeling was then extrapolated to pregnant women based on pregnancy- related physiological parameters in Mobi® Simulator. The mean fold error (MFE) and geometric mean fold error (GMFE) methods were used to compare the differences between predicted and observed values of PK parameters (Cmax, tmax, AUC0-∞) to evaluate the accuracy of PBPK modeling. The developed PBPK modeling successfully predicted the TDF disposition in the non-pregnant population, wherein the MFE average and GMFE of all predicted PK parameters were within a 1.5-fold error range, and more than 96.30% of the predicted drug concentration values were within a 2-fold error range of the measured values. After the extrapolation of these models to the third trimester of pregnancy, the scaling anatomy/physiology and hepatic intrinsic clearance made the pregnant population PBPK modeling meet the standard requirement of 0.5 < MFE and GMFE value < 2. It was more appropriate to simulate the in vivo process of low-dose TDF in pregnant women. The non-pregnant population PBPK modeling of TDF established in our study can be extrapolated to pregnant women. Our study provides a reference for realizing clinical personalized medication for pregnant women.

Prediction for optimal dosage of pazopanib under various clinical situations using physiologically based pharmacokinetic modeling.

This study aimed to apply a physiologically based pharmacokinetic (PBPK) model to predict optimal dosing regimens of pazopanib (PAZ) for safe and effective administration when co-administered with CYP3A4 inhibitors, acid-reducing agents, food, and administered in patients with hepatic impairment. Here, we have successfully developed the population PBPK model and the predicted PK variables by this model matched well with the clinically observed data. Most ratios of prediction to observation were between 0.5 and 2.0. Suitable dosage modifications of PAZ have been identified using the PBPK simulations in various situations, i.e., 200 mg once daily (OD) or 100 mg twice daily (BID) when co-administered with the two CYP3A4 inhibitors, 200 mg BID when simultaneously administered with food or 800 mg OD when avoiding food uptake simultaneously. Additionally, the PBPK model also suggested that dosing does not need to be adjusted when co-administered with esomeprazole and administration in patients with wild hepatic impairment. Furthermore, the PBPK model also suggested that PAZ is not recommended to be administered in patients with severe hepatic impairment. In summary, the present PBPK model can determine the optimal dosing adjustment recommendations in multiple clinical uses, which cannot be achieved by only focusing on AUC linear change of PK.

Physiologically based pharmacokinetic (PBPK) modeling of piroxicam with regard to CYP2C9 genetic polymorphism.

Piroxicam is a non-steroidal anti-inflammatory drug used to alleviate symptoms of osteoarthritis and rheumatoid arthritis. CYP2C9 genetic polymorphism significantly influences the pharmacokinetics of piroxicam. The objective of this study was to develop and validate the piroxicam physiologically based pharmacokinetic (PBPK) model related to CYP2C9 genetic polymorphism. PK-Sim® version 10.0 was used for the PBPK modeling. The PBPK model was evaluated by predicted and observed plasma concentration-time profiles, fold errors of predicted to observed pharmacokinetic parameters, and a goodness-of-fit plot. The turnover number (kcat) of CYP2C9 was adjusted to capture the pharmacokinetics of piroxicam in different CYP2C9 genotypes. The population PBPK model overall accurately described and predicted the plasma concentration-time profiles in different CYP2C9 genotypes. In our simulations, predicted AUCinf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.83-, 2.07-, and 6.43-fold higher than CYP2C9*1/*1 genotype, respectively. All fold error values for AUC, Cmax, and t1/2 were included in the acceptance criterion with the ranges of 0.57-1.59, 0.63-1.39, and 0.65-1.51, respectively. The range of fold error values for predicted versus observed plasma concentrations was 0.11-3.13. 93.9% of fold error values were within the two-fold range. Average fold error, absolute average fold error, and root mean square error were 0.93, 1.27, and 0.72, respectively. Our model accurately captured the pharmacokinetic alterations of piroxicam according to CYP2C9 genetic polymorphism.

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.

Current use and future potential of (physiologically based) pharmacokinetic modelling of radiopharmaceuticals: a review.

Rationale: Physiologically based pharmacokinetic (PBPK) and population pharmacokinetic (PK) modelling approaches are widely accepted in non-radiopharmaceutical drug development and research, while there is no major role for these approaches in radiopharmaceutical development yet. In this review, a literature search was performed to specify different research purposes and questions that have previously been answered using both PBPK and population PK modelling for radiopharmaceuticals. Methods: The literature search was performed using the databases PubMed and Embase. Wide search terms included radiopharmaceutical, tracer, radioactivity, physiologically based pharmacokinetic model, PBPK, population pharmacokinetic model and nonlinear mixed-effects model. Results: Eight articles and twenty articles were included for this review based on this literature search for population PK modelling and PBPK modelling, respectively. Included population PK analyses showed to have an added value to develop predictive models for a population and to describe individual variability sources. Main purposes of PBPK models appeared related to optimizing treatment (planning), or more specifically: to find the optimal combination of peptide amount and radioactivity, to optimize treatment planning by reducing the number of measurements, to individualize treatment, to get insights in differences between pre-therapeutic and therapeutic scans or to understand inter-patient differences. Other main research subjects were regarding radiopharmaceutical comparisons, selecting ligands based on their peptide characteristics and gaining a better understanding of drug-drug interactions. Conclusions: The use of PK modelling approaches in radiopharmaceutical research remains scarce, but can be expanded to obtain a better understanding of PK and whole-body distribution of radiopharmaceuticals in general. PK modelling of radiopharmaceuticals has great potential for the nearby future and could contribute to the evolving research of radiopharmaceuticals.

Diversity and Inclusion in Drug Development: Rethinking Intrinsic and Extrinsic Factors with Patient Centricity.

Diversity and Inclusion in Drug Development: Rethinking Intrinsic and Extrinsic Factors with Patient Centricity.

Quantitative Characterization of Population-Wide Tissue- and Metabolite-Specific Variability in Perchloroethylene Toxicokinetics in Male Mice.

Quantification of interindividual variability is a continuing challenge in risk assessment, particularly for compounds with complex metabolism and multi-organ toxicity. Toxicokinetic variability for perchloroethylene (perc) was previously characterized across 3 mouse strains and in 1 mouse strain with various degrees of liver steatosis. To further characterize the role of genetic variability in toxicokinetics of perc, we applied Bayesian population physiologically based pharmacokinetic (PBPK) modeling to the data on perc and metabolites in blood/plasma and tissues of male mice from 45 inbred strains from the Collaborative Cross (CC) mouse population. After identifying the most influential PBPK parameters based on global sensitivity analysis, we fit the model with a hierarchical Bayesian population analysis using Markov chain Monte Carlo simulation. We found that the data from 3 commonly used strains were not representative of the full range of variability in perc and metabolite blood/plasma and tissue concentrations across the CC population. Using interstrain variability as a surrogate for human interindividual variability, we calculated dose-dependent, chemical-, and tissue-specific toxicokinetic variability factors (TKVFs) as candidate science-based replacements for the default uncertainty factor for human toxicokinetic variability of 100.5. We found that toxicokinetic variability factors for glutathione conjugation metabolites of perc showed the greatest variability, often exceeding the default, whereas those for oxidative metabolites and perc itself were generally less than the default. Overall, we demonstrate how a combination of a population-based mouse model such as the CC with Bayesian population PBPK modeling can reduce uncertainty in human toxicokinetic variability and increase accuracy and precision in quantitative risk assessment.

Development of a Virtual Chinese Pediatric Population Physiological Model Targeting Specific Metabolism and Kidney Elimination Pathways.

Background: Physiologically based pharmacokinetic (PBPK) modeling and simulating may be a powerful tool in predicting drug behaviors in specific populations. It is a mathematical model that relates the pharmacokinetic (PK) profile of a compound with human anatomical characteristics, physiological characteristics, and biochemical parameters. Predictions using PBPK models offer a promising way to guide drug development and can be used to optimize clinical dosing regimens. However, PK data of new drugs in the pediatric population are too limited to guide clinical therapy, which may lead to frequent adverse events or insufficient efficacy for pediatric patients, particularly in neonates and infants. Objective: The objective of this study was to establish a virtual Chinese pediatric population based on the physiological parameters of Chinese children that could be utilized in PBPK models. Methods: A Chinese pediatric PBPK model was developed in Simcyp Simulator by collecting published Chinese pediatric physiological and anthropometric data to use as system parameters. This pediatric population model was then evaluated in the Chinese pediatric population by predicting the pharmacokinetic characteristics of four probe drugs: theophylline (major CYP1A2 substrate), fentanyl (major CYP3A4 substrate), vancomycin, and ceftazidime (renal-eliminated). Results: The predicted maximum concentration (Cmax), area under the curve of concentration-time (AUC), and clearance (CL) for theophylline (CYP1A2 metabolism pathway) and fentanyl (CYP3A4 metabolism pathway) were within two folds of the observed data. For drugs mainly eliminated by renal clearance (vancomycin and ceftazidime) in the Chinese pediatric population, the ratio of prediction to observation for major PK parameters was within a 2-fold error range. Conclusion: The model is a supplement to the previous Chinese population PBPK model. We anticipate the model to be a better representative of the pediatric Chinese population for drugs PK, offering greater clinical precision for medication given to the pediatric population, ultimately advancing clinical development of pediatric drugs. We can refine this model further by collecting more physiological parameters of Chinese children.

Bayesian population physiologically-based pharmacokinetic model for robustness evaluation of withdrawal time in tilapia aquaculture administrated to florfenicol

The antimicrobial residues of aquacultural production is a growing public concern, leading to reexamine the method for establishing robust withdrawal time and ensuring food safety. Our study aims to develop the optimizing population physiologically-based pharmacokinetic (PBPK) model for assessing florfenicol residues in the tilapia tissues, and for evaluating the robustness of the withdrawal time (WT). Fitting with published pharmacokinetic profiles that experimented under temperatures of 22 and 28 °C, a PBPK model was constructed by applying with the Bayesian Markov chain Monte Carol (MCMC) algorithm to estimate WTs under different physiological, environmental and dosing scenarios. Results show that the MCMC algorithm improves the estimates of uncertainty and variability of PBPK-related parameters, and optimizes the simulation of the PBPK model. It is noteworthy that posterior sets generated from temperature-associated datasets to be respectively used for simulating residues under corresponding temperature conditions. Simulating the residues under regulated regimen and overdosing scenarios for Taiwan, the estimated WTs were 12–16 days at 22 °C and 9–12 days at 28 °C, while for the USA, the estimated WTs were 14–18 and 11–14 days, respectively. Comparison with the regulated WT of 15 days, results indicate that the current WT has well robustness and resilience in the environment of higher temperatures. The optimal Bayesian population PBPK model provides effective analysis for determining WTs under scenario-specific conditions. It is a new insight into the increasing body of literature on developing the Bayesian-PBPK model and has practical implications for improving the regulation of food safety.

Use of normalized prediction distribution errors for assessing population physiologically-based pharmacokinetic model adequacy.

Currently employed methods for qualifying population physiologically-based pharmacokinetic (Pop-PBPK) model predictions of continuous outcomes (e.g., concentration-time data) fail to account for within-subject correlations and the presence of residual error. In this study, we propose a new method for evaluating Pop-PBPK model predictions that account for such features. The approach focuses on deriving Pop-PBPK-specific normalized prediction distribution errors (NPDE), a metric that is commonly used for population pharmacokinetic model validation. We describe specific methodological steps for computing NPDE for Pop-PBPK models and define three measures for evaluating model performance: mean of NPDE, goodness-of-fit plots, and the magnitude of residual error. Utility of the proposed evaluation approach was demonstrated using two simulation-based study designs (positive and negative control studies) as well as pharmacokinetic data from a real-world clinical trial. For the positive-control simulation study, where observations and model simulations were generated under the same Pop-PBPK model, the NPDE-based approach denoted a congruency between model predictions and observed data (mean of NPDE = - 0.01). In contrast, for the negative-control simulation study, where model simulations and observed data were generated under different Pop-PBPK models, the NPDE-based method asserted that model simulations and observed data were incongruent (mean of NPDE = - 0.29). When employed to evaluate a previously developed clindamycin PBPK model against prospectively collected plasma concentration data from 29 children, the NPDE-based method qualified the model predictions as successful (mean of NPDE = 0). However, when pediatric subpopulations (e.g., infants) were evaluated, the approach revealed potential biases that should be explored.

Physiological parameter values for physiologically based pharmacokinetic models in food-producing animals. Part I: Cattle and swine.

Physiologically based pharmacokinetic (PBPK) models for chemicals in food animals are a useful tool in estimating chemical tissue residues and withdrawal intervals. Physiological parameters such as organ weights and blood flows are an important component of a PBPK model. The objective of this study was to compile PBPK‐related physiological parameter data in food animals, including cattle and swine. Comprehensive literature searches were performed in PubMed, Google Scholar, ScienceDirect, and ProQuest. Relevant literature was reviewed and tables of relevant parameters such as relative organ weights (% of body weight) and relative blood flows (% of cardiac output) were compiled for different production classes of cattle and swine. The mean and standard deviation of each parameter were calculated to characterize their variability and uncertainty and to allow investigators to conduct population PBPK analysis via Monte Carlo simulations. Regression equations using weight or age were created for parameters having sufficient data. These compiled data provide a comprehensive physiological parameter database for developing PBPK models of chemicals in cattle and swine to support animal‐derived food safety assessment. This work also provides a basis to compile data in other food animal species, including goats, sheep, chickens, and turkeys.

Integration of Food Animal Residue Avoidance Databank (FARAD) empirical methods for drug withdrawal interval determination with a mechanistic population-based interactive physiologically based pharmacokinetic (iPBPK) modeling platform: example for flunixin meglumine administration.

Violative chemical residues in animal-derived food products affect food safety globally and have impact on the trade of international agricultural products. The Food Animal Residue Avoidance Databank program has been developing scientific tools to provide appropriate withdrawal interval (WDI) estimations after extralabel drug use in food animals for the past three decades. One of the tools is physiologically based pharmacokinetic (PBPK) modeling, which is a mechanistic-based approach that can be used to predict tissue residues and WDIs. However, PBPK models are complicated and difficult to use by non-modelers. Therefore, a user-friendly PBPK modeling framework is needed to move this field forward. Flunixin was one of the top five violative drug residues identified in the United States from 2010 to 2016. The objective of this study was to establish a web-based user-friendly framework for the development of new PBPK models for drugs administered to food animals. Specifically, a new PBPK model for both cattle and swine after administration of flunixin meglumine was developed. Population analysis using Monte Carlo simulations was incorporated into the model to predict WDIs following extralabel administration of flunixin meglumine. The population PBPK model was converted to a web-based interactive PBPK (iPBPK) framework to facilitate its application. This iPBPK framework serves as a proof-of-concept for further improvements in the future and it can be applied to develop new models for other drugs in other food animal species, thereby facilitating the application of PBPK modeling in WDI estimation and food safety assessment.

Pitfalls of using numerical predictive checks for population physiologically-based pharmacokinetic model evaluation.

Comparisons between observed data and model simulations represent a critical component for establishing confidence in population physiologically-based pharmacokinetic (Pop-PBPK) models. Numerical predictive checks (NPC) that assess the proportion of observed data that correspond to Pop-PBPK model prediction intervals (PIs) are frequently used to qualify such models. We evaluated the effects of three components on the performance of NPC for qualifying Pop-PBPK model concentration-time predictions: (1) correlations (multiple samples per subject), (2) residual error, and (3) discrepancies in the distribution of demographics between observed and virtual subjects. Using a simulation-based study design, we artificially created observed pharmacokinetic (PK) datasets and compared them to model simulations generated under the same Pop-PBPK model. Observed datasets containing uncorrelated and correlated observations (± residual error) were formulated using different random-sampling techniques. In addition, we created observed datasets where the distribution of subject body weights differed from that of the virtual population used to generate model simulations. NPC for each observed dataset were computed based on the Pop-PBPK model's 90% PI. NPC were associated with inflated type-I-error rates (> 0.10) for observed datasets that contained correlated observations, residual error, or both. Additionally, the performance of NPC were sensitive to the demographic distribution of observed subjects. Acceptable use of NPC was only demonstrated for the idealistic case where observed data were uncorrelated, free of residual error, and the demographic distribution of virtual subjects matched that of observed subjects. Considering the restricted applicability of NPC for Pop-PBPK model evaluation, their use in this context should be interpreted with caution.