Supplementing clinical lactation studies with PBPK modeling to inform drug therapy in lactating mothers: Prediction of primaquine exposure as a case example.

Evaluating the safety of primaquine (PQ) during breastfeeding requires an understanding of its pharmacokinetics (PKs) in breast milk and its exposure in the breastfed infant. Physiologically-based PK (PBPK) modeling is primed to assess the complex interplay of factors affecting the exposure of PQ in both the mother and the nursing infant. A published PBPK model for PQ describing the metabolism by monoamine oxidase A (MAO-A; 90% contribution) and cytochrome P450 2D6 (CYP2D6; 10%) in adults was applied to predict the exposure of PQ in mothers and their breastfeeding infants. Plasma exposures following oral daily dosing of 0.5 mg/kg in the nursing mothers in a clinical lactation study were accurately captured, including the observed ranges. Reported infant daily doses based on milk data from the clinical study were used to predict the exposure of PQ in breastfeeding infants greater than or equal to 28 days. On average, the predicted exposures were less than or equal to 0.13% of the mothers. Furthermore, in simulations involving neonates less than 28 days, PQ exposures remain less than 0.16% of the mothers. Assuming that MAO-A increases slowly with age, the predicted relative exposure of PQ remains low in neonates (<0.46%). Thus, the findings of our study support the recommendation made by the authors who reported the results of the clinical lactation study, that is, that when put into context of safety data currently available in children, PQ should not be withheld in lactating women as it is unlikely to cause adverse events in breastfeeding infants greater than or equal to 28 days old.

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

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

Safety of elexacaftor/tezacaftor/ivacaftor dose reduction: Mechanistic exploration through physiologically based pharmacokinetic modeling and a clinical case series.

INTRODUCTION: Elexacaftor/tezacaftor/ivacaftor (ETI) treatment is associated with significant improvement in lung function in people with cystic fibrosis (pwCF); however, some patients experience adverse effects (AEs) including hepatotoxicity. One potential strategy is dose reduction in ETI with the goal of maintaining therapeutic efficacy while resolving AEs. We report our experience of dose reduction in individuals who experienced AEs following ETI therapy. We provide mechanistic support for ETI dose reduction by exploring predicted lung exposures and underlying pharmacokinetics-pharmacodynamics (PK-PD) relationships. METHOD: Adults prescribed ETI who underwent dose reduction due to the AEs were included in this case series, and their percent predicted forced expiratory volume in 1 s (ppFEV1 ) and self-reported respiratory symptoms were collected. The full physiologically based pharmacokinetic (PBPK) models of ETI were developed incorporating physiological information and drug-dependent parameters. The models were validated against available pharmacokinetic and dose-response relationship data. The models were then used to predict lung concentrations of ETI at steady-state. RESULTS: Fifteen patients underwent dose reduction in ETI due to AEs. Clinical stability without significant changes in ppFEV1 after dose reduction was observed in all patients. Resolution or improvement of AEs occurred in 13 of the 15 cases. The model-predicted lung concentrations of reduced dose ETI exceeded the reported half maximal effective concentration (EC50 ) from measurement of in vitro chloride transport, providing a hypothesis as to why therapeutic efficacy was maintained. CONCLUSION: Albeit in a small number of patients, this study provides evidence that reduced ETI doses in pwCF who have experienced AEs may be effective. The PBPK models enable exploration of a mechanistic basis for this finding by simulating target tissue concentrations of ETI that can be compared with drug efficacy in vitro.

Physiologically-Based Pharmacokinetic-Led Guidance for Patients With Cystic Fibrosis Taking Elexacaftor-Tezacaftor-Ivacaftor With Nirmatrelvir-Ritonavir for the Treatment of COVID-19.

Cystic fibrosis transmembrane conductance regulator (CFTR) modulating therapies, including elexacaftor-tezacaftor-ivacaftor, are primarily eliminated through cytochrome P450 (CYP) 3A-mediated metabolism. This creates a therapeutic challenge to the treatment of coronavirus disease 2019 (COVID-19) with nirmatrelvir-ritonavir in people with cystic fibrosis (CF) due to the potential for significant drug-drug interactions (DDIs). However, the population with CF is more at risk of serious illness following COVID-19 infection and hence it is important to manage the DDI risk and provide treatment options. CYP3A-mediated DDI of elexacaftor-tezacaftor-ivacaftor was evaluated using a physiologically-based pharmacokinetic modeling approach. Modeling was performed incorporating physiological information and drug-dependent parameters of elexacaftor-tezacaftor-ivacaftor to predict the effect of ritonavir (the CYP3A inhibiting component of the combination) on the pharmacokinetics of elexacaftor-tezacaftor-ivacaftor. The elexacaftor-tezacaftor-ivacaftor models were verified using independent clinical pharmacokinetic and DDI data of elexacaftor-tezacaftor-ivacaftor with a range of CYP3A modulators. When ritonavir was administered on Days 1 through 5, the predicted area under the curve (AUC) ratio of ivacaftor (the most sensitive CYP3A substrate) on Day 6 was 9.31, indicating that its metabolism was strongly inhibited. Based on the predicted DDI, the dose of elexacaftor-tezacaftor-ivacaftor should be reduced when coadministered with nirmatrelvir-ritonavir to elexacaftor 200 mg-tezacaftor 100 mg-ivacaftor 150 mg on Days 1 and 5, with delayed resumption of full-dose elexacaftor-tezacaftor-ivacaftor on Day 9, considering the residual inhibitory effect of ritonavir as a mechanism-based inhibitor. The simulation predicts a regimen of elexacaftor-tezacaftor-ivacaftor administered concomitantly with nirmatrelvir-ritonavir in people with CF that will likely decrease the impact of the drug interaction.

Impact of Disease on Plasma and Lung Exposure of Chloroquine, Hydroxychloroquine and Azithromycin: Application of PBPK Modeling.

We use a mechanistic lung model to demonstrate that accumulation of chloroquine (CQ), hydroxychloroquine (HCQ), and azithromycin (AZ) in the lungs is sensitive to changes in lung pH, a parameter that can be affected in patients with coronavirus disease 2019 (COVID-19). A reduction in pH from 6.7 to 6 in the lungs, as observed in respiratory disease, led to 20-fold, 4.0-fold, and 2.7-fold increases in lung exposure of CQ, HCQ, and AZ, respectively. Simulations indicated that the relatively high concentrations of CQ and HCQ in lung tissue were sustained long after administration of the drugs had stopped. Patients with COVID-19 often present with kidney failure. Our simulations indicate that renal impairment (plus lung pH reduction) caused 30-fold, 8.0-fold, and 3.4-fold increases in lung exposures for CQ, HCQ, and AZ, respectively, with relatively small accompanying increases (20 to 30%) in systemic exposure. Although a number of different dosage regimens were assessed, the purpose of our study was not to provide recommendations for a dosing strategy, but to demonstrate the utility of a physiologically-based pharmacokinetic modeling approach to estimate lung concentrations. This, used in conjunction with robust in vitro and clinical data, can help in the assessment of COVID-19 therapeutics going forward.

Prediction of Drug-Drug Interactions Arising from CYP3A induction Using a Physiologically Based Dynamic Model.

Using physiologically based pharmacokinetic modeling, we predicted the magnitude of drug-drug interactions (DDIs) for studies with rifampicin and seven CYP3A4 probe substrates administered i.v. (10 studies) or orally (19 studies). The results showed a tendency to underpredict the DDI magnitude when the victim drug was administered orally. Possible sources of inaccuracy were investigated systematically to determine the most appropriate model refinement. When the maximal fold induction (Indmax) for rifampicin was increased (from 8 to 16) in both the liver and the gut, or when the Indmax was increased in the gut but not in liver, there was a decrease in bias and increased precision compared with the base model (Indmax = 8) [geometric mean fold error (GMFE) 2.12 vs. 1.48 and 1.77, respectively]. Induction parameters (mRNA and activity), determined for rifampicin, carbamazepine, phenytoin, and phenobarbital in hepatocytes from four donors, were then used to evaluate use of the refined rifampicin model for calibration. Calibration of mRNA and activity data for other inducers using the refined rifampicin model led to more accurate DDI predictions compared with the initial model (activity GMFE 1.49 vs. 1.68; mRNA GMFE 1.35 vs. 1.46), suggesting that robust in vivo reference values can be used to overcome interdonor and laboratory-to-laboratory variability. Use of uncalibrated data also performed well (GMFE 1.39 and 1.44 for activity and mRNA). As a result of experimental variability (i.e., in donors and protocols), it is prudent to fully characterize in vitro induction with prototypical inducers to give an understanding of how that particular system extrapolates to the in vivo situation when using an uncalibrated approach.

Incorporating population variability and susceptible subpopulations into dosimetry for high-throughput toxicity testing.

Momentum is growing worldwide to use in vitro high-throughput screening (HTS) to evaluate human health effects of chemicals. However, the integration of dosimetry into HTS assays and incorporation of population variability will be essential before its application in a risk assessment context. Previously, we employed in vitro hepatic metabolic clearance and plasma protein binding data with in vitro in vivo extrapolation (IVIVE) modeling to estimate oral equivalent doses, or daily oral chemical doses required to achieve steady-state blood concentrations (Css) equivalent to media concentrations having a defined effect in an in vitro HTS assay. In this study, hepatic clearance rates of selected ToxCast chemicals were measured in vitro for 13 cytochrome P450 and five uridine 5'-diphospho-glucuronysyltransferase isozymes using recombinantly expressed enzymes. The isozyme-specific clearance rates were then incorporated into an IVIVE model that captures known differences in isozyme expression across several life stages and ethnic populations. Comparison of the median Css for a healthy population against the median or the upper 95th percentile for more sensitive populations revealed differences of 1.3- to 4.3-fold or 3.1- to 13.1-fold, respectively. Such values may be used to derive chemical-specific human toxicokinetic adjustment factors. The IVIVE model was also used to estimate subpopulation-specific oral equivalent doses that were directly compared with subpopulation-specific exposure estimates. This study successfully combines isozyme and physiologic differences to quantitate subpopulation pharmacokinetic variability. Incorporation of these values with dosimetry and in vitro bioactivities provides a viable approach that could be employed within a high-throughput risk assessment framework.

Application of a systems approach to the bottom-up assessment of pharmacokinetics in obese patients: expected variations in clearance.

BACKGROUND AND OBJECTIVES: The maintenance dose of a drug is dependent on drug clearance, and thus any biochemical and physiological changes in obesity that affect parameters such as cardiac output, renal function, expression of drug-metabolizing enzymes and protein binding may result in altered clearance compared with that observed in normal-weight subjects (corrected or uncorrected for body weight). Because of the increasing worldwide incidence of obesity, there is a need for more information regarding the optimal dosing of drug therapy to be made available to prescribers. This is usually provided via clinical studies in obese people; however, such studies are not available for all drugs that might be used in obese subjects. Incorporation of the relevant physiological and biochemical changes into predictive bottom-up pharmacokinetic models in order to optimize dosage regimens may offer a logical way forward for the cases where no clinical data exist. The aims of the current report are to apply such a 'systems approach' to identify the likelihood of observing variations in the clearance of drugs in obesity and morbid obesity for a set of compounds for which clinical data, as well as the necessary in vitro information, are available, and to provide a framework for assessing other drugs in the future. METHODS: The population-specific changes in demographic, physiological and biochemical parameters that are known to be relevant to obese and morbidly obese subjects were collated and incorporated into two separate population libraries. These libraries, together with mechanistic in vitro-in vivo extrapolations (IVIVE) within the Simcyp Population-based Simulator™, were used to predict the clearance of oral alprazolam, oral caffeine, oral chlorzoxazone, oral ciclosporin, intravenous and oral midazolam, intravenous phenytoin, oral theophylline and oral triazolam. The design of the simulated studies was matched as closely as possible with that of the clinical studies. Outcome was measured by the predicted ratio of the clearance of the drug in obese and lean subjects ± its 90% confidence interval, compared with observed values. The overall statistical measures of the performance of the model to detect differences in compound clearance between obese and lean populations were investigated by measuring sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV). A power calculation was carried out to investigate the impact of the sample size on the overall outcome of clinical studies. RESULTS: The model was successful in predicting clearance in obese subjects, with the degree to which simulations could mimic the outcome of in vivo studies being greater than 60% for six of the eight drugs. A clear difference in the clearance of chlorzoxazone was correctly picked up via simulation. The overall statistical measures of the performance of the Simcyp Simulator were 100% sensitivity, 66% specificity, 60% PPV and 100% NPV. Studies designed on the basis of the ratio of the absolute values required substantial numbers of participants in order to detect a significant difference, except for phenytoin and chlorzoxazone, where the ratios of the weight-normalized clearances generally showed statistically significant differences with a smaller number of subjects. CONCLUSION: Extension of a mechanistic predictive pharmacokinetic model to accommodate physiological and biochemical changes associated with obesity and morbid obesity allowed prediction of changes in drug clearance on the basis of in vitro data, with reasonable accuracy across a range of compounds that are metabolized by different enzymes. Prediction of the effects of obesity on drug clearance, normalized by various body size scalars, is of potential value in the design of clinical studies during drug development and in the introduction of dosage adjustments that are likely to be needed in clinical practice.

Intracellular and plasma pharmacokinetics of efavirenz in HIV-infected individuals.

OBJECTIVES: The site of action of efavirenz is inside HIV-infected cells. Measurement of intracellular (IC) concentrations of efavirenz may therefore provide further understanding of therapeutic failure, especially where virological rebound occurs despite adequate plasma levels, and a lack of detectable viral resistance. Here, we determined IC and plasma pharmacokinetics of efavirenz and their relationship with plasma protein binding and P-glycoprotein (P-gp, an active drug efflux transporter) expression. PATIENTS AND METHODS: Venous blood samples from 10 HIV-infected patients receiving efavirenz (600 mg once a day plus two nucleoside reverse transcriptase inhibitors) were collected over the 24 h dosing interval. Plasma and peripheral blood mononuclear cells (PBMCs) were isolated. Plasma protein bound and unbound efavirenz were separated using ultrafiltration. IC (or cell-associated), total plasma and unbound plasma efavirenz levels were quantified using HPLC-UV. P-gp expression was measured by flow cytometry. Area under the concentration-time curves (AUC0-24) were then calculated using non-compartmental analyses and the IC accumulation expressed as a ratio of IC to plasma AUC0-24. RESULTS: The median (range) % unbound and IC accumulation ratio was 0.6% (0.4-1.5%) and 1.3 (0.7-3.3), respectively. There was a linear relationship between IC and total AUC0-24 (r2= 0.59, P = 0.01) but not unbound AUC0-24 (r2= 0.13, P = 0.75). An inverse correlation between IC AUC0-24 and % unbound was observed (r2= 41, P = 0.05). There was no relationship between IC AUC0-24 and P-gp expression on the cell surface (r2< 0.01, P = 0.98). CONCLUSIONS: There was a direct relationship between % bound efavirenz in plasma and IC accumulation implying that the IC accumulation of efavirenz is related to binding to IC proteins or other cellular constituents. Studies investigating the unbound concentration of antiretrovirals inside the cell are now required.

Intracellular and plasma pharmacokinetics of nevirapine in human immunodeficiency virus-infected individuals.

BACKGROUND AND OBJECTIVE: Plasma concentrations of nevirapine have been linked to human immunodeficiency virus (HIV) treatment outcome. However, because the site of action of nevirapine is within HIV-infected cells, intracellular concentrations may better relate to antiviral exposure. Investigation of factors that alter the intracellular pharmacokinetics of nevirapine may also aid in our understanding of therapeutic failure. Our objective was to determine intracellular (or cell-associated) nevirapine concentrations over the full dosing interval and to relate protein binding and P-glycoprotein (P-gp) expression to intracellular exposure. METHODS: Plasma and peripheral blood mononuclear cells were isolated from blood samples taken from 10 HIV-infected patients at 0, 2, 4, 8, and 12 hours after dosing. Intracellular and plasma (total and unbound) concentrations were determined by liquid chromatography-tandem mass spectrometry, and the ratios of intracellular to total plasma exposure (area under the concentration-time curves) were calculated. P-gp expression was measured by flow cytometry. RESULTS: The median intracellular accumulation ratio was 0.005 (range, 0.001-0.054) and remained unchanged over the dosing interval. There was an association between higher plasma concentrations and lower cellular concentrations of nevirapine (total r(2) = 0.62, P = .007). There was no relationship between percent unbound nevirapine and intracellular nevirapine. There was a correlation between higher plasma nevirapine exposure and higher P-gp expression (r(2) = 0.77, P = .03), whereas intracellular nevirapine exposure decreased with higher P-gp expression (r(2) = 0.62, P = .01). CONCLUSIONS: The intracellular accumulation of nevirapine was low, did not change over the dosing interval, and was not related to protein binding. In this small study, cells with higher P-gp expression had lower cellular concentrations of nevirapine. Further studies are required to explore the influx and efflux transporter profile of this drug.

The relationship between nevirapine plasma concentrations and abnormal liver function tests.

Abnormal liver function tests are frequently observed in HIV-infected individuals receiving nevirapine (NVP). Here we investigate the relationship between total and unbound plasma concentrations of NVP and the liver enzymes alanine aminotransferase (ALT) and gamma-glutamyl transferase (gammaGT). HIV-infected individuals [n = 85, 22 female, 34 hepatitis C or B virus (HCV or HBV(+))] receiving NVP (200 mg bd; median duration 66 weeks, range 3-189) and two nucleoside reverse transcriptase inhibitors (NRTIs) were enrolled into this study. Blood samples were taken at C(trough) (12 hr postdose) for measurement of NVP and liver function tests (ALT and gammaGT). Plasma protein bound and unbound drug was separated using ultrafiltration and NVP concentrations quantified using HPLC-MS/MS. A linear relationship was observed between total and unbound NVP C(trough) (r(2) = 0.77, p < 0.0001). Patients with elevated ALT (>37 IU/liter; n = 31) had higher NVP unbound C(trough) than those with ALT within the normal range (median 2268 vs. 1694 ng/ml, p = 0.04) but there was no difference in total C(trough). Logistic regression revealed no association between higher NVP C(trough) and ALT elevations. Significantly higher NVP total and unbound C(trough) were observed in patients with increased gammaGT (>40 IU/liter; n = 63; total 6747 vs. 4530 ng/ml, p = 0.001; unbound 2113 vs. 1557 ng/ml, p = 0.03). Significantly higher unbound NVP C(trough) was observed in HCV/HBV(+) (median 2275 vs. 1726 ng/ml, p = 0.02) and on bivariate analysis, higher NVP C(trough) was associated with HCV/HBV coinfection (chi(2) = 4.228; p = 0.04). Overall we found no strong association between NVP concentrations and hepatotoxicity. Although in this study NVP was well tolerated in HCV/HBV coinfected patients, higher plasma concentrations were observed.

Time-dependent changes in HIV nucleoside analogue phosphorylation and the effect of hydroxyurea.

BACKGROUND: Nucleoside analogues are activated to their triphosphates, which compete with endogenous deoxynucleoside triphosphate (dNTP) pools to inhibit HIV reverse transcriptase. Hydroxyurea has been administered with nucleoside analogues to modulate intracellular dNTP pools and thus the ratio of drug triphosphate:endogenous triphosphate. OBJECTIVES: To examine changes in drug activation over time and investigate the effects of hydroxyurea on intracellular phosphorylation of antiretroviral nucleoside analogues. PATIENTS: A total of 229 HIV-infected individuals receiving abacavir, lamivudine and zidovudine were randomly assigned to receive or not nevirapine and hydroxyurea. Twenty-four patients were recruited to an observational substudy measuring intracellular drug triphosphate and dNTP concentrations at 0, 2, 6, 12, 24 and 48 weeks. METHODS: Drugs were extracted from isolated peripheral blood mononuclear cells before analysis of endogenous dNTP and drug triphosphates by primer extension assays. RESULTS: Twenty-two out of 24 patients were followed to completion of the substudy. Hydroxyurea had no demonstrable effect on endogenous dNTP or drug triphosphate levels at any timepoint. However, the ratio of zidovudine triphosphate to endogenous deoxythymidine triphosphate was significantly increased with hydroxyurea. A significant decrease in lamivudine triphosphate (3TCTP) and the 3TCTP:endogenous deoxycytidine triphosphate ratio was seen over 48 weeks. In five patients who failed therapy in the first 24 weeks, significantly reduced 3TCTP was seen. CONCLUSION: Hydroxyurea does not affect measurable pools of endogenous nucleosides in vivo. Decreased lamivudine phosphorylation over time may provide a novel pharmacological explanation for the mechanism of resistance to this drug.