Physiologically Based Pharmacokinetic Research Articles

Utility of PBPK Modeling in Predicting and Characterizing Clinical Drug Interactions.

Physiologically-based pharmacokinetic (PBPK) modeling is a mechanistic dynamic modeling approach that can be used to predict or retrospectively describe changes in drug exposure due to drug-drug interactions. With advancements in commercially available PBPK software, PBPK DDI modeling has become a mainstream approach from early drug discovery through to late stage drug development and is often utilized to support regulatory packages for new drug applications. This minireview will briefly describe the approaches to predicting DDI utilizing PBPK and static modeling approaches, the basic model structures and features inherent to PBPK DDI models and key examples where PBPK DDI models have been used to describe complex DDI mechanisms. Future directions aimed at using PBPK models to characterize transporter-mediated DDI, predict DDI in special populations and assess the DDI potential of protein therapeutics will be discussed. A summary of the 209 PBPK DDI examples published to date in 2023 will be provided. Overall, current data and trends suggest a continued role for PBPK models in the characterization and prediction of DDI for therapeutic molecules. Significance Statement PBPK models have been a key tool in the characterization of various pharmacokinetic phenomenon, including drug-drug interactions. This minireview will highlight recent advancements and publications around PBPK DDI modeling, an important area of drug discovery and development research in light of the increasing prevalence of polypharmacology in clinical settings.

Mechanistic Modeling of Intrathecal Chemotherapy Pharmacokinetics in the Human Central Nervous System.

The pharmacokinetics of intrathecally administered antibody or small-molecule drugs in the human central nervous system (CNS) remains poorly understood. This study aimed to provide mechanistic and quantitative perspectives on the CNS pharmacokinetics of intrathecal chemotherapy, by using a physiologically based pharmacokinetic (PBPK) modeling approach. A novel CNS PBPK model platform was developed and verified, which accounted for the human CNS general anatomy and physiologic processes governing drug distribution and disposition. The model was used to predict CNS pharmacokinetics of antibody (trastuzumab) and small-molecule drugs (methotrexate, abemaciclib, tucatinib) following intraventricular injection or intraventricular 24-hour infusion, and to assess the key determinants of drug penetration into the deep brain parenchyma. Intraventricularly administered antibody and small-molecule drugs exhibited distinct temporal and spatial distribution and disposition in human CNS. Both antibody and small-molecule drugs achieved supratherapeutic or therapeutic concentrations in the cerebrospinal fluid (CSF) compartments and adjacent brain tissue. While intrathecal small-molecule drugs penetrated the deep brain parenchyma to a negligible extent, intrathecal antibodies may achieve therapeutic concentrations in the deep brain parenchyma. Intraventricular 24-hour infusion enabled prolonged CNS exposure to therapeutically relevant concentrations while avoiding excessively high and potentially neurotoxic drug concentrations. CNS PBPK modeling, in line with available clinical efficacy data, confirms the therapeutic value of intrathecal chemotherapy with antibody or small-molecule drugs for treating neoplastic meningitis and warrants further clinical investigation of intrathecal antibody drugs to treat brain parenchyma tumors. Compared with intraventricular injection, intraventricular 24-hour infusion may mitigate neurotoxicity while retaining potential efficacy.

Evaluation of the drug-drug interaction potential of brigatinib using a physiologically-based pharmacokinetic modeling approach.

Brigatinib is an oral anaplastic lymphoma kinase (ALK) inhibitor approved for the treatment of ALK-positive metastatic non-small cell lung cancer. Invitro studies indicated that brigatinib is primarily metabolized by CYP2C8 and CYP3A4 and inhibits P-gp, BCRP, OCT1, MATE1, and MATE2K. Clinical drug-drug interaction (DDI) studies with the strong CYP3A inhibitor itraconazole or the strong CYP3A inducer rifampin demonstrated that CYP3A-mediated metabolism was the primary contributor to overall brigatinib clearance in humans. A physiologically-based pharmacokinetic (PBPK) model for brigatinib was developed to predict potential DDIs, including the effect of moderate CYP3A inhibitors or inducers on brigatinib pharmacokinetics (PK) and the effect of brigatinib on the PK of transporter substrates. The developed model was able to predict clinical DDIs with itraconazole (area under the plasma concentration-time curve from time 0 to infinity [AUC∞] ratio [with/without itraconazole]: predicted 1.86; observed 2.01) and rifampin (AUC∞ ratio [with/without rifampin]: predicted 0.16; observed 0.20). Simulations using the developed model predicted that moderate CYP3A inhibitors (e.g., verapamil and diltiazem) may increase brigatinib AUC∞ by ~40%, whereas moderate CYP3A inducers (e.g., efavirenz) may decrease brigatinib AUC∞ by ~50%. Simulations of potential transporter-mediated DDIs predicted that brigatinib may increase systemic exposures (AUC∞) of P-gp substrates (e.g., digoxin and dabigatran) by 15%-43% and MATE1 substrates (e.g., metformin) by up to 29%; however, negligible effects were predicted on BCRP-mediated efflux and OCT1-mediated uptake. The PBPK analysis results informed dosing recommendations for patients receiving moderate CYP3A inhibitors (40% brigatinib dose reduction) or inducers (up to 100% increase in brigatinib dose) during treatment, as reflected in the brigatinib prescribing information.

Applications of Advanced Natural Language Processing for Clinical Pharmacology.

Natural language processing (NLP) is a branch of artificial intelligence, which combines computational linguistics, machine learning, and deep learning models to process human language. Although there is a surge in NLP usage across various industries in recent years, NLP has not been widely evaluated and utilized to support drug development. To demonstrate how advanced NLP can expedite the extraction and analyses of information to help address clinical pharmacology questions, inform clinical trial designs, and support drug development, three use cases are described in this article: (1) dose optimization strategy in oncology, (2) common covariates on pharmacokinetic (PK) parameters in oncology, and (3) physiologically-based PK (PBPK) analyses for regulatory review and product label. The NLP workflow includes (1) preparation of source files, (2) NLP model building, and (3) automation of data extraction. The Clinical Pharmacology and Biopharmaceutics Summary Basis of Approval (SBA) documents, US package inserts (USPI), and approval letters from the US Food and Drug Administration (FDA) were used as our source data. As demonstrated in the three example use cases, advanced NLP can expedite the extraction and analyses of large amounts of information from regulatory review documents to help address important clinical pharmacology questions. Although this has not been adopted widely, integrating advanced NLP into the clinical pharmacology workflow can increase efficiency in extracting impactful information to advance drug development.

Personalized Chronomodulated 5-Fluorouracil Treatment: A Physiologically-Based Pharmacokinetic Precision Dosing Approach for Optimizing Cancer Therapy.

The discovery of circadian clock genes greatly amplified the study of diurnal variations impacting cancer therapy, transforming it into a rapidly growing field of research. Especially, use of chronomodulated treatment with 5-fluorouracil (5-FU) has gained significance. Studies indicate high interindividual variability (IIV) in diurnal variations in dihydropyrimidine dehydrogenase (DPD) activity - a key enzyme for 5-FU metabolism. However, the influence of individual DPD chronotypes on chronomodulated therapy remains unclear and warrants further investigation. To optimize precision dosing of chronomodulated 5-FU, this study aims to: (i) build physiologically-based pharmacokinetic (PBPK) models for 5-FU, uracil, and their metabolites, (ii) assess the impact of diurnal variation on DPD activity, (iii) estimate individual DPD chronotypes, and (iv) personalize chronomodulated 5-FU infusion rates based on a patient's DPD chronotype. Whole-body PBPK models were developed with PK-Sim(R) and MoBi(R). Sinusoidal functions were used to incorporate variations in enzyme activity and chronomodulated infusion rates as well as to estimate individual DPD chronotypes from DPYD mRNA expression or DPD enzymatic activity. Four whole-body PBPK models for 5-FU, uracil, and their metabolites were established utilizing data from 41 5-FU and 10 publicly available uracil studies. IIV in DPD chronotypes was assessed and personalized chronomodulated administrations were developed to achieve (i) comparable 5-FU peak plasma concentrations, (ii) comparable 5-FU exposure, and (iii) constant 5-FU plasma levels via "noise cancellation" chronomodulated infusion. The developed PBPK models capture the extent of diurnal variations in DPD activity and can help investigate individualized chronomodulated 5-FU therapy through testing alternative personalized dosing strategies.

Physiologically based pharmacokinetic modeling to predict the clinical effect of azole antifungal agents as CYP3A inhibitors on azelnidipine pharmacokinetics

In this study, a physiologically based pharmacokinetic (PBPK) model of the cytochrome P450 3A (CYP3A) substrate azelnidipine was developed using in vitro and clinical data to predict the effects of azole antifungals on azelnidipine pharmacokinetics. Modeling and simulations were conducted using the Simcyp™ PBPK simulator. The azelnidipine model consisted of a full PBPK model and a first-order absorption model. CYP3A was assumed as the only azelnidipine elimination route, and CYP3A clearance was optimized using the pharmacokinetic profile of single-dose 5-mg azelnidipine in healthy participants. The model reproduced the results of a clinical drug–drug interaction study and met validation criteria. PBPK model simulations using azole antifungals (itraconazole, voriconazole, posaconazole, fluconazole, fosfluconazole) and azelnidipine or midazolam (CYP3A index substrate) were performed. Increases in the simulated area under the plasma concentration–time curve from time zero extrapolated to infinity with inhibitors were comparable between azelnidipine (range, 2.11–6.47) and midazolam (range, 2.26–9.22), demonstrating that azelnidipine is a sensitive CYP3A substrate. Increased azelnidipine plasma concentrations are expected when co-administered with azole antifungals, potentially affecting azelnidipine safety. These findings support the avoidance of azole antifungals in patients taking azelnidipine and demonstrate the utility of PBPK modeling to inform appropriate drug use.

Physiologically based radiopharmacokinetic (PBRPK) modeling to simulate and analyze radiopharmaceutical therapies: studies of non-linearities, multi-bolus injections, and albumin binding

BackgroundWe aimed to develop a publicly shared computational physiologically based pharmacokinetic (PBPK) model to reliably simulate and analyze radiopharmaceutical therapies (RPTs), including probing of hot-cold ligand competitions as well as alternative injection scenarios and drug designs, towards optimal therapies.ResultsTo handle the complexity of PBPK models (over 150 differential equations), a scalable modeling notation called the “reaction graph” is introduced, enabling easy inclusion of various interactions. We refer to this as physiologically based radiopharmacokinetic (PBRPK) modeling, fine-tuned specifically for radiopharmaceuticals. As three important applications, we used our PBRPK model to (1) study the effect of competition between hot and cold species on delivered doses to tumors and organs at risk. In addition, (2) we evaluated an alternative paradigm of utilizing multi-bolus injections in RPTs instead of prevalent single injections. Finally, (3) we used PBRPK modeling to study the impact of varying albumin-binding affinities by ligands, and the implications for RPTs. We found that competition between labeled and unlabeled ligands can lead to non-linear relations between injected activity and the delivered dose to a particular organ, in the sense that doubling the injected activity does not necessarily result in a doubled dose delivered to a particular organ (a false intuition from external beam radiotherapy). In addition, we observed that fractionating injections can lead to a higher payload of dose delivery to organs, though not a differential dose delivery to the tumor. By contrast, we found out that increased albumin-binding affinities of the injected ligands can lead to such a differential effect in delivering more doses to tumors, and this can be attributed to several factors that PBRPK modeling allows us to probe.ConclusionsAdvanced computational PBRPK modeling enables simulation and analysis of a variety of intervention and drug design scenarios, towards more optimal delivery of RPTs.

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.

Prediction of combination effect of quinidine on the pharmacokinetics of tipepidine using a physiologically based pharmacokinetic model

Tipepidine, an antitussive drug, has been reported to have central pharmacological effects and can be expected to be safely repositioned as treatment for psychiatric disorders. Since tipepidine requires three doses per day, development of a once-daily medication would be highly beneficial. Previously, we reported that combination use with quinidine, a CYP2D6 inhibitor, prolongs the half-life of tipepidine in chimeric mice with humanised liver. In this study, to predict this combination effect in humans, a physiologically based pharmacokinetic (PBPK) model was developed, and quantitative simulation was conducted. The simulation results indicated that concomitant administration of tipepidine with quinidine increased the predicted C max, AUC, and t 1/2 of tipepidine in the Japanese population by 3.4-, 6.6-, and 2.4-fold, respectively. Furthermore, to compare with another approach that aims to prolong the half-life, the PK profile of tipepidine administered in hypothetical extended-release form was simulated. Extended-release form was predicted to be more influenced by CYP2D6 genotype than combination with quinidine, and the predicted plasma exposure was markedly increased in poor metabolizers, potentially leading to adverse effects. In conclusion, quantitative simulation using the PBPK model suggests the feasibility of the safe repositioning of tipepidine as a once-daily medication in combination with quinidine.

Effects of H2-Receptor Antagonists on the Exposure of Dacomitinib.

Dacomitinib is an irreversible epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor indicated for the treatment of patients with advanced non-small-cell lung cancer (NSCLC) and EGFR-activating mutations. Proton-pump inhibitors decreased dacomitinib exposure. This analysis summarizes the effect of Histamine-2 receptor antagonists (H2RAs) on dacomitinib exposure. A within-patient comparison of the steady-state trough concentrations (Ctrough,ss) of dacomitinib and its active metabolite and active moiety with and without concomitant use of H2RAs was conducted using a linear mixed effects model with pooled data from 11 clinical studies in patients with NSCLC. An oral absorption physiologically based pharmacokinetic (PBPK) model was constructed and verified using clinical pharmacokinetic (PK) data after a single dose of dacomitinib in healthy volunteers to estimate the effect of gastric pH altered by an H2RA on dacomitinib's PKs. The adjusted geometric mean of the dacomitinib Ctrough,ss of the dacomitinib parent, metabolite and active moiety following co-administration with an H2RA was approximately 86%, 104% and 100% relative to that following dacomitinib 45 mg administration without an H2RA (p > 0.05). The PBPK modeling showed negligible change in dacomitinib maximum concentration (Cmax) and area under the drug concentration-time curve (AUC) over 0-24 h after H2RA administration when compared with those administered dacomitinib alone. Co-administration of an H2RA with dacomitinib is not expected to have any clinically relevant effect on dacomitinib exposure.

Physiologically-based pharmacokinetic modeling of prominent oral contraceptive agents and applications in drug-drug interactions.

Considerable interest remains across the pharmaceutical industry and regulatory landscape in capabilities to model oral contraceptives (OCs), whether combined (COCs) with ethinyl estradiol (EE) or progestin-only pill. Acceptance of COC drug-drug interaction (DDI) assessment using physiologically-based pharmacokinetic (PBPK) is often limited to the estrogen component (EE), requiring further verification, with extrapolation from EE to progestins discouraged. There is a paucity of published progestin component PBPK models to support the regulatory DDI guidance for industry to evaluate a new chemical entity's (NCE's) DDI potential with COCs. Guidance recommends a clinical interaction study to be considered if an investigational drug is a weak or moderate inducer, or a moderate/strong inhibitor, of CYP3A4. Therefore, availability of validated OC PBPK models within one software platform, will be useful in predicting the DDI potential with NCEs earlier in the clinical development. Thus, this work was focused on developing and validating PBPK models for progestins, DNG, DRSP, LNG, and NET, within Simcyp, and assessing the DDI potential with known CYP3A4 inhibitors (e.g., ketoconazole) and inducers (e.g., rifampicin) with published clinical data. In addition, this work demonstrated confidence in the Simcyp EE model for regulatory and clinical applications by extensive verification in 70+ clinical PK and CYP3A4 interaction studies. The results provide greater capability to prospectively model clinical CYP3A4 DDI with COCs using Simcyp PBPK to interrogate the regulatory decision-tree to contextualize the potential interaction by known perpetrators and NCEs, enabling model-informed decision making, clinical study designs, and delivering potential alternative COC options for women of childbearing potential.

PBPK Modeling Approach to Predict the Behavior of Drugs Cleared by Metabolism in Pregnant Subjects and Fetuses.

This study aimed to develop a physiologically based pharmacokinetic (PBPK) model that simulates metabolically cleared compounds' pharmacokinetics (PK) in pregnant subjects and fetuses. This model accounts for the differences in tissue sizes, blood flow rates, enzyme expression levels, plasma protein binding, and other physiological factors affecting the drugs' PK in both the pregnant woman and the fetus. The PBPKPlus™ module in GastroPlus® was used to model the PK of metoprolol, midazolam, and metronidazole for both non-pregnant and pregnant groups. For each of the three compounds, the model was first developed and validated against PK data in healthy non-pregnant volunteers and then applied to predict the PK in the pregnant groups. The model accurately described the PK in both the non-pregnant and pregnant groups and explained well the differences in the plasma concentration due to pregnancy. When available, the fetal plasma concentration, placenta, and fetal tissue concentrations were also predicted reasonably well at different stages of pregnancy. The work described the use of a PBPK approach for drug development and demonstrates the ability to predict differences in PK in pregnant subjects and fetal exposure for metabolically cleared compounds.

Physiologically based pharmacokinetic modeling of ritonavir-oxycodone drug interactions and its implication for dosing strategy

The concomitant administration of ritonavir and oxycodone may significantly increase the plasma concentrations of oxycodone. This study was aimed to simulate DDI between ritonavir and oxycodone, a widely used opioid, and to formulate dosing protocols for oxycodone by using physiologically based pharmacokinetic (PBPK) modeling. We developed a ritonavir PBPK model incorporating induction and competitive and time-dependent inhibition of CYP3A4 and competitive inhibition of CYP2D6. The ritonavir model was evaluated with observed clinical pharmacokinetic data and validated for its CYP3A4 inhibition potency. We then used the model to simulate drug interactions between oxycodone and ritonavir under various dosing scenarios. The developed model captured the pharmacokinetic characteristics of ritonavir from clinical studies. The model also accurately predicts exposure changes of midazolam, triazolam, and oxycodone in the presence of ritonavir. According to model simulations, the steady-state maximum, minimum and average concentrations of oxycodone increased by up to 166% after co-administration with ritonavir, and the total exposure increased by approximately 120%. To achieve similar steady-state concentrations, halving the dose with an unchanged dosing interval or doubling the dosing interval with an unaltered single dose should be practical for oxycodone, whether formulated in uncoated or controlled-release tablets during long-term co-medication with ritonavir. The results revealed exposure-related risks of oxycodone-ritonavir interactions that have not been studied clinically and emphasized PBPK as a workable method to direct judicious dosage.

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

Assessment of Quinidine-Induced Torsades de Pointes Risks Using a Whole-Body Physiologically Based Pharmacokinetic Model Linked to Cardiac Ionic Current Inhibition.

The lethality of torsades de pointes (TdP) by drugs is one of main reasons that some drugs were withdrawn from the market. In order to assess drug-induced TdP risks, a model of cardiac ionic current suppression in human ventricular myocytes (ToR-ORd model), combined with the maximum effective free therapeutic plasma concentration or the maximum effective free therapeutic myocyte concentration was often used, with the latter proved to be more relevant and more accurate. We aimed to develop a whole-body physiologically-based pharmacokinetic (PBPK) model, incorporated with a human cardiomyocyte pharmacodynamic (PD) model, to provide a comprehensive assessment of drug-induced TdP risks in normal and specific scenarios. Quinidine served as an example to validate the PBPK-PD model via predicting plasma quinidine concentrations and quinidine-induced changes in QT interval (ΔQTc). The predicted plasma quinidine concentrations and ΔQTc values following oral administration or intravenous administration of quinidine were comparable to clinic observations. Visual predictive checks showed that most of the observed plasma concentrations and ΔQTc values fell within the 5th and 95th percentiles of simulations. The validated PBPK-PD model was further applied to assess the TdP risks using frequencies of early afterdepolarization and long-QT syndrome occurrence in 4 scenarios, such as therapeutic dose, supra-therapeutic dose, alkalosis, and hyperkalemia in 200 human subjects. In conclusion, the developed PBPK-PD model may be applied to predict the quinidine pharmacokinetics and quinidine-induced TdP risks in healthy subjects, but also simulate quinidine-induced TdP risks under disease conditions, such as hypokalemia and alkalosis.

Building a Predictive PBPK Model for Human OATP Substrates: a Strategic Framework for Early Evaluation of Clinical Pharmacokinetic Variations Using Pitavastatin as an Example.

To select a drug candidate for clinical development, accurately and promptly predicting human pharmacokinetic (PK) profiles, assessing drug-drug interactions (DDIs), and anticipating potential PK variations in disease populations are crucial steps in drug discovery. The complexity of predicting human PK significantly increases when hepatic transporters are involved in drug clearance (CL) and volume of distribution (Vss). A strategic framework is developed here, utilizing pitavastatin as an example. The framework includes the construction of a monkey physiologically-based PK (PBPK) model, model calibration to obtain scaling factors (SF) of in vitro-in vivo extrapolation (IVIVE) for various clearance parameters, human model development and validation, and assessment of DDIs and PK variations in disease populations. Through incorporating in vitro human parameters and calibrated SFs from the monkey model of 3.45, 0.14, and 1.17 for CLint,active, CLint,passive, and CLint,bile, respectively, and together with the relative fraction transported by individual transporters obtained from in vitro studies and the optimized Ki values for OATP inhibition, the model reasonably captured observed pitavastatin PK profiles, DDIs and PK variations in human subjects carrying genetic polymorphisms, i.e., AUC within 20%. Lastly, when applying the functional reduction based on measured OATP1B biomarkers, the model adequately predicted PK changes in the hepatic impairment population. The present study presents a strategic framework for early-stage drug development, enabling the prediction of PK profiles and assessment of PK variations in scenarios like DDIs, genetic polymorphism, and hepatic impairment-related disease states, specifically focusing on OATP substrates.

Local Identifiability Analysis, Parameter Subset Selection and Verification for a Minimal Brain PBPK Model.

Physiologically-based pharmacokinetic (PBPK) modeling is important for studying drug delivery in the central nervous system, including determining antibody exposure, predicting chemical concentrations at target locations, and ensuring accurate dosages. The complexity of PBPK models, involving many variables and parameters, requires a consideration of parameter identifiability; i.e., which parameters can be uniquely determined from data for a specified set of concentrations. We introduce the use of a local sensitivity-based parameter subset selection algorithm in the context of a minimal PBPK (mPBPK) model of the brain for antibody therapeutics. This algorithm is augmented by verification techniques, based on response distributions and energy statistics, to provide a systematic and robust technique to determine identifiable parameter subsets in a PBPK model across a specified time domain of interest. The accuracy of our approach is evaluated for three key concentrations in the mPBPK model for plasma, brain interstitial fluid and brain cerebrospinal fluid. The determination of accurate identifiable parameter subsets is important for model reduction and uncertainty quantification for PBPK models.

Development and application of neonatal physiology-based pharmacokinetic models of amikacin and fosfomycin to assess pharmacodynamic target attainment.

Antimicrobial resistance increasingly complicates neonatal sepsis in a global context. Fosfomycin and amikacin are two agents being tested in an ongoing multicenter neonatal sepsis trial. Although neonatal pharmacokinetics (PKs) have been described for these drugs, the physiological variability within neonatal populations makes population PKs in this group uncertain. Physiologically-based pharmacokinetic (PBPK) models were developed in Simcyp for fosfomycin and amikacin sequentially for adult, pediatric, and neonatal populations, with visual and quantitative validation compared to observed data at each stage. Simulations were performed using the final validated neonatal models to determine drug exposures for each drug across a demographic range, with probability of target attainment (PTA) assessments. Successfully validated neonatal PBPK models were developed for both fosfomycin and amikacin. PTA analysis demonstrated high probability of target attainment for amikacin 15 mg/kg i.v. q24h and fosfomycin 100 mg/kg (in neonates aged 0-7 days) or 150 mg/kg (in neonates aged 7-28 days) i.v. q12h for Enterobacterales with fosfomycin and amikacin minimum inhibitory concentrations at the adult breakpoints. Repeat analysis in premature populations demonstrated the same result. PTA analysis for a proposed combination fosfomycin-amikacin target was also performed. The simulated regimens, tested in a neonatal sepsis trial, are likely to be adequate for neonates across different postnatal ages and gestational age. This work demonstrates a template for determining target attainment for antimicrobials (alone or in combination) in special populations without sufficient available PK data to otherwise assess with traditional pharmacometric methods.

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.