Applied Concepts in PBPK Modeling: How to Build a PBPK/PD Model

Assessment of the CYP3A4 Induction Potential by Carbamazepine: Insights from Two Clinical DDI Studies and PBPK Modeling

In the past, rifampicin was well-established as strong index CYP3A inducer in clinical drug-drug interaction (DDI) studies. However, due to identified potentially genotoxic nitrosamine impurities, it should not any longer be used in healthy volunteer studies. Available clinical data suggest carbamazepine as an alternative to rifampicin as strong index CYP3A4 inducer in clinical DDI studies. Further, physiologically-based pharmacokinetic (PBPK) modeling is a tool with increasing importance to support the DDI risk assessment of drugs during drug development. CYP3A4 induction properties and the safety profile of carbamazepine were investigated in two open-label, fixed sequence, crossover clinical pharmacology studies in healthy volunteers using midazolam as a sensitive index CYP3A4 substrate. Carbamazepine was up-titrated from 100 mg twice daily (b.i.d.) to 200 mg b.i.d., and to a final dose of 300 mg b.i.d. for 10 consecutive days. Mean area under plasma concentration-time curve from zero to infinity (AUC(0-∞)) of midazolam consistently decreased by 71.8% (ratio: 0.282, 90% confidence interval (CI): 0.235-0.340) and 67.7% (ratio: 0.323, 90% CI: 0.256-0.407) in study 1 and study 2, respectively. The effect was adequately described by an internally developed PBPK model for carbamazepine which has been made freely available to the scientific community. Further, carbamazepine was safe and well-tolerated in the investigated dosing regimen in healthy participants. The results demonstrated that the presented design is appropriate for the use of carbamazepine as alternative inducer to rifampicin in DDI studies acknowledging its CYP3A4 inductive potency and safety profile.

Bioequivalence risk assessment of oral formulations containing racemic ibuprofen through a chiral physiologically based pharmacokinetic model of ibuprofen enantiomers

The characterization of the time course of ibuprofen enantiomers can be useful in the selection of the most sensitive analyte in bioequivalence studies. Physiologically based pharmacokinetic (PBPK) modelling and simulation represents the most efficient methodology to virtually assess bioequivalence outcomes. In this work, we aim to develop and verify a PBPK model for ibuprofen enantiomers administered as a racemic mixture with different immediate release dosage forms to anticipate bioequivalence outcomes based on different particle size distributions. A PBPK model incorporating stereoselectivity and non-linearity in plasma protein binding and metabolism as well as R-to-S unidirectional inversion has been developed in Simcyp®. A dataset composed of 11 Phase I clinical trials with 54 scenarios (27 per enantiomer) and 14,452 observations (7129 for R-ibuprofen and 7323 for S-ibuprofen) was used. Prediction errors for AUC0-t and Cmax for both enantiomers fell within the 0.8-1.25 range in 50/54 (93 %) and 42/54 (78 %) of scenarios, respectively. Outstanding model performance, with 10/10 (100 %) of Cmax and 9/10 (90 %) of AUC0-t within the 0.9-1.1 range, was demonstrated for oral suspensions, which strongly supported its use for bioequivalence risk assessment. The deterministic bioequivalence risk assessment has revealed R-ibuprofen as the most sensitive analyte to detect differences in particle size distribution for oral suspensions containing 400 mg of racemic ibuprofen, suggesting that achiral bioanalytical methods would increase type II error and declare non-bioequivalence for formulations that are bioequivalent for the eutomer.

The Effect of Inflammatory Bowel Disease and Irritable Bowel Syndrome on Pravastatin Oral Bioavailability: In vivo and in silico evaluation using bottom-up wbPBPK modeling

The common disorders irritable bowel syndrome (IBS) and inflammatory bowel disease (IBD) can modify the drugs' pharmacokinetics via their induced pathophysiological changes. This work aimed to investigate the impact of these two diseases on pravastatin oral bioavailability. Rat models for IBS and IBD were used to experimentally test the effects of IBS and IBD on pravastatin pharmacokinetics. Then, the observations made in rats were extrapolated to humans using a mechanistic whole-body physiologically-based pharmacokinetic (wbPBPK) model. The rat in vivo studies done herein showed that IBS and IBD decreased serum albumin (> 11% for both), decreased PRV binding in plasma, and increased pravastatin absolute oral bioavailability (0.17 and 0.53 compared to 0.01) which increased plasma, muscle, and liver exposure. However, the wbPBPK model predicted muscle concentration was much lower than the pravastatin toxicity thresholds for myotoxicity and rhabdomyolysis. Overall, IBS and IBD can significantly increase pravastatin oral bioavailability which can be due to a combination of increased pravastatin intestinal permeability and decreased pravastatin gastric degradation resulting in higher exposure. This is the first study in the literature investigating the effects of IBS and IBD on pravastatin pharmacokinetics. The high interpatient variability in pravastatin concentrations as induced by IBD and IBS can be reduced by oral administration of pravastatin using enteric-coated tablets. Such disease (IBS and IBD)-drug interaction can have more drastic consequences for narrow therapeutic index drugs prone to gastric degradation, especially for drugs with low intestinal permeability.

Use of quantitative in vitro to in vivo extrapolation (QIVIVE) for the assessment of non-combustible next-generation product aerosols

With the use of in vitro new approach methodologies (NAMs) for the assessment of non-combustible next-generation nicotine delivery products, new extrapolation methods will also be required to interpret and contextualize the physiological relevance of these results. Quantitative in vitro to in vivo extrapolation (QIVIVE) can translate in vitro concentrations into in-life exposures with physiologically-based pharmacokinetic (PBPK) modelling and provide estimates of the likelihood of harmful effects from expected exposures. A major challenge for evaluating inhalation toxicology is an accurate assessment of the delivered dose to the surface of the cells and the internalized dose. To estimate this, we ran the multiple-path particle dosimetry (MPPD) model to characterize particle deposition in the respiratory tract and developed a PBPK model for nicotine that was validated with human clinical trial data for cigarettes. Finally, we estimated a Human Equivalent Concentration (HEC) and predicted plasma concentrations based on the minimum effective concentration (MEC) derived after acute exposure of BEAS-2B cells to cigarette smoke (1R6F), or heated tobacco product (HTP) aerosol at the air liquid interface (ALI). The MPPD-PBPK model predicted the in vivo data from clinical studies within a factor of two, indicating good agreement as noted by WHO International Programme on Chemical Safety (2010) guidance. We then used QIVIVE to derive the exposure concentration (HEC) that matched the estimated in vitro deposition point of departure (POD) (MEC cigarette = 0.38 puffs or 11.6 µg nicotine, HTP = 22.9 puffs or 125.6 µg nicotine) and subsequently derived the equivalent human plasma concentrations. Results indicate that for the 1R6F cigarette, inhaling 1/6th of a stick would be required to induce the same effects observed in vitro, in vivo. Whereas, for HTP it would be necessary to consume 3 sticks simultaneously to induce in vivo the effects observed in vitro. This data further demonstrates the reduced physiological potency potential of HTP aerosol compared to cigarette smoke. The QIVIVE approach demonstrates great promise in assisting human health risk assessments, however, further optimization and standardization are required for the substantiation of a meaningful contribution to tobacco harm reduction by alternative nicotine delivery products.

Research Landscape of Physiologically Based Pharmacokinetic Model Utilization in Different Fields: A Bibliometric Analysis (1999-2023)

PURPOSE

The physiologically based pharmacokinetic (PBPK) modeling has received increasing attention owing to its excellent predictive abilities. However, there has been no bibliometric analysis about PBPK modeling. This research aimed to summarize the research development and hot points in PBPK model utilization overall through bibliometric analysis.

METHODS

We searched for publications related to the PBPK modeling from 1999 to 2023 in the Web of Science Core Collection (WoSCC) database. The Microsoft Office Excel, CiteSpace and VOSviewers were used to perform the analyses.

RESULTS

A total of 4,649 records from 1999 to 2023 were identified, and the largest number of publications focused in the period 2018-2023. The United States was the leading country, and the Environmental Protection Agency (EPA) was the leading institution. The journal Drug Metabolism and Disposition published and co-cited the most articles. Drug-drug interactions, special populations, and new drug development are the main topics in this research field.

CONCLUSION

We first visualize the research landscape and hotspots of the PBPK modeling through bibliometric methods. Our study provides a better understanding for researchers, especially beginners about the dynamization of PBPK modeling and presents the relevant trend in the future.

Maximum likelihood estimation of renal transporter ontogeny profiles for pediatric PBPK modeling

Optimal treatment of infants with many renally cleared drugs must account for maturational differences in renal transporter (RT) activity. Pediatric physiologically-based pharmacokinetic (PBPK) models may incorporate RT activity, but this requires ontogeny profiles for RT activity in children, especially neonates, to predict drug disposition. Therefore, RT expression measurements from human kidney postmortem cortical tissue samples were normalized to represent a fraction of mature RT activity. Using these data, maximum likelihood estimated the distributions of RT activity across the pediatric age spectrum, including preterm and term neonates. PBPK models of four RT substrates (acyclovir, ciprofloxacin, furosemide, and meropenem) were evaluated with and without ontogeny profiles using average fold error (AFE), absolute average fold error (AAFE), and proportion of observations within the 5-95% prediction interval. Novel maximum likelihood profiles estimated ontogeny distributions for the following RT: OAT1, OAT3, OCT2, P-gp, URAT1, BCRP, MATE1, MRP2, MRP4, and MATE-2 K. Profiles for OAT3, P-gp, and MATE1 improved infant furosemide and neonate meropenem PBPK model AFE from 0.08 to 0.70 and 0.53 to 1.34 and model AAFE from 12.08 to 1.44 and 2.09 to 1.36, respectively, and improved the percent of data within the 5-95% prediction interval from 48% to 98% for neonatal ciprofloxacin simulations, respectively. Even after accounting for other critical population-specific maturational differences, novel RT ontogeny profiles substantially improved neonatal PBPK model performance, providing validated estimates of maturational differences in RT activity for optimal dosing in children.

Docetaxel, cyclophosphamide, and epirubicin: application of PBPK modeling to gain new insights for drug-drug interactions

The new adjuvant chemotherapy of docetaxel, epirubicin, and cyclophosphamide has been recommended for treating breast cancer. It is necessary to investigate the potential drug-drug Interactions (DDIs) since they have a narrow therapeutic window in which slight differences in exposure might result in significant differences in treatment efficacy and tolerability. To guide clinical rational drug use, this study aimed to evaluate the DDI potentials of docetaxel, cyclophosphamide, and epirubicin in cancer patients using physiologically based pharmacokinetic (PBPK) models. The GastroPlus™ was used to develop the PBPK models, which were refined and validated with observed data. The established PBPK models accurately described the pharmacokinetics (PKs) of three drugs in cancer patients, and the predicted-to-observed ratios of all the PK parameters met the acceptance criterion. The PBPK model predicted no significant changes in plasma concentrations of these drugs during co-administration, which was consistent with the observed clinical phenomenon. Besides, the verified PBPK models were then used to predict the effect of other Cytochrome P450 3A4 (CYP3A4) inhibitors/inducers on these drug exposures. In the DDI simulation, strong CYP3A4 modulators changed the exposure of three drugs by 0.71-1.61 fold. Therefore, patients receiving these drugs in combination with strong CYP3A4 inhibitors should be monitored regularly to prevent adverse reactions. Furthermore, co-administration of docetaxel, cyclophosphamide, or epirubicin with strong CYP3A4 inducers should be avoided. In conclusion, the PBPK models can be used to further investigate the DDI potential of each drug and to develop dosage recommendations for concurrent usage by additional perpetrators or victims.

Predictions of Bedaquiline Central Nervous System Exposure in Patients with Tuberculosis Meningitis Using Physiologically based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The use of bedaquiline as a treatment option for drug-resistant tuberculosis meningitis (TBM) is of interest to address the increased prevalence of resistance to first-line antibiotics. To this end, we describe a whole-body physiologically based pharmacokinetic (PBPK) model for bedaquiline to predict central nervous system (CNS) exposure.

METHODS

A whole-body PBPK model was developed for bedaquiline and its metabolite, M2. The model included compartments for brain and cerebrospinal fluid (CSF). Model predictions were evaluated by comparison to plasma PK time profiles following different dosing regimens and sparse CSF concentrations data from patients. Simulations were then conducted to compare CNS and lung exposures to plasma exposure at clinically relevant dosing schedules.

RESULTS

The model appropriately described the observed plasma and CSF bedaquiline and M2 concentrations from patients with pulmonary tuberculosis (TB). The model predicted a high impact of tissue binding on target site drug concentrations in CNS. Predicted unbound exposures within brain interstitial exposures were comparable with unbound vascular plasma and unbound lung exposures. However, unbound brain intracellular exposures were predicted to be 7% of unbound vascular plasma and unbound lung intracellular exposures.

CONCLUSIONS

The whole-body PBPK model for bedaquiline and M2 predicted unbound concentrations in brain to be significantly lower than the unbound concentrations in the lung at clinically relevant doses. Our findings suggest that bedaquiline may result in relatively inferior efficacy against drug-resistant TBM when compared with efficacy against drug-resistant pulmonary TB.

Physiologically based pharmacokinetic modeling to predict maternal pharmacokinetics and fetal carbamazepine exposure during pregnancy

Carbamazepine is an antiepileptic drug commonly used in pregnant women, during which the physiological changes may affect its efficacy. The aim of this study was to establish a physiologically based pharmacokinetic (PBPK) model of carbamazepine and its active metabolite carbamazepine-10,11-epoxide, and simulate maternal and fetal pharmacokinetic changes of carbamazepine and carbamazepine-10,11-epoxide in different trimesters and propose dose adjustment. We established pregnancy PBPK models for carbamazepine and carbamazepine-10,11-epoxide in PK-Sim® and Mobi® and validated the models with observed data from clinical reports. The placental transfer parameters obtained using different methods were also imported into the model and compared with the observed data to establish and validate fetal pharmacokinetic curves. The simulated results showed that mean steady-state trough plasma concentration of carbamazepine decreased by 27, 43.1, and 52 % during the first, second, and third trimesters, respectively. Therefore, to achieve an optimum therapeutic concentration, administering at least 1.4, 1.8, and 2.1 times the baseline dose of carbamazepine in the first, second, and third trimesters, respectively can be used as a dose reference. In conclusion, this study established and validated a pregnancy PBPK model of carbamazepine and carbamazepine-10,11-epoxide to assess exposure in pregnant women and fetuses, which provided a reference for the dosage adjustment of carbamazepine during pregnancy.

Physiologically Based Pharmacokinetic Modeling of Nilotinib for Drug-Drug Interactions, Pediatric Patients, and Pregnancy and Lactation

Nilotinib is a second-generation BCR-ABL tyrosine kinase inhibitor for the treatment of Philadelphia chromosome-positive chronic myeloid leukemia in both adult and pediatric patients. The pharmacokinetics (PK) of nilotinib in specific populations such as pregnant and lactating people remain poorly understood. Therefore, the objectives of the current study were to develop a physiologically based pharmacokinetic (PBPK) model to predict nilotinib PK in virtual drug-drug interaction (DDI) studies, as well as in pediatric, pregnant, and lactating populations. The nilotinib PBPK model was built in PK-Sim, which is part of the free and open-source software Open Systems Pharmacology. The observed clinical data for the validation of the nilotinib models were obtained from the literature. The model reasonably predicted nilotinib concentrations in the adult population; the DDIs between nilotinib and rifampin or ketoconazole in the adult population; and the PK in the pediatric, pregnant, and lactating populations, although in the latter 2 populations plasma concentrations were slightly underestimated. The ratio of predicted versus observed PK parameters for the adult model ranged from 0.71 to 1.11 for area under the concentration-time curve and 0.55 to 0.95 for maximum concentration. For the DDI, the predicted area under the concentration-time curve ratio and maximum concentration ratio fell within the Guest criterion. The current study demonstrated the utility of using PBPK modeling to understand the mechanistic basis of PK differences between adults and specific populations, such as pediatrics, and pregnant and lactating individuals, indicating that this technology can potentially inform or optimize dosing conditions in specific populations.

Physiologically-based pharmacokinetic model for evaluating gender-specific exposures of N-nitrosodimethylamine (NDMA)

N-nitrosodimethylamine (NDMA) is classified as a human carcinogen and could be produced by both natural and industrial processes. Although its toxicity and histopathology have been well-studied in animal species, there is insufficient data on the blood and tissue exposures that can be correlated with the toxicity of NDMA. The purpose of this study was to evaluate gender-specific pharmacokinetics/toxicokinetics (PKs/TKs), tissue distribution, and excretion after the oral administration of three different doses of NDMA in rats using a physiologically-based pharmacokinetic (PBPK) model. The major target tissues for developing the PBPK model and evaluating dose metrics of NDMA included blood, gastrointestinal (GI) tract, liver, kidney, lung, heart, and brain. The predictive performance of the model was validated using sensitivity analysis, (average) fold error, and visual inspection of observations versus predictions. Then, a Monte Carlo simulation was performed to describe the magnitudes of inter-individual variability and uncertainty of the single model predictions. The developed PBPK model was applied for the exposure simulation of daily oral NDMA to estimate blood concentration ranges affecting health effects following acute-duration (≤ 14 days), intermediate-duration (15-364 days), and chronic-duration (≥ 365 days) intakes. The results of the study could be used as a scientific basis for interpreting the correlation between in vivo exposures and toxicological effects of NDMA.

Physiologically Based Pharmacokinetic (PBPK) Modeling to Predict CYP3A-Mediated Drug Interaction between Saxagliptin and Nicardipine: Bridging Rat-to-Human Extrapolation

The aim of this study was to predict the cytochrome P450 3A (CYP3A)-mediated drug-drug interactions (DDIs) between saxagliptin and nicardipine using a physiologically based pharmacokinetic (PBPK) model. Initially, in silico and in vitro parameters were gathered from experiments or the literature to construct PBPK models for each drug in rats. These models were integrated to predict the DDIs between saxagliptin, metabolized via CYP3A2, and nicardipine, exhibiting CYP3A inhibitory activity. The rat DDI PBPK model was completed by optimizing parameters using experimental rat plasma concentrations after co-administration of both drugs. Following co-administration in Sprague-Dawley rats, saxagliptin plasma concentration significantly increased, resulting in a 2.60-fold rise in AUC, accurately predicted by the rat PBPK model. Subsequently, the workflow of the rat PBPK model was applied to humans, creating a model capable of predicting DDIs between the two drugs in humans. Simulation from the human PBPK model indicated that nicardipine co-administration in humans resulted in a nearly unchanged AUC of saxagliptin, with an approximate 1.05-fold change, indicating no clinically significant changes and revealing a lack of direct translation of animal interaction results to humans. The animal-to-human PBPK model extrapolation used in this study could enhance the reliability of predicting drug interactions in clinical settings where DDI studies are challenging.

Physiologically based pharmacokinetic modeling of apixaban to predict exposure in populations with hepatic and renal impairment and elderly populations

BACKGROUND

Apixaban is a factor Xa inhibitor with a limited therapeutic index that belongs to the family of oral direct anticoagulants. The pharmacokinetic (PK) behavior of apixaban may be altered in elderly populations and populations with renal or hepatic impairment, necessitating dosage adjustments.

METHODS

This study was conducted to examine how the physiologically based pharmacokinetic (PBPK) model describes the PKs of apixaban in adult and elderly populations and to determine the PKs of apixaban in elderly populations with renal and hepatic impairment. After PBPK models were constructed using the reported physicochemical properties of apixaban and clinical data, they were validated using data from clinical studies involving various dose ranges. Comparing predicted and observed blood concentration data and PK parameters was utilized to evaluate the model's fit performance.

RESULTS

Doses should be reduced to approximately 70% of the healthy adult population for the healthy elderly population to achieve the same PK exposure; approximately 88%, 71%, and 89% of that for the elderly populations with mild, moderate, and severe renal impairment, respectively; and approximately 96%, 81%, and 58% of that for the Child Pugh-A, Child Pugh-B, and Child Pugh-C hepatic impairment elderly populations, respectively to achieve the same PK exposure.

CONCLUSION

The findings indicate that the renal and hepatic function might be considered for apixaban therapy in Chinese elderly patients and the PBPK model can be used to optimize dosage regimens for specific populations.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

Integrating the gut microbiome and pharmacology

The gut microbiome harbors trillions of organisms that contribute to human health and disease. These bacteria can also affect the properties of medical drugs used to treat these diseases, and drugs, in turn, can reshape the microbiome. Research addressing interdependent microbiome-host-drug interactions thus has broad impact. In this Review, we discuss these interactions from the perspective of drug bioavailability, absorption, metabolism, excretion, toxicity, and drug-mediated microbiome modulation. We survey approaches that aim to uncover the mechanisms underlying these effects and opportunities to translate this knowledge into new strategies to improve the development, administration, and monitoring of medical drugs.

Applied physiologically-based pharmacokinetic modeling to assess uridine diphosphate-glucuronosyltransferase-mediated drug-drug interactions for Vericiguat

Vericiguat (Verquvo; US: Merck, other countries: Bayer) is a novel drug for the treatment of chronic heart failure. Preclinical studies have demonstrated that the primary route of metabolism for vericiguat is glucuronidation, mainly catalyzed by uridine diphosphate-glucuronosyltransferase (UGT)1A9 and to a lesser extent UGT1A1. Whereas a drug-drug interaction (DDI) study of the UGT1A9 inhibitor mefenamic acid showed a 20% exposure increase, the effect of UGT1A1 inhibitors has not been assessed clinically. This modeling study describes a physiologically-based pharmacokinetic (PBPK) approach to complement the clinical DDI liability assessment and support prescription labeling. A PBPK model of vericiguat was developed based on in vitro and clinical data, verified against data from the mefenamic acid DDI study, and applied to assess the UGT1A1 DDI liability by running an in silico DDI study with the UGT1A1 inhibitor atazanavir. A minor effect with an area under the plasma concentration-time curve (AUC) ratio of 1.12 and a peak plasma concentration ratio of 1.04 was predicted, which indicates that there is no clinically relevant DDI interaction anticipated. Additionally, the effect of potential genetic polymorphisms of UGT1A1 and UGT1A9 was evaluated, which showed that an average modest increase of up to 1.7-fold in AUC may be expected in the case of concomitantly reduced UGT1A1 and UGT1A9 activity for subpopulations expressing non-wild-type variants for both isoforms. This study is a first cornerstone to qualify the PK-Sim platform for use of UGT-mediated DDI predictions, including PBPK models of perpetrators, such as mefenamic acid and atazanavir, and sensitive UGT substrates, such as dapagliflozin and raltegravir.

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

Background

Candesartan cilexetil is a widely used angiotensin II receptor blocker with minimal adverse effects and high tolerability for the treatment of hypertension. Candesartan is administered orally as the prodrug candesartan cilexetil, which is wholly and swiftly converted to the active metabolite candesartan by carboxylesterase during absorption in the intestinal tract. In populations with renal or hepatic impairment, candesartan's pharmacokinetic (PK) behavior may be altered, necessitating dosage adjustments.

Objectives

This study was conducted to examine how the physiologically based PK (PBPK) model characterizes the PKs of candesartan in adult and geriatric populations and to predict the PKs of candesartan in elderly populations with renal and hepatic impairment.

Design

After developing PBPK models using the reported physicochemical properties of candesartan and clinical data, these models were validated using data from clinical investigations involving various dose ranges.

Methods

Comparing predicted and observed blood concentration data and PK parameters was used to assess the fit performance of the models.

Results

Doses should be reduced to approximately 94% of Chinese healthy adults for the Chinese healthy elderly population; approximately 92%, 68%, and 64% of that of the Chinese healthy adult dose in elderly populations with mild, moderate, and severe renal impairment, respectively; and approximately 72%, 71%, and 52% of that of the Chinese healthy adult dose in elderly populations with Child-Pugh-A, Child-Pugh-B, and Child-Pugh-C hepatic impairment, respectively.

Conclusion

The results suggest that the PBPK model of candesartan can be utilized to optimize dosage regimens for special populations.

A Physiologically Based Pharmacokinetic Model to Predict Systemic Ondansetron Concentration in Liver Cirrhosis Patients

INTRODUCTION

Ondansetron is a drug that is routinely prescribed for the management of nausea and vomiting associated with cancer, radiation therapy, and surgical operations. It is mainly metabolized in the liver, and it might accumulate in patients with hepatic impairment and lead to unwanted adverse events.

METHODS

A physiologically based pharmacokinetic (PBPK) model was developed to predict the exposure of ondansetron in healthy and liver cirrhosis populations. The population-based PBPK simulator PK-Sim was utilized for simulating ondansetron exposure in healthy and liver cirrhosis populations.

RESULTS

The developed model successfully described the pharmacokinetics of ondansetron in healthy and liver cirrhosis populations. The predicted area under the curve, maximum systemic concentration, and clearance were within the allowed twofold range. The exposure of ondansetron in the population of Child-Pugh class C has doubled in comparison to Child-Pugh class A. The dose has to be adjusted for liver cirrhosis patients to ensure comparable exposure to a healthy population.

CONCLUSION

In this study, the developed PBPK model has described the pharmacokinetics of ondansetron successfully. The PBPK model has been successfully evaluated to be used as a tool for dose adjustments in liver cirrhosis patients.

Physiologically based pharmacokinetic (PBPK) modeling to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

Combining pharmacokinetic and electrophysiological models for early prediction of drug-induced arrhythmogenicity

BACKGROUND AND OBJECTIVE

In silico methods are gaining attention for predicting drug-induced Torsade de Pointes (TdP) in different stages of drug development. However, many computational models tended not to account for inter-individual response variability due to demographic covariates, such as sex, or physiologic covariates, such as renal function, which may be crucial when predicting TdP. This study aims to compare the effects of drugs in male and female populations with normal and impaired renal function using in silico methods.

METHODS

Pharmacokinetic models considering sex and renal function as covariates were implemented from data published in pharmacokinetic studies. Drug effects were simulated using an electrophysiologically calibrated population of cellular models of 300 males and 300 females. The population of models was built by modifying the endocardial action potential model published by O'Hara et al. (2011) according to the experimentally measured gene expression levels of 12 ion channels.

RESULTS

Fifteen pharmacokinetic models for CiPA drugs were implemented and validated in this study. Eight pharmacokinetic models included the effect of renal function and four the effect of sex. The mean difference in action potential duration (APD) between male and female populations was 24.9 ms (p<0.05). Our simulations indicated that women with impaired renal function were particularly susceptible to drug-induced arrhythmias, whereas healthy men were less prone to TdP. Differences between patient groups were more pronounced for high TdP-risk drugs. The proposed in silico tool also revealed that individuals with impaired renal function, electrophysiologically simulated with hyperkalemia (extracellular potassium concentration [K+]o = 7 mM) exhibited less pronounced APD prolongation than individuals with normal potassium levels. The pharmacokinetic/electrophysiological framework was used to determine the maximum safe dose of dofetilide in different patient groups. As a proof of concept, 3D simulations were also run for dofetilide obtaining QT prolongation in accordance with previously reported clinical values.

CONCLUSIONS

This study presents a novel methodology that combines pharmacokinetic and electrophysiological models to incorporate the effects of sex and renal function into in silico drug simulations and highlights their impact on TdP-risk assessment. Furthermore, it may also help inform maximum dose regimens that ensure TdP-related safety in a specific sub-population of patients.

Comparison of the predictive accuracy of the physiologically based pharmacokinetic (PBPK) model and population pharmacokinetic (PPK) model of vancomycin in Japanese patients with MRSA infection

INTRODUCTION

The latest therapeutic drug monitoring guidelines for vancomycin (VCM) recommend that area under the concentration-time curve is estimated based on model-informed precision dosing and used to evaluate efficacy and safety. Therefore, we predicted VCM concentrations in individual methicillin-resistant Staphylococcus aureus-infected patients using existing a physiologically based pharmacokinetic (PBPK) model and 1- and 2-compartment population pharmacokinetic (PPK) models and confirmed and verified the accuracy of the PBPK model in estimating VCM concentrations with the PPK model.

METHODS

The subjects of the study are 20 patients, and the predicted concentrations were evaluated by comparing the observed and predicted trough and peak values of VCM concentrations for individual patients.

RESULTS

The results showed good correlation between the observed and predicted trough and peak concentrations of VCM was observed generally in the PBPK model, R2 values of 0.72, 0.62, and 0.40 with trough values of 0.49, 0.40, and 0.34 with peak values for PBPK model, 1-compartment, and 2-compartment model, respectively.

CONCLUSIONS

Although the performance of the PBPK model is not as predictive as the PPK model, generally similar predictive trends were obtained, suggesting that it may be a valuable tool for rapid and accurate prediction of AUC for VCM.

Increased sinusoidal export of drug glucuronides is a compensative mechanism in liver cirrhosis of mice

Rationale: Liver cirrhosis is known to affect drug pharmacokinetics, but the functional assessment of the underlying pathophysiological alterations in drug metabolism is difficult. Methods: Cirrhosis in mice was induced by repeated treatment with carbon tetrachloride for 12 months. A cocktail of six drugs was administered, and parent compounds as well as phase I and II metabolites were quantified in blood, bile, and urine in a time-dependent manner. Pharmacokinetics were modeled in relation to the altered expression of metabolizing enzymes. In discrepancy with computational predictions, a strong increase of glucuronides in blood was observed in cirrhotic mice compared to vehicle controls. Results: The deviation between experimental findings and computational simulations observed by analyzing different hypotheses could be explained by increased sinusoidal export and corresponded to increased expression of export carriers (Abcc3 and Abcc4). Formation of phase I metabolites and clearance of the parent compounds were surprisingly robust in cirrhosis, although the phase I enzymes critical for the metabolism of the administered drugs in healthy mice, Cyp1a2 and Cyp2c29, were downregulated in cirrhotic livers. RNA-sequencing revealed the upregulation of numerous other phase I metabolizing enzymes which may compensate for the lost CYP isoenzymes. Comparison of genome-wide data of cirrhotic mouse and human liver tissue revealed similar features of expression changes, including increased sinusoidal export and reduced uptake carriers. Conclusion: Liver cirrhosis leads to increased blood concentrations of glucuronides because of increased export from hepatocytes into the sinusoidal blood. Although individual metabolic pathways are massively altered in cirrhosis, the overall clearance of the parent compounds was relatively robust due to compensatory mechanisms.

Physiologically-Based Pharmacokinetic Modeling for Drug Dosing in Pediatric Patients: A Tutorial for a Pragmatic Approach in Clinical Care

It is well-accepted that off-label drug dosing recommendations for pediatric patients should be based on the best available evidence. However, the available traditional evidence is often low. To bridge this gap, physiologically-based pharmacokinetic (PBPK) modeling is a scientifically well-founded tool that can be used to enable model-informed dosing (MID) recommendations in children in clinical practice. In this tutorial, we provide a pragmatic, PBPK-based pediatric modeling workflow. For this approach to be successfully implemented in pediatric clinical practice, a thorough understanding of the model assumptions and limitations is required. More importantly, careful evaluation of an MID approach within the context of overall benefits and the potential risks is crucial. The tutorial is aimed to help modelers, researchers, and clinicians, to effectively use PBPK simulations to support pediatric drug dosing.

A physiologically based model of bile acid metabolism in mice

Bile acid (BA) metabolism is a complex system that includes a wide variety of primary and secondary, as well as conjugated and unconjugated BAs that undergo continuous enterohepatic circulation (EHC). Alterations in both composition and dynamics of BAs have been associated with various diseases. However, a mechanistic understanding of the relationship between altered BA metabolism and related diseases is lacking. Computational modeling may support functional analyses of the physiological processes involved in the EHC of BAs along the gut-liver axis. In this study, we developed a physiologically based model of murine BA metabolism describing synthesis, hepatic and microbial transformations, systemic distribution, excretion, and EHC of BAs at the whole-body level. For model development, BA metabolism of specific pathogen-free (SPF) mice was characterized in vivo by measuring BA levels and composition in various organs, expression of transporters along the gut, and cecal microbiota composition. We found significantly different BA levels between male and female mice that could only be explained by adjusted expression of the hepatic enzymes and transporters in the model. Of note, this finding was in agreement with experimental observations. The model for SPF mice could also describe equivalent experimental data in germ-free mice by specifically switching off microbial activity in the intestine. The here presented model can therefore facilitate and guide functional analyses of BA metabolism in mice, e.g., the effect of pathophysiological alterations on BA metabolism and translation of results from mouse studies to a clinically relevant context through cross-species extrapolation.

Development and Validation of Physiologically Based Pharmacokinetic Model of Levetiracetam to Predict Exposure and Dose Optimization in Pediatrics

Levetiracetam (Lev) is an antiepileptic drug that has been increasingly used in the epilepsy pediatric population in recent years, but its pharmacokinetic behavior in pediatric population needs to be characterized clearly. Clinical trials for the pediatric drug remain difficult to conduct due to ethical and practical factors. The purpose of this study was to use the physiologically based pharmacokinetic (PBPK) model to predict changes in plasma exposure of Lev in pediatric patients and to provide recommendations for dose adjustment. A PBPK model of Lev in adults was developed using PK-Sim® software and extrapolated to the entire age range of the pediatric population. The model was evaluated using clinical pharmacokinetic data. The results showed the good fit between predictions and observations of the adult and pediatric models. The recommended doses for neonates, infants and children are 0.78, 1.67 and 1.22 times that of adults, respectively. Moreover, at the same dose, plasma exposure in adolescents was similar to that of adults. The PBPK models of Lev for adults and pediatrics were successfully developed and validated to provide a reference for the rational administration of drugs in the pediatric population.

Mathematical Modeling of In Vivo Alpha Particle Generators and Chelator Stability

Background: Targeted α particle therapy using long-lived in vivo α particle generators is cytotoxic to target tissues. However, the redistribution of released radioactive daughters through the circulation should be considered. A mathematical model was developed to describe the physicochemical kinetics of 212Pb-labeled pharmaceuticals and its radioactive daughters. Materials and Methods: A bolus of 212Pb-labeled pharmaceuticals injected in a developed compartmental model was simulated. The contributions of chelated and free radionuclides to the total released energy were investigated for different dissociation fractions of 212Bi for different chelators, for example, 36% for DOTA. The compartmental model was applied to describe a 212Bi retention study and to assess the stability of the 212Bi-1,4,7,10-tetrakis(carbamoylmethyl)-1,4,7,10-tetraazacyclododecane (212Bi-DOTAM) complex after β- decay of 212Pb. Results: The simulation of the injection showed that α emissions contribute 75% to the total released energy, mostly from 212Po (72%). The simulation of the 212Bi retention study showed that (16 ± 5)% of 212Bi atoms dissociate from the 212Bi-DOTAM complexes. The fractions of energies released by free radionuclides were 21% and 38% for DOTAM and DOTA chelators, respectively. Conclusion: The developed α particle generator model allows for simulating the radioactive kinetics of labeled and unlabeled pharmaceuticals being released from the chelating system due to a preceding disintegration.

Development of a physiologically based toxicokinetic model for lead in pregnant women: The role of bone tissue in the maternal and fetal internal exposure

Epidemiological studies have shown associations between prenatal exposure to lead (Pb) and neurodevelopmental effects in young children. Prenatal exposure is generally characterized by measuring the concentration in the umbilical cord at delivery or in the maternal blood during pregnancy. To assess internal Pb exposure during prenatal life, we developed a pregnancy physiologically based pharmacokinetic (p-PBPK) model that to simulates Pb levels in blood and target tissues in the fetus, especially during critical periods for brain development. An existing Pb PBPK model was adapted to pregnant women and fetuses. Using data from literature, both the additional maternal bone remodeling, that causes Pb release into the blood, and the Pb placental transfers were estimated by Bayesian inference. Additional maternal bone remodeling was estimated to start at 21.6 weeks. Placental transfers were estimated between 4.6 and 283 L.day-1 at delivery with high interindividual variability. Once calibrated, the p-PBPK model was used to simulate fetal exposure to Pb. Internal fetal exposure greatly varies over the pregnancy with two peaks of Pb levels in blood and brain at the end of the 1st and 3rd trimesters. Sensitivity analysis shows that the fetal blood lead levels are affected by the maternal burden of bone Pb via maternal bone remodeling and by fetal bone formation at different pregnancy stages. Coupling the p-PBPK model with an effect model such as an adverse outcome pathway could help to predict the effects on children's neurodevelopment.

Uncertainty assessment of proarrhythmia predictions derived from multi-level in silico models

In silico methods can be used for an early assessment of arrhythmogenic properties of drug candidates. However, their use for decision-making is conditioned by the possibility to estimate the predictions' uncertainty. This work describes our efforts to develop uncertainty quantification methods for the predictions produced by multi-level proarrhythmia models. In silico models used in this field usually start with experimental or predicted IC50 values that describe drug-induced ion channel blockade. Using such inputs, an electrophysiological model computes how the ion channel inhibition, exerted by a drug in a certain concentration, translates to an altered shape and duration of the action potential in cardiac cells, which can be represented as arrhythmogenic risk biomarkers such as the APD90. Using this framework, we identify the main sources of aleatory and epistemic uncertainties and propose a method based on probabilistic simulations that replaces single-point estimates predicted using multiple input values, including the IC50s and the electrophysiological parameters, by distributions of values. Two selected variability types associated with these inputs are then propagated through the multi-level model to estimate their impact on the uncertainty levels in the output, expressed by means of intervals. The proposed approach yields single predictions of arrhythmogenic risk biomarkers together with value intervals, providing a more comprehensive and realistic description of drug effects on a human population. The methodology was tested by predicting arrhythmogenic biomarkers on a series of twelve well-characterised marketed drugs, belonging to different arrhythmogenic risk classes.

Physiologically based pharmacokinetic (PBPK) modelling of oral drug absorption in older adults - an AGePOP review

The older population consisting of persons aged 65 years or older is the fastest-growing population group and also the major consumer of pharmaceutical products. Due to the heterogenous ageing process, this age group shows high interindividual variability in the dose-exposure-response relationship and, thus, a prediction of drug safety and efficacy is challenging. Although physiologically based pharmacokinetic (PBPK) modelling is a well-established tool to inform and confirm drug dosing strategies during drug development for special population groups, age-related changes in absorption are poorly accounted for in current PBPK models. The purpose of this review is to summarise the current state-of-knowledge in terms of physiological changes with increasing age that can influence the oral absorption of dosage forms. The capacity of common PBPK platforms to incorporate these changes and describe the older population is also discussed, as well as the implications of extrinsic factors such as drug-drug interactions associated with polypharmacy on the model development process. The future potential of this field will rely on addressing the gaps identified in this article, which can subsequently supplement in-vitro and in-vivo data for more robust decision-making on the adequacy of the formulation for use in older adults and inform pharmacotherapy.

Investigating Tacrolimus Disposition in Paediatric Patients with a Physiologically Based Pharmacokinetic Model Incorporating CYP3A4 Ontogeny, Mechanistic Absorption and Red Blood Cell Binding

Tacrolimus is a crucial immunosuppressant for organ transplant patients, requiring therapeutic drug monitoring due to its variable exposure after oral intake. Physiologically based pharmacokinetic (PBPK) modelling has provided insights into tacrolimus disposition in adults but has limited application in paediatrics. This study investigated age dependency in tacrolimus exposure at the levels of absorption, metabolism, and distribution. Based on the literature data, a PBPK model was developed to predict tacrolimus exposure in adults after intravenous and oral administration. This model was then extrapolated to the paediatric population, using a unique reference dataset of kidney transplant patients. Selecting adequate ontogeny profiles for hepatic and intestinal CYP3A4 appeared critical to using the model in children. The best model performance was achieved by using the Upreti ontogeny in both the liver and intestines. To mechanistically evaluate the impact of absorption on tacrolimus exposure, biorelevant in vitro solubility and dissolution data were obtained. A relatively fast and complete release of tacrolimus from its amorphous formulation was observed when mimicking adult or paediatric dissolution conditions (dose, fluid volume). In both the adult and paediatric PBPK models, the in vitro dissolution profiles could be adequately substituted by diffusion-layer-based dissolution modelling. At the level of distribution, sensitivity analysis suggested that differences in blood plasma partitioning of tacrolimus may contribute to the variability in exposure in paediatric patients.

Physiologically Based Pharmacokinetic Modeling of Extracellular Vesicles

Extracellular vesicles (EVs) are lipid membrane bound-cell-derived structures that are a key player in intercellular communication and facilitate numerous cellular functions such as tumor growth, metastasis, immunosuppression, and angiogenesis. They can be used as a drug delivery platform because they can protect drugs from degradation and target specific cells or tissues. With the advancement in the technologies and methods in EV research, EV-therapeutics are one of the fast-growing domains in the human health sector. Therapeutic translation of EVs in clinics requires assessing the quality, safety, and efficacy of the EVs, in which pharmacokinetics is very crucial. We report here the application of physiologically based pharmacokinetic (PBPK) modeling as a principal tool for the prediction of absorption, distribution, metabolism, and excretion of EVs. To create a PBPK model of EVs, researchers would need to gather data on the size, shape, and composition of the EVs, as well as the physiological processes that affect their behavior in the body. The PBPK model would then be used to predict the pharmacokinetics of drugs delivered via EVs, such as the rate at which the drug is absorbed and distributed throughout the body, the rate at which it is metabolized and eliminated, and the maximum concentration of the drug in the body. This information can be used to optimize the design of EV-based drug delivery systems, including the size and composition of the EVs, the route of administration, and the dose of the drug. There has not been any dedicated review article that describes the PBPK modeling of EV. This review provides an overview of the absorption, distribution, metabolism, and excretion (ADME) phenomena of EVs. In addition, we will briefly describe the different computer-based modeling approaches that may help in the future of EV-based therapeutic research.

Revisiting the in-vitro and in-vivo considerations for in-silico modelling of complex injectable drug products

Complex injectable drug products (CIDPs) have often been developed to modulate the pharmacokinetics along with efficacy for therapeutic agents used for remediation of chronic disorders. The effective development of CIDPs has exhibited complex kinetics associated with multiphasic drug release from the prepared formulations. Consequently, predictability of pharmacokinetic modelling for such CIDPs has been difficult and there is need for advanced complex computational models for the establishment of accurate prediction models for in-vitro-in-vivo correlation (IVIVC). The computational modelling aims at supplementing the existing knowledge with mathematical equations to develop formulation strategies for generation of predictable and discriminatory IVIVC. Such an approach would help in reduction of the burden of effect of hidden factors on preclinical to clinical translations. Computational tools like physiologically based pharmacokinetics (PBPK) modelling have combined physicochemical and physiological properties along with IVIVC characteristics of clinically used formulations. Such techniques have helped in prediction and understanding of variability in pharmacodynamic parameters of potential generic products to clinically used formulations like Doxil®, Ambisome®, Abraxane® in healthy and diseased population using mathematical equations. The current review highlights the important formulation characteristics, in-vitro, preclinical in-vivo aspects which need to be considered while developing a stimulatory predictive PBPK model in establishment of an IVIVC and in-vitro-in-vivo relationship (IVIVR).

Physiologically-based pharmacokinetic modeling for single and multiple dosing regimens of ceftriaxone in healthy and chronic kidney disease populations: a tool for model-informed precision dosing

Introduction: Ceftriaxone is one of commonly prescribed beta-lactam antibiotics with several label and off-label clinical indications. A high fraction of administered dose of ceftriaxone is excreted renally in an unchanged form, and it may accumulate significantly in patients with impaired renal functions, which may lead to toxicity. Methods: In this study, we employed a physiologically-based pharmacokinetic (PBPK) modeling, as a tool for precision dosing, to predict the biological exposure of ceftriaxone in a virtually-constructed healthy and chronic kidney disease patient populations, with subsequent dosing optimizations. We started developing the model by integrating the physicochemical properties of the drug with biological system information in a PBPK software platform. A PBPK model in an adult healthy population was developed and evaluated visually and numerically with respect to experimental pharmacokinetic data. The model performance was evaluated based on the fold error criteria of the predicted and reported values for different pharmacokinetic parameters. Then, the model was applied to predict drug exposure in CKD patient populations with various degrees of severity. Results: The developed PBPK model was able to precisely describe the pharmacokinetic behavior of ceftriaxone in adult healthy population and in mild, moderate, and severe CKD patient populations. Decreasing the dose by approximately 25% in mild and 50% in moderate to severe renal disease provided a comparable exposure to the healthy population. Based on the simulation of multiple dosing regimens in severe CKD population, it has been found that accumulation of 2 g every 24 h is lower than the accumulation of 1 g every 12 h dosing regimen. Discussion: In this study, the observed concentration time profiles and pharmacokinetic parameters for ceftriaxone were successfully reproduced by the developed PBPK model and it has been shown that PBPK modeling can be used as a tool for precision dosing to suggest treatment regimens in population with renal impairment.

Modelling Based Approaches to Support Generic Drug Regulatory Submissions-Practical Considerations and Case Studies

Model informed drug development (MiDD) is useful to predict in vivo exposure of drugs during various stages of the drug development process. This approach employs a variety of quantitative tools to assess the risks during the drug development process. One important tool in the MiDD tool kit is the Physiologically Based Pharmacokinetic Modelling (PBPK). This tool is extensively used to reduce the development cost and to accelerate the access of medicines to the patients. In this work, we provide an overview of PBPK modelling approaches in the generic drug development process, with a special emphasis on the bio-waiver applications. We describe herein approaches and common pitfalls while submitting model based justifications as a response to the regulatory deficiencies during the generic drug development process. With some in-house case studies, we have attempted to provide a clear path for PBPK model based justifications for bio-waivers. With this review, the gap between theoretical knowledge and practical application of modelling and simulation tools for generic drug product development could be potentially reduced.

Generic Workflow to Predict Medicine Concentrations in Human Milk Using Physiologically-Based Pharmacokinetic (PBPK) Modelling-A Contribution from the ConcePTION Project

Women commonly take medication during lactation. Currently, there is little information about the exposure-related safety of maternal medicines for breastfed infants. The aim was to explore the performance of a generic physiologically-based pharmacokinetic (PBPK) model to predict concentrations in human milk for ten physiochemically diverse medicines. First, PBPK models were developed for "non-lactating" adult individuals in PK-Sim/MoBi v9.1 (Open Systems Pharmacology). The PBPK models predicted the area-under-the-curve (AUC) and maximum concentrations (C_max) in plasma within a two-fold error. Next, the PBPK models were extended to include lactation physiology. Plasma and human milk concentrations were simulated for a three-months postpartum population, and the corresponding AUC-based milk-to-plasma (M/P) ratios and relative infant doses were calculated. The lactation PBPK models resulted in reasonable predictions for eight medicines, while an overprediction of human milk concentrations and M/P ratios (>2-fold) was observed for two medicines. From a safety perspective, none of the models resulted in underpredictions of observed human milk concentrations. The present effort resulted in a generic workflow to predict medicine concentrations in human milk. This generic PBPK model represents an important step towards an evidence-based safety assessment of maternal medication during lactation, applicable in an early drug development stage.

Parameters to Consider for Successful Medication Use in Older Adults - an AGePOP Review

Older adults are the main users of medicine and due to their multimorbidity are often faced/confronted with a complex medication management. This review article provides a brief overview on aspects of medication management such as maintaining a stock of the required medicine, understanding and following the instructions for use, coping with the primary and secondary packaging as well as preparation prior to use. The main focus however is on the drug intake itself and provides an overview about the current understanding of real life dosing conditions of older adults and geriatric patients. Furthermore, it elaborates the acceptability of dosage forms, in particular solid oral dosage forms as they represent the majority of dosage forms taken by these patient populations. An improved understanding of the needs of older adults and geriatric patients, their acceptability of various dosage forms and the circumstances under which they manage their medications, will make the design of more patient-centric drug products possible.

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

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

All Roads Lead to Rome: Enhancing the Probability of Target Attainment with Different Pharmacokinetic/Pharmacodynamic Modelling Approaches

In light of rising antimicrobial resistance and a decreasing number of antibiotics with novel modes of action, it is of utmost importance to accelerate development of novel treatment options. One aspect of acceleration is to understand pharmacokinetics (PK) and pharmacodynamics (PD) of drugs and to assess the probability of target attainment (PTA). Several in vitro and in vivo methods are deployed to determine these parameters, such as time-kill-curves, hollow-fiber infection models or animal models. However, to date the use of in silico methods to predict PK/PD and PTA is increasing. Since there is not just one way to perform the in silico analysis, we embarked on reviewing for which indications and how PK and PK/PD models as well as PTA analysis has been used to contribute to the understanding of the PK and PD of a drug. Therefore, we examined four recent examples in more detail, namely ceftazidime-avibactam, omadacycline, gepotidacin and zoliflodacin as well as cefiderocol. Whereas the first two compound classes mainly relied on the ‘classical’ development path and PK/PD was only deployed after approval, cefiderocol highly profited from in silico techniques that led to its approval. Finally, this review shall highlight current developments and possibilities to accelerate drug development, especially for anti-infectives.

Artificial intelligence and machine learning disciplines with the potential to improve the nanotoxicology and nanomedicine fields: a comprehensive review

The use of nanomaterials in medicine depends largely on nanotoxicological evaluation in order to ensure safe application on living organisms. Artificial intelligence (AI) and machine learning (MI) can be used to analyze and interpret large amounts of data in the field of toxicology, such as data from toxicological databases and high-content image-based screening data. Physiologically based pharmacokinetic (PBPK) models and nano-quantitative structure-activity relationship (QSAR) models can be used to predict the behavior and toxic effects of nanomaterials, respectively. PBPK and Nano-QSAR are prominent ML tool for harmful event analysis that is used to understand the mechanisms by which chemical compounds can cause toxic effects, while toxicogenomics is the study of the genetic basis of toxic responses in living organisms. Despite the potential of these methods, there are still many challenges and uncertainties that need to be addressed in the field. In this review, we provide an overview of artificial intelligence (AI) and machine learning (ML) techniques in nanomedicine and nanotoxicology to better understand the potential toxic effects of these materials at the nanoscale.

Process and System Clearances in Pharmacokinetic Models: Our Basic Clearance Concepts Are Correct

Clearances are important parameters in pharmacokinetic (PK) models. All clearances in PK models are either process clearances that include diffusion, transport, and metabolism clearances or system clearances that include organ and systemic clearance. Clearance and volume of distribution are two independent parameters that characterize drug disposition in both individual compartments and systems of compartments. In this minireview, we show that systemic and organ clearances are net clearances that can be easily derived by partition analysis. When drugs are eliminated from the central compartment by first-order processes, systemic clearance is constant. When drugs are eliminated from a peripheral compartment, instantaneous systemic clearance will vary with time. However, average clearance and clearance at steady state will be constant and will equal dose divided by area under the curve. We show that peripheral elimination will not have a large impact on most pharmacokinetic analyses and that standard models of organ and systemic clearance are useful and appropriate. SIGNIFICANCE STATEMENT: There are two basic kinds of clearances used in pharmacokinetic models, process and system clearances. We show that organ and systemic clearances are net clearances with blood or plasma as the driving concentration. For linear pharmacokinetics, clearance is constant for elimination from the central compartment but varies with time for peripheral elimination. Despite the different kinds of clearance parameters and models, standard clearance models and concepts remain valid.

Considering developmental neurotoxicity in vitro data for human health risk assessment using physiologically-based kinetic modeling: deltamethrin case study

Developmental neurotoxicity (DNT) is a potential hazard of chemicals. Recently, an in vitro testing battery (DNT IVB) was established to complement existing rodent in vivo approaches. Deltamethrin (DLT), a pyrethroid with a well-characterized neurotoxic mode of action, has been selected as a reference chemical to evaluate the performance of the DNT IVB. The present study provides context for evaluating the relevance of these DNT IVB results for the human health risk assessment of DLT by estimating potential human fetal brain concentrations after maternal exposure to DLT. We developed a physiologically based kinetic (PBK) model for rats which was then translated to humans considering realistic in vivo exposure conditions (acceptable daily intake [ADI] for DLT). To address existing uncertainties, we designed case studies considering the most relevant drivers of DLT uptake and distribution. Calculated human fetal brain concentrations were then compared with the lowest benchmark concentration achieved in the DNT IVB. The developed rat PBK model was validated on in vivo rat toxicokinetic data of DLT over a broad range of doses. The uncertainty based case study evaluation confirmed that repeated exposure to DLT at an ADI level would likely result in human fetal brain concentrations far below the in vitro benchmark. The presented results indicate that DLT concentrations in the human fetal brain are highly unlikely to reach concentrations associated with in vitro findings under realistic exposure conditions. Therefore, the new in vitro DNT results are considered to have no impact on the current risk assessment approach.

PBPK-based dose finding for sildenafil in pregnant women for antenatal treatment of congenital diaphragmatic hernia

Sildenafil is a potent vasodilator and phosphodiesterase type five inhibitor, commercially known as Revatio® and approved for the treatment of pulmonary arterial hypertension. Maternal administration of sildenafil during pregnancy is being evaluated for antenatal treatment of several conditions, including the prevention of pulmonary hypertension in fetuses with congenital diaphragmatic hernia. However, determination of a safe and effective maternal dose to achieve adequate fetal exposure to sildenafil remains challenging, as pregnancy almost always is an exclusion criterion in clinical studies. Physiologically-based pharmacokinetic (PBPK) modelling offers an attractive approach for dose finding in this specific population. The aim of this study is to exploit physiologically-based pharmacokinetic modelling to predict the required maternal dose to achieve therapeutic fetal exposure for the treatment congenital diaphragmatic hernia. A full-PBPK model was developed for sildenafil and N-desmethyl-sildenafil using the Simcyp simulator V21 platform, and verified in adult reference individuals, as well as in pregnant women, taking into account maternal and fetal physiology, along with factors known to determine hepatic disposition of sildenafil. Clinical pharmacokinetic data in mother and fetus were previously obtained in the RIDSTRESS study and were used for model verification purposes. Subsequent simulations were performed relying either on measured values for fetal fraction unbound (fu = 0.108) or on values predicted by the simulator (fu = 0.044). Adequate doses were predicted according to the efficacy target of 15 ng/mL (or 38 ng/mL) and safety target of 166 ng/mL (or 409 ng/mL), assuming measured (or predicted) fu values, respectively. Considering simulated median profiles for average steady state sildenafil concentrations, dosing regimens of 130 mg/day or 150 mg/day (administered as t.i.d.), were within the therapeutic window, assuming either measured or predicted fu values, respectively. For safety reasons, dosing should be initiated at 130 mg/day, under therapeutic drug monitoring. Additional experimental measurements should be performed to confirm accurate fetal (and maternal) values for fu. Additional characterization of pharmacodynamics in this specific population is required and may lead to further optimization of the dosing regimen.

Quantitative Characterization of Clinically Relevant Drug-Metabolizing Enzymes and Transporters in Rat Liver and Intestinal Segments for Applications in PBPK Modeling

Rats are extensively used as a preclinical model for assessing drug pharmacokinetics (PK) and tissue distribution; however, successful translation of the rat data requires information on the differences in drug metabolism and transport mechanisms between rats and humans. To partly fill this knowledge gap, we quantified clinically relevant drug-metabolizing enzymes and transporters (DMETs) in the liver and different intestinal segments of Sprague-Dawley rats. The levels of DMET proteins in rats were quantified using the global proteomics-based total protein approach (TPA) and targeted proteomics. The abundance of the major DMET proteins was largely comparable using quantitative global and targeted proteomics. However, global proteomics-based TPA was able to detect and quantify a comprehensive list of 66 DMET proteins in the liver and 37 DMET proteins in the intestinal segments of SD rats without the need for peptide standards. Cytochrome P450 (Cyp) and UDP-glycosyltransferase (Ugt) enzymes were mainly detected in the liver with the abundance ranging from 8 to 6502 and 74 to 2558 pmol/g tissue. P-gp abundance was higher in the intestine (124.1 pmol/g) as compared to that in the liver (26.6 pmol/g) using the targeted analysis. Breast cancer resistance protein (Bcrp) was most abundant in the intestinal segments, whereas organic anion transporting polypeptides (Oatp) 1a1, 1a4, 1b2, and 2a1 and multidrug resistance proteins (Mrp) 2 and 6 were predominantly detected in the liver. To demonstrate the utility of these data, we modeled digoxin PK by integrating protein abundance of P-gp and Cyp3a2 into a physiologically based PK (PBPK) model constructed using PK-Sim software. The model was able to reliably predict the systemic as well as tissue concentrations of digoxin in rats. These findings suggest that proteomics-informed PBPK models in preclinical species can allow mechanistic PK predictions in animal models including tissue drug concentrations.

Gene expression databases for physiologically based pharmacokinetic modeling of humans and animal species

In drug research, developing a sound understanding of the key mechanistic drivers of pharmacokinetics (PK) for new molecular entities is essential for human PK and dose predictions. Here, characterizing the absorption, distribution, metabolism, and excretion (ADME) processes is crucial for a mechanistic understanding of the drug-target and drug-body interactions. Sufficient knowledge on ADME processes enables reliable interspecies and human PK estimations beyond allometric scaling. The physiologically based PK (PBPK) modeling framework allows the explicit consideration of organ-specific ADME processes. The sum of all passive and active ADME processes results in the observed plasma PK. Gene expression information can be used as surrogate for protein abundance and activity within PBPK models. The absolute and relative expression of ADME genes can differ between species and strains. This is affecting both, the PK and pharmacodynamics and is therefore posing a challenge for the extrapolation from preclinical findings to humans. We developed an automated workflow that generates whole-body gene expression databases for humans and other species relevant in drug development, animal health, nutritional sciences, and toxicology. Solely, bulk RNA-seq data curated and provided by the Swiss Institute of Bioinformatics from healthy, normal, and untreated primary tissue samples were considered as an unbiased reference of normal gene expression. The databases are interoperable with the Open Systems Pharmacology Suite (PK-Sim and MoBi) and enable seamless access to a central source of curated cross-species gene expression data. This will increase data transparency, increase reliability and reproducibility of PBPK model simulations, and accelerate mechanistic PBPK model development in the future.

Mathematical modeling of regulatory networks of intracellular processes - Aims and selected methods

Regulatory networks structure and signaling pathways dynamics are uncovered in time- and resource consuming experimental work. However, it is increasingly supported by modeling, analytical and computational techniques as well as discrete mathematics and artificial intelligence applied to to extract knowledge from existing databases. This review is focused on mathematical modeling used to analyze dynamics and robustness of these networks. This paper presents a review of selected modeling methods that facilitate advances in molecular biology.

Quantitative systems pharmacology of the eye: Tools and data for ocular QSP

Good eyesight belongs to the most-valued attributes of health, and diseases of the eye are a significant healthcare burden. Case numbers are expected to further increase in the next decades due to an aging society. The development of drugs in ophthalmology, however, is difficult due to limited accessibility of the eye, in terms of drug administration and in terms of sampling of tissues for drug pharmacokinetics (PKs) and pharmacodynamics (PDs). Ocular quantitative systems pharmacology models provide the opportunity to describe the distribution of drugs in the eye as well as the resulting drug-response in specific segments of the eye. In particular, ocular physiologically-based PK (PBPK) models are necessary to describe drug concentration levels in different regions of the eye. Further, ocular effect models using molecular data from specific cellular systems are needed to develop dose-response correlations. We here describe the current status of PK/PBPK as well as PD models for the eyes and discuss cellular systems, data repositories, as well as animal models in ophthalmology. The application of the various concepts is highlighted for the development of new treatments for postoperative fibrosis after glaucoma surgery.

Nanoparticle biodistribution coefficients: A quantitative approach for understanding the tissue distribution of nanoparticles

The objective of this manuscript is to provide quantitative insights into the tissue distribution of nanoparticles. Published pharmacokinetics of nanoparticles in plasma, tumor and 13 different tissues of mice were collected from literature. A total of 2018 datasets were analyzed and biodistribution of graphene oxide, lipid, polymeric, silica, iron oxide and gold nanoparticles in different tissues was quantitatively characterized using Nanoparticle Biodistribution Coefficients (NBC). It was observed that typically after intravenous administration most of the nanoparticles are accumulated in the liver (NBC = 17.56 %ID/g) and spleen (NBC = 12.1 %ID/g), while other tissues received less than 5 %ID/g. NBC values for kidney, lungs, heart, bones, brain, stomach, intestine, pancreas, skin, muscle and tumor were found to be 3.1 %ID/g, 2.8 %ID/g, 1.8 %ID/g, 0.9 %ID/g, 0.3 %ID/g, 1.2 %ID/g, 1.8 %ID/g, 1.2 %ID/g, 1.0 %ID/g, 0.6 %ID/g and 3.4 %ID/g, respectively. Significant variability in nanoparticle distribution was observed in certain organs such as liver, spleen and lungs. A large fraction of this variability could be explained by accounting for the differences in nanoparticle physicochemical properties such as size and material. A critical overview of published nanoparticle physiologically-based pharmacokinetic (PBPK) models is provided, and limitations in our current knowledge about in vitro and in vivo pharmacokinetics of nanoparticles that restrict the development of robust PBPK models is also discussed. It is hypothesized that robust quantitative assessment of whole-body pharmacokinetics of nanoparticles and development of mathematical models that can predict their disposition can improve the probability of successful clinical translation of these modalities.

In vitro to in vivo acetaminophen hepatotoxicity extrapolation using classical schemes, pharmacodynamic models and a multiscale spatial-temporal liver twin

In vitro to in vivo extrapolation represents a critical challenge in toxicology. In this paper we explore extrapolation strategies for acetaminophen (APAP) based on mechanistic models, comparing classical (CL) homogeneous compartment pharmacodynamic (PD) models and a spatial-temporal (ST), multiscale digital twin model resolving liver microarchitecture at cellular resolution. The models integrate consensus detoxification reactions in each individual hepatocyte. We study the consequences of the two model types on the extrapolation and show in which cases these models perform better than the classical extrapolation strategy that is based either on the maximal drug concentration (Cmax) or the area under the pharmacokinetic curve (AUC) of the drug blood concentration. We find that an CL-model based on a well-mixed blood compartment is sufficient to correctly predict the in vivo toxicity from in vitro data. However, the ST-model that integrates more experimental information requires a change of at least one parameter to obtain the same prediction, indicating that spatial compartmentalization may indeed be an important factor.

An Overview of Physiologically-Based Pharmacokinetic Models for Forensic Science

A physiologically-based pharmacokinetic (PBPK) model represents the structural components of the body with physiologically relevant compartments connected via blood flow rates described by mathematical equations to determine drug disposition. PBPK models are used in the pharmaceutical sector for drug development, precision medicine, and the chemical industry to predict safe levels of exposure during the registration of chemical substances. However, one area of application where PBPK models have been scarcely used is forensic science. In this review, we give an overview of PBPK models successfully developed for several illicit drugs and environmental chemicals that could be applied for forensic interpretation, highlighting the gaps, uncertainties, and limitations.

Use of physiological based pharmacokinetic modeling for cross-species prediction of pharmacokinetic and tissue distribution profiles of a novel niclosamide prodrug

Introduction: Niclosamide (Nc) is an FDA-approved anthelmintic drug that was recently identified in a drug repurposing screening to possess antiviral activity against SARS-CoV-2. However, due to the low solubility and permeability of Nc, its in vivo efficacy was limited by its poor oral absorption. Method: The current study evaluated a novel prodrug of Nc (PDN; NCATS-SM4705) in improving in vivo exposure of Nc and predicted pharmacokinetic profiles of PDN and Nc across different species. ADME properties of the prodrug were determined in humans, hamsters, and mice, while the pharmacokinetics (PK) of PDN were obtained in mice and hamsters. Concentrations of PDN and Nc in plasma and tissue homogenates were measured by UPLC-MS/MS. A physiologically based pharmacokinetic (PBPK) model was developed based on physicochemical properties, pharmacokinetic and tissue distribution data in mice, validated by the PK profiles in hamsters and applied to predict pharmacokinetic profiles in humans. Results: Following intravenous and oral administration of PDN in mice, the total plasma clearance (CL_p) and volume of distribution at steady-state (Vd_ss) were 0.061-0.063 L/h and 0.28-0.31 L, respectively. PDN was converted to Nc in both liver and blood, improving the systemic exposure of Nc in mice and hamsters after oral administration. The PBPK model developed for PDN and in vivo formed Nc could adequately simulate plasma and tissue concentration-time profiles in mice and plasma profiles in hamsters. The predicted human CL_p/F and Vd_ss/F after an oral dose were 2.1 L/h/kg and 15 L/kg for the prodrug respectively. The predicted Nc concentrations in human plasma and lung suggest that a TID dose of 300 mg PDN would provide Nc lung concentrations at 8- to 60-fold higher than in vitro IC₅₀ against SARS-CoV-2 reported in cell assays. Conclusion: In conclusion, the novel prodrug PDN can be efficiently converted to Nc in vivo and improves the systemic exposure of Nc in mice after oral administration. The developed PBPK model adequately depicts the mouse and hamster pharmacokinetic and tissue distribution profiles and highlights its potential application in the prediction of human pharmacokinetic profiles.