Understanding CYP3A4 and P-gp mediated drug-drug interactions through PBPK modeling - Case example of pralsetinib

Pralsetinib, a potent and selective inhibitor of oncogenic RET fusion and RET mutant proteins, is a substrate of the drug metabolizing enzyme CYP3A4 and a substrate of the efflux transporter P-gp based on in vitro data. Therefore, its pharmacokinetics (PKs) may be affected by co-administration of potent CYP3A4 inhibitors and inducers, P-gp inhibitors, and combined CYP3A4 and P-gp inhibitors. With the frequent overlap between CYP3A4 and P-gp substrates/inhibitors, pralsetinib is a challenging and representative example of the need to more quantitatively characterize transporter-enzyme interplay. A physiologically-based PK (PBPK) model for pralsetinib was developed to understand the victim drug-drug interaction (DDI) risk for pralsetinib. The key parameters driving the magnitude of pralsetinib DDIs, the P-gp intrinsic clearance and the fraction metabolized by CYP3A4, were determined from PBPK simulations that best captured observed DDIs from three clinical studies. Sensitivity analyses and scenario simulations were also conducted to ensure these key parameters were determined with sound mechanistic rationale based on current knowledge, including the worst-case scenarios. The verified pralsetinib PBPK model was then applied to predict the effect of other inhibitors and inducers on the PKs of pralsetinib. This work highlights the challenges in understanding DDIs when enzyme-transporter interplay occurs, and demonstrates an important strategy for differentiating enzyme/transporter contributions to enable PBPK predictions for untested scenarios and to inform labeling.

Machine Learning-Based Quantification of Patient Factors Impacting Remission in Patients With Ulcerative Colitis: Insights from Etrolizumab Phase III Clinical Trials

Etrolizumab, an investigational anti-β7 integrin monoclonal antibody, has undergone evaluation for safety and efficacy in phase III clinical trials on patients with moderate to severe ulcerative colitis (UC). Etrolizumab was terminated because mixed efficacy results were shown in the induction and maintenance phase in patients with UC. In this post hoc analysis, we characterized the impact of explanatory variables on the probability of remission using XGBoost machine learning (ML) models alongside with the SHapley Additive exPlanations framework for explainability. We used patient-level data encompassing demographics, physiology, disease history, clinical questionnaires, histology, serum biomarkers, and etrolizumab drug exposure to develop ML models aimed at predicting remission. Baseline covariates and early etrolizumab exposure at week 4 in the induction phase were utilized to develop an induction ML model, whereas covariates from the end of the induction phase and early etrolizumab exposure at week 4 in the maintenance phase were used to develop a maintenance ML model. Both the induction and maintenance ML models exhibited good predictive performance, achieving an area under the receiver operating characteristic curve (AUROC) of 0.74 ± 0.03 and 0.75 ± 0.06 (mean ± SD), respectively. Compared with placebo, the highest tertile of etrolizumab exposure contributed to 15.0% (95% confidence interval (CI): 9.7-19.9) and 17.0% (95% CI: 8.1-26.4) increases in remission probability in the induction and maintenance phases, respectively. Additionally, the key covariates that predicted remission were CRP, MAdCAM-1, and stool frequency for the induction phase and white blood cells, fecal calprotectin and age for the maintenance phase. These findings hold significant implications for establishing stratification factors in the design of future clinical trials.

Population pharmacokinetics and pharmacodynamics of efmarodocokin alfa (IL-22Fc)

Efmarodocokin alfa (IL-22Fc) is a fusion protein of human IL-22 linked to the crystallizable fragment (Fc) of human IgG4. It has been tested in multiple indications including inflammatory bowel disease (IBD). The purposes of the present analyses were to describe the population pharmacokinetics (PK) of efmarodocokin alfa and perform pharmacodynamic (PD) analysis on the longitudinal changes of the PD biomarker REG3A after efmarodocokin alfa treatment as well as identify covariates that affect efmarodocokin alfa PK and REG3A PD. The data used for this analysis included 182 subjects treated with efmarodocokin alfa in two clinical studies. The population PK and PD analyses were conducted sequentially. Efmarodocokin alfa concentration-time data were analyzed using a nonlinear mixed-effects modeling approach, and an indirect response model was adopted to describe the REG3A PD data with efmarodocokin alfa serum concentration linked to the increase in REG3A. The analysis software used were NONMEM and R. A 3-compartment model with linear elimination best described the PK of efmarodocokin alfa. The estimated population-typical value for clearance (CL) was 1.12 L/day, and volume of central compartment was 6.15 L. Efmarodocokin alfa CL increased with higher baseline body weight, C-reactive protein, and CL was 27.6% higher in IBD patients compared to healthy subjects. The indirect response PD model adequately described the longitudinal changes of REG3A after efmarodocokin alfa treatment. A popPK and PD model for efmarodocokin alfa and REG3A was developed and covariates affecting the PK and PD were identified.

Exposure-Response Relationships for Pralsetinib in Patients with RET-Altered Thyroid Cancer or RET Fusion-Positive Nonsmall Cell Lung Cancer

Pralsetinib is a highly potent oral kinase inhibitor of oncogenic RET (rearranged during transfection) fusions and mutations. Pralsetinib received approval from the United States Food and Drug Administration for the treatment of patients with metastatic RET fusion-positive non-small cell lung cancer (NSCLC), and received accelerated approval for the treatment of patients with RET fusion-positive thyroid cancer. Exposure-response (ER) analyses of efficacy were performed separately in patients with thyroid cancer and in patients with NSCLC, but data for all patients were pooled for the safety analysis. ER models were developed with time-varying exposure; the effect of covariates was also examined. For patients with NSCLC, a higher starting dose was associated with improved progression-free survival (PFS), but this improvement did not correlate with a higher exposure overall. Significant covariates included sex and baseline Eastern Cooperative Oncology Group (ECOG) score. For patients with thyroid cancer, a higher exposure was associated with improved PFS. Significant covariates included prior systemic cancer therapy and ECOG score. For safety, higher exposure was associated with a greater risk of grade ≥3 anemia, pneumonia, and lymphopenia. Patients with an ECOG score of ≥1 had an increased risk of grade ≥3 pneumonia. Non-White patients had a lower risk of grade ≥3 lymphopenia. ER analysis revealed that higher pralsetinib exposure was associated with improved PFS in thyroid cancer, but not in NSCLC. However, a higher starting dose (ie, 400 vs ≤300 mg daily) was correlated with better PFS for all indications. Higher exposure was also associated with an increased risk of grade ≥3 adverse events (AEs); however, the overall incidence of these events was acceptably low (≤20%). This analysis supports the use of a 400 mg starting dose of pralsetinib, allowing for dose reduction in the event of AEs.

A disease model predicting placebo response and remission status of patients with ulcerative colitis using modified Mayo score

The Mayo Clinical Score is used in clinical trials to describe the clinical status of patients with ulcerative colitis (UC). It comprises four subscores: rectal bleeding (RB), stool frequency (SF), physician's global assessment, and endoscopy (ENDO). According to recent US Food and Drug Administration guidelines (Ulcerative colitis: developing drugs for treatment, Guidance Document, https://www.fda.gov/regulatory-information/s. 2022), clinical response and remission should be based on modified Mayo Score (mMS) relying on RB, SF, and ENDO. Typically, ENDO is performed at the beginning and end of each phase, whereas RB and SF are more frequently available. Item response theory (IRT) models allow the shared information to be used for prediction of all subscores at each observation time; therefore, it leverages information from RB and SF to predict ENDO. A UC disease IRT model was developed based on four etrolizumab phase III studies to describe the longitudinal mMS subscores, placebo response, and remission at the end of induction and maintenance. For each subscore, a bounded integer model was developed. The placebo response was characterized by a mono-exponential function acting on all mMS subscores similarly. The final model reliably predicted longitudinal mMS data. In addition, remission was well-predicted by the model, with only 5% overprediction at the end of induction and 3% underprediction at the end of maintenance. External evaluation of the final model using placebo arms from five different studies indicated adequate performance for both longitudinal mMS subscores and remission status. These results suggest utility of the current disease model for informed decision making in UC clinical development, such as assisting future clinical trial designs and evaluations.

External validation of a tumor growth inhibition-overall survival model in non-small-cell lung cancer based on atezolizumab studies using alectinib data

BACKGROUND

A modeling framework was previously developed to simulate overall survival (OS) using tumor growth inhibition (TGI) data from six randomized phase 2/3 atezolizumab monotherapy or combination studies in non-small-cell lung cancer (NSCLC). We aimed to externally validate this framework to simulate OS in patients with treatment-naive advanced anaplastic lymphoma kinase (ALK)-positive NSCLC in the alectinib ALEX study.

METHODS

TGI metrics were estimated from a biexponential model using longitudinal tumor size data from a Phase 3 study evaluating alectinib compared with crizotinib in patients with treatment-naive ALK-positive advanced NSCLC. Baseline prognostic factors and TGI metric estimates were used to predict OS.

RESULTS

286 patients were evaluable (at least baseline and one post-baseline tumor size measurements) out of 303 (94%) followed for up to 5 years (cut-off: 29 November 2019). The tumor growth rate estimate and baseline prognostic factors (inflammatory status, tumor burden, Eastern Cooperative Oncology Group performance status, race, line of therapy, and sex) were used to simulate OS in ALEX study. Observed survival distributions for alectinib and crizotinib were within model 95% prediction intervals (PI) for approximately 2 years. Predicted hazard ratio (HR) between alectinib and crizotinib was in agreement with the observed HR (predicted HR 0.612, 95% PI 0.480-0.770 vs. 0.625 observed HR).

CONCLUSION

The TGI-OS model based on unselected or PD-L1 selected NSCLC patients included in atezolizumab trials is externally validated to predict treatment effect (HR) in a biomarker-selected (ALK-positive) population included in alectinib ALEX trial suggesting that TGI-OS models may be treatment independent.

Population-Based Pharmacodynamic Modeling of Omalizumab in Pediatric Patients with Moderate to Severe Persistent Inadequately Controlled Allergic Asthma

Omalizumab is the first approved anti-immunoglobulin E (IgE) agent for the treatment of moderate to severe persistent inadequately controlled allergic asthma in adults and adolescents (≥ 12 years old). In 2016, it was approved in pediatric patients (6-11 years old). The objective of this study was to quantitatively characterize the relationship between serum free IgE and pulmonary function (as measured by forced expiratory volume in 1 s [FEV1]) in pediatrics using a population-based pharmacodynamic model. Data collected during the steroid-stable period (first 24 weeks) of an omalizumab trial with pediatric asthma patients (Study IA05) were used to build the pediatric IgE-FEV1 model. The previously developed population IgE-FEV1 model in adults/adolescents was adapted to characterize the FEV1 and IgE relationship in pediatrics with different magnitude and onset of response. The pediatric IgE-FEV1 model adequately characterized the IgE-FEV1 relationship in pediatrics, particularly at the extremes of the observed body weights (i.e., ≤ 30 kg) and IgE values at screening (i.e., > 700 IU/mL). The estimated sigmoidal free IgE-FEV1 curves were similar in shape and maximum effect, but the estimated free IgE concentration leading to 50% maximum effect (IC50) in pediatric patients (39.4, 95% confidence interval [CI] 24.3-63.9 ng/mL) was higher than estimated in adults (19.8, 95% CI 15.1-24.5 ng/mL). The model further confirmed that the current omalizumab dosing rationale based on the mean target free IgE level of 25 ng/ml was appropriate. The pediatric model can be used to predict population FEV1 response for omalizumab when combined with an omalizumab pharmacokinetic-IgE model.

Machine Learning for Exposure-Response Analysis: Methodological Considerations and Confirmation of Their Importance via Computational Experimentations

Exposure-response (E-R) is a key aspect of pharmacometrics analysis that supports drug dose selection. Currently, there is a lack of understanding of the technical considerations necessary for drawing unbiased estimates from data. Due to recent advances in machine learning (ML) explainability methods, ML has garnered significant interest for causal inference. To this end, we used simulated datasets with known E-R "ground truth" to generate a set of good practices for the development of ML models required to avoid introducing biases when performing causal inference. These practices include the use of causal diagrams to enable the careful consideration of model variables by which to obtain desired E-R relationship insights, keeping a strict separation of data for model-training and for inference generation to avoid biases, hyperparameter tuning to improve the reliability of models, and estimating proper confidence intervals around inferences using a bootstrap sampling with replacement strategy. We computationally confirm the benefits of the proposed ML workflow by using a simulated dataset with nonlinear and non-monotonic exposure-response relationships.

Population Pharmacokinetics of Tenecteplase in Patients With Acute Myocardial Infarction and Application to Patients With Acute Ischemic Stroke

The pharmacokinetics (PK) of tenecteplase in patients with acute ischemic stroke has not been extensively studied. This study aimed to describe PK characteristics of tenecteplase in patients with acute myocardial infarction (AMI) using a population PK approach and to assess applicability of the findings to patients with acute ischemic stroke by means of external validation. A population PK model was developed using nonlinear mixed-effects modeling based on the phase II TIMI 10B study in patients with AMI (785 PK observations from 103 patients). The statistical and clinical impact of selected covariates on PK parameters were evaluated by a stepwise covariate modeling procedure and simulations, respectively. The performance of the final model was evaluated for patients with acute ischemic stroke using summary statistics of tenecteplase concentrations of 75 patients from investigator-initiated study N1811s. Tenecteplase PK was well described by a 2-compartment linear model, incorporating allometric scaling of clearance and volume parameters and weight-normalized creatinine clearance on clearance. Simulations showed that the identified covariates (weight and creatinine clearance) were of limited influence on exposure at the intended dosing regimen for patients with acute ischemic stroke. The model overpredicted mean tenecteplase plasma concentrations from N1811s by 39%, but 72% of the distribution from N1811s was within the 90% prediction interval of the model predictions. The PK characteristics of tenecteplase in patients with AMI were well described by the final model. Simulations from the model indicated that no specific dose recommendations based on covariates are warranted for patients with AMI.

Exposure‐response relationships of etrolizumab in patients with moderately‐to‐severely active ulcerative colitis

Etrolizumab is an IgG1‐humanized monoclonal anti‐β7 integrin antibody. Phase III trials with induction and/or maintenance phases were conducted in patients with moderately‐to‐severely active ulcerative colitis (UC) who were either previously treated with tumor necrosis factor (TNF) inhibitors (HICKORY) or were TNF inhibitor naïve (HIBISCUS I/II, LAUREL, and GARDENIA). A total of eight exposure‐response analyses were conducted for two clinical outcomes (remission and endoscopic improvement) at the end of induction for studies HIBISCUS I/II (combined) and HICKORY and at the end of maintenance for studies HICKORY and LAUREL. Trough concentration at week 4 (C_trough,wk4) of induction was selected as the exposure metric. Exposure‐response (ER) modeling was conducted using logistic regression. A full covariate model was used to examine the impact of covariates on clinical outcomes. Linear models with a single intercept for placebo and active treatments adequately described the data for all eight analyses. The etrolizumab exposure‐response slope was significant (p < 0.05) for seven of the eight analyses. Baseline Mayo Clinic Score (MCS) was the only statistically significant covariate that impacted induction remission and endoscopic improvement. No statistically significant covariate was identified to impact maintenance outcomes except for baseline fecal calprotectin on endoscopic improvement for LAUREL study. A statistically significant positive ER relationship was identified for most of the clinical outcomes tested, reflecting a better treatment effect in patients with UC with higher etrolizumab C_trough,wk4 of induction. Baseline MCS was the only other significant covariate impacting induction efficacy. Besides C_trough,wk4 of induction, no consistent covariate was identified to impact maintenance efficacy.

Precision dosing of intravenous busulfan in pediatric hematopoietic stem cell transplantation: Results from a multicenter population pharmacokinetic study

Busulfan (Bu) is a common component of conditioning regimens before hematopoietic stem cell transplantation (HSCT) and is known for high interpatient pharmacokinetic (PK) variability. This study aimed to develop and externally validate a multicentric, population PK (PopPK) model for intravenous Bu in pediatric patients before HSCT to first study the influence of glutathione‐s‐transferase A1 (GSTA1) polymorphisms on Bu's PK in a large multicentric pediatric population while accounting for fludarabine (Flu) coadministration and, second, to establish an individualized, model‐based, first‐dose recommendation for intravenous Bu that can be widely used in pediatric patients. The model was built using data from 302 patients from five transplantation centers who received a Bu‐based conditioning regimen. External model validation used data from 100 patients. The relationship between body weight and Bu clearance (CL) was best described by an age‐dependent allometric scaling of a body weight model. A stepwise covariate analysis identified Day 1 of Bu conditioning, GSTA1 metabolic groups based on GSTA1 polymorphisms, and Flu coadministration as significant covariates influencing Bu CL. The final model adequately predicted Bu first‐dose CL in the external cohort, with 81% of predicted area under the curves within the therapeutic window. The final model showed minimal bias (mean prediction error, −0.5%; 95% confidence interval [CI], −3.1% to 2.0%) and acceptable precision (mean absolute prediction error percentage, 18.7%; 95% CI, 17.0%–20.5%) in Bu CL prediction for dosing. This multicentric PopPK study confirmed the influence of GSTA1 polymorphisms and Flu coadministration on Bu CL. The developed model accurately predicted Bu CL and first doses in an external cohort of pediatric patients.

Population pharmacokinetic and exposure-response analysis of apremilast in Japanese subjects with moderate to severe psoriasis

Apremilast is an orally available phosphodiesterase 4 inhibitor used for the treatment of moderate to severe psoriasis. The aims of this analysis were to develop a population pharmacokinetic (PPK) model of apremilast based on observed data from phase 1 studies combined with clinical trial data from subjects with moderate to severe psoriasis, and to develop exposure–response (E‐R) models to determine whether Japanese subjects with moderate to severe psoriasis achieve response to apremilast treatment similar to that observed in non‐Japanese, predominantly Caucasian subjects with moderate to severe psoriasis. The PPK model demonstrated that apremilast plasma concentrations and overall apparent clearance rate were comparable between the Japanese and Caucasian subgroups. The E‐R analyses of ≥75% or ≥50% improvement from baseline in Psoriasis Area and Severity Index score and achievement of static Physician Global Assessment score of 0 (clear) or 1 (almost clear) at week 16 indicated that apremilast treatment in Japanese subjects approached the maximal effect with response rates comparable to those in predominantly Caucasian subjects. Overall, the analyses confirm that the approved apremilast 30 mg b.i.d. dose is appropriate for Japanese subjects with moderate to severe psoriasis, with an efficacy profile similar to that previously observed in Caucasian subjects.

Population Pharmacokinetics and Exposure-Response Relationship of Luspatercept, an Erythroid Maturation Agent, in Anemic Patients With β-Thalassemia

β‐Thalassemia is an inherited blood disorder resulting from defects in hemoglobin production, leading to premature death of red blood cells (RBCs) or their precursors. Patients with transfusion‐dependent β‐thalassemia often need lifelong regular RBC transfusions to maintain adequate hemoglobin levels. Frequent transfusions may lead to iron overload and organ damage. Thus, there is a large unmet need for alternative therapies. Luspatercept, a first‐in‐class erythroid maturation agent, is the first approved therapy in the United States for the treatment of anemia in adult patients with β‐thalassemia who require regular RBC transfusions. The population pharmacokinetics and exposure‐response relationship of luspatercept were evaluated in 285 patients with β‐thalassemia. Luspatercept displayed linear and time‐invariant pharmacokinetics when administered subcutaneously once every 3 weeks. Body weight was the only clinically relevant covariate of luspatercept clearance, favoring weight‐based dosing. Magnitude and frequency of hemoglobin increase, if not influenced by RBC transfusions, was positively correlated with luspatercept area under the serum concentration‐time curve (AUC), 0.2‐1.25 mg/kg, whereas a significant reduction in RBC units transfused was observed in frequently transfused patients. The probability of achieving ≥33% or ≥50% reduction in RBC transfusion burden was similar across the time‐averaged AUC (0.6‐1.25 mg/kg), with the 1 mg/kg starting dose sufficient for most early responders (71%‐80%). Increasing luspatercept AUC (0.2‐1.25 mg/kg) did not increase incidence or severity of treatment‐emergent adverse events. These results provide a positive benefit‐risk profile for the recommended luspatercept doses (1‐1.25 mg/kg) in treating adult patients with β‐thalassemia who require regular RBC transfusions.

Population Pharmacokinetics and Exposure-Response of Luspatercept, an Erythroid Maturation Agent, in Anemic Patients With Myelodysplastic Syndromes

Luspatercept is a recombinant fusion protein that enhances late‐stage erythroid maturation. This report describes the population pharmacokinetics and exposure–response relationship of luspatercept in 260 patients with anemia due to myelodysplastic syndromes. Luspatercept displayed linear and time‐invariant pharmacokinetics over a dose range of 0.125–1.75 mg/kg administered subcutaneously once every 3 weeks. Body weight was the only clinically relevant covariate of luspatercept exposure, supporting the weight‐based dosing. The probability of achieving transfusion independence ≥ 8 weeks increased with time‐averaged luspatercept serum exposure, reaching the plateau at doses 1.0–1.75 mg/kg. The probability of achieving multiple efficacy end points increased with slower luspatercept clearance, independent of effects of luspatercept exposure or disease characteristics. The probability of experiencing severe treatment‐emergent adverse events decreased with increasing luspatercept exposure, especially during long‐term treatment. These results provide a positive benefit–risk profile for the titration‐to‐response dose regimen (1.0–1.75 mg/kg) recommended for this population.

Population Pharmacokinetics and Exposure-Response Analysis of nab-Paclitaxel in Pediatric Patients With Recurrent or Refractory Solid Tumors

Pediatric malignancies are most commonly of primary central nervous system or hematopoietic origin. The main reason for cancer death in pediatrics is refractory and relapsed disease, and improved therapeutic options are needed in the pediatric population. Nanoparticle albumin-bound (nab)-paclitaxel (Abraxane) is a human albumin-stabilized formulation of paclitaxel and was designed to improve the chemotherapeutic effects of paclitaxel and to reduce toxicities. Although nab-paclitaxel pharmacokinetics (PK) has been extensively studied in adults, no information is available on its PK in children. ABI-007-PST-001 was the first nab-paclitaxel clinical trial conducted in pediatrics, and the current analysis is the first study of nab-paclitaxel PK in pediatrics. Our analyses suggested that ontogeny and maturation play a role in nab-paclitaxel PK disposition, as demonstrated by the finding that both blood clearance and volume of distribution increased from younger to older pediatric age groups and from pediatrics to adults. A 3-compartment population PK (PPK) model with saturable elimination was developed to describe the paclitaxel whole blood concentrations in pediatrics. The PPK model was customized by estimating the allometric function on PK parameters to take into account the ontogeny/maturation of patients. PPK estimates are consistent with the fast and deep distribution of paclitaxel that was previously observed in adults. Finally, the exposure-safety analysis showed an increased probability of drug-related adverse events (>grade 2) in cycle 1 and the first cycle of neutropenia (>grade 2) associated with higher doses. However, there is no statistically significant association between exposures (measured by area under the concentration-time curve) and the probabilities of either safety event.

Pharmacokinetics, Pharmacodynamics, and Exposure-Response of Lanadelumab for Hereditary Angioedema

Hereditary angioedema (HAE) with C1 inhibitor deficiency is a rare disorder characterized by unpredictable, potentially life‐threatening recurrent angioedema attacks. Lanadelumab is a fully human monoclonal antibody with selective binding to active plasma kallikrein, and prevents the formation of cleaved high molecular weight kininogen (cHMWK) and bradykinin, thereby preventing HAE attacks. The clinical pharmacology of lanadelumab was characterized following subcutaneous administration in 257 subjects (24 healthy subjects and 233 patients with HAE). The pharmacokinetics of lanadelumab were described using a one‐compartment model with first‐order rate of absorption and linear clearance, showing slow absorption and a long half‐life (14.8 days). A covariate analysis retained body weight and health status on apparent clearance (CL/F) and body weight on volume of distribution (V/F). Population estimates of CL/F and V/F were 0.0249 L/hour (0.586 L/day) and 12.8 L, respectively. An indirect‐response Imax model showed 53.7% maximum suppression in cHMWK formation with a low potential for interactions with concomitant medications (analgesic, anti‐inflammatory, and antirheumatic medications). A 300 mg dose administered Q2W was associated with a mean steady‐state minimum concentration (C_min,ss; 25.4 μg/mL) that was ~ 4.5‐fold higher than the half‐maximal inhibitory concentration for cHMWK reduction (5.71 μg/mL). Exposure‐response analyses suggest that 300 mg Q2W dosing was associated with a significantly reduced HAE attack rate, prolonged time to first attack after treatment initiation, and lower need for concomitant medications. The response was comparable across patient body weight groups. Findings from this analysis support the dosing rationale for lanadelumab to prevent attacks in patients with HAE.

Population Pharmacokinetics and Exposure Response Analysis of Pomalidomide in Subjects With Relapsed or Refractory Multiple Myeloma From the Novel Combination Treatment of Pomalidomide, Bortezomib, and Low-Dose Dexamethasone

Multiple myeloma is an incurable progressive neoplastic disease that accounts for 10% of all hematologic malignancies. Even though significant progress has been made in the treatment of newly diagnosed multiple myeloma, the disease follows a relapsing course in the majority of patients, and there is a need for more effective therapeutic options for the treatment of relapsed or refractory multiple myeloma. CC-4047-MM-005 and CC-4047-MM-007 were phase 1 and 3 studies to evaluate the novel combination of pomalidomide, bortezomib, and low-dose dexamethasone for the treatment of patients with relapsed or refractory multiple myeloma who have already received lenalidomide-based treatments early. This analysis was performed to characterize the population pharmacokinetics (PK) of pomalidomide from the combination treatment and to examine exposure-response relationships. Our analysis showed that pomalidomide concentration-time profiles from the combination treatment were adequately described with a 1-compartment PK model, with first-order absorption and elimination and pomalidomide exhibiting linear and time-invariant PK with moderate variability from the combination treatment. Except for the body surface area, none of the tested covariates had an effect on pomalidomide PK. Although body surface area was identified as a statistically significant covariate of pomalidomide PK, the impact was not deemed clinically relevant. A flat exposure-response curve was observed, consistent with a near-saturated drug effect at the tested exposure range suggesting an appropriately recommended clinical dose of 4 mg of pomalidomide for the combination treatment. Finally, pomalidomide exposure was not associated with higher probabilities of dose interruption during cycle 1 or dose reduction during the treatment period.

Population Pharmacokinetics and Pharmacodynamics of Ontamalimab (SHP647), a Fully Human Monoclonal Antibody Against Mucosal Addressin Cell Adhesion Molecule-1 (MAdCAM-1), in Patients With Ulcerative Colitis or Crohn's Disease

Ontamalimab (SHP647) is a fully human, immunoglobulin G₂, antihuman mucosal addressin cell adhesion molecule‐1 (MAdCAM‐1) monoclonal antibody being developed for the treatment of ulcerative colitis (UC) and Crohn's disease (CD). A population pharmacokinetic/pharmacodynamic (PK/PD) analysis was conducted using clinical phase 2 study data to evaluate the PK and PD of ontamalimab following subcutaneous administrations of 7.5, 22.5, 75, and 225 mg every 4 weeks in patients with moderate to severe UC or CD. A total of 440 patients with UC (n = 249; 56.6%) or CD (n = 191; 43.4%) were included in the analysis. A 2‐compartment model with parallel linear and nonlinear elimination adequately characterized concentration‐time profiles of ontamalimab. The apparent clearance and volume of distribution were 0.0127 L/h (0.305 L/day) and 6.53 L, respectively. Apparent clearance and volume of distribution were mainly dependent on baseline albumin and body weight, respectively. No differences in the PK properties of ontamalimab were observed between patients with UC or CD. The presence of antidrug antibodies did not impact the PK of ontamalimab. Nonlinear elimination occurred at very low concentrations and was unlikely to contribute to the elimination half‐life under steady‐state conditions. A linear PK/PD model described the relationship between ontamalimab and free MAdCAM‐1. Minimum concentrations of ontamalimab at steady state following 75 mg every 4 weeks were associated with >95% suppression of circulating free MAdCAM‐1. The PK/PD properties characterized support phase 3 testing in UC and CD.

Incorporation of GSTA1 genetic variations into a population pharmacokinetic model for IV busulfan in paediatric hematopoietic stem cell transplantation

AIMS

The aim of this study is to develop a population pharmacokinetic (PopPK) model for intravenous busulfan in children that incorporates variants of GSTA1, gene coding for the main enzyme in busulfan metabolism.

METHODS

Busulfan concentration-time data was collected from 112 children and adolescents (median 5.4 years old, range: 0.1-20) who received intravenous busulfan during the conditioning regimen prior to stem cell transplantation. Weight, sex, baseline disease (malignant vs. non-malignant), age, conditioning regimen and GSTA1 diplotypes were evaluated as covariates of pharmacokinetic parameters by using nonlinear mixed effects analysis. The ability to achieve the target AUC_24h (3600-6000 μM min^-1 ) was assessed by estimating the first dose based on the present PopPK model and by comparing the results with other available models in children.

RESULTS

A one-compartment model with first-order elimination best described the data. Allometric scaling of weight and a factor of busulfan metabolism maturation were included in the base model. GSTA1 diplotypes were found to be a significant covariate of busulfan clearance, which was 7% faster in rapid metabolizers and 12% slower in poor metabolizers, in comparison with normal ones. Busulfan doses calculated using the parameters of the proposed PopPK model were estimated to achieve the target AUC in 85.2% of the cases (95% CI 78.7-91.7%).

CONCLUSION

This is the first PopPK for busulfan that successfully incorporated GSTA1 genotype in a paediatric population. Its use may contribute to better prediction of busulfan exposure in children and adolescents since the first dose, by tailoring the dose according to the individual metabolic capacity.

Dose-Exposure Simulation for Piperacillin-Tazobactam Dosing Strategies in Infants and Young Children

BACKGROUND

Extended piperacillin-tazobactam (TZP) infusions have been associated with favourable outcomes. There are currently no pediatric dosing recommendations.

OBJECTIVES

To determine appropriate TZP dosing strategies in children 2 months - 6 years according to age and different minimal inhibitory concentrations (MICs).

METHODS

Age and weight were simulated for 1000 children. Post-hoc pharmacokinetic parameter estimates were generated using published clearance and volume of distribution data. For different dosing regimens, we estimated the probability of target attainment (PTA) over a range of MICs from 4 to 128 mg/L. The pharmacodynamic (PD) target was defined as free piperacillin concentrations above the MIC for ≥ 50% of the dosing interval. A PTA ≥ 90% was defined as optimal.

RESULTS

PTA decreased as MIC and age increased. In all age groups, standard dosing regimens (240-300 mg/kg/day, 0.5h infusions) failed to reach PTAs ≥ 90% at MICs ≥ 16 mg/L. Standard 0.5h infusions reached PTAs ≥ 90% at MICs up to 8 mg/L in infants > 2 to 6m. No 0.5h infusion reached PTAs ≥ 90% for MICs ≥ 4 mg/L in children > 6m. While none of the tested regimens were optimal at MICs > 16 mg/L in children > 6m, 100 mg/kg/dose every 6h as a 3h infusion reached PD target at MICs of 32 mg/L in infants > 2 to 6m.

CONCLUSIONS

Up to MICs of 16 mg/L, 90 mg/kg/dose every 8h as a 2h infusion in infants > 2 to 6m and 100 mg/kg/dose every 8h as a 4h infusion in children > 6m-6y achieved PTAs ≥ 90%.

The Utility of a Population Approach in Drug-Drug Interaction Assessments: A Simulation Evaluation

This study aims at evaluating the utility of the population pharmacokinetics approach in therapeutic protein drug-drug-interaction (DDI) assessment. Simulations were conducted for 2 representative victim drugs, methotrexate and trastuzumab, using a parallel-group design with and without the interaction drug. The effect of a perpetrator on the exposure of the victim drug is described as the ratio of clearance/apparent clearance of the victim drug given with or without the perpetrator. The power of DDI assessment was calculated as the percentage of runs with 90% confidence interval of the estimated DDI effect within 80% to 125% for the scenarios of no DDI, benchmarked with the noncompartmental approach with intensive sampling. The impact of the number of subjects, the number of sampling points per subject, sampling time error, and model misspecification on the power of DDI determination were evaluated. Results showed that with equal numbers of subjects in each arm, the population pharmacokinetics approach with sparse sampling may need about the same or a higher number of subjects compared to a noncompartmental approach in order to achieve similar power. Increasing the number of subjects, even if only in the study drug alone arm, can increase the power. Sampling or dosing time error had notable impacts on the power for methotrexate but not for trastuzumab. Model misspecification had no notable impacts on the power for trastuzumab. Overall, the population pharmacokinetics approach with sparse sampling built in phase 2/3 studies allows appropriate DDI assessment with adequate study design and analysis and can be considered as an alternative to dedicated DDI studies.

Population pharmacokinetics and exposure-response of trametinib, a MEK inhibitor, in patients with BRAF V600 mutation-positive melanoma

PURPOSE

To characterize the pharmacokinetics of oral trametinib, a first in class MEK inhibitor, identify covariates, and describe the relationship between exposure and clinical effects in patients with BRAF V600 metastatic melanoma.

EXPERIMENTAL DESIGN

Trametinib concentrations obtained in three clinical studies were included in the population pharmacokinetic analysis. Trametinib 2 mg once daily was administered in the Phase 2 and 3 studies. The impact of exposure [trough (C min) or average concentration] on response rates and progression-free survival (PFS) was examined.

RESULTS

Plasma concentrations (n = 3120) obtained in 493 patients were described using a two-compartment model. Trametinib oral clearance was lower in women relative to men (1.26-fold) and increased with body weight. There was no significant effect of age, mild or moderate renal impairment, or mild hepatic impairment on oral clearance. Between-subject variability was low (24 %). The number of responders was consistent across median exposure range, although tended to be lower at trough concentration <10 ng/mL. Disease stage was found to be a significant predictor of response with a lower response rate in patients with disease stage of M1c. Lactate dehydrogenase was significant in the analysis of PFS. Patients with observed C min above the median had longer PFS than those below median based on Phase 2 study (median 10.6 ng/mL), while the effect of exposure was not statistically significant in the Phase 3 study (median 13.6 ng/mL).

CONCLUSIONS

No dosage adjustments are required with any of the covariates tested. Clinical efficacy was associated with trametinib trough concentrations greater than 10 ng/mL.

Population pharmacokinetics and Bayesian estimation of tacrolimus exposure in paediatric liver transplant recipients

AIMS

The objectives of this study were to develop a population pharmacokinetic (PopPK) model for tacrolimus in paediatric liver transplant patients and determine optimal sampling strategies to estimate tacrolimus exposure accurately.

METHODS

Twelve hour intensive pharmacokinetic profiles from 30 patients (age 0.4-18.4 years) receiving tacrolimus orally were analysed. The PopPK model explored the following covariates: weight, age, sex, type of transplant, age of liver donor, liver function tests, albumin, haematocrit, drug interactions, drug formulation and time post-transplantation. Optimal sampling strategies were developed and validated with jackknife.

RESULTS

A two-compartment model with first-order absorption and elimination and lag time described the data. Weight was included on all pharmacokinetic parameters. Typical apparent clearance and central volume of distribution were 12.1 l h(-1) and 31.3 l, respectively. The PopPK approach led to the development of optimal sampling strategies, which allowed estimation of tacrolimus pharmacokinetics and area under the concentration–time curve (AUC) on the basis of practical sampling schedules (three or four sampling times within 4 h) with clinically acceptable prediction error limit. The mean bias and precision of the Bayesian vs. reference (trapezoidal) AUCs ranged from -2.8 to -1.9% and from 7.4 to 12.5%, respectively.

CONCLUSIONS

The PopPK of tacrolimus and empirical Bayesian estimates represent an accurate and convenient method to predict tacrolimus AUC(0-12) in paediatric liver transplant recipients, despite high between-subject variability in pharmacokinetics and patient demographics. The developed optimal sampling strategies will allow the undertaking of prospective trials to define the tacrolimus AUC-based therapeutic window and dosing guidelines in this population.

Conversion from twice- to once-daily tacrolimus in pediatric kidney recipients: a pharmacokinetic and bioequivalence study

BACKGROUND

The objectives of this study were to investigate pharmacokinetic and pharmacogenetic parameters during the conversion on a 1:1 (mg:mg) basis from a twice-daily (Prograf) to once-daily (Advagraf) tacrolimus formulation in pediatric kidney transplant recipients.

METHODS

Twenty-four-hour pharmacokinetic profiles were analyzed before and after conversion in 19 stable renal transplant recipients (age 7-19 years). Tacrolimus pharmacokinetic parameters [area under the concentration-time curve (AUC0-24), minimum whole-blood concentration (Cmin), maximum whole-blood concentration (Cmax), and time to achieve maximum whole-blood concentration (tmax)] were compared between Tac formulations and between CYP3A5 and MDR1 genotypes after dose normalization.

RESULTS

Both AUC0-24 and Cmin decreased after conversion (223.3 to 197.5 ng.h/ml and 6.5 to 5.6 ng/ml; p = 0.03 and 0.01, respectively). However, the ratio of the least square means (LSM) for AUC0-24 was 90.8 %, with 90 % CI limits of 85.3 to 96.7 %, falling within bioequivalence limits. The CYP3A5 genotype influences the dose-normalized Cmin with the twice-daily formulation only.

CONCLUSIONS

Both tacrolimus formulations are bioequivalent in pediatric renal recipients. However, we observed a decrease in AUC0-24 and Cmin after the conversion, requiring close pharmacokinetic monitoring during the conversion period.

Pharmacokinetic and pharmacodynamic modeling to determine the dose of ST-246 to protect against smallpox in humans

Although smallpox has been eradicated, the United States government considers it a "material threat" and has funded the discovery and development of potential therapeutic compounds. As reported here, the human efficacious dose for one of these compounds, ST-246, was determined using efficacy studies in nonhuman primates (NHPs), together with pharmacokinetic and pharmacodynamic analysis that predicted the appropriate dose and exposure levels to provide therapeutic benefit in humans. The efficacy analysis combined the data from studies conducted at three separate facilities that evaluated treatment following infection with a closely related virus, monkeypox virus (MPXV), in a total of 96 NHPs. The effect of infection on ST-246 pharmacokinetics in NHPs was applied to humans using population pharmacokinetic models. Exposure at the selected human dose of 600 mg is more than 4-fold higher than the lowest efficacious dose in NHPs and is predicted to provide protection to more than 95% of the population.

Safety and population pharmacokinetic analysis of intravenous acetaminophen in neonates, infants, children, and adolescents with pain or Fever

OBJECTIVES

The administration of acetaminophen via the oral and rectal routes may be contraindicated in specific clinical settings. Intravenous administration provides an alternative route for fever reduction and analgesia. This phase 1 study of intravenous acetaminophen (Ofirmev, Cadence Pharmaceuticals, Inc., San Diego, CA) in inpatient pediatric patients with pain or fever requiring intravenous therapy was designed to assess the safety and pharmacokinetics of repeated doses over 48 hours.

METHODS

Neonates (full-term to 28 days) received either 12.5 mg/kg every 6 hours or 15 mg/kg every 8 hours. Infants (29 days to <2 years), children (2 to <12 years) and adolescents (≥12 years) received either 12.5 mg/kg every 4 hours or 15 mg/kg every 6 hours. Both noncompartmental and population nonlinear mixed-effects modeling approaches were used. Urinary metabolite data were analyzed, and safety and tolerability were assessed.

RESULTS

Pharmacokinetic parameters of acetaminophen were estimated using a two-compartment disposition model with weight allometrically expressed on clearances and central and peripheral volumes of distribution (Vds). Postnatal age, with a maturation function, was a significant covariate on clearance. Total systemic normalized clearance was 18.4 L/hr per 70 kg, with a plateau reached at approximately 2 years. Total central and peripheral Vds of acetaminophen were 16 and 59.5 L/70 kg, respectively. The drug was well tolerated based on the incidence of adverse events. The primary and minor pathways of elimination were acetaminophen glucuronidation, sulfation, and glutathione conjugate metabolites across all age groups.

CONCLUSIONS

Intravenous acetaminophen in infants, children, and adolescents was well tolerated and achieved plasma concentrations similar to those achieved with labeled 15 mg/kg body weight doses by oral or rectal administration.