Pharmacokinetic Modeling Approach Research Articles

Predict the Drug-Drug Interaction of a Novel PI3Kα/δ Inhibitor, TQ-B3525, and Its Two Metabolites Using Physiologically Based Pharmacokinetic Modeling.

A novel dual PI3K α/δ inhibitor, TQ-B3525, has been developed for the targeted treatment of lymphoma and solid tumors. TQ-B3525 is primarily metabolized by CYP3A4 and FOM3, while also serving as a substrate for the P-glycoprotein transporter. The aim of this study was to anticipate the drug-drug interaction (DDI) of TQ-B3525 and its two metabolites with CYP3A4 enzyme potent inducer (rifampicin) and CYP3A4/P-gp inhibitor (itraconazole) utilizing a physiologically based pharmacokinetic (PBPK) modeling approach. Clinical data from healthy and cancer patient adults were employed to construct and evaluate the PBPK model for TQ-B3525, M3, and M8-3. Models involving rifampicin combined with midazolam, itraconazole combined with midazolam or digoxin were utilized to showcase the robustness of evaluating DDI effects. The simulated drug exposure of TQ-B3525, M3, and M8-3 in healthy and patient adults were consistent with clinical data, and the mean fold error values were within the acceptable ranges. The simulated results of positive substrates correspond to those reported in the literature. Co-administration with rifampicin reduces Cmax and AUC of TQ-B3525 to 76.1% and 46.0%, while increasing the levels of M3 and M8-3. With itraconazole, Cmax and AUC of TQ-B3525 rise to 131% and 204%, but decrease substantially for M3 and M8-3. PBPK model simulation results showed that the systemic exposure of TQ-B3525 was significantly affected when co-administered with CYP3A4/P-gp inducers and inhibitors. This indicates that the combination with strong inducers and inhibitors should be carefully avoided or adjust the dosage of TQ-B3525 in clinic.

A physiologically-based pharmacokinetic modeling approach for dosing amiodarone in children on ECMO.

Extracorporeal membrane oxygenation (ECMO) is a cardiopulmonary bypass device commonly used to treat cardiac arrest in children. The American Heart Association guidelines for cardiopulmonary resuscitation (CPR) and emergency cardiovascular care recommend using amiodarone as a first-line agent to treat ventricular arrhythmias in children with cardiac arrest. However, there are no dosing recommendations for amiodarone to treat ventricular arrhythmias in pediatric patients on ECMO. Amiodarone has a high propensity for adsorption to the ECMO components due to its physicochemical properties leading to altered pharmacokinetics (PK) in ECMO patients. The change in amiodarone PK due to interaction with ECMO components may result in a difference in optimal dosing in patients on ECMO when compared with non-ECMO patients. To address this clinical knowledge gap, a physiologically-based pharmacokinetic model of amiodarone was developed in adults and scaled to children, followed by the addition of an ECMO compartment. The pediatric model included ontogeny functions of cytochrome P450 (CYP450) enzyme maturation across various age groups. The ECMO compartment was parameterized using the adsorption data of amiodarone obtained from exvivo studies. Model predictions captured observed concentrations of amiodarone in pediatric patients with ECMO well with an average fold error between 0.5 and 2. Model simulations support an amiodarone intravenous (i.v) bolus dose of 22 mg/kg (neonates), 13 mg/kg (infants), 8 mg/kg (children), and 6 mg/kg (adolescents). This PBPK modeling approach can be applied to explore the dosing of other drugs used in children on ECMO.

Model-informed drug development of envafolimab, a subcutaneously injectable PD-L1 antibody, in patients with advanced solid tumors.

Envafolimab is the first and only globally approved subcutaneously injectable PD-L1 antibody for the treatment of instability-high (MSI-H) or DNA mismatch repair deficient (dMMR) advanced solid tumors in adults, including those with advanced colorectal cancer that has progressed after treatment with a fluoropyrimidine, oxaliplatin, and irinotecan. The aim of this investigation was to examine the pharmacokinetic and exposure-response (E-R) profile of envafolimab in patients with solid tumors to support the approval of fixed and alternative dose regimens. In this study, a population pharmacokinetic (PopPK) modeling approach will be employed to quantitatively evaluate intrinsic and extrinsic covariates. Additionally, PopPK-estimated exposure parameters were used to evaluate E-R relationship for safety and efficacy to provide a theoretical basis for recommending optimal treatment regimens. Simulations were performed on the dosing regimens of body weight-based regimen of 2.50mg/kg QW, fixed dose 150mg QW, and 300mg Q2W for the selection of alternative dosing regimens. Data from 4 clinical studies (NCT02827968, NCT03101488, NCT03248843, and NCT03667170) were utilized. The PopPK dataset comprised 182 patients with 1810 evaluable envafolimab concentration records. Finally, a one-compartment model incorporating first-order absorption, first-order linear elimination, and time-dependent elimination according to an Emax function was found to accurately describe the concentration-time data of envafolimab in patients with advanced solid tumors. Creatinine clearance and country were identified as statistically significant factors affecting clearance, but had limited clinical significance. A relative flat exposure-response relationship was observed between early measures of safety and efficacy to verify that no dose adjustment is required. Simulation results indicated that 2.50mg/kg QW, 150mg QW, and 300mg Q2W regimen yield similar steady-state exposure. No statistically significant difference was observed between weight-based and fixed dose regimens. Model-based simulation supports the adoption of a 150mg weekly or 300mg biweekly dosing regimen of envafolimab in the solid tumor population, as these schedules effectively balance survival benefits and safety risks.

Determining tissue distribution of the oral antileishmanial agent miltefosine: a physiologically-based pharmacokinetic modeling approach.

Miltefosine (MTS) is the only approved oral drug for treating leishmaniasis caused by intracellular Leishmania parasites that localize in macrophages of the liver, spleen, skin, bone marrow, and lymph nodes. MTS is extensively distributed in tissues and has prolonged elimination half-lives due to its high plasma protein binding, slow metabolic clearance, and minimal urinary excretion. Thus, understanding and predicting the tissue distribution of MTS help assess therapeutic and toxicologic outcomes of MTS, especially in special populations, e.g., pediatrics. In this study, a whole-body physiologically-based pharmacokinetic (PBPK) model of MTS was built on mice and extrapolated to rats and humans. MTS plasma and tissue concentration data obtained by intravenous and oral administration to mice were fitted simultaneously to estimate model parameters. The resulting high tissue-to-plasma partition coefficient values corroborate extensive distribution in all major organs except the bone marrow. Sensitivity analysis suggests that plasma exposure is most susceptible to changes in fraction unbound in plasma. The murine oral-PBPK model was further validated by assessing overlay of simulations with plasma and tissue profiles obtained from an independent study. Subsequently, the murine PBPK model was extrapolated to rats and humans based on species-specific physiological and drug-related parameters, as well as allometrically scaled parameters. Fold errors for pharmacokinetic parameters were within acceptable range in both extrapolated models, except for a slight underprediction in the human plasma exposure. These animal and human PBPK models are expected to provide reliable estimates of MTS tissue distribution and assist dose regimen optimization in special populations.

Investigating clinical pharmacokinetics of brivaracetam by using a pharmacokinetic modeling approach

The development of technology and the processing speed of computing machines have facilitated the evaluation of advanced pharmacokinetic (PK) models, making modeling processes simple and faster. The present model aims to analyze the PK of brivaracetam (BRV) in healthy and diseased populations. A comprehensive literature review was conducted to incorporate the BRV plasma concentration data and its input parameters into PK-Sim software, leading to the creation of intravenous (IV) and oral models for both populations. The developed physiologically based pharmacokinetic (PBPK) model of BRV was then assessed using the visual predictive checks, mean observed/predicted ratios (Robs/pre), and average fold error for PK parameters including the maximum systemic concentration (Cmax), the area under the curve at time 0 to t (AUC0–∞), and drug clearance (CL). The PBPK model of BRV demonstrated that mean Robs/pre ratios of the PK parameters remained within the acceptable limits when assessed against a twofold error margin. Furthermore, model predictions were carried out to assess how AUC0–∞ is affected following the administration of BRV in individuals with varying degrees of liver cirrhosis, ranging from different child–pugh (CP) scores like A, B, and C. Moreover, dose adjustments were recommended by considering the variations in Cmax and CL in various kidney disease stages (mild to severe).

Rerouting cardiovascular management following gastric bypass surgery: Dose optimization of carvedilol using population‐based analysis

AimsA population‐based pharmacokinetic (PK) modeling approach (PopPK) was used to investigate the impact of Roux‐en‐Y gastric bypass (RYGB) on the PK of (R)‐ and (S)‐carvedilol. We aimed to optimize carvedilol dosing for these patients utilizing a pharmacokinetic/pharmacodynamic (PK/PD) link model.MethodsPopPK models were developed utilizing data from 52 subjects, including nonobese, obese, and post‐ RYGB patients who received rac‐ carvedilol orally. Covariate analysis included anthropometric and laboratory data, history of RYGB surgery, CYP2D6 and CYP3A4 in vivo activity, and relative intestinal abundance of major drug‐ metabolizing enzymes and transporters. A direct effect inhibitory Emax pharmacodynamic model was linked to the PK model of (S)‐ carvedilol to simulate the changes in exercise‐ induced heart rate.ResultsA 2‐compartmental model with linear elimination and parallel first‐order absorptions best described (S)‐carvedilol PK. RYGB led to a twofold reduction in relative oral bioavailability compared to nonoperated subjects, along with delayed absorption of both enantiomers. The intestinal ABCC2 mRNA expression increases the time to reach the maximum plasma concentration. The reduced exposure (AUC) of (S)‐carvedilol post‐RYGB corresponded to a 33% decrease in the predicted area under the effect curve (AUEC) for the 24‐hour β‐blocker response. Simulation results suggested that a 50‐mg daily dose in post‐RYGB patients achieved comparable AUC and AUEC to 25‐mg dose in nonoperated subjects.ConclusionIntegrated PK/PD modeling indicated that standard dosage regimens for nonoperated subjects do not provide equivalent β‐blocking activity in RYGB patients. This study highlights the importance of personalized dosing strategies to attain desired therapeutic outcomes in this patient cohort.

Evaluating gender effect in the generic bioequivalence studies by physiologically based pharmacokinetic modeling - A case study of dextromethorphan modified release tablets.

The United States Food and Drug Administration guidelines for the bioequivalence (BE) testing of the generic drug products suggests that there should be an equal proportion of male and female population in the BE study. Despite this requirement, many generic drug companies do not maintain the suggested proportion of female population in their studies. Several socio-economic and cultural factors lead to lower participation of the females in the BE studies. More recently, the regulatory agencies across the globe are requesting the generic drug companies to demonstrate the performance of their drug products in the under-represented sex via additional studies. In this work, we describe the case of Dextromethorphan modified release tablets where the gender effect on the product performance was evaluated by physiologically based pharmacokinetic (PBPK) modeling approach. We have compared the drug product's performance by population simulations considering four different scenarios. The data from all-male population (from in house Pharmacokinetic [PK] BE studies) was considered as a reference and other scenarios were compared against the all-male population data. In the first scenario, we made a comparison between all-male (100% male) vs all-female (100% female) population. Second scenario was as per agency's requirements-equal proportion of male and female in the BE study. As an extreme scenario, 100% male vs 30:70 male:female was considered (higher females than males in the BE studies). Finally, as a more realistic scenario, 100% male versus 70:30 male:female was considered (lower females than males in the BE studies). Population PK followed by virtual BE was employed to demonstrate the similarity/differences in the drug product performance between the sexes. This approach can be potentially utilized to seek BE study waivers thus saving cost and accelerating the entry of the generic products to the market.

Model-Informed Clinical Development of 6-Monthly Injection of Paliperidone Palmitate in Patients with Schizophrenia: Dosing Strategies Guided by Population Pharmacokinetic Modeling and Simulation (Part II).

Paliperidone palmitate 6-month (PP6M) intramuscular (IM) injection is the longest-acting treatment available for patients with schizophrenia. A population pharmacokinetic (popPK) modeling and simulation approach was deployed to inform dosing strategies for PP6M. The extensive analysis database included 15,932 paliperidone samples from 700 patients receiving gluteal paliperidone palmitate 3-month (PP3M) or PP6M injections in the double-blind phase of a phase-3 noninferiority study (NCT03345342). Exposure parameters for paliperidone appeared to increase dose-proportionally within each dosing schedule (PP3M/PP6M). The range of paliperidone exposures after IM administration of PP6M overlaps with that of corresponding doses of oral paliperidone extended release, PP 1-month (PP1M), and PP3M. Model-based simulations were performed to investigate paliperidone exposures in different PP6M dosing scenarios and relevant subpopulations. A dosing window of ≤ 2 weeks earlier and ≤ 3 weeks later than the target 6-month interval for maintenance treatment with PP6M dosing maintains paliperidone exposures at levels that are not expected to substantially impact its safety and efficacy. For missed-dose scenarios, tailored re-initiation regimens are proposed that should be applied before resuming PP6M maintenance dosing. Regarding subpopulations, PP6M 700 mg eq. is the highest dose recommended in mild renal-impairment patients; the paliperidone pharmacokinetics after PP6M administration is not affected by sex, body mass index, or age in a clinically meaningful way. Paliperidone concentration-time profiles after PP6M and PP3M dosing were adequately described by the popPK model. Model-based simulation results provide guidance for clinicians on initiating PP6M therapy, transitioning between paliperidone formulations, the dosing windows to use for maintenance dosing, and managing missed PP6M doses.

Impact of Obesity on Hepatic Drug Clearance: What are the Influential Variables?

Drug clearance in obese subjects varies widely among different drugs and across subjects with different severity of obesity. This study investigates correlations between plasma clearance (CLp) and drug- and patient-related characteristics in obese subjects, and evaluates the systematic accuracy of common weight-based dosing methods. A physiologically-based pharmacokinetic (PBPK) modeling approach that uses recent information on obesity-related changes in physiology was used to simulate CLp for a normal-weight subject (body mass index [BMI] = 20) and subjects with various severities of obesity (BMI 25–60) for hypothetical hepatically cleared drugs with a wide range of properties. Influential variables for CLp change were investigated. For each drug and obese subject, the exponent that yields perfect allometric scaling of CLp from normal-weight subjects was assessed. Among all variables, BMI and relative changes in enzyme activity resulting from obesity proved highly correlated with obesity-related CLp changes. Drugs bound to α1-acid glycoprotein (AAG) had lower CLp changes compared to drugs bound to human serum albumin (HSA). Lower extraction ratios (ER) corresponded to higher CLp changes compared to higher ER. The allometric exponent for perfect scaling ranged from -3.84 to 3.34 illustrating that none of the scaling methods performed well in all situations. While all three dosing methods are generally systematically accurate for drugs with unchanged or up to 50% increased enzyme activity in subjects with a BMI below 30 kg/m2, in any of the other cases, information on the different drug properties and severity of obesity is required to select an appropriate dosing method for individuals with obesity.Graphical

Optimization of the betamethasone and dexamethasone dosing regimen during pregnancy: a combined placenta perfusion and pregnancy physiologically based pharmacokinetic modeling approach

BackgroundAntenatal betamethasone and dexamethasone are prescribed to women who are at high risk of premature birth to prevent neonatal respiratory distress syndrome. The current treatment regimens, effective to prevent neonatal RDS, may be suboptimal. Recently, concerns have been raised regarding possible adverse long-term neurological outcomes due to high fetal drug exposures. Data from non-human primates and sheep suggest maintaining a fetal plasma concentration above 1 ng/mL for 48 hours to retain efficacy, while avoiding undesirable high fetal plasma levels. ObjectiveWe aimed to re-evaluate the current betamethasone and dexamethasone dosing strategies to assess estimated fetal exposure and provide new dosing proposals that meet the efficacy target but avoid excessive peak exposures. Study designA pregnancy physiologically-based pharmacokinetic model was used to predict fetal drug exposures. To allow prediction of the extent of betamethasone and dexamethasone exposure in the fetus, placenta perfusion experiments were conducted to determine placental transfer. Placental transfer rates were integrated in the PBPK model to predict fetal exposure and model performance was verified using published maternal and fetal PK data. The verified pregnancy physiologically-based pharmacokinetic models were then used to simulate alternative dosing regimens to establish a model-informed dose. ResultsEx vivo data showed that both drugs extensively cross the placenta. For betamethasone 15.7 ± 1.7% and for dexamethasone 14.4 ± 1.5% of the initial maternal perfusate concentration reached the fetal circulations at the end of the 3-hour perfusion period. Pregnancy physiologically-based pharmacokinetic models that include these ex vivo-derived placental transfer rates, accurately predicted maternal and fetal exposures resulting from current dosing regimens. The dose simulations suggest that for betamethasone intramuscular a dose reduction from 2 dosages 11.4 mg, 24 hours apart, to 4 dosages 1.425 mg, 12 hours apart would avoid excessive peak exposures and still meet the fetal response threshold. For dexamethasone, the dose may be reduced from four times 6 mg every 12 hours to 8 times 1.5 mg every 6 hours. ConclusionA combined placenta perfusion and pregnancy physiologically-based pharmacokinetic modeling approach adequately predicted both maternal and fetal drug exposures of two antenatal corticosteroids. Strikingly, our PBPK simulations suggest that drug doses might be reduced drastically to still meet earlier proposed efficacy targets and minimize peak exposures. We propose the provided model-informed dosing regimens are used to support further discussion on an updated antenatal corticosteroid scheme and design of clinical trials to confirm the effectiveness and safety of lower doses.

Physiologically-based pharmacokinetic modelling guided dose evaluations of nirmatrelvir/ritonavir in renal impairment for the management of COVID-19.

We aimed to address factors contributing to the pharmacokinetic changes of nirmatrelvir/ritonavir in renal impaired (RI) patients and recommend dosing adjustment via a physiologically-based pharmacokinetic (PBPK) modelling approach. A PBPK model of nirmatrelvir/ritonavir was developed via Simcyp® Simulator. Sensitivity analysis of the influence of hepatic CYP3A4 intrinsic clearance and abundance, as well as hepatic non-CYP3A4 metabolism (other human liver microsomes [HLM] CLint) was performed to evaluate the effects of RI on oral clearance of nirmatrelvir. Other HLM CLint, the most sensitive parameter, was adjusted, and the simulated plasma concentration profiles of nirmatrelvir in severe RI subjects were within the therapeutic index of 292-10 000 ng/mL for dosing regimens of loading doses of 300/100 mg followed by 150/100 mg or 75/100 mg twice daily of nirmatrelvir/ritonavir. Considering that nirmatrelvir is available as a 150 mg tablet, we recommend 300/100 mg followed by 150/100 mg twice daily as the dosing regimen to be investigated in severe RI.

Predictions of tissue concentrations of myclobutanil, oxyfluorfen, and pronamide in rat and human after oral exposures via GastroPlusTM physiologically based pharmacokinetic modelling

ABSTRACT Heritage agrochemicals like myclobutanil, oxyfluorfen, and pronamide, are extensively used in agriculture, with well-established studies on their animal toxicity. Yet, human toxicity assessment relies on conventional human risk assessment approaches including the utilization of animal-based ADME (Absorption, Distribution, Metabolism, and Excretion) data. In recent years, Physiologically Based Pharmacokinetic (PBPK) modelling approaches have played an increasing role in human risk assessment of many chemicals including agrochemicals. This study addresses the absence of PBPK-type data for myclobutanil, oxyfluorfen, and pronamide by generating in vitro data for key input PBPK parameters (Caco-2 permeability, rat plasma binding, rat blood to plasma ratio, and rat liver microsomal half-life), followed by generation of PBPK models for these three chemicals via the GastroPlusTM software. Incorporating these experimental input parameters into PBPK models, the prediction accuracy of plasma AUC (area under curve) was significantly improved. Validation against rat oral administration data demonstrated substantial enhancement. Steady-state plasma concentrations (Css) of pronamide aligned well with published data using measured PBPK parameters. Following validation, parent-based tissue concentrations for these agrochemicals were predicted in humans and rats after single or 30-day repeat exposure of 10 mg/kg/day. These predicted concentrations contribute valuable information for future human toxicity risk assessments of these agrochemicals.

Understanding adefovir pharmacokinetics as a component of a transporter phenotyping cocktail

PurposeAdefovir (as dipivoxil) was selected as a probe drug in a previous transporter cocktail phenotyping study to assess renal organic anion transporter 1 (OAT1), with renal clearance (CLR) as the primary parameter describing renal elimination. An approximately 20% higher systemic exposure of adefovir was observed when combined with other cocktail components (metformin, sitagliptin, pitavastatin, and digoxin) compared to sole administration. The present evaluation applied a population pharmacokinetic (popPK) modeling approach to describe adefovir pharmacokinetics as a cocktail component in more detail.MethodsData from 24 healthy subjects were reanalyzed. After establishing a base model, covariate effects, including the impact of co-administered drugs, were assessed using forward inclusion then backward elimination.ResultsA one-compartment model with first-order absorption (including lag time) and a combination of nonlinear renal and linear nonrenal elimination best described the data. A significantly higher apparent bioavailability (73.6% vs. 59.0%) and a lower apparent absorption rate constant (2.29 h−1 vs. 5.18 h−1) were identified in the combined period compared to the sole administration period, while no difference was seen in renal elimination. The population estimate for the Michaelis-Menten constant (Km) of the nonlinear renal elimination was 170 nmol/L, exceeding the observed range of adefovir plasma maximum concentration, while the maximum rate (Vmax) of nonlinear renal elimination was 2.40 µmol/h at the median absolute estimated glomerular filtration rate of 105 mL/min.ConclusionThe popPK modeling approach indicated that the co-administration primarily affected the apparent absorption and/or prodrug conversion of adefovir dipivoxil, resulting in the minor drug-drug interaction observed for adefovir as a victim. However, renal elimination remained unaffected. The high Km value suggests that assessing renal OAT1 activity by CLR has no relevant misspecification error with the cocktail doses used.

Simulated pharmacokinetics of inhaled caffeine and melatonin from existing products indicate the lack of dosimetric considerations

Numerous commercially available inhalable products claim to improve sleep-wake cycle-related target indications by delivering a wide variety of chemicals like caffeine and melatonin. The resulting exposure-responses from inhaling different doses are unknown and obtaining early understanding of resulting pharmacokinetics is beneficial. This study applied a physiologically based pharmacokinetic modeling approach to predict the inhalation pharmacokinetics of caffeine and melatonin for different target indications related to the sleep-wake cycle. The model predicted rapid systemic delivery of caffeine and melatonin based on airway regional deposition of inhaled aerosol. A low inhaled dose of 1 mg of caffeine resulted in a 72.3-times lower plasma maximal concentration and was predicted to not improve cognitive performance task outcomes compared to oral consumption of coffee containing 80 mg of caffeine. Conversely, 2-mg oral and inhaled doses of melatonin under recommended directions of use result in more than 25.1- and 645-times higher plasma concentrations compared to endogenous melatonin, respectively. The recommended doses for inhalation products for potential improvement in the target indications vary widely. Additional research is needed to evaluate the human pharmacokinetics, efficacy, and safety of inhaled products. Given the lack of assessments, inhaled caffeine and melatonin must be consumed with caution as the toxicological concerns are not known and could outweigh the potential beneficial effects.

Assessing the relative contribution of CYP3A-and P-gp-mediated pathways to the overall disposition and drug-drug interaction of dabigatran etexilate using a comprehensive mechanistic physiological-based pharmacokinetic model.

Dabigatran etexilate (DABE) is a clinical probe substrate for studying drug-drug interaction (DDI) through an intestinal P-glycoprotein (P-gp). A recent in vitro study, however, has suggested a potentially significant involvement of CYP3A-mediated oxidative metabolism of DABE and its intermediate monoester BIBR0951 in DDI following microdose administration of DABE. In this study, the relative significance of CYP3A- and P-gp-mediated pathways to the overall disposition of DABE has been explored using mechanistic physiologically based pharmacokinetic (PBPK) modeling approach. The developed PBPK model linked DABE with its 2 intermediate (BIBR0951 and BIBR1087) and active (dabigatran, DAB) metabolites, and with all relevant drug-specific properties known to date included. The model was successfully qualified against several datasets of DABE single/multiple dose pharmacokinetics and DDIs with CYP3A/P-gp inhibitors. Simulations using the qualified model supported that the intestinal CYP3A-mediated oxidation of BIBR0951, and not the gut P-gp-mediated efflux of DABE, was a key contributing factor to an observed difference in the DDI magnitude following the micro-versus therapeutic doses of DABE with clarithromycin. Both the saturable CYP3A-mediated metabolism of BIBR0951 and the solubility-limited DABE absorption contributed to the relatively modest nonlinearity in DAB exposure observed with increasing doses of DABE. Furthermore, the results suggested a limited role of the gut P-gp, but an appreciable, albeit small, contribution of gut CYP3A in mediating the DDIs following the therapeutic dose of DABE with dual CYP3A/P-gp inhibitors. Thus, a possibility exists for a varying extent of CYP3A involvement when using DABE as a clinical probe in the DDI assessment, across DABE dose levels.

Physiologically Based Pharmacokinetic (PBPK) Model Predictions of Disease Mediated Changes in Drug Disposition in Patients with Nonalcoholic Fatty Liver Disease (NAFLD).

This study was designed to verify a virtual population representing patients with nonalcoholic fatty liver disease (NAFLD) to support the implementation of a physiologically based pharmacokinetic (PBPK) modeling approach for prediction of disease-related changes in drug pharmacokinetics. A virtual NAFLD patient population was developed in GastroPlus (v.9.8.2) by accounting for pathophysiological changes associated with the disease and proteomics-informed alterations in the abundance of metabolizing enzymes and transporters pertinent to drug disposition. The NAFLD population model was verified using exemplar drugs where elimination is influenced predominantly by cytochrome P450 (CYP) enzymes (chlorzoxazone, caffeine, midazolam, pioglitazone) or by transporters (rosuvastatin, 11C-metformin, morphine and the glucuronide metabolite of morphine). PBPK model predictions of plasma concentrations of all the selected drugs and hepatic radioactivity levels of 11C-metformin were consistent with the clinically-observed data. Importantly, the PBPK simulations using the virtual NAFLD population model provided reliable estimates of the extent of changes in key pharmacokinetic parameters for the exemplar drugs, with mean predicted ratios (NAFLD patients divided by healthy individuals) within 0.80- to 1.25-fold of the clinically-reported values, except for midazolam (prediction-fold difference of 0.72). A virtual NAFLD population model within the PBPK framework was successfully developed with good predictive capability of estimating disease-related changes in drug pharmacokinetics. This supports the use of a PBPK modeling approach for prediction of the pharmacokinetics of new investigational or repurposed drugs in patients with NAFLD and may help inform dose adjustments for drugs commonly used to treat comorbidities in this patient population.

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.

Determining the optimal pharmacokinetic modelling and simplified quantification method of [18F]AlF-P16-093 for patients with primary prostate cancer (PPCa).

This paper discusses the optimization of pharmacokinetic modelling and alternate simplified quantification method for [18F]AlF-P16-093, a novel tracer for in vivo imaging of prostate cancer. Dynamic PET/CT scans were conducted on eight primary prostate cancer patients, followed by a whole-body scan at 60 min post-injection. Time-activity curves (TACs) were obtained by drawing volumes of interest for primary prostatic and metastatic lesions. Optimal kinetic modelling involved evaluating three compartmental models (1T2K, 2T3K, and 2T4K) accounting for fractional blood volume (Vb). The simplified quantification method was then determined based on the correlation between the static uptake measure and total distribution volume (Vt) obtained from the optimal pharmacokinetic analysis. In total, 17 intraprostatic lesions, 10 lymph nodes, and 36 osseous metastases were evaluated. Visually, the contrast of the tumor increased and showed the steepest incline within the first few minutes, whereas background activity decreased over time. Full pharmacokinetic analysis revealed that a reversible two-compartmental (2T4K) model is the preferred kinetic model for the given tracer. The kinetic parameters K1, k3, Vb, and Vt were all significantly higher in lesions when compared with normal tissue (P < 0.01). Several simplified protocols were tested for approximating comprehensive dynamic quantification in tumors, with image-based SURmean (the ratio of tumor SUVmean to blood SUVmean) within the 28-34 min window found to be sufficient for approximating the total distribution Vt values (R2 = 0.949, P < 0.01). Both Vt and SURmean correlated significantly with the total serum prostate-specific antigen (tPSA) levels (P < 0.01). This study introduced an optimized pharmacokinetic modelling approach and a simplified acquisition method for [18F]AlF-P16-093, a novel PSMA-targeted radioligand, highlighting the feasibility of utilizing one static PET imaging (between 30 and 60 min) for the diagnosis of prostate cancer. Note that the image-derived input function in this study may not reflect the true corrected plasma input function, therefore the interpretation of the associated kinetic parameter estimates should be done with caution.

Advances and Challenges in Modeling Cannabidiol Pharmacokinetics and Hepatotoxicity.

Cannabidiol (CBD) is a pharmacologically active metabolite of cannabis that is US Food and Drug Administration approved to treat seizures associated with Lennox-Gastaut syndrome, Dravet syndrome, and tuberous sclerosis complex in children aged 1 year and older. During clinical trials, CBD caused dose-dependent hepatocellular toxicity at therapeutic doses. The risk for toxicity was increased in patients taking valproate, another hepatotoxic antiepileptic drug, through an unknown mechanism. With the growing popularity of CBD in the consumer market, an improved understanding of the safety risks associated with CBD is needed to ensure public health. This review details current efforts to describe CBD pharmacokinetics and mechanisms of hepatotoxicity using both pharmacokinetic models and in vitro models of the liver. In addition, current evidence and knowledge gaps related to intracellular mechanisms of CBD-induced hepatotoxicity are described. The authors propose future directions that combine systems-based models with markers of CBD-induced hepatotoxicity to understand how CBD pharmacokinetics may influence the adverse effect profile and risk of liver injury for those taking CBD. SIGNIFICANCE STATEMENT: This review describes current pharmacokinetic modeling approaches to capture the metabolic clearance and safety profile of cannabidiol (CBD). CBD is an increasingly popular natural product and US Food and Drug Administration-approved antiepileptic drug known to cause clinically significant enzyme-mediated drug interactions and hepatotoxicity at therapeutic doses. CBD metabolism, pharmacokinetics, and putative mechanisms of CBD-induced liver injury are summarized from available preclinical data to inform future modeling efforts for understanding CBD toxicity.