
==== Front
Cancer Chemother Pharmacol
Cancer Chemother Pharmacol
Cancer Chemotherapy and Pharmacology
0344-5704
1432-0843
Springer Berlin Heidelberg Berlin/Heidelberg

38878207
4679
10.1007/s00280-024-04679-z
Original Article
Population pharmacokinetics of TLD-1, a novel liposomal doxorubicin, in a phase I trial
Mc Laughlin Anna M. 12
Hess Dagmar 3
Michelet Robin 1
Colombo Ilaria 4
Haefliger Simon 5
Bastian Sara 6
Rabaglio Manuela 5
Schwitter Michael 7
Fischer Stefanie 3
Eckhardt Katrin 8
Hayoz Stefanie 8
Kopp Christoph 8
Klose Marian 12
Sessa Cristiana 4
Stathis Anastasios 910
Halbherr Stefan 11
Huisinga Wilhelm 12
Joerger Markus markus.joerger@kssg.ch

3
Kloft Charlotte charlotte.kloft@fu-berlin.de

1
1 https://ror.org/046ak2485 grid.14095.39 0000 0001 2185 5786 Department of Clinical Pharmacy and Biochemistry, Institute of Pharmacy, Freie Universitaet Berlin, Kelchstr. 31, 12169 Berlin, Germany
2 https://ror.org/03bnmw459 grid.11348.3f 0000 0001 0942 1117 Graduate Research Training Program PharMetrX, Freie Universitaet Berlin/University of Potsdam, Berlin/Potsdam, Germany
3 https://ror.org/00gpmb873 grid.413349.8 0000 0001 2294 4705 Department of Medical Oncology and Haematology, Cantonal Hospital St. Gallen, Rorschacher Strasse 95, 9007 St. Gallen, Switzerland
4 Department of Medical Oncology, Oncology Institute of Southern Switzerland, EOC, Bellinzona, Switzerland
5 grid.5734.5 0000 0001 0726 5157 Department of Medical Oncology, Inselspital Bern University Hospital, University of Bern, Bern, Switzerland
6 https://ror.org/04wpn1218 grid.452286.f 0000 0004 0511 3514 Department of Medical Oncology, Kantonsspital Graubünden, Chur, Switzerland
7 https://ror.org/04wpn1218 grid.452286.f 0000 0004 0511 3514 Oncology/Hematology, Kantonsspital Graubünden, Chur, Switzerland
8 https://ror.org/04rtrpb08 grid.476782.8 0000 0001 1955 3199 Coordinating Center, Swiss Group for Clinical Cancer Research, Bern, Switzerland
9 Oncology Institute of Southern Switzerland, EOC, Bellinzona, Switzerland
10 https://ror.org/03c4atk17 grid.29078.34 0000 0001 2203 2861 Faculty of Biomedical Sciences, Universita della Svizzera Italiana, Lugano, Switzerland
11 Innomedica Switzerland AG, Bern, Switzerland
12 https://ror.org/03bnmw459 grid.11348.3f 0000 0001 0942 1117 Institute of Mathematics, University of Potsdam, Potsdam, Germany
15 6 2024
15 6 2024
2024
94 3 349360
4 2 2024
18 5 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
Study objectives

TLD-1 is a novel pegylated liposomal doxorubicin (PLD) formulation aiming to optimise the PLD efficacy-toxicity ratio. We aimed to characterise TLD-1’s population pharmacokinetics using non-compartmental analysis and nonlinear mixed-effects modelling.

Methods

The PK of TLD-1 was analysed by performing a non-compartmental analysis of longitudinal doxorubicin plasma concentration measurements obtained from a clinical trial in 30 patients with advanced solid tumours across a 4.5-fold dose range. Furthermore, a joint parent-metabolite PK model of doxorubicinentrapped, doxorubicinfree, and metabolite doxorubicinol was developed. Interindividual and interoccasion variability around the typical PK parameters and potential covariates to explain parts of this variability were explored.

Results

Medians ± standard deviations of dose-normalised doxorubicinentrapped+free Cmax and AUC0−∞ were 0.342 ± 0.134 mg/L and 40.1 ± 18.9 mg·h/L, respectively. The median half-life (95 h) was 23.5 h longer than the half-life of currently marketed PLD. The novel joint parent-metabolite model comprised a one-compartment model with linear release (doxorubicinentrapped), a two-compartment model with linear elimination (doxorubicinfree), and a one-compartment model with linear elimination for doxorubicinol. Body surface area on the volumes of distribution for free doxorubicin was the only significant covariate.

Conclusion

The population PK of TLD-1, including its release and main metabolite, were successfully characterised using non-compartmental and compartmental analyses. Based on its long half-life, TLD-1 presents a promising candidate for further clinical development. The PK characteristics form the basis to investigate TLD-1 exposure-response (i.e., clinical efficacy) and exposure-toxicity relationships in the future. Once such relationships have been established, the developed population PK model can be further used in model-informed precision dosing strategies.

Clinical trial registration

ClinicalTrials.gov–NCT03387917–January 2, 2018

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-024-04679-z.

Keywords

Nanoparticles
Liposomes
Doxorubicin
Pharmacokinetics
Pharmacometrics
Nonlinear mixed-effects model
InnoMedica, Swiss State Secretary for Education, Research and Innovation (SERI), Swiss Cancer Research Foundation (SCS), Swiss Cancer League (SCL)Freie Universität Berlin (1008)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Efficacy and safety of doxorubicin

Doxorubicin is a well-established and highly efficacious drug used for the treatment of several tumour entities such as breast cancer, ovarian cancer, Kaposi’s sarcoma, lymphoma, and multiple myeloma [1–3]. The mechanism of action of this anthracycline is a combination of free radical formation, cellular membrane interaction, topoisomerase II inhibition, and DNA intercalation, all ultimately leading to apoptosis [1, 3, 4]. Unfortunately, only a small fraction of doxorubicin accumulates in the tumour, resulting in high drug exposure in healthy tissue and associated toxicities [5]. Degradation of doxorubicin in the blood stream and absorption and accumulation of the formed reactive oxygen species into the heart tissue leads to cardiomyocyte apoptosis and irreversible cardiac damage [2]. This cardiotoxicity of doxorubicin is further amplified by the cardiotoxic effect of doxorubicinol, which is the primary circulating metabolite of doxorubicin [2, 6, 7]. Due to its irreversible cardiac damage, doxorubicin is only given until a pre-determined cumulative lifetime dose is reached, often resulting in early termination of the otherwise efficacious treatment [2, 5].

Benefits, challenges, and pharmacokinetics of pegylated liposomal doxorubicin

Entrapping doxorubicin into polyethylene glycol modified (PEGylated) liposomes (such as in Caelyx® (Europe) /Doxil® (US), hereinafter referred to as Caelyx®) [8] largely reduces the observed cardiotoxicity [5] by increasing drug accumulation and release at the tumour site by approximately 10-fold [9]. Due to their large molecular size, the liposomes exploit the enhanced permeability and retention effect (EPR) [10, 11] by only extravasating into tissue with increased vascular permeability, such as tumours [12]. Thus, this targeted drug delivery strategy strongly increases the antitumour effect while simultaneously reducing the adverse effects in healthy tissue. Caelyx® is approved for the treatment of several tumour entities, such as ovarian cancer, breast cancer, myeloma, and Kaposi-Sarcoma [8]. However, approximately 50% of patients treated with Caelyx® suffer from palmar-plantar erythrodysaesthesia (PPE) [5, 13], a painful inflammation of the palms of the hands and soles of the feet. Furthermore, approximately one in four patients experiences mucositis [5].

The clinical pharmacokinetics (PK) of PEGylated doxorubicin are significantly different from the PK of free doxorubicin [3, 4, 10, 14] and characterised by a small volume of distribution [8, 15], a low clearance, a long half-life of 50–80 h, and an approximately 300-fold higher area under the concentration-time curve (AUC) [10]. The drug effect of PEGylated liposomal doxorubicin is associated with the liposome-released unbound (free) concentration of the parent drug doxorubicin and to a lesser extent with its main metabolite doxorubicinol in the tumour cell; however, the PK of liposome-entrapped doxorubicin and its interplay with the two free species has not been well studied [10, 20]. Due to their long circulation time, it is hypothesised that the doxorubicin-containing liposomes extravasate into and accumulate in the skin at the pressure points of the hands and the feet, and that this accumulation is a key factor in the development of PPE [10, 16, 17].

TLD-1, a novel pegylated liposomal doxorubicin

To maximise the efficacy and minimise the toxicity of PEGylated liposomal doxorubicin, TLD-1, a novel PEGylated liposomal doxorubicin formulation [18], is currently under investigation. Compared to Caelyx®, TLD-1 consists of smaller, uniform, and more stable liposomes with an average diameter of 36 nm (average diameter of Caelyx®: 70 nm) [18]. TLD-1 is being developed to optimise the antitumour activity-toxicity ratio of PEGylated liposomal doxorubicin, and preclinical studies suggest a potential for an improved safety profile, including the lack of PPE in animal models. In the dose escalation part of the first-in-human phase I clinical trial including 12 patients (NCT03387917), grade 2/3 cumulative PPE was observed in four patients [19]. Moreover, two patients experienced grade 2 mucositis and further two patients experienced grade 2 rash [19].

Objectives

In this work, we aimed to characterise the PK of TLD-1 using non-compartmental analysis and nonlinear mixed-effects modeling to jointly evaluate the pharmacokinetics of entrapped doxorubicin, free doxorubicin, and the main metabolite doxorubicinol in patients with advanced solid tumours. Typical PK values along with interindividual and interoccasion variability (i.e., between-patient and between-cycle variability) for the key kinetic processes, such as liposomal release and doxorubicin elimination, were estimated. Additionally, patient characteristics (covariates) influencing pharmacokinetic parameters were explored.

Materials and methods

Clinical study

The open-label, single-arm, multicentre, first-in-human phase I TLD-1 dose-escalation trial SAKK 65/16 (NCT03387917) in patients with advanced solid tumours was conducted at four phase I centres in Switzerland. The compound was administered via intravenous infusion in 21 days cycles for a maximum of 6 cycles for patients previously treated with anthracyclines or 9 cycles for patients previously not treated with anthracycline until disease progression, unacceptable toxicity, or withdrawal of consent. The infusion durations were 60 min for dose levels 1–6 (10–40 mg/m2) and 90 min for dose level 7 (45 mg/m2). Dose escalation followed an accelerated titration design until first occurrence of a dose-limiting toxicity. Afterwards, a continual reassessment method using cohorts of three was applied [21, 22]. Individual doses were based on one of seven dose levels and the individual body surface area (BSA) (Table 1). Upon treatment of 12 patients, dose level 7 was identified as the tentative maximum tolerated dose and nine additional patients were treated at this dose level. Due to several late appearing cumulative toxicities in this expansion cohort, nine additional patients were treated at dose level 6 (40 mg/m2).Table 1 Dose levels 1–7 of TLD-1 used in the first-in-human phase I clinical trial SAKK 65/16

Dose level	Dose [mg] per m2 body surface area	
1	10	
2	16	
3	23	
4	30	
5	35	
6	40	
7	45	

Concentrations of total doxorubicin (doxorubicinentrapped+free), unencapsulated doxorubicin (doxorubicinfree), and the main metabolite doxorubicinol were measured by Swiss BioQuant AG (Reinach, Switzerland) using a validated liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) method as described below. Samples were collected in the first two cycles at the following pre-defined time points: prior to infusion (t = 0), mid-infusion (0.5 h for dose level 1–6 and 0.75 h for dose level 7), end of infusion (1 h for dose levels 1–6 and 1.5 h for dose level 7), as well as, 0.5 h, 1 h, 3 h, 5 h, and 7 h after end of infusion, at 24 h, 48 h (only in cycle 1), 168 h (day 8), and 336 h (day 15).

Bioanalytical method and assay performance

For the analysis of doxorubicinfree, doxorubicinol, and doxorubicinentrapped+free, two validated bioanalytical assays were used.

Sample preparation for the quantification of doxorubicinfree and doxorubicinol was performed in an ice-bath and under light protected conditions. To an aliquot of 50 µL matrix, 100 µL of PBS buffer with 1.0% BSA containing the internal standard was added. After gentle shaking and storage for 5 min, the samples were transferred to an ultrafiltration tube and filtrated for 15 min at approximately 10,000g. The temperature of the centrifuge was set to 8 ℉C. An aliquot of 50 µL of the filtrate was transferred to an Eppendorf tube and precipitated with 100 µL of acetonitrile. After vortex mixing, the samples were centrifuged for 5 min at approximately 50,000 g. The temperature of the centrifuge was set to 8 °C. An aliquot of the supernatant was transferred to an autosampler vial for subsequent HPLC-MS/MS analysis.

Sample preparation for the quantification of doxorubicinentrapped+free was done in an ice-bath and under light protected conditions. To an aliquot of 50.0 µL human plasma, 200 µL of acetonitrile containing the internal standard were added. After mixing, the samples were centrifuged for 10 min at approximately 50,000 g. The temperature of the centrifuge was set to 8 °C. An aliquot of the supernatant was transferred to an autosampler tube for subsequent HPLC-MS/MS analysis.

For the quantification of doxorubicinentrapped+free, sample analysis was done by column separation using reversed-phase liquid chromatography followed by detection with triple-stage quadrupole MS/MS in the selected reaction monitoring mode. Chromatography was performed by gradient elution using acidified water and acetonitrile (ACN; 5–95%).

For the quantification of doxorubicinfree and doxorubicinol, sample analysis was done by column separation using reversed-phase liquid chromatography followed by detection with triple-stage quadrupole MS/MS in the selected reaction monitoring mode. Chromatography was performed by gradient elution using acidified water and acetonitrile (ACN; 5–95%). On-line solid phase extraction with a reversed-phase trapping column was used to further purify and concentrate the sample prior to MS/MS quantification.

For doxorubicinfree, the descriptive statistics of the QC of the batch in the calibration range of 2.00–2000 ng/mL showed an inter-batch precision of 4.8–10.9%, whereas the inter-batch accuracy was in the range of 99.0–101.3% of the nominal concentration.

For doxorubicinentrapped+free, the descriptive statistics of the QC of the batch in the calibration range of 20.0–20,000 ng/mL showed an inter-batch precision of 3.2–6.3% whereas the inter-batch accuracy was in the range of 100.0–104.3% of the nominal concentration.

For doxorubicinol, the descriptive statistics of the QC of the batch in the calibration range of 0.500–500 ng/mL showed an inter-batch precision of 5.3–12.1%, whereas the inter-batch accuracy was in the range of 98.0–99.0% of the nominal concentration.

Analysis dataset generation

For each PK sampling timepoint, the concentration of doxorubicinentrapped was calculated by subtracting the measured concentration of doxorubicinfree from the concentration of doxorubicinentrapped+free. In total, 1870 concentration measurements were available (n = 624 for total doxorubicin, n = 623 for doxorubicinfree, and n = 623 for metabolite doxorubicinol). For the one sample with missing doxorubicinfree and doxorubicinol concentrations, it was also not possible to derive the corresponding doxorubicinentrapped concentration by subtracting the concentration of doxorubicinfree from the concentration of doxorubicinentrapped+free. Concentration measurements below the lower limit of quantification (BLOQ) (n = 314, all doxorubicinol) were removed from the analysis dataset. Based on in-house data showing that usually > 99% of doxorubicin was entrapped in liposomes in the final product, it was assumed that 100% of the doxorubicin was entrapped at time of infusion. Patient characteristics age, body weight, body height, and BSA, and clinical chemistry parameters serum creatinine, alanine aminotransferase, aspartate amino transferase, alkaline phosphatase, bilirubin, alkaline phosphatase, serum creatinine and the estimated glomerular filtration rate according to the CKD-EPI formula [23] were included in the dataset and available for testing as potential covariates during model development. Moreover, additional body size descriptors lean body weight [24] and body mass index (BMI) were calculated and available for testing as potential covariates.

Non-compartmental analysis

Dose-normalised maximum concentration (Cmax), area under the concentration-time curve from t = 0 until infinity (AUC0-∞), and the terminal half-life of total doxorubicin (doxorubicinentrapped+free) were calculated for every patient and cycle using R Statistical Software [25] and the package pkr. Dose proportionality was assessed by inspecting dose-normalised Cmax and AUC0−∞ vs. dose for trends. Due to the unbalanced number of patients in each dose level, no additional statistical tests were performed to assess dose proportionality [26], however, the possibility of a nonlinear clearance was further explored in the nonlinear mixed-effects analysis. A possible cycle-dependent clearance was investigated by inspecting the ratios of individual dose-normalised AUC0−∞ in cycles 1 and 2 for trends.

Nonlinear mixed-effects pharmacokinetic model

A nonlinear mixed-effects modelling approach [27, 28] was chosen to develop a joint parent (entrapped-free)-metabolite model for TLD-1 consisting of three submodels: the structural submodel, characterising the typical concentration-time profile of doxorubicinentrapped, doxorubicinfree, and metabolite doxorubicinol using ordinary differential equations, the stochastic submodel, characterising different levels of variability around PK parameters of the structural model and concentrations, and the covariate submodel, aiming to identify patient characteristics that explain parts of this variability.

Model development was performed using NONMEM® Version 7.4 (Icon Development Solutions, Ellicott City, MD, USA). Model parameters were estimated using the First-Order Conditional Estimation with Interaction (FOCE+I) algorithm implemented in NONMEM. Relative standard errors (RSE) were obtained using the $COVARIANCE function. For the final model, sampling importance resampling (SIR) was additionally performed to assess parameter precision [29, 30]. SIR was chosen over the more common non-parametric bootstrap procedure since non-parametric bootstraps can result in incorrect confidence intervals when applied to small and heterogeneous datasets [31]. In contrast to the non-parametric bootstrap, SIR does not rely on resampling individuals in new datasets and is therefore more suitable for smaller datasets [31, 32].

Structural submodel

Based on the reported small volume of distribution of doxorubicinentrapped approximating the plasma volume [8] and previous models for PEGylated liposomal doxorubicin [15, 33, 34], a one-compartment model was assumed for doxorubicinentrapped. Several elimination pathways, using (i) linear [15, 34], (ii) nonlinear [8], and (iii) parallel linear and nonlinear processes were investigated for doxorubicinentrapped to characterise the release of doxorubicin from the liposomes and the degradation of the liposomes by the reticuloendothelial system (RES) [33, 34]. Based on previously published PK models for doxorubicinfree, two- [4, 34] and three- [35–37] compartment models were investigated. Linear and nonlinear formation processes from doxorubicinfree to the metabolite doxorubicinol were investigated and, based on previous PK analyses, one- [4, 35, 36] and two- [37] compartment models were investigated for doxorubicinol. The respective best submodel was selected for each doxorubicin species based on PK parameter value plausibility, model fit, and parameter precision (RSE ≤ 30%).

Stochastic submodel

Interindividual variability (IIV) parameters were investigated on structural PK parameters using exponential functions and only retained if the inclusion improved model fit (as indicated by a decrease in the objective function value, an indicator of model fit, and in improved goodness-of-fit plots), parameter precision was adequate (RSE ≤ 50%), and inclusion did not lead to model overparameterisation as indicated by a condition number > 1000 [38]. Correlations between parameters were assessed and, if the correlation exceeded 0.8, a “shared-η approach” was used [28]. Interoccasion variability (IOV), for which the start of each cycle represented a new occasion, was next investigated for all parameters using exponential functions. Final IOV parameter selection followed the same criteria as the IIV parameter selection. To characterise the residual unexplained variability (RUV), additive, proportional, combined additive proportional, and log-transformed both sides approaches were investigated.

Covariate submodel

Potential patient characteristics to be implemented as covariates in the model were pre-selected based on plausibility, previous reports [36, 39], and availability in the dataset. Exploratory graphical analyses were used to assess potential trends between covariate values and PK parameter estimates. Continuous covariates (e.g., body surface area) were normalised to the respective median value of the study population and implemented using power relationships. Categorical covariates (e.g., sex) were implemented using fractional change models. Stepwise covariate modelling [40] using significance criteria of changes in the objective function value of 3.84 (α = 0.05, df = 1) for the forward inclusion and 7.88 (α = 0.005, df = 1) for the backward elimination were applied for final covariate selection. Furthermore, only precisely estimated (RSE: ≤ 30%) covariate effects were retained in the final model.

PK model evaluation

The final population PK model was evaluated using goodness-of-fit plots and plots showing observed concentrations overlaid with typical and individual model predictions. Systematic bias was assessed by plotting conditional weighted residuals (CWRES) vs. time and typical predictions. To assess predictive performance, a prediction-corrected visual predictive check (pcVPC, n = 2000 simulations) [41] was performed.

Results

Clinical data

Among the 30 patients included in the trial, the most frequent tumour types were breast cancer (43.3%), ovarian cancer (20.0%), and gastrointestinal cancer (3.3%). The remaining 33.3% of patients had other solid tumour types. In line with the high frequency of breast cancer and ovarian cancer, 80% of the patients were female. Median age and BSA were 67.5 years (range: 38–83 years) and 1.75 m2 (range: 1.44–2.44 m2), respectively. Median body mass index (BMI) and lean body weight were 24.7 kg/m2 (range: 16.5–42.2 kg/m2) and 48.5 kg (range: 38.7–77.6 kg), respectively. Clinical chemistry parameters were mostly within the normal range. The raw concentration-time data of all patients are shown in Supplementary Fig. 1.

Non-compartmental analysis results

For the non-compartmental analysis, two PK concentration-time profiles were not available, as no second cycle dose had been administered in two individuals (one in dose level 6 and one in dose level 7). Median dose-normalised AUC0-∞ and Cmax of doxorubicinentrapped+free were 40.1 h/L (range: 16.8–79.5 h/L) and 0.342 L−1 (range: 0.196– 0.859 L−1), respectively. The median half-life of doxorubicinentrapped+free was 95 h with a large variability (range: 46–213 h). There was no clear relationship of increasing dose-normalised AUC0−∞ vs. dose and thus no clear sign of dose nonlinearity (Fig. 1).Fig. 1 Dose-normalised AUC0-inf (panal a, left) and dose-normalised Cmax (panel b, right) vs. individual doses in the two first cycles of the 30 patients treated with TLD-1. AUC0-inf AUC from t = 0 to infinity; Cmax maximum concentration

However, due to the limited interpretability based on the unbalanced number of patients in each dose level, a potential nonlinear clearance was considered during the following compartmental model development. AUC0-∞ and half-life data for both cycles were available for 28 of the 30 patients. There was no trend for a cycle-dependent clearance, as indicated by the ratios of dose-normalised AUC0−∞ in cycle 2/cycle 1 being randomly scattered around 1 (Supplementary Fig. 2).

Nonlinear mixed-effects pharmacokinetic model

Structural model

In the joint population PK model, the PK of doxorubicinentrapped was best characterised by a one-compartment structural model with central volume of distribution V1 (3.39 L, RSE: 7%) and linear clearance CL1 (0.0271 L/h, RSE: 10%) to the doxorubicinfree compartment. As the clearance parameter for the elimination pathway aiming to capture the removal of liposomes by the RES shrank to zero, it was not retained in the model. For un-encapsulated doxorubicin (doxorubicinfree), a two-compartment model with central volume of distribution V2 (0.531 L at a BSA of 1.75 m2, RSE: 16%), peripheral volume of distribution V3 (61.3 L at a BSA of 1.75 m2, RSE: 19%), intercompartmental clearance Q (0.136 L/h, RSE: 18%), and linear clearance for the metabolism to doxorubicinol CL2 (0.450 L/h, RSE: 11%) best characterised the concentration-time profile. For doxorubicinol, a one-compartmental model with central volume of distribution V4 (8152 L, RSE: 12%) and linear clearance (CL4: 74.6 L/h, RSE: 7%) was sufficient. The final structural model is shown in Fig. 2.Fig. 2 Schematic structure of the joint parent-metabolite PK model of entrapped doxorubicin, free doxorubicin and doxorubicinol. Abbreviations: V1: volume of distribution of entrapped doxorubicin; V2: central volume of distribution of free doxorubicin; V3: peripheral volume of distribution of free doxorubicin; V4: volume of distribution of doxorubicinol; CL1: release clearance of the entrapped doxorubicin; CL2: clearance of free doxorubicin for the metabolism to doxorubicinol; CL4: clearance of doxorubicinol; QDoxo,f: intercompartmental clearance between the central and peripheral compartment of free doxorubicin; BSA: body surface area

Stochastic model

IIV parameters for the structural parameters V1, CL1, CL2, V2, and CL4 fulfilled the requirements for model inclusion and were thus retained in the model. As the correlation between the IIV parameter estimates for V1 and CL1 was high (r = 0.92), a shared-η approach was implemented. IOV was modelled for parameters CL1 (14.4% CV, RSE: 12%), V1 (8.85% CV, RSE: 10%), CL2 (22.4% CV, RSE: 12%) and V2 (126% CV, RSE: 10%) and significantly improved model fit. After implementation of IOV on V2, the IIV parameter estimate for V2 became imprecise and was thus removed from the model without worsening model fit. Implementation of covariances between the other IIV parameters did not improve model fit and was thus not included in the model. The final IIV parameter estimates for CL1, V1, CL2, and CL4 were 45.1% CV (RSE: 11%), 28.2% CV (RSE: 11%), 34.2% CV (RSE 12%), and 15.1% CV (RSE: 10%), respectively. A log-transformed both sides approach with an additive component in the log-domain with separate and uncorrelated parameters for each model species (parameter estimates for doxorubicinentrapped, doxorubicinfree, and doxorubicinol: 19.6% CV (RSE: 5%), 64.2% CV (RSE: 5%), and 65.0% CV (RSE: 7%), respectively) best characterised the residual unexplained variability.

Covariate model

BSA was identified as covariate on the central and peripheral volumes of distribution of doxorubicinfree V2 and V3, respectively (V2_BSA: 4.47, RSE: 19% and V3_BSA: 11.5, RSE: 18%, Supplementary Fig. 3), significantly improving model fit. A trend of increasing volume of distribution of doxorubicinentrapped with increasing BSA was observed as well, however, implementation of this covariate did not significantly improve model fit. Replacing BSA with other body size descriptors, including body weight or lean body weight, did not improve model fit. No other covariates were identified.

Model evaluation

The model predictions captured the observed concentrations well (Fig. 3a, b). Furthermore, no systematic bias was identified as indicated by random distributions of conditional weighted residuals vs. time and conditional weighted residuals vs. typical predictions around zero (Fig. 3c, d). Plots overlaying measured concentrations with typical and individual predictions showed a good concordance of predictions and measured concentrations across dose levels and individuals (Supplementary Fig. 4–6). For some measured concentrations, doxorubicinfree Cmax was underpredicted (Supplementary Fig. 5). The final parameter relative standard errors acquired by SIR were low (≤ 19%, Table 2) and the pcVPC showed good predictive model performance (Fig. 4).Fig. 3 Goodness-of-fit plots of (panel a, upper left) population predictions vs. measured concentrations (“observations”) and (panel b, upper right) individual predictions vs. measured concentrations as well as (panel c, lower left) CWRES vs. Time and (panel d, lower right) CWRES vs. population predictions. Blue lines: trend lines. Red points: entrapped doxorubicin. Blue points: free doxorubicin. Purple points: metabolite doxorubicinol. Abbreviations: CWRES: Conditional weighted residuals

Table 2 Parameter estimates for the final joint parent-metabolite population pharmacokinetic model of liposomal doxorubicin, free doxorubicin and doxorubicinol using the clinical dataset (n = 30 patients)

Parameter [unit]	Parameter description	Estimate	RSE, %	
V1 [L]	Volume of distribution of entrapped doxorubicin	3.39	7	
CL1 [L/h]	Release clearance of entrapped doxorubicin	0.0271	11	
V2 [L]	Baseline volume of distribution of free doxorubicin	0.531	16	
V2_BSA [-]b	Exponent of the power covariate model of BSA on V2	4.47	19	
Q [L/h]	Intercompartmental clearance of free doxorubicin	0.136	18	
V3 [L]	Peripheral volume of distribution of free doxorubicin	61.3	19	
V3_BSA [-]b	Exponent of the power covariate model of BSA on V3	11.5	18	
CL2 [L/h]	Clearance of free doxorubicin to doxorubicinol	0.450	11	
V4 [L]	Volume of distribution of doxorubicinol	8152	12	
CL4 [L/h]	Clearance of doxorubicinol	74.6	7	
θshared	Shared η scale factor for V1	0.643	10	
IIV CL1	Interindividual variability in the release of entrapped doxorubicin	45.1% CV	11	
IIV V1a	Interindividual variability in the volume of distribution of entrapped doxorubicin	28.2% CV	11	
IIV CL2	Interindividual variability in the clearance of free doxorubicin for metabolism to doxorubinol	34.2% CV	11	
IIV CL4	Interindividual variability in the clearance of doxorubinol	15.1% CV	10	
IOV CL1	Interoccasion variability for CL1	14.4% CV	12	
IOV CL2	Interoccasion variability for CL2	22.4% CV	12	
IOV V1	Interoccasion variability for V1	8.85% CV	10	
IOV V2	Interoccasion variability for V2	126% CV	10	
RUV Entrapped doxorubicin	Residual unexplained variability in the observed concentrations of entrapped doxorubicin	19.6% CV	5	
RUV Free doxorubicin	Residual unexplained variability in the observed concentrations of free doxorubicin	64.2% CV	5	
RUV Doxorubicinol	Residual unexplained variability in the observed concentrations of doxorubicinol	65.0% CV	7	
aCalculated using a shared-η approach [27] with ωCL1·θshared2; bimplemented as power covariate model, normalised to the median BSA of 1.75 m2 using the equation V2∙BSA1.75V2_BSA and V3∙BSA1.75V3_BSA, respectively; CV: coefficient of variation; IIV: interindividual variability; IOV: interoccasion variability; RUV: residual unexplained variability; RSE: relative standard error = (standard error/estimate) 100

Fig. 4 Prediction-corrected visual predictive check for the joint parent-metabolite TLD-1 model characterising entrapped doxorubicin (panel a, upper left), free doxorubicin (panel b, upper right), and doxorubicinol (panel c, lower left). Black solid lines: observed median concentrations; black dashed lines: 10th and 90th percentiles of the observed concentrations; blue solid lines: simulated median concentrations; blue dashed lines: 10th and 90th percentiles of the simulated concentrations; blue shaded areas: 95% confidence intervals around the predicted median concentrations, purple shaded areas: 95% confidence intervals around the 10th and the 90th predicted percentiles, respectively. Black open circles: observed concentrations

Discussion

TLD-1 is a novel pegylated liposomal doxorubicin currently in clinical development. Doxorubicin is a widely used antineoplastic drug for the treatment of solid and haematological malignancies, however, its cardiotoxicity is limiting its long-time use [2, 5]. Entrapping it into PEGylated liposomes, such as in Caelyx®, virtually eliminates most of its irreversible cardiotoxicity and favourably modifies its PK characteristics. However, the widespread occurrence of PPE during treatment with Caelyx® warrants further efforts to improve the safety of PEGylated liposomal doxorubicin formulations. Additionally, the PK of pegylated liposomal formulations is poorly understood in detail, since it has never been thoroughly studied in dose-escalation studies spanning a large dose range [10, 20, 42].

In this work, we comprehensively characterised the population PK of TLD-1 by non-compartmental analysis and by developing a novel joint parent-metabolite NLME model comprising the three relevant species doxorubicinentrapped, doxorubicinfree, and metabolite doxorubicinol, based on densely sampled patients over a 4.5-fold dose range. Our model predictions captured the measured concentrations well, however, for 10 patients, the initial peak doxorubicinfree concentrations were underestimated. The reason for this is yet to be elucidated – it could be hypothesised that for these patients, a fraction of doxorubicin had already been released from the liposomes at the time of infusion. This would lead to an unexpectedly high peak in doxorubicinfree under the assumption that all doxorubicin was entrapped at the time of infusion. However, this hypothesis is not supported by in-house manufacturer data showing that > 99% of doxorubicin is entrapped in the liposomes in the drug product. Furthermore, during model development, we attempted to extend our model based on this hypothesis by estimating the fractions of TLD-1 dose being entrapped and free at time of infusion. However, we estimated a fraction of 99% to be entrapped at the time of infusion, supporting independent stability data from the manufacturer. Furthermore, the described model extension did not improve model predictions. As it additionally inflated the model run time and parameter estimate imprecision, we did not include it in the final model by assuming all drug being encapsulated at the time of infusion. Additional exploratory analysis did not reveal any correlations between any of the patient characteristics available in our dataset and the doxorubicinfree concentrations. Further research should thus focus on the occurrence of unexpectedly high initial concentrations of doxorubicinfree. Our model parameter estimates, e.g., the estimate for the volume of distribution of doxorubicinentrapped V1: 3.39 L, were plausible and in line with previously reported volumes of distribution of pegylated liposomal doxorubicin (median: 3.90 L, range: 2.10–10.0 L) [10].

Furthermore, we identified and quantified different levels of variability such as IIV, IOV and RUV. Both IIV and IOV were estimated to be moderate for all parameters (≤42.6% CV) except for the IOV on the volume of distribution for doxorubicinfree V2 (IOV V2: 125% CV). A possible explanation for this high variability between cycles could be the fast elimination of doxorubicin compared to the relatively slow release of doxorubicin from the liposomes (0.450 L/h vs. 0.0271 L/h). Thus, the distribution time of doxorubicinfree between release and elimination is small and a robust estimation of its volume of distribution challenging. The estimated release rate of doxorubicin from the liposomes (TLD-1: CL1: 0.0271 L/h, corresponding to a leakage half-life of 86.7 h at the estimated liposomal volume of distribution of 3.39 L) was lower than the previously published leakage half-life of doxorubicin from Caelyx® liposomes (118.4 h) [17]. In general, a longer leakage half-life has been associated with better efficacy [17]. Interestingly, the median half-life of total doxorubicin was longer in TLD-1 compared to Caelyx® (95 h vs. 71.5 h [8]). The increased half-life could be due to the novel liposome manufacturing process used for TLD-1, ensuring the localisation of PEG only on the outer layer of the liposomes and leading to more effective protection of the liposomes from the RES. Significant clinical correlations between longer half-lives, smaller clearance, and higher dose-normalised AUC with longer survival have been observed in a clinical study on the PK of mitomycin-entrapped liposomes [43]. Thus, a correlation of half-life with efficacy could also be investigated for TLD-1 in the future. Moreover, the available half-life data could aid in the continued development of the compound by optimising the dosing interval(s).

Our thorough PK characterisation of a seven dose level, 4.5-fold range dataset, including the free drug and the main metabolite at the population (typical parameter values) as well as on the individual level (considering the different variability components) allows to explore potential links of plasma concentrations with clinical outcome and toxicity data of the phase I trial [44] next. These future investigations shall focus on exploratory exposure-response and exposure-toxicity relationships and, most importantly, clinical correlations between the estimated PK parameters and outcome. As the IIV for the doxorubicin release efficiency CL1 was moderate (45.1% CV) and the IOV lower (14.4% CV), predicting the probability for treatment success by calculating individual CL1 based on individual PK samples should be investigated in the future. If successful, this could be then applied in individualised dosing [45]. Our model can additionally be used in various model-informed drug development applications, such as in clinical trial simulations or in optimal design analysis, informing future clinical study designs.

In summary, TLD-1 is a new compound aiming to further improve efficacy and reduce toxicity of pegylated liposomal doxorubicin. In this work, we presented the thorough PK data analysis of TLD-1, which is currently in continuing clinical development. The developed joint parent-metabolite NLME model can now be integrated with recently published efficacy and toxicity data [44] to explore potential exposure-response relationships.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (DOCX 2202 KB)

Acknowledgements

The authors thank all the investigators, patients, and their families for their participation in the study. The authors thank the High-Performance Computing Service of ZEDAT at Freie Universitaet Berlin (http://www.zedat.fuberlin.de/Compute) for computing time. An interim analysis of the first 21 patients has been included in the doctoral thesis of A.M.L.

Author contributions

A.M.L, C.K. conceptualised the work. A.M.L. curated the data. A.M.L, R.M.; perfomed the formal analysis. D.H., I.C., S.H., S.B., M.R., M.S., S.F., K.E., S.Hay., C.Ko., C.S., A.S., M.J. performed the investigation. A.M.L, R.M., W.H., M.J., C.K. discussed the methodology. A.M.L, M.K. visualised. D.H., I.C., S.H.., S.B., M.R., M.S., S.F., A.S., M.J., C.K. provided resources. W.H., C.K. supervised; A.M.L. wrote the first draft and all authors reviewed the manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. This trial was supported by InnoMedica and the Swiss State Secretary for Education, Research and Innovation (SERI), Swiss Cancer Research Foundation (SCS) and Swiss Cancer League (SCL).

Data availability

The datasets generated during and/or analysed for the presented study are available from the corresponding authors on reasonable request.

Declarations

Conflict of interest

C.K. and W.H. report grants from an industry consortium (AbbVie Deutschland GmbH & Co. KG, Astra Zeneca Ltd., Boehringer Ingelheim Pharma GmbH & Co. KG, Grünenthal GmbH, F. Hoffmann-La Roche Ltd., Merck KGaA, Novo Nordisk and Sanofi) for the PharMetrX PhD program. C.K. reports grants for the Innovative Medicines Initiative-Joint Undertaking (“DDMoRe”), Diurnal Ltd., the Federal Ministry of Education and Research within the Joint Programming Initiative on Antimicrobial Resistance Initiative (JPIAMR) and from the European Commission within the Horizon 2020 framework program (“FAIR”). A.M.L. is a current employee of Pharmetheus AB and a paid consultant to multiple pharmaceutical companies. I.C. provided advisory/expert opinion for GSK, Novartis, Astra Zeneca, and MSD and recceived travel grants from Tesaro and institutional grants for clinical trials (Principal Investigator): MSD, Bayer, Oasmia. A.T. received institutional research funding from Innomedica, MEI Pharma, Merck, Bayer, Roche, Novartis, Pfizer, ADC Therapeutics, and Eli Lilly, and consulting fees from Bayer, Eli Lilly, Roche, and Novartis. M.J. is investigators in clinical trials for AstraZeneca, Basilea Pharmaceutica, Bayer, BMS, Daiichi Sankyo, Immunophotonics, Innomedia, Janssen, Lilly, MSD, Novartis, Pfizer, Pharmamar, Roche, Sanofi, Takeda, and received travel grants from BSM, Roche, MSD. S.H. provided adivsory/expert opinion for Bayer, Novartis, Lilly, AstraZeneca, and MSD. A.S. received institutional funding for clinical trials for AbbVie, ADC Therapeutics, Amgen, AstraZeneca, Bayer, Cellestia, Incyte, Loxo Oncology, Merck MSD, Novartis, Pfizer, Philogen and Roche, provided paid consultancy services for Debiopharm, Janssen, AstraZeneca, Incyte, Eli Lilly, Novartis, Roche, and Lox Oncology, and received travel grants from Incyte and AstraZeneca. The other authors declare that they have no conflict of interest.

Ethical approval

The institutional review board/ethics committees of participating centers approved the trial. The trial followed the ethical principles of the Declaration of Helsinki, the International Conference on Harmonization Guideline for Good Clinical Practice, and local regulations.

Consent to participate

All patients provided written informed consent before enrollment.

Consent for publication

Not applicable.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Anna M. Mc Laughlin and Dagmar Hess shared first authorship.

Markus Joerger and Charlotte Kloft shared senior authorship.
==== Refs
References

1. Tacar O Sriamornsak P Dass CR Doxorubicin: an update on anticancer molecular action, toxicity and novel drug delivery systems J Pharm Pharmacol 2013 65 157 170 10.1111/j.2042-7158.2012.01567.x 23278683
Tacar O, Sriamornsak P, Dass CR (2013) Doxorubicin: an update on anticancer molecular action, toxicity and novel drug delivery systems. J Pharm Pharmacol 65:157–170. 10.1111/j.2042-7158.2012.01567.x23278683
2. Rahman AM Yusuf SW Ewer MS Anthracycline-induced cardiotoxicity and the cardiac-sparing effect of liposomal formulation Int J Nanomed 2007 2 567 583
Rahman AM, Yusuf SW, Ewer MS (2007) Anthracycline-induced cardiotoxicity and the cardiac-sparing effect of liposomal formulation. Int J Nanomed 2:567–583
3. Speth PAJ van Hoesel QGCM Haanen C Clinical pharmacokinetics of doxorubicin Clin Pharmacokinet 1988 14 287 310 10.2165/00003088-198814050-00002 3293868
Speth PAJ, van Hoesel QGCM, Haanen C (1988) Clinical pharmacokinetics of doxorubicin. Clin Pharmacokinet 14:287–310. 10.2165/00003088-198814050-000023293868
4. Joerger M Huitema ADR Meenhorst PL Pharmacokinetics of low-dose doxorubicin and metabolites in patients with AIDS-related Kaposi sarcoma Cancer Chemother Pharmacol 2005 55 488 496 10.1007/s00280-004-0900-4 15726371
Joerger M, Huitema ADR, Meenhorst PL et al (2005) Pharmacokinetics of low-dose doxorubicin and metabolites in patients with AIDS-related Kaposi sarcoma. Cancer Chemother Pharmacol 55:488–496. 10.1007/s00280-004-0900-415726371
5. O’Brien MER Wigler N Inbar M Reduced cardiotoxicity and comparable efficacy in a phase III trial of pegylated liposomal doxorubicin HCl (CAELYXTM/Doxil®) versus conventional doxorubicin for first-line treatment of metastatic breast cancer Ann Oncol 2004 15 440 449 10.1093/annonc/mdh097 14998846
O’Brien MER, Wigler N, Inbar M et al (2004) Reduced cardiotoxicity and comparable efficacy in a phase III trial of pegylated liposomal doxorubicin HCl (CAELYXTM/Doxil®) versus conventional doxorubicin for first-line treatment of metastatic breast cancer. Ann Oncol 15:440–449. 10.1093/annonc/mdh09714998846
6. Olson RD Mushlin PS Brenner DE Doxorubicin cardiotoxicity may be caused by its metabolite, doxorubicinol Proc Natl Acad Sci USA 1988 85 3585 3589 10.1073/pnas.85.10.3585 2897122
Olson RD, Mushlin PS, Brenner DE et al (1988) Doxorubicin cardiotoxicity may be caused by its metabolite, doxorubicinol. Proc Natl Acad Sci USA 85:3585–3589. 10.1073/pnas.85.10.35852897122
7. Boucek RJ Kunkel EM Graham TP Doxorubicinol, the metabolite of doxorubucin, is more cardiotoxic than doxorubicin Pediatr Res 1987 21 187A 187A 10.1203/00006450-198704010-00127
Boucek RJ, Kunkel EM, Graham TP et al (1987) Doxorubicinol, the metabolite of doxorubucin, is more cardiotoxic than doxorubicin. Pediatr Res 21:187A-187A. 10.1203/00006450-198704010-00127
8. Committee for Medicinal Products for Human Use (CHMP) European Medicines Agency. Caelyx summary of product characteristics. https://www.ema.europa.eu/en/documents/product-information/caelyx-pegylated-liposomal-epar-product-information_en.pdf. Accessed 20 Nov 2020
9. Symon Z Peyser A Tzemach D Selective delivery of doxorubicin to patients with breast carcinoma metastases by stealth liposomes Cancer 1999 86 72 78 10.1002/(SICI)1097-0142(19990701)86:1<72::AID-CNCR12>3.0.CO;2-1 10391566
Symon Z, Peyser A, Tzemach D et al (1999) Selective delivery of doxorubicin to patients with breast carcinoma metastases by stealth liposomes. Cancer 86:72–78. 10.1002/(SICI)1097-0142(19990701)86:1%3c72::AID-CNCR12%3e3.0.CO;2-110391566
10. Gabizon A Shmeeda H Barenholz Y Pharmacokinetics of pegylated liposomal doxorubicin: review of animal and human studies Clin Pharmacokinet 2003 42 419 436 10.2165/00003088-200342050-00002 12739982
Gabizon A, Shmeeda H, Barenholz Y (2003) Pharmacokinetics of pegylated liposomal doxorubicin: review of animal and human studies. Clin Pharmacokinet 42:419–436. 10.2165/00003088-200342050-0000212739982
11. Maeda H Wu J Sawa T Tumor vascular permeability and the EPR effect in macromolecular therapeutics: a review J Control Release 2000 65 271 284 10.1016/S0168-3659(99)00248-5 10699287
Maeda H, Wu J, Sawa T et al (2000) Tumor vascular permeability and the EPR effect in macromolecular therapeutics: a review. J Control Release 65:271–284. 10.1016/S0168-3659(99)00248-510699287
12. Rosenblum D Joshi N Tao W Progress and challenges towards targeted delivery of cancer therapeutics Nat Commun 2018 9 1410 10.1038/s41467-018-03705-y 29650952
Rosenblum D, Joshi N, Tao W et al (2018) Progress and challenges towards targeted delivery of cancer therapeutics. Nat Commun 9:1410. 10.1038/s41467-018-03705-y29650952
13. Gordon AN Fleagle JT Guthrie D Recurrent epithelial ovarian carcinoma: a randomized phase III study of pegylated liposomal doxorubicin versus topotecan J Clin Oncol 2001 19 3312 3322 10.1200/JCO.2001.19.14.3312 11454878
Gordon AN, Fleagle JT, Guthrie D et al (2001) Recurrent epithelial ovarian carcinoma: a randomized phase III study of pegylated liposomal doxorubicin versus topotecan. J Clin Oncol 19:3312–3322. 10.1200/JCO.2001.19.14.331211454878
14. Amantea MA Forrest A Northfelt DW Mamelok R Population pharmacokinetics and pharmacodynamics of pegylated-liposomal doxorubicin in patients with AIDS-related Kaposi’s sarcoma Clin Pharmacol Ther 1997 61 301 311 10.1016/S0009-9236(97)90162-4 9084455
Amantea MA, Forrest A, Northfelt DW, Mamelok R (1997) Population pharmacokinetics and pharmacodynamics of pegylated-liposomal doxorubicin in patients with AIDS-related Kaposi’s sarcoma. Clin Pharmacol Ther 61:301–311. 10.1016/S0009-9236(97)90162-49084455
15. Xu L Wang W Sheng YC Zheng QS Pharmacokinetics and its relation to toxicity of pegylated-liposomal doxorubicin in chinese patients with breast tumours J Clin Pharm Ther 2010 35 593 601 10.1111/j.1365-2710.2009.01128.x 20831683
Xu L, Wang W, Sheng YC, Zheng QS (2010) Pharmacokinetics and its relation to toxicity of pegylated-liposomal doxorubicin in chinese patients with breast tumours. J Clin Pharm Ther 35:593–601. 10.1111/j.1365-2710.2009.01128.x20831683
16. Yokomichi N Nagasawa T Coler-Reilly A Pathogenesis of hand-foot syndrome induced by PEG-modified liposomal doxorubicin Hum Cell 2013 26 8 18 10.1007/s13577-012-0057-0 23386177
Yokomichi N, Nagasawa T, Coler-Reilly A et al (2013) Pathogenesis of hand-foot syndrome induced by PEG-modified liposomal doxorubicin. Hum Cell 26:8–18. 10.1007/s13577-012-0057-023386177
17. Charrois GJR Allen TM Drug release rate influences the pharmacokinetics, biodistribution, therapeutic activity, and toxicity of pegylated liposomal doxorubicin formulations in murine breast cancer Biochim Biophys Acta Biomembr 2004 1663 167 177 10.1016/j.bbamem.2004.03.006
Charrois GJR, Allen TM (2004) Drug release rate influences the pharmacokinetics, biodistribution, therapeutic activity, and toxicity of pegylated liposomal doxorubicin formulations in murine breast cancer. Biochim Biophys Acta Biomembr 1663:167–177. 10.1016/j.bbamem.2004.03.006
18. Innomedica Talidox Brochure. https://relaunch.innomedica.com/wp-content/uploads/2021/04/Talidox.pdf. Accessed 4 Jan 2022
19. Hess D Colombo I Haefliger S 575P TLD-1, a novel liposomal doxorubicin, in patients (pts) with advanced solid tumours: dose escalation part of a multicenter open-label phase I trial (SAKK 65/16) Ann Oncol 2020 31 S490 10.1016/j.annonc.2020.08.689
Hess D, Colombo I, Haefliger S et al (2020) 575P TLD-1, a novel liposomal doxorubicin, in patients (pts) with advanced solid tumours: dose escalation part of a multicenter open-label phase I trial (SAKK 65/16). Ann Oncol 31:S490. 10.1016/j.annonc.2020.08.689
20. Methaneethorn J Tengcharoen K Leelakanok N AlEjielat R Population pharmacokinetics of doxorubicin: a systematic review Asia Pac J Clin Oncol 2023 19 9 26 10.1111/ajco.13776 35415961
Methaneethorn J, Tengcharoen K, Leelakanok N, AlEjielat R (2023) Population pharmacokinetics of doxorubicin: a systematic review. Asia Pac J Clin Oncol 19:9–26. 10.1111/ajco.1377635415961
21. Wheeler GM Mander AP Bedding A How to design a dose-finding study using the continual reassessment method BMC Med Res Methodol 2019 19 1 15 10.1186/s12874-018-0638-z 30611213
Wheeler GM, Mander AP, Bedding A et al (2019) How to design a dose-finding study using the continual reassessment method. BMC Med Res Methodol 19:1–15. 10.1186/s12874-018-0638-z30611213
22. Garrett-Mayer E The continual reassessment method for dose-finding studies: a tutorial Clin Trials 2006 3 57 71 10.1191/1740774506cn134oa 16539090
Garrett-Mayer E (2006) The continual reassessment method for dose-finding studies: a tutorial. Clin Trials 3:57–71. 10.1191/1740774506cn134oa16539090
23. Levey AS Stevens LA Schmid CH A new equation to estimate glomerular filtration rate Ann Intern Med 2009 150 604 612 10.7326/0003-4819-150-9-200905050-00006 19414839
Levey AS, Stevens LA, Schmid CH et al (2009) A new equation to estimate glomerular filtration rate. Ann Intern Med 150:604–612. 10.7326/0003-4819-150-9-200905050-0000619414839
24. Janmahasatian S Duffull SB Ash S Quantification of lean bodyweight Clin Pharmacokinet 2005 44 1051 1065 10.2165/00003088-200544100-00004 16176118
Janmahasatian S, Duffull SB, Ash S et al (2005) Quantification of lean bodyweight. Clin Pharmacokinet 44:1051–1065. 10.2165/00003088-200544100-0000416176118
25. R Core Team (2022) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/.
26. Calvo E Zafar H Goetz A Analysis of dose proportionality testing methods in phase I clinical trials of anticancer agents Cancer Res 2005 65 973 974
Calvo E, Zafar H, Goetz A et al (2005) Analysis of dose proportionality testing methods in phase I clinical trials of anticancer agents. Cancer Res 65:973–974
27. Owen JS Fiedler-Kelly J Introduction to population pharmacokinetic/pharmacodynamic analysis with nonlinear mixed effects models 2014 1 Hoboken John Wiley and Sons Ltd
Owen JS, Fiedler-Kelly J (2014) Introduction to population pharmacokinetic/pharmacodynamic analysis with nonlinear mixed effects models, 1st edn. John Wiley and Sons Ltd, Hoboken
28. Mould DR Upton RN Basic concepts in population modeling, simulation, and model-based drug development–Part 2: introduction to pharmacokinetic modeling methods CPT Pharmacometrics Syst Pharmacol 2013 2 e38 10.1038/psp.2013.14 23887688
Mould DR, Upton RN (2013) Basic concepts in population modeling, simulation, and model-based drug development–part 2: introduction to pharmacokinetic modeling methods. CPT Pharmacometrics Syst Pharmacol 2:e38. 10.1038/psp.2013.1423887688
29. Dosne AG Bergstrand M Karlsson MO An automated sampling importance resampling procedure for estimating parameter uncertainty J Pharmacokinet Pharmacodyn 2017 44 509 520 10.1007/s10928-017-9542-0 28887735
Dosne AG, Bergstrand M, Karlsson MO (2017) An automated sampling importance resampling procedure for estimating parameter uncertainty. J Pharmacokinet Pharmacodyn 44:509–520. 10.1007/s10928-017-9542-028887735
30. Dosne AG Bergstrand M Harling K Karlsson MO Improving the estimation of parameter uncertainty distributions in nonlinear mixed effects models using sampling importance resampling J Pharmacokinet Pharmacodyn 2016 43 583 596 10.1007/s10928-016-9487-8 27730482
Dosne AG, Bergstrand M, Harling K, Karlsson MO (2016) Improving the estimation of parameter uncertainty distributions in nonlinear mixed effects models using sampling importance resampling. J Pharmacokinet Pharmacodyn 43:583–596. 10.1007/s10928-016-9487-827730482
31. Dosne AG Niebecker R Karlsson MO dOFV distributions: a new diagnostic for the adequacy of parameter uncertainty in nonlinear mixed-effects models applied to the bootstrap J Pharmacokinet Pharmacodyn 2016 43 597 608 10.1007/s10928-016-9496-7 27730481
Dosne AG, Niebecker R, Karlsson MO (2016) dOFV distributions: a new diagnostic for the adequacy of parameter uncertainty in nonlinear mixed-effects models applied to the bootstrap. J Pharmacokinet Pharmacodyn 43:597–608. 10.1007/s10928-016-9496-727730481
32. Broeker A Wicha SG Assessing parameter uncertainty in small-n pharmacometric analyses: value of the log-likelihood profiling-based sampling importance resampling (LLP-SIR) technique J Pharmacokinet Pharmacodyn 2020 47 219 228 10.1007/s10928-020-09682-4 32248328
Broeker A, Wicha SG (2020) Assessing parameter uncertainty in small-n pharmacometric analyses: value of the log-likelihood profiling-based sampling importance resampling (LLP-SIR) technique. J Pharmacokinet Pharmacodyn 47:219–228. 10.1007/s10928-020-09682-432248328
33. Harashima H Tsuchihashi M Iida S Pharmacokinetic/pharmacodynamic modeling of antitumor agents encapsulated into liposomes Adv Drug Deliv Rev 1999 40 39 61 10.1016/S0169-409X(99)00039-3 10837779
Harashima H, Tsuchihashi M, Iida S et al (1999) Pharmacokinetic/pharmacodynamic modeling of antitumor agents encapsulated into liposomes. Adv Drug Deliv Rev 40:39–61. 10.1016/S0169-409X(99)00039-310837779
34. Hsu L Investigation of the discriminatory ability of pharmacokinetic metrics for the bioequivalence assessment of PEGylated liposomal doxorubicin Pharm Res 2018 35 106 10.1007/s11095-018-2387-4 29564719
Hsu L (2018) Investigation of the discriminatory ability of pharmacokinetic metrics for the bioequivalence assessment of PEGylated liposomal doxorubicin. Pharm Res 35:106. 10.1007/s11095-018-2387-429564719
35. Völler S Boos J Krischke M Age-dependent pharmacokinetics of doxorubicin in children with cancer Clin Pharmacokinet 2015 54 1139 1149 10.1007/s40262-015-0272-4 25925711
Völler S, Boos J, Krischke M et al (2015) Age-dependent pharmacokinetics of doxorubicin in children with cancer. Clin Pharmacokinet 54:1139–1149. 10.1007/s40262-015-0272-425925711
36. Kontny NE Würthwein G Joachim B Population pharmacokinetics of doxorubicin: establishment of a NONMEM model for adults and children older than 3 years Cancer Chemother Pharmacol 2013 71 749 763 10.1007/s00280-013-2069-1 23314734
Kontny NE, Würthwein G, Joachim B et al (2013) Population pharmacokinetics of doxorubicin: establishment of a NONMEM model for adults and children older than 3 years. Cancer Chemother Pharmacol 71:749–763. 10.1007/s00280-013-2069-123314734
37. García MJ FernándezdeGatta MD Martín A Population pharmacokinetics of doxorubicin and doxorubicinol in patients diagnosed with non-Hodgkin’s lymphoma Br J Clin Pharmacol 2016 82 1517 1527 10.1111/bcp.13070 27447545
García MJ, FernándezdeGatta MD, Martín A et al (2016) Population pharmacokinetics of doxorubicin and doxorubicinol in patients diagnosed with non-Hodgkin’s lymphoma. Br J Clin Pharmacol 82:1517–1527. 10.1111/bcp.1307027447545
38. Bonate PL Strougo A Desai A Guidelines for the quality control of population pharmacokinetic–pharmacodynamic analyses: an industry perspective AAPS J 2012 14 749 758 10.1208/s12248-012-9387-9 22826033
Bonate PL, Strougo A, Desai A et al (2012) Guidelines for the quality control of population pharmacokinetic–pharmacodynamic analyses: an industry perspective. AAPS J 14:749–758. 10.1208/s12248-012-9387-922826033
39. Joerger M Huitema ADR Richel DJ Population pharmacokinetics and pharmacodynamics of doxorubicin and cyclophosphamide in breast cancer patients Clin Pharmacokinet 2007 46 1051 1068 10.2165/00003088-200746120-00005 18027989
Joerger M, Huitema ADR, Richel DJ et al (2007) Population pharmacokinetics and pharmacodynamics of doxorubicin and cyclophosphamide in breast cancer patients. Clin Pharmacokinet 46:1051–1068. 10.2165/00003088-200746120-0000518027989
40. Jonsson EN Karlsson MO Automated covariate model building within NONMEM Pharm Res 1998 15 1463 1468 10.1023/A:1011970125687 9755901
Jonsson EN, Karlsson MO (1998) Automated covariate model building within NONMEM. Pharm Res 15:1463–14689755901
41. Bergstrand M Hooker AC Wallin JE Karlsson MO Prediction-corrected visual predictive checks for diagnosing nonlinear mixed-effects models AAPS J 2011 13 143 151 10.1208/s12248-011-9255-z 21302010
Bergstrand M, Hooker AC, Wallin JE, Karlsson MO (2011) Prediction-corrected visual predictive checks for diagnosing nonlinear mixed-effects models. AAPS J 13:143–151. 10.1208/s12248-011-9255-z21302010
42. Xu G Yang D He C Population pharmacokinetics and toxicity correlation analysis of free and liposome-encapsulated doxorubicin in Chinese patients with advanced breast cancer Cancer Chemother Pharmacol 2023 92 181 192 10.1007/s00280-023-04559-y 37378676
Xu G, Yang D, He C et al (2023) Population pharmacokinetics and toxicity correlation analysis of free and liposome-encapsulated doxorubicin in Chinese patients with advanced breast cancer. Cancer Chemother Pharmacol 92:181–192. 10.1007/s00280-023-04559-y37378676
43. Gabizon AA Tahover E Golan T Pharmacokinetics of mitomycin-c lipidic prodrug entrapped in liposomes and clinical correlations in metastatic colorectal cancer patients Invest New Drugs 2020 38 1411 1420 10.1007/s10637-020-00897-3 31955309
Gabizon AA, Tahover E, Golan T et al (2020) Pharmacokinetics of mitomycin-c lipidic prodrug entrapped in liposomes and clinical correlations in metastatic colorectal cancer patients. Invest New Drugs 38:1411–1420. 10.1007/s10637-020-00897-331955309
44. Colombo I Koster KL Holer L TLD-1, a novel liposomal doxorubicin, in patients with advanced solid tumors: Dose escalation and expansion part of a multicenter open-label phase I trial (SAKK 65/16) Eur J Cancer 2024 201 113588 10.1016/j.ejca.2024.113588 38377773
Colombo I, Koster KL, Holer L et al (2024) TLD-1, a novel liposomal doxorubicin, in patients with advanced solid tumors: dose escalation and expansion part of a multicenter open-label phase I trial (SAKK 65/16). Eur J Cancer 201:113588. 10.1016/j.ejca.2024.11358838377773
45. Kluwe F Michelet R Mueller-Schoell A Perspectives on model-informed precision dosing in the digital health era: challenges, opportunities, and recommendations Clin Pharmacol Ther 2021 109 29 36 10.1002/cpt.2049 33068303
Kluwe F, Michelet R, Mueller-Schoell A et al (2021) Perspectives on model-informed precision dosing in the digital health era: challenges, opportunities, and recommendations. Clin Pharmacol Ther 109:29–36. 10.1002/cpt.204933068303
