Population pharmacokinetics–pharmacodynamics of sunitinib in pediatric patients with solid tumors

Purpose The safety profile of sunitinib in children, including the impact of sunitinib exposure on safety endpoints, was assessed using population pharmacokinetic (PK) and pharmacokinetic–pharmacodynamic (PK–PD) models. Methods Data were from two clinical studies in 59 children with solid tumors (age range 2–21 years, 28 male/31 female, body weight range 16.2–100 kg, body surface are [BSA] range 0.7–2.1 m2). Analysis of covariates that affected PK and PD parameters was conducted using a nonlinear mixed-effects model. Safety and tolerability endpoints were absolute neutrophil count, hepatic transaminases, diastolic blood pressure, hemoglobin, lymphocyte count, platelet count, white blood cell count, hand-foot syndrome, fatigue, nausea, intracranial hemorrhage, and vomiting. Results The models well described the time courses of concentrations of sunitinib and its primary active metabolite SU012662, as well as safety and tolerability endpoints. In PK models for sunitinib and SU012662, BSA was the only covariate that statistically significantly affected apparent clearance (CL/F) and apparent central volume of distribution (Vc/F). Higher BSA was associated with greater CL/F and Vc/F. No statistically significant covariates were identified in the PK–PD models. For safety endpoints that had a sufficient number of adverse events, a higher probability of adverse events was associated with higher average plasma sunitinib concentrations. Conclusion In PK models, BSA was the only covariate that affected major PK parameters of sunitinib and SU012662. Based on analysis of safety and tolerability endpoints, the PK–PD relationships were mainly driven by sunitinib plasma exposures and were not affected by age, sex, respective baseline safety endpoint values, baseline Eastern Cooperative Oncology Group performance status, or body size. Trial registration ClinicalTrials.gov: NCT00387920 (registered October 13, 2006), NCT01462695 (registered October 31, 2011). Electronic supplementary material The online version of this article (10.1007/s00280-020-04106-z) contains supplementary material, which is available to authorized users.

Population pharmacokinetic (PK) models have been previously developed in healthy adults or adults with solid tumors to assess the PK of sunitinib and its active metabolite SU012662 that examined covariates that might explain the variability in exposure of sunitinib and SU012662 and to make predictions on their efficacy and safety [13][14][15][16]. Based on regulatory proposals that promote the use of model-informed drug development in pediatric patients [17,18], population PK and pharmacokinetic-pharmacodynamic (PK-PD) models in these patients may also help to explore further the safety profile of sunitinib and, in particular, the impact of sunitinib exposure on safety endpoints. The objectives of the current study were to develop a population PK model for sunitinib and SU012662, identify covariates that account for the inter-individual variability in the PK of sunitinib and SU012662, and develop sequential PK-PD and/or exposure-response models with respect to key safety and tolerability endpoints using sunitinib PK model post hoc predictions.

Patient population
Data were pooled from two Phase I/II clinical trials in pediatric patients with solid tumors (studies ADVL0612 and ACNS1021), including predominantly high-grade glioma, ependymoma, brain stem glioma, or sarcoma. Details of these trials have been published previously [8][9][10]. The original trials were approved by the institutional review boards at all participating centers and The National Cancer Institute Pediatric Central institutional review board also approved study ACNS1021. All participants in the original trials or their parent/legal guardian signed a document of informed consent and assent was obtained as appropriate. A detailed description of the eligibility and exclusion criteria can be found in the Online Resource Methods. Patients were between 18 months and 22 years of age and received sunitinib at a starting dose of 15 or 20 mg/m 2 on a schedule of 4 weeks on treatment followed by 2 weeks off treatment (schedule 4/2) [8][9][10]. In total, data for analysis were available from 59 patients (ADVL0612, n = 35; ACNS1021, n = 24).

Bioanalytical methods
Sample collections (1.5 mL) for PK analysis were completed at pre-specified visits. For study ADVL0612, the series of samples for PK analysis was conducted at 0, 1, 2, 4, 6, 8-10, and 24-48 h post-dose on day 1 of cycle one and pre-dose (trough) assessments were made on days 7, 14, 21, and 28 of cycle one. For study ACNS1021, the series of samples for PK analysis were taken at 2, 4, 6-8, and 24 ± 1 h post-dose on day 1 of cycle one. In addition, pre-dose (trough) assessments were made on days 7, 14, 21, and 28 of cycle one, and days 1 and 28 of cycle two. The majority of the safety assessment data that were to be used for PK-PD modeling were captured at each study visit, but some safety assessments were performed less frequently per the specific requirements of each study protocol. Plasma samples were analyzed for the determination of sunitinib and SU012662 concentrations using a sensitive, specific, and validated liquid chromatography with tandem mass spectrometry assay (BASi, West Lafayette, IN), as previously published [8][9][10]. For study ADVL0612, calibration standard responses met acceptance criteria over the range of 1-200 ng/mL for sunitinib and 1-100 ng/mL for SU012662, using a quadratic weighted (l/ concentration 2 ) regression. The lower limit of quantification (LLOQ) for both sunitinib and SU012662 was 1 ng/mL. The between-day assay accuracy, expressed as percent relative error, for quality control (QC) concentrations, ranged from − 3.3-1.3% for sunitinib and − 1-3.8% for SU012662 for the low, medium, high, and diluted QC samples. Assay precision, expressed as the between-day percent coefficient of variation of the mean estimated concentrations of QC samples, was ≤ 5.1% for the low (3.00 ng/mL), medium (100 ng/ , and diluted (9.00 ng/mL after dilution) concentrations of sunitinib, and was ≤ 13.7% for the low (0.300 ng/ mL), medium (9.00 ng/mL), high (15.0 ng/mL), and diluted (3.00 ng/mL after dilution) concentrations of SU012662.

Model development
The PK and PK-PD modeling approaches were described in detail in the analysis plan before initiating the analyses.
Regulatory guidance and quality control were taken into account. A systematic multistep approach to model development consisted of base model development, random effects model development, full model development, final model development, assessment of model adequacy (goodnessof-fit), and assessment of model predictive performance (validation). Analysis was conducted using nonlinear mixed-effects modeling methodology, as implemented in NONMEM (v7.1.2, University of California at San Francisco, CA), and using the first-order conditional estimation method with interaction to estimate all parameters. During model development, the goodness-of-fit of different models to the data were evaluated using change in the objective function, visual inspection of different scatterplots, precision of the parameter estimates, and decreases in inter-individual and residual variability. The base model consisted of a two-compartment model with first-order absorption and lag time (t lag ) to fit sunitinib and SU012662 concentrations. The type of base models used for the PK-PD modeling portion were transit compartments in series with feedback loop models or indirect response models, as used previously [16]. Sequential PK-PD models for the safety and tolerability endpoints were built using the final PK model-predicted sunitinib concentrations. SU012662 data were not included in this process because a previous study showed that inclusion of the predicted SU012662 concentrations did not improve the model objective function value (OFV) [16]. During the modeling portion, safety endpoints were used when there was a sufficient number of patients with PK and PD data (i.e., at least two patients with both PK and PD data, either both baseline and post-baseline for continuous safety endpoints, or post-baseline for categorical safety endpoints). The safety and tolerability endpoints assessed were related to the most common safety events observed with sunitinib use [5,6], which were absolute neutrophil count (ANC), alanine transaminase (ALT), aspartate transaminase (AST), diastolic blood pressure (BP), hemoglobin, lymphocyte count, platelet count, white blood cell (WBC) count, hand-foot syndrome, fatigue, nausea, intracranial hemorrhage, and vomiting. For categorical safety endpoints, PK-PD modeling was not used. Instead, the relationships between the average daily plasma exposures (i.e., the accumulated dose divided by the time after the first dose, divided by the individual sunitinib clearance × 1000) up to the time of earliest worst grade Common Terminology Criteria for Adverse Events (CTCAE) and the incidence rate were explored by looking at different CTCAE grade incidence rates at the lower and upper halves for the mean plasma exposures, as well as using ordered logistic regression approaches, where there was a sufficient number of events (i.e., more than one event). Interindividual variability in the PK or PK-PD parameters was modeled using multiplicative exponential random effects of the form θ i = θ × e ηi , where θ is the typical or central value of the  Based on prior experience, a group of potential covariates was examined with respect to the PK and PK-PD parameters (Table 1). Covariates were tested for significance in a stepwise manner using a stepwise covariate model (SCM) building procedure and statistical criteria of α = 0.01 for the forward inclusion step. The full model was then subjected to a backward elimination step with a statistical criterion of α = 0.001. In the SCM approach, linear and power functions were tested for any continuous covariates on PK parameters, with the SCM model reporting the best one selected (i.e., the one with the lower OFV). Baseline body weight and body surface area (BSA) were examined to ensure that the body size measure with the largest effect was included in the final model, and other measures related to body size that were highly correlated with body weight or BSA, such as lean body weight, body mass index, and height, were not included in the covariate analysis. BSA was calculated based on the method of DuBois and DuBois [i.e., BSA = 0.20247 × height (m) 0.725 × weight (kg) 0.425 ]. No adaptation was required because no patients had undergone prior amputation. In addition, for both sunitinib and SU012662, baseline bodyweight was tested in a separate SCM run where BSA was replaced with bodyweight and the final model with the lower OFV was selected.

Model validation
For the validation of the base and final models, visual predictive check (VPC) techniques comprising 1000 simulations were carried out, and the median and upper and lower bounds of the 95% prediction interval (PI) for PK or PD profiles were compared against the observed data median and confidence intervals (CIs). The number of observations not within the 95% PI was to remain within 5% of the total number of observations, and the mean prediction profile was expected to follow the observed mean profile. In additon, the 90% PIs for the median and the lower and upper bounds of the 95% PIs were identified to ensure that they included the observed median and 95% CI bounds. Bootstrapping techniques (1000 boostrap datasets) were applied to generate the nominal 95% CIs around the point estimates and to confirm the 95% CIs generated by the base or final models based on asymptotic standard errors generated from the NONMEM covariance step.

Results
The final dataset comprised 365 sunitinib and 340 SU012662 post-baseline measurable plasma observations from the 59 patients who received sunitinib [57 of 59 patients (96.6%) provided ≥ 2 samples for PK analysis]. Twenty-four of 389 (6.2%) samples were below the limit of quantification for sunitinib, and 49 of 389 (12.6%) for SU012662, and were not included in the analysis. Baseline patient characteristics overall and by age group (2-5, 6-11, 12-17, and 18-21 years) are shown in Table 2. Online Resource Fig. S1 shows the concentration-time profiles for each study using time after first dose.

Sunitinib and SU012662 base and final PK models
The sunitinib and SU012662 base models comprised a twocompartment model with first-order absorption and elimination rates, consistent with previous work. Based on previous observations [13], a conversion of 21% of sunitinib to SU012662 was assumed to bring the magnitude of the parameters for SU012662 to a more physiologically relevant level. In both the sunitinib and SU012662 base models, the diagnostic plots were satisfactory. The effects of extreme outliers on the population PK parameter estimates and on the diagnostic plots were tested in both models, and no outliers met the criteria for exclusion from the datasets [i.e., |conditional weighted residual (CWRES) |> 6]. Bootstrap results were consistent with the population parameter estimates, indicating that the base models were stable and the parameter estimates represented the final datasets adequately.
Using SCM, the effects of different covariates on apparent clearance (CL/F) and apparent central volume of distribution (Vc/F) in both models were examined. For sunitinib, the estimated typical values for CL/F and Vc/F were 24.1 L/h and 1070 L, respectively, and for SU012662 the values were estimated to be 10.9 L/h and 1030 L, respectively. For sunitinib, the effect of BSA on CL/F was statistically significant (P < 0.001) using a linear function: CL/F = 24.1 L/h · [1 + 0.557 (BSA-1.47)], and for SU012662, the effect of BSA on CL/F was statistically significant (P < 0.001) using a power function: CL/F = 10.9 L/h · (BSA/1.47) 0.843 . For sunitinib, the effect of BSA on Vc/F was also statistically significant (P < 0.001) using a power function: Vc/F = 1070 L · (BSA/1.47) 1.47 , and for SU012662, the effect of BSA was Fig. 1 Final model prediction and variance-corrected visual predictive check plots for sunitinib and SU012662 plasma concentrations, and safety and tolerability endpoints, up to 3000 h post-dose. Visual predictive check plots for a sunitinib, b SU012662, c ANC, d ALT, e AST, f diastolic BP, g hemoglobin, h lymphocyte count, i platelet count, and j WBC count. Blue circles represent the observed data. Red lines represent the median (solid line) and 2.5th and 97.5th percentiles (dashed lines) of the observed data. Black lines represent the median (solid line) and 2.5th and 97.5th percentiles (dashed lines) of the simulated data. The 95% CIs for simulated median and each percentile are shown by pink and blue shaded areas, respectively. ALT alanine transaminase, ANC absolute neutrophil count, AST aspartate transaminase, BP blood pressure, WBC white blood cell ◂ also statistically significant on Vc/F (P < 0.001) using a power function: Vc/F = 1030 L · (BSA/1.47) 1.72 . Therefore, for sunitinib and SU012662, higher BSA was associated with greater CL/F and Vc/F. Baseline bodyweight was also tested in a separate SCM run where BSA was replaced with bodyweight. However, the final model with BSA was the one with the lower OFV and hence was selected. The effects on PK parameters of the other covariates tested (Table 1) were not found to be statistically significant for sunitinib or SU012662.
Goodness-of-fit diagnostic plots for the observed versus individual-predicted or population-predicted values, as well as CWRES versus time or population-predicted values, in the final PK models for sunitinib and SU012662 plasma concentrations demonstrated the adequacy of the models (Online Resource Fig. S2A and S2B, S3A and S3B, S4A and S4B, and S5A and S5B). Moreover, the prediction-and variance-corrected VPC plots of sunitinib and SU012662 plasma concentration demonstrated the similarity between the predicted and the observed data during the first 12 h post-dose (Online Resource Fig. S6) and up to 3000 h postdose (Fig. 1a, b). A summary of sunitinib and SU012662 PK parameters from the final models and following bootstrapping are shown in Table 3. Bootstrap results were consistent with the population parameter estimates, indicating that the final models were stable and that the population parameter estimates represented the dataset adequately. For sunitinib, inclusion of baseline BSA as a covariate into the final model reduced interindividual variability on CL/F and Vc/F by 25.9 and 78.0%, respectively, and for SU012662, interindividual variability on CL/F and Vc/F was reduced by 28.7 and 53.9%, respectively.

PK-PD modeling of safety and tolerability endpoints
PK-PD models were based solely on sunitinib data, as described previously [16]. For ANC, hemoglobin, lymphocyte count, platelet count, and WBC count, a sequential transit compartments in series with feedback loop PK-PD model was used as the base model. For ANC, lymphocyte count, platelet count, and WBC count, the model included a maximum sunitinib effect (E max ) on the proliferation rate constant (k prol ) in the stem cell compartment, whereas for hemoglobin, the model included an effect first-order rate constant (k PD ) effect on the k prol constant in the stem cell compartment. For ALT, AST, and diastolic BP, a sequential indirect response PK-PD model was used as the base model. For ALT and AST, the models included a k PD effect on the output/elimination rate constant (k out ), and for diastolic BP,  the model included a k PD effect on the input rate constant (k in ). Each of these models has been used previously to describe the time course of each endpoint following sunitinib dosing [16]. Following SCM analysis for each endpoint, no covariates (Table 1) were found to statistically significantly (P > 0.001) affect the concentration at half maximum effect (EC 50 ) or k PD ; therefore, the final model used was the base model. Goodness-of-fit plots for the observed values versus individual-predicted or population-predicted values, as well as CRWES values versus time or predicted values, in the final PK-PD models for the continuous safety endpoints (ANC, ALT, AST, diastolic BP, hemoglobin, lymphocyte count, platelet count, and WBC count) demonstrated the adequacy of the models (Online Resource Fig. S2C-J, S3C-J, S4C-J, and S5C-J). Moreover, the predictionand variance-corrected VPC plots of each safety endpoint demonstrated the similarity between the predicted and the observed data (Fig. 1c-j), further supporting that the PK-PD models described the time course of each of the safety and tolerability endpoints adequately. Summaries of the PK-PD parameters from each of the final models for the continuous safety endpoints are shown in Table 4. For each of these endpoints, the mean from the model bootstrap runs was consistent with those from the model run, indicating that the population parameter estimates from the models adequately represented the datasets.

Discussion
In this population analysis of pooled data from pediatric patients with solid tumors administered sunitinib at 15-20 mg/m 2 on schedule 4/2, the PK of sunitinib and SU012662 were well described using a two-compartment model with first-order absorption and t lag , using nonlinear mixed-effect modeling approaches. Covariate analysis identified BSA as the only statistically significant (P < 0.001) covariate for CL/F and Vc/F in the final models for both sunitinib and SU012662. Higher BSA was associated with greater CL/F and Vc/F, meaning lower sunitinib and SU012662 exposure. As part of the PK-PD analyses, the PK-PD models well described the time course of each safety and tolerability endpoint. Furthermore, no covariates were identified as statistically significant on either EC 50 or k PD , indicating that the PK-PD relationship did not appear to be affected by body size or any other baseline characteristic, including age, sex, respective baseline safety endpoint values, or baseline Eastern Cooperative Oncology group performance status. Use of exposure-response logistic regression analysis showed that for each of the safety endpoints that had sufficient number of events to conduct and achieve a successful analysis, there was a higher probability of events with higher average sunitinib plasma concentrations. For hemoglobin decreased, hand-foot syndrome, and intracranial hemorrhage, there were too few events to conduct or achieve a successful logistic regression analysis.
To the best of our knowledge, this is the first population PK-PD analysis of sunitinib in a pediatric patient group. Therefore, these results cannot be directly compared with previous findings in other pediatric populations. However, several similar studies have been conducted in adult patient populations. In a study by Houk et al. [14] in adults with advanced solid tumors, including metastatic renal cell carcinoma (RCC) and gastrointestinal stromal tumor (GIST), there was a higher probability of fatigue and diastolic hypertension, and lower ANC, with greater sunitinib exposure (25-150 mg/day or every other day). In a second study by Lindauer et al. [15] in healthy adult volunteers, PK-PD models were successfully built for systolic and diastolic BP. In a third study by Khosravan et al. [16] in adult patients with advanced RCC or GIST, similar safety PK-PD relationships were identified as those shown in the current study in pediatric patients.
This study had some limitations. The dataset was relatively small, which sometimes led to large η-shrinkage values in EC 50 and k PD for some of the safety and tolerability endpoints. Therefore, failure to identify any statistically significant covariates on EC 50 and k PD should be interpreted with caution. In addition, the small sample size and the sampling schedule may have contributed to the wide range of bootstrap results for sunitinib Vp/F and SU012662 Q/F in the final PK models. In the PK-PD modeling and exposure-response analysis portions, there were too few events for some safety endpoints to explore fully the relationship between the incidence of AEs and sunitinib exposure. Finally, patients in the original clinical trials (studies ADVL0612 and ACNS1021) predominantly had high-grade glioma, ependymoma, brain stem glioma, or sarcoma [8][9][10]. Therefore, no patients in the current study had indications for which sunitinib is approved (e.g., RCC or GIST) albeit in adults [5,6].
In conclusion, the PK of sunitinib and SU012662 was described using a two-compartment model with first-order absorption and t lag , and BSA was the only covariate that statistically significantly affected CL/F and Vc/F. Higher BSA was associated with lower sunitinib and SU012662 exposure. No covariates statistically significantly affected parameters in the PK-PD models. For each of the safety endpoints with a sufficient incidence of events to conduct or achieve a successful analysis/model, there was a higher probability of AEs with higher average plasma sunitinib concentrations. The exposure-response relationships of safety endpoints of sunitinib in pediatric patients with solid tumors were mainly driven by sunitinib plasma exposures and were not affected by age, sex, respective baseline safety endpoint values, baseline Eastern Cooperative Oncology Group performance status, or body size. Acknowledgements Pfizer would like to thank the Children's Oncology Group for conducting studies ADVL0612 and ACNS1021 and for allowing Pfizer to use their data. Medical writing support was provided by David Cope, PhD, of Engage Scientific Solutions and funded by Pfizer.
Author contributions Conception and design-EW, RK; development of methodology-EW, RK; acquisition of data-EW, SGD, CW, RK; analysis and interpretation of data-EW, RK; writing, review, and/or revision of the manuscript-EW, SGD, CW, RK; clinical involvement in the original trials-SGD, CW.
Funding This study was sponsored by Pfizer.
Data availability Upon request, and subject to certain criteria, conditions and exceptions (see https ://www.pfize r.com/scien ce/clini cal-trial s/ trial -data-and-resul ts for more information), Pfizer will provide access to individual de-identified participant data from Pfizer-sponsored global interventional clinical studies conducted for medicines, vaccines and medical devices (1) for indications that have been approved in the US and/or EU or (2) in programs that have been terminated (i.e., development for all indications has been discontinued). Pfizer will also consider requests for the protocol, data dictionary, and statistical analysis plan. Data may be requested from Pfizer trials 24 months after study completion. The de-identified participant data will be made available to researchers whose proposals meet the research criteria and other conditions, and for which an exception does not apply, via a secure portal. To gain access, data requestors must enter into a data access agreement with Pfizer.

Compliance with ethical standards
Conflict of interest Erjian Wang and Reza Khosravan are employees of Pfizer and hold stock or stock options with Pfizer. Steven G. DuBois has received prior travel expenses from Loxo Oncology, Roche and Salarius, and prior consulting fees from Loxo Oncology. Cynthia Wetmore has received research support and prior travel expenses from Eli Lilly & Co.
Ethical approval The original trials were approved by the institutional review boards at all participating centers. The National Cancer Institute Pediatric Central institutional review board also approved study ACNS1021.
Consent to participate All participants in the original trials or their parent/legal guardian signed a document of informed consent and assent was obtained as appropriate.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.