Venetoclax and Cobimetinib in Relapsed/Refractory AML: A Phase 1b Trial.

BACKGROUND: Therapies for relapsed/refractory acute myeloid leukemia remain limited and outcomes poor, especially amongst patients who are ineligible for cytotoxic chemotherapy or targeted therapies. PATIENTS AND METHODS: This phase 1b trial evaluated venetoclax, a B-cell lymphoma-2 (BCL-2) inhibitor, plus cobimetinib, a MEK1/2 inhibitor, in patients with relapsed/refractory acute myeloid leukemia, ineligible for cytotoxic chemotherapy. Two-dimensional dose-escalation was performed for venetoclax dosed daily, and for cobimetinib dosed on days 1-21 of each 28-day cycle. RESULTS: Thirty patients (median [range] age: 71.5 years [60-84]) received venetoclax-cobimetinib. The most common adverse events (AEs; in ≥40.0% of patients) were diarrhea (80.0%), nausea (60.0%), vomiting (40.0%), febrile neutropenia (40.0%), and fatigue (40.0%). Overall, 66.7% and 23.3% of patients experienced AEs leading to dose modification/interruption or treatment withdrawal, respectively. The composite complete remission (CRc) rate (complete remission [CR] + CR with incomplete blood count recovery + CR with incomplete platelet recovery) was 15.6%; antileukemic response rate (CRc + morphologic leukemia-free state/partial remission) was 18.8%. For the recommended phase 2 dose (venetoclax: 600 mg; cobimetinib: 40 mg), CRc and antileukemic response rates were both 12.5%. Failure to achieve an antileukemic response was associated with elevated baseline phosphorylated ERK and MCL-1 levels, but not BCL-xL. Baseline mutations in ≥1 signaling gene or TP53 were noted in nonresponders and emerged on treatment. Pharmacodynamic biomarkers revealed inconsistent, transient inhibition of the mitogen-activated protein kinase (MAPK) pathway. CONCLUSION: Venetoclax-cobimetinib showed limited preliminary efficacy similar to single-agent venetoclax, but with added toxicity. Our findings will inform future trials of BCL-2/MAPK pathway inhibitor combinations.

"Digital twins elucidate critical role of Tscm in clinical persistence of TCR-engineered cell therapy".

Despite recent progress in adoptive T cell therapy for cancer, understanding and predicting the kinetics of infused T cells remains a challenge. Multiple factors can impact the distribution, expansion, and decay or persistence of infused T cells in patients. We have developed a novel quantitative systems pharmacology (QSP) model of TCR-transgenic T cell therapy in patients with solid tumors to describe the kinetics of endogenous T cells and multiple memory subsets of engineered T cells after infusion. These T cells undergo lymphodepletion, proliferation, trafficking, differentiation, and apoptosis in blood, lymph nodes, tumor site, and other peripheral tissues. Using the model, we generated patient-matched digital twins that recapitulate the circulating T cell kinetics reported from a clinical trial of TCR-engineered T cells targeting E7 in patients with metastatic HPV-associated epithelial cancers. Analyses of key parameters influencing cell kinetics and differences among digital twins identify stem cell-like memory T cells (Tscm) cells as an important determinant of both expansion and persistence and suggest that Tscm-related differences contribute significantly to the observed variability in cellular kinetics among patients. We simulated in silico clinical trials using digital twins and predict that Tscm enrichment in the infused product improves persistence of the engineered T cells and could enable administration of a lower dose. Finally, we verified the broader relevance of the QSP model, the digital twins, and findings on the importance of Tscm enrichment by predicting kinetics for two patients with pancreatic cancer treated with KRAS G12D targeting T cell therapy. This work offers insight into the key role of Tscm biology on T cell kinetics and provides a quantitative framework to evaluate cellular kinetics for future efforts in the development and clinical application of TCR-engineered T cell therapies.

Venetoclax in biomarker-selected multiple myeloma patients: Impact of exposure on clinical efficacy and safety.

Venetoclax, a potent BCL-2 inhibitor, is currently under development for treatment of t(11;14) Multiple myeloma (MM). The objective of this research was to investigate the exposure-response relationships of venetoclax for a phase 1/2 study evaluating venetoclax monotherapy or in combination with dexamethasone in relapsed or refractory MM. A total of 117 patients receiving venetoclax at 300, 600, 800, 900, or 1200 mg were included in the analysis. The impact of venetoclax exposures on efficacy (objective response rate [ORR], progression-free survival [PFS] and overall survival [OS]) as well as safety (treatment-emergent adverse effects (grade ≥3) of neutropenia, infection, and any grade of serious treatment-emergent adverse effects) was evaluated. In the t(11;14)-positive subpopulation, venetoclax exposure relationships to PFS and OS indicated a trend of longer PFS and OS with higher exposures. Moreover, logistic regression analyses for clinical response (ORR and ≥VGPR rate) demonstrated a statistically significant (p < 0.05) relationship with exposure. Evaluation of the exposure-safety relationships demonstrated a lack of a relationship between venetoclax exposures (AUCavg ) and grade ≥3 infections, grade ≥3 neutropenia, grade ≥3 treatment-emergent adverse events or any grade serious treatment-emergent adverse events. These findings support further study of venetoclax at 800 mg QD dose in combination with dexamethasone in the t(11;14)-positive patient population where increased efficacy was observed without an increase in safety events.Clinical Trial: NCT01794520 registered 20 February 2013.

Venetoclax pharmacokinetics in participants with end-stage renal disease undergoing haemodialysis.

AIMS: Renal insufficiency is a common comorbidity in patients with haematological malignancies. This study aimed to assess how end-stage renal disease (ESRD) might affect the pharmacokinetics of venetoclax, a Bcl-2 inhibitor, in participants with ESRD undergoing haemodialysis. METHODS: Venetoclax was administered as a single 100-mg dose to 6 female participants with ESRD (estimated glomerular filtration rate <15 mL/min) both prior to haemodialysis and between haemodialysis days and 7 healthy female participants with normal renal function (estimated glomerular filtration rate >90 mL/min). Intensive pharmacokinetic and protein binding samples were collected from all participants. Arterial and venous samples were collected from ESRD participants during haemodialysis to assess the effect of haemodialysis on venetoclax pharmacokinetics. Pharmacokinetic parameters were estimated using noncompartmental methods. RESULTS: There was no difference in plasma venetoclax concentrations between arterial and venous samples, suggesting that haemodialysis did not affect the pharmacokinetics of venetoclax. The fraction unbound (fu ) of venetoclax was ~2-fold higher for participants with ESRD compared to participants with normal renal function. The unbound maximum plasma concentration and area under the plasma concentration-time curve from time 0 to 48 h were comparable between ESRD and normal function groups. The mean half-life ranged from 10.4 to 12.2 h across groups, demonstrating that ESRD did not affect the half-life of venetoclax. No new safety signals were observed during this study. CONCLUSION: ESRD and dialysis do not alter unbound venetoclax plasma concentrations. No pharmacokinetics driven dose adjustment is needed for patients with renal insufficiency.

Effect of isotretinoin on CYP2D6 and CYP3A activity in patients with severe acne.

AIMS: Previously, retinoids have decreased CYP2D6 mRNA expression in vitro and induced CYP3A4 in vitro and in vivo. This study aimed to determine whether isotretinoin administration changes CYP2D6 and CYP3A activities in patients with severe acne. METHODS: Thirty-three patients (22 females and 11 males, 23.5 ± 6.0 years old) expected to receive isotretinoin treatment completed the study. All participants were genotyped for CYP2D6 and CYP3A5. Participants received dextromethorphan (DM) 30 mg orally as a dual-probe substrate of CYP2D6 and CYP3A activity at two study timepoints: pre-isotretinoin treatment and with isotretinoin for at least 1 week. The concentrations of isotretinoin, DM and their metabolites were measured in 2-h postdose plasma samples and in cumulative 0-4-h urine collections using liquid chromatography-mass spectrometry. RESULTS: In CYP2D6 extensive metabolizers, the urinary dextrorphan (DX)/DM metabolic ratio (MR) (CYP2D6 activity marker) was numerically, but not significantly, lower with isotretinoin administration compared to pre-isotretinoin (geometric mean ratio [GMR] [90% confidence interval (CI)] 0.78 [0.55, 1.11]). The urinary 3-hydroxymorphinan (3HM)/DX MR (CYP3A activity marker) was increased (GMR 1.18 [1.03, 1.35]) and the urinary DX-O-glucuronide/DX MR (proposed UGT2B marker) was increased (GMR 1.22 [1.06, 1.39]) with isotretinoin administration compared to pre-isotretinoin. CONCLUSIONS: Administration of isotretinoin did not significantly reduce CYP2D6 activity in extensive metabolizers, suggesting that the predicted downregulation of CYP2D6 based on in vitro data does not translate into humans. We observed a modest increase in CYP3A activity (predominantly CYP3A4) with isotretinoin treatment. The data also suggest that DX glucuronidation is increased following isotretinoin administration.

Predicting complex kidney drug handling using a physiologically-based pharmacokinetic model informed by biomarker-estimated secretory clearance and blood flow.

Kidney function-adjusted drug dosing is currently based solely on the estimated glomerular filtration rate (GFR), however, kidney drug handling is accomplished by a combination of filtration, tubular secretion, and re-absorption. Mechanistic physiologically-based pharmacokinetic (PBPK) models recapitulate anatomic compartments to predict elimination from estimated perfusion, filtration, secretion, and re-absorption, but clinical applications are limited by a lack of empiric individual-level measurements of these functions. We adapted and validated a PBPK model to predict drug clearance from individual biomarker-based estimates of kidney perfusion and secretory clearance. We estimated organic anion transporter-mediated secretion via kynurenic acid clearance and kidney blood flow (KBF) via isovalerylglycine clearance in human participants, incorporating these measurements with GFR into the model to predict kidney drug clearance. We compared measured and model-predicted clearances of administered tenofovir and oseltamivir, which are cleared by both filtration and secretion. There were 27 outpatients (age 55 ± 15 years, mean iohexol-GFR [iGFR] 76 ± 31 mL/min/1.73 m2 ) in this drug clearance study. The mean observed and mechanistic model-predicted tenofovir clearances were 169 ± 102 mL/min and 163 ± 80 mL/min, respectively; estimated mean error of the mechanistic model was 37.1 mL/min (95% confidence interval [CI]: 24-52.9), compared to a mean error of 41.8 mL/min (95% CI: 25-61.6) from regression model. The mean observed and model-predicted oseltamivir carboxylate clearances were 183 ± 104 mL/min and 179 ± 89 mL/min, respectively; estimated mean error of the mechanistic model was 42.9 mL/min (95% CI: 29.7-56.4), versus error of 48.1 mL/min (95% CI: 31.2-67.3) from the regression model. Individualized estimates of tubular secretion and KBF improved the accuracy of PBPK model-predicted tenofovir and oseltamivir kidney clearances, suggesting the potential for biomarker-informed measures of kidney function to refine personalized drug dosing.

Incorporating Uremic Solute-mediated Inhibition of OAT1/3 Improves PBPK Prediction of Tenofovir Renal and Systemic Disposition in Patients with Severe Kidney Disease.

BACKGROUND: Dose modification of renally secreted drugs in patients with chronic kidney disease (CKD) has relied on serum creatinine concentration as a biomarker to estimate glomerular filtration (GFR) under the assumption that filtration and secretion decline in parallel. A discrepancy between actual renal clearance and predicted renal clearance based on GFR alone is observed in severe CKD patients with tenofovir, a compound secreted by renal OAT1/3. Uremic solutes that inhibit OAT1/3 may play a role in this divergence. METHODS: To examine the impact of transporter inhibition by uremic solutes on tenofovir renal clearance, we determined the inhibitory potential of uremic solutes hippuric acid, indoxyl sulfate, and p-cresol sulfate. The inhibition parameters (IC50) were incorporated into a previously validated mechanistic kidney model; simulated renal clearance and plasma PK profile were compared to data from clinical studies. RESULTS: Without the incorporation of uremic solute inhibition, the PBPK model failed to capture the observed data with an absolute average fold error (AAFE) > 2. However, when the inhibition of renal uptake transporters and uptake transporters in the slow distribution tissues were included, the AAFE value was within the pre-defined twofold model acceptance criterion, demonstrating successful model extrapolation to CKD patients. CONCLUSION: A PBPK model that incorporates inhibition by uremic solutes has potential to better predict renal clearance and systemic disposition of secreted drugs in patients with CKD. Ongoing research is warranted to determine if the model can be expanded to include other OAT1/3 substrate drugs and to evaluate how these findings can be translated to clinical guidance for drug selection and dose optimization in patients with CKD.

CYP2D6 Activity Is Correlated with Changes in Plasma Concentrations of Taurocholic Acid during Pregnancy and Postpartum in CYP2D6 Extensive Metabolizers.

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of >20% of marketed drugs. CYP2D6 expression and activity exhibit high interindividual variability and is induced during pregnancy. The farnesoid X receptor (FXR) is a transcriptional regulator of CYP2D6 that is activated by bile acids. In pregnancy, elevated plasma bile acid concentrations are associated with maternal and fetal risks. However, modest changes in bile acid concentrations may occur during healthy pregnancy, thereby altering FXR signaling. A previous study demonstrated that hepatic tissue concentrations of bile acids positively correlated with the hepatic mRNA expression of CYP2D6. This study sought to characterize the plasma bile acid metabolome in healthy women (n = 47) during midpregnancy (25-28 weeks gestation) and ≥3 months postpartum and to determine if plasma bile acids correlate with CYP2D6 activity. It is hypothesized that during pregnancy, plasma bile acids would favor less hydrophobic bile acids (cholic acid vs. chenodeoxycholic acid) and that plasma concentrations of cholic acid and its conjugates would positively correlate with the urinary ratio of dextrorphan/dextromethorphan. At 25-28 weeks gestation, taurine-conjugated bile acids comprised 23% of the quantified serum bile acids compared with 7% ≥3 months postpartum. Taurocholic acid positively associated with the urinary ratio of dextrorphan/dextromethorphan, a biomarker of CYP2D6 activity. Collectively, these results confirm that the bile acid plasma metabolome differs between pregnancy and postpartum and provide evidence that taurocholic acid may impact CYP2D6 activity during pregnancy. SIGNIFICANCE STATEMENT: Bile acid homeostasis is altered in pregnancy, and plasma concentrations of taurocholic acid positively correlate with CYP2D6 activity. Differences between plasma and/or tissue concentrations of farnesoid X receptor ligands such as bile acids may contribute to the high interindividual variability in CYP2D6 expression and activity.

Venetoclax and idasanutlin in relapsed/refractory AML: a nonrandomized, open-label phase 1b trial.

This phase 1b trial (NCT02670044) evaluated venetoclax-idasanutlin in patients with relapsed/refractory (R/R) acute myeloid leukemia (AML) ineligible for cytotoxic chemotherapy. Two-dimensional dose escalation (DE, n = 50) was performed for venetoclax daily with idasanutlin on days 1 to 5 in 28-day cycles, followed by dosing schedule optimization (n = 6) to evaluate reduced venetoclax schedules (21-/14-day dosing). Common adverse events (occurring in ≥40% of patients) included diarrhea (87.3% of patients), nausea (74.5%), vomiting (52.7%), hypokalemia (50.9%), and febrile neutropenia (45.5%). During DE, across all doses, composite complete remission (CRc; CR + CR with incomplete blood count recovery + CR with incomplete platelet count recovery) rate was 26.0% and morphologic leukemia-free state (MLFS) rate was 12%. For anticipated recommended phase 2 doses (venetoclax 600 mg + idasanutlin 150 mg; venetoclax 600 mg + idasanutlin 200 mg), the combined CRc rate was 34.3% and the MLFS rate was 14.3%. Pretreatment IDH1/2 and RUNX1 mutations were associated with higher CRc rates (50.0% and 45.0%, respectively). CRc rate in patients with TP53 mutations was 20.0%, with responses noted among those with co-occurring IDH and RUNX1 mutations. In 12 out of 36 evaluable patients, 25 emergent TP53 mutations were observed; 22 were present at baseline with low TP53 variant allele frequency (median 0.0095% [range, 0.0006-0.4]). Venetoclax-idasanutlin showed manageable safety and encouraging efficacy in unfit patients with R/R AML. IDH1/2 and RUNX1 mutations were associated with venetoclax-idasanutlin sensitivity, even in some patients with co-occurring TP53 mutations; most emergent TP53 clones were preexisting. Our findings will aid ongoing/future trials of BCL-2/MDM2 inhibitor combinations. This trial was registered at www.clinicaltrials.gov as #NCT02670044.

Impact of Pregnancy and Vitamin A Supplementation on CYP2D6 Activity.

The mechanism of cytochrome P450 2D6 (CYP2D6) induction during pregnancy has not been evaluated in humans. This study assessed the changes in CYP2D6 and CYP3A activities during pregnancy and postpartum, and the effect of vitamin A administration on CYP2D6 activity. Forty-seven pregnant CYP2D6 extensive metabolizers (with CYP2D6 activity scores of 1 to 2) received dextromethorphan (DM) 30 mg orally as a single dose during 3 study windows (at 25 to 28 weeks of gestation, study day 1; at 28 to 32 weeks of gestation, study day 2; and at ≥3 months postpartum, study day 3). Participants were randomly assigned to groups with no supplemental vitamin A (control) or with supplemental vitamin A (10 000 IU/day orally for 3 to 4 weeks) after study day 1. Concentrations of DM and its metabolites, dextrorphan (DX) and 3-hydroxymorphinan (3HM), were determined from a 2-hour post-dose plasma sample and cumulative 4-hour urine sample using liquid chromatography-mass spectrometry. Change in CYP2D6 activity was assessed using DX/DM plasma and urine metabolic ratios. The activity change in CYP3A was also assessed using the 3HM/DM urine metabolic ratio. The DX/DM urine ratio was significantly higher (43%) in pregnancy compared with postpartum (P = .03), indicating increased CYP2D6 activity. The DX/DM plasma ratio was substantially higher in the participants, with an activity score of 1.0 during pregnancy (P = .04) compared with postpartum. The 3HM/DM urinary ratio was significantly higher (92%) during pregnancy, reflecting increased CYP3A activity (P = .02). Vitamin A supplementation did not change CYP2D6 activity during pregnancy; however, plasma all-trans retinoic acid (atRA) concentrations were positively correlated with increased CYP2D6 activity during pregnancy and postpartum. Further research is needed to elucidate the mechanisms of increased CYP2D6 activity during pregnancy.

Development of a pediatric physiologically-based pharmacokinetic model to support recommended dosing of atezolizumab in children with solid tumors.

Background: Atezolizumab has been studied in multiple indications for both pediatric and adult patient populations. Generally, clinical studies enrolling pediatric patients may not collect sufficient pharmacokinetic data to characterize the drug exposure and disposition because of operational, ethical, and logistical challenges including burden to children and blood sample volume limitations. Therefore, mechanistic modeling and simulation may serve as a tool to predict and understand the drug exposure in pediatric patients. Objective: To use mechanistic physiologically-based pharmacokinetic (PBPK) modeling to predict atezolizumab exposure at a dose of 15 mg/kg (max 1,200 mg) in pediatric patients to support dose rationalization and label recommendations. Methods: A minimal mechanistic PBPK model was used which incorporated age-dependent changes in physiology and biochemistry that are related to atezolizumab disposition such as endogenous IgG concentration and lymph flow. The PBPK model was developed using both in vitro data and clinically observed data in adults and was verified across dose levels obtained from a phase I and multiple phase III studies in both pediatric patients and adults. The verified model was then used to generate PK predictions for pediatric and adult subjects ranging from 2- to 29-year-old. Results: Individualized verification in children and in adults showed that the simulated concentrations of atezolizumab were comparable (76% within two-fold and 90% within three-fold, respectively) to the observed data with no bias for either over- or under-prediction. Applying the verified model, the predicted exposure metrics including Cmin, Cmax, and AUCtau were consistent between pediatric and adult patients with a geometric mean of pediatric exposure metrics between 0.8- to 1.25-fold of the values in adults. Conclusion: The results show that a 15 mg/kg (max 1,200 mg) atezolizumab dose administered intravenously in pediatric patients provides comparable atezolizumab exposure to a dose of 1,200 mg in adults. This suggests that a dose of 15 mg/kg will provide adequate and effective atezolizumab exposure in pediatric patients from 2- to 18-year-old.

Population Pharmacokinetics and Exposure-Response Analyses for Venetoclax in Combination with R-CHOP in Relapsed/Refractory and Previously Untreated Patients with Diffuse Large B Cell Lymphoma.

INTRODUCTION: Outcomes remain poor in patients with diffuse large B cell lymphoma (DLBCL) who overexpress BCL-2 protein. We present population pharmacokinetics (PopPK) and exposure-response (ER) analyses for venetoclax (a selective BCL-2 inhibitor) administered with rituximab-cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) in patients with relapsed/refractory (R/R) and previously untreated (1L) non-Hodgkin lymphoma (NHL) from the phase 1b/2 CAVALLI study, to confirm dose selection for future studies. METHODS: Analyses included 216 patients with R/R or 1L NHL treated for eight 21-day cycles with 400-800 mg venetoclax (cycle 1: days 4-10; cycles 2-8: days 1-10) in combination with R for eight cycles and CHOP for 6-8 cycles. A legacy PopPK model for venetoclax was used to describe the observed data and provide post hoc PK parameters. Venetoclax steady-state exposure (AUCss) was used to predict clinical efficacy, safety, or tolerability. To isolate the effect of venetoclax, ER analyses referenced data from the R-CHOP arm of a historical control study, GOYA, in 1L DLBCL. RESULTS: There was no significant association between venetoclax AUCss and progression-free survival or complete response either for all-comers or the BCL-2-immunohistochemistry-positive subpopulation. No statistically significant trends were observed with venetoclax AUCss and the key grade ≥ 3 adverse events and serious adverse events. Similar dose intensities were observed for venetoclax and R-CHOP components across venetoclax exposures, suggesting venetoclax did not impact delivery of the R-CHOP backbone. CONCLUSIONS: The PopPK and ER analyses, in addition to the positive benefit-risk observed in the clinical data, support the selection of 800 mg venetoclax given with R-CHOP for future studies in BCL-2-immunohistochemistry-positive patients with 1L DLBCL. TRIAL REGISTRATION: ClinicalTrials.gov Identifier NCT02055820.

Clinical Pharmacology Perspectives for Adoptive Cell Therapies in Oncology.

Adoptive cell therapies (ACTs) have shown transformative efficacy in oncology with five US Food and Drug Administration (FDA) approvals for chimeric antigen receptor (CAR) T-cell therapies in hematological malignancies, and promising activity for T cell receptor T-cell therapies in both liquid and solid tumors. Clinical pharmacology can play a pivotal role in optimizing ACTs, aided by modeling and simulation toolboxes and deep understanding of the underlying biological and immunological processes. Close collaboration and multilevel data integration across functions, including chemistry, manufacturing, and control, biomarkers, bioanalytical, and clinical science and safety teams will be critical to ACT development. As ACT is comprised of alive, polyfunctional, and heterogeneous immune cells, its overall physicochemical and pharmacological property is vastly different from other platforms/modalities, such as small molecule and protein therapeutics. In this review, we first describe the unique kinetics of T cells and the appropriate bioanalytical strategies to characterize cellular kinetics. We then assess the distinct aspects of clinical pharmacology for ACTs in comparison to traditional small molecule and protein therapeutics. Additionally, we provide a review for the five FDA-approved CAR T-cell therapies and summarize their properties, cellular kinetic characteristics, dose-exposure-response relationship, and potential baseline factors/variables in product, patient, and regimen that may affect the safety and efficacy. Finally, we probe into existing empirical and mechanistic quantitative techniques to understand how various modeling and simulation approaches can support clinical pharmacology strategy and propose key considerations to be incorporated and explored in future models.

Objective Identification of Cannabis Use Levels in Clinical Populations Is Critical for Detecting Pharmacological Outcomes.

Introduction: Cannabis is widely used for recreational and medical purposes, but its therapeutic efficacy remains unresolved for many applications as data from retrospective studies show dramatic discrepancy. We hypothesized that false self-reporting of cannabis use and lack of differentiation of heavy users from light or occasional users contribute to the conflicting outcomes. Objective: The goal of this study was to develop an objective biomarker of cannabis use and test how application of such biomarker impacts clinical study outcomes and dose-response measures. Methods and Analysis: Population pharmacokinetic (PK) models of (-)-trans-Δ9-tetrahydrocannabinol (THC) and its metabolites 11-hydroxy-Δ9-tetrahydrocannabinol (11-OH-THC) and 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (11-COOH-THC) were developed based on published studies reporting cannabinoid disposition in individual subjects following intravenous administration or smoking of cannabis. Plasma 11-COOH-THC concentration distributions in different cannabis user groups smoking cannabis were generated via Monte Carlo simulations, and plasma concentration cutoff values of 11-COOH-THC were developed to differentiate light and heavy daily cannabis users in clinical studies. The developed cutoff value was then applied to a retrospective study that assessed the impact of cannabis use on T cell activation in subjects with HIV who self-reported as either nonuser or daily user of cannabis. Results: The developed population PK models established plasma 11-COOH-THC concentration of 73.1 µg/L as a cutoff value to identify heavy daily users, with a positive predictive value of 80% in a mixed population of equal proportions of once daily and three times a day users. The stratification allowed detection of changes in T cell activation in heavy users which was not detected based on self-reporting or detectability of plasma cannabinoids. A proof-of-concept power analysis demonstrated that implementation of such cutoff value greatly increases study power and sensitivity to detect pharmacological effects of cannabis use. Conclusions: This study shows that the use of plasma 11-COOH-THC concentration cutoff value as an objective measure to classify cannabis use in target populations is critical for study sensitivity and specificity and provides much needed clarity for addressing dose-response relationships and therapeutic effects of cannabis.

Bridging the gap between in silico and in vivo by modeling opioid disposition in a kidney proximal tubule microphysiological system.

Opioid overdose, dependence, and addiction are a major public health crisis. Patients with chronic kidney disease (CKD) are at high risk of opioid overdose, therefore novel methods that provide accurate prediction of renal clearance (CLr) and systemic disposition of opioids in CKD patients can facilitate the optimization of therapeutic regimens. The present study aimed to predict renal clearance and systemic disposition of morphine and its active metabolite morphine-6-glucuronide (M6G) in CKD patients using a vascularized human proximal tubule microphysiological system (VPT-MPS) coupled with a parent-metabolite full body physiologically-based pharmacokinetic (PBPK) model. The VPT-MPS, populated with a human umbilical vein endothelial cell (HUVEC) channel and an adjacent human primary proximal tubular epithelial cells (PTEC) channel, successfully demonstrated secretory transport of morphine and M6G from the HUVEC channel into the PTEC channel. The in vitro data generated by VPT-MPS were incorporated into a mechanistic kidney model and parent-metabolite full body PBPK model to predict CLr and systemic disposition of morphine and M6G, resulting in successful prediction of CLr and the plasma concentration-time profiles in both healthy subjects and CKD patients. A microphysiological system together with mathematical modeling successfully predicted renal clearance and systemic disposition of opioids in CKD patients and healthy subjects.

Altered vitamin A metabolism in human liver slices corresponds to fibrogenesis.

All-trans-retinoic acid (atRA), the active metabolite of vitamin A, has antifibrogenic properties in vitro and in animal models. Liver vitamin A homeostasis is maintained by cell-specific enzymatic activities including storage in hepatic stellate cells (HSCs), secretion into circulation from hepatocytes, and formation and clearance of atRA. During chronic liver injury, HSC activation is associated with a decrease in liver retinyl esters and retinol concentrations. atRA is synthesized through two enzymatic steps from retinol, but it is unknown if the loss of retinoid stores is associated with changes in atRA formation and which cell types contribute to the metabolic changes. The aim of this study was to determine if the vitamin A metabolic flux is perturbed in acute liver injury, and if changes in atRA concentrations are associated with HSC activation and collagen expression. At basal levels, HSC and Kupffer cells expressed key genes involved in vitamin A metabolism, whereas after acute liver injury, complex changes to the metabolic flux were observed in liver slices. These changes include a reproducible spike in atRA tissue concentrations, decreased retinyl ester and atRA formation rate, and time-dependent changes to the expression of metabolizing enzymes. Kinetic simulations suggested that oxidoreductases are important in determining retinoid metabolic flux after liver injury. These early changes precede HSC activation and upregulation of profibrogenic gene expression, which were inversely correlated with atRA tissue concentrations, suggesting that HSC and Kupffer cells are key cells involved in changes to vitamin A metabolic flux and signaling after liver injury. Study Highlights WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC? Vitamin A is metabolized in the liver for storage as retinyl esters in hepatic stellate cell (HSCs) or to all-trans-retinoic acid (atRA), an active metabolite with antifibrogenic properties. Following chronic liver injury, vitamin A metabolic flux is perturbed, and HSC activation leads to diminished retinoid stores. WHAT QUESTION DID THIS STUDY ADDRESS? Do changes in the expression of vitamin A metabolizing enzymes explain changes in atRA concentrations and the regulation of fibrosis following acute liver injury? WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE? In healthy liver, both HSC and Kupffer cells may mediate vitamin A homeostasis. Following acute liver injury, complex changes in metabolizing enzyme expression/activity alter the metabolic flux of retinoids, resulting in a transient peak in atRA concentrations. The atRA concentrations are inversely correlated with profibrogenic gene expression, HSC activation, and collagen deposition. HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE? Improved understanding of altered vitamin A metabolic flux in acute liver injury may provide insight into cell-specific contributions to vitamin A loss and lead to novel interventions in liver fibrosis.

Knockout of Cyp26a1 and Cyp26b1 during postnatal life causes reduced lifespan, dermatitis, splenomegaly, and systemic inflammation in mice.

All-trans-retinoic acid (atRA), the active metabolite of vitamin A, is an essential signaling molecule in all chordates. Global knockouts of the atRA clearing enzymes Cyp26a1 or Cyp26b1 are embryonic lethal. In adult rodents, inhibition of Cyp26a1 and Cyp26b1 increases atRA concentrations and signaling. However, postnatal knockout of Cyp26a1 does not cause a severe phenotype. We hypothesized that Cyp26b1 is the main atRA clearing Cyp in postnatal mammals. This hypothesis was tested by generating tamoxifen-inducible knockout mouse models of Cyp26b1 alone or with Cyp26a1. Both mouse models showed dermatitis, blepharitis, and splenomegaly. Histology showed infiltration of inflammatory cells including neutrophils and T lymphocytes into the skin and hyperkeratosis/hyperplasia of the nonglandular stomach. The mice lacking both Cyp26a1 and Cyp26b1 also had a reduced lifespan, failed to gain weight, and showed fat atrophy. There were significant changes in vitamin A homeostasis. Postnatal knockout of Cyp26b1 resulted in increased atRA concentrations in the skin while the postnatal knockout of both Cyp26a1 and Cyp26b1 resulted in increased atRA concentrations in the liver, serum, skin, spleen, and intestines. This study demonstrates the paramount role of Cyp26b1 in regulating retinoid homeostasis in postnatal life.

Novel Mechanistic PBPK Model to Predict Renal Clearance in Varying Stages of CKD by Incorporating Tubular Adaptation and Dynamic Passive Reabsorption.

Chronic kidney disease (CKD) has significant effects on renal clearance (CLr ) of drugs. Physiologically-based pharmacokinetic (PBPK) models have been used to predict CKD effects on transporter-mediated renal active secretion and CLr for hydrophilic nonpermeable compounds. However, no studies have shown systematic PBPK modeling of renal passive reabsorption or CLr for hydrophobic permeable drugs in CKD. The goal of this study was to expand our previously developed and verified mechanistic kidney model to develop a universal model to predict changes in CLr in CKD for permeable and nonpermeable drugs that accounts for the dramatic nonlinear effect of CKD on renal passive reabsorption of permeable drugs. The developed model incorporates physiologically-based tubular changes of reduced water reabsorption/increased tubular flow rate per remaining functional nephron in CKD. The final adaptive kidney model successfully (absolute fold error (AFE) all < 2) predicted renal passive reabsorption and CLr for 20 permeable and nonpermeable test compounds across the stages of CKD. In contrast, use of proportional glomerular filtration rate reduction approach without addressing tubular adaptation processes in CKD to predict CLr generated unacceptable CLr predictions (AFE = 2.61-7.35) for permeable compounds in severe CKD. Finally, the adaptive kidney model accurately predicted CLr of para-amino-hippuric acid and memantine, two secreted compounds, in CKD, suggesting successful integration of active secretion into the model, along with passive reabsorption. In conclusion, the developed adaptive kidney model enables mechanistic predictions of in vivo CLr through CKD progression without any empirical scaling factors and can be used for CLr predictions prior to assessment of drug disposition in renal impairment.

Mechanistic PBPK Modeling of Urine pH Effect on Renal and Systemic Disposition of Methamphetamine and Amphetamine.

The effect of urine pH on renal excretion and systemic disposition has been observed for many drugs and metabolites. When urine pH is altered, tubular ionization, passive reabsorption, renal clearance, and systemic exposure of drugs and metabolites may all change dramatically, raising clinically significant concerns. Surprisingly, the urine pH effect on drug disposition is not routinely explored in humans, and regulatory agencies have neither developed guidance on this issue nor required industry to conduct pertinent human trials. In this study, we hypothesized that physiologically based pharmacokinetic (PBPK) modeling could be used as a cost-effective method to examine potential urine pH effect on drug and metabolite disposition. Our previously developed and verified mechanistic kidney model was integrated with a full-body PBPK model to simulate renal clearance and area under the plasma concentration-time curve (AUC) with varying urine pH statuses using methamphetamine and amphetamine as model compounds. We first developed and verified drug models for methamphetamine and amphetamine under normal urine pH condition [absolute average fold error (AAFE) < 1.25 at study level]. Then, acidic and alkaline urine scenarios were simulated. Our simulation results show that the renal excretion and plasma concentration-time profiles for methamphetamine and amphetamine could be recapitulated under different urine pH (AAFE < 2 at individual level). The methamphetamine-amphetamine parent-metabolite full-body PBPK model also successfully simulated amphetamine plasma concentration-time profiles (AAFE < 1.25 at study level) and amphetamine/methamphetamine urinary concentration ratios (AAFE < 2 at individual level) after dosing methamphetamine. This demonstrates that our mechanistic PBPK model can predict urine pH effect on systemic and urinary disposition of drugs and metabolites. SIGNIFICANCE STATEMENT: Our study shows that integrating mechanistic kidney model with full-body physiologically based pharmacokinetic model can predict the magnitude of alteration in renal excretion and area under the plasma concentration-time curve (AUC) of drugs and metabolites when urine pH is changed. This provides a cost-effective method to evaluate the likelihood of renal and systemic disposition changes due to varying urine pH. This is important because multiple drugs and diseases can alter urine pH, leading to quantitatively and clinically significant changes in drug and metabolite disposition that may require adjustment of therapy.

Sampling Site Has a Critical Impact on Physiologically Based Pharmacokinetic Modeling.

It has been shown that arterial (central) and venous (peripheral) plasma drug concentrations can be very different. While pharmacokinetic studies typically measure drug concentrations from the peripheral vein such as the arm vein, physiologically based pharmacokinetic (PBPK) models generally output simulated concentrations from the central venous compartment that physiologically represents the right atrium, a merge of the superior and inferior vena cava. In this study, a physiologically based peripheral forearm sampling site model was developed and verified using nicotine, ketamine, lidocaine, and fentanyl as model drugs. This verified model allows output of simulated peripheral venous concentrations that can be meaningfully compared with observed pharmacokinetic data from the arm vein. The generalized effect of PBPK model sampling site on simulation output was investigated. Drugs and metabolites with large volumes of distribution showed considerable concentration discrepancy between the simulated central venous compartment and the peripheral arm vein after intravenous or oral administration, resulting in significant differences in values for C max and time taken to reach C max (t max ) In addition, the simulated central venous metabolite profile showed an unexpected profile that was not observed in the peripheral arm vein. Using fentanyl as a model compound, we show that using the wrong sampling site in PBPK models can lead to biased model evaluation and subsequent erroneous model parameter optimization. Such an error in model parameters along with the discrepant sampling site could dramatically mislead the pharmacokinetic prediction in unstudied clinical scenarios, affecting the assessment of drug safety and efficacy. Overall, this study shows that PBPK model publications should specify the model sampling sites and match them with those employed in clinical studies. SIGNIFICANCE STATEMENT: Our study shows that sampling from the central venous compartment (right atrium) during physiologically based pharmacokinetic model development gives rise to biased model evaluation and erroneous model parameterization when observed data are collected from the peripheral arm vein. This can lead to a clinically significant error in predictions of plasma concentration-time profiles in unstudied scenarios. To address this error, we developed and verified a novel peripheral sampling site model to simulate arm vein drug concentrations that can be applied to different drug dosing scenarios.