Yes / Adequate drug distribution through tumours is essential for treatment to be effective. Palbociclib is a cyclin-dependent kinase (CDK) 4/6 inhibitor approved for use in patients with hormone receptor (HR) positive, HER2 negative metastatic breast cancer (BC). It has unusual physicochemical properties, which may significantly influence its distribution in tumour tissue. We studied the penetration and distribution of palbociclib in vitro, including the use of multicellular three-dimensional models and mathematical modelling. MCF-7 and DLD-1 cell lines were grown as single cell suspensions (SCS) and spheroids; palbociclib uptake and efflux were studied using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Intracellular concentrations of palbociclib for MCF-7 SCS (Cmax 3.22 µM) and spheroids (Cmax 2.91 µM) were 32 and 29 fold higher and in DLD-1, 13 and 7 fold higher, respectively than the media concentration (0.1 µM). Total palbociclib uptake was lower in DLD-1 cells than MCF-7 cells both in SCS and in spheroids. Both uptake and efflux of palbociclib were slower in spheroids than SCS. These data were used to develop a mathematical model of palbociclib transport that quantifies key parameters determining drug penetration and distribution. The model reproduced qualitatively most features of the experimental data and distinguished between SCS and spheroids, providing additional support for hypotheses derived from the experimental data. Mathematical modelling has the potential for translating in vitro data into clinically relevant estimates of tumour drug concentrations. / Grant for Translational Research and a grant from Leeds NHS Charitable Trustees.

Objetivo: O objetivo do presente estudo foi descrever a relação entre as concentrações plasmáticas totais e livres tumorais do etoposídeo (ETO) e a inibição do crescimento do tumor observada em ratos Wistar portadores de tumor Walker- 256 (W256) utilizando a modelagem farmacocinética/farmacodinâmica (PK/PD). Métodos: Os procedimentos com animais foram aprovados no CEUA/UFRGS sob o número 22302. Os experimentos de farmacocinética foram realizados para determinar concentrações plasmáticas e livres em duas regiões do tumor sólido W256 através de microdiálise. Após a administração do ETO nas doses de 10 ou 20 mg/kg i.v. bolus em ratos Wistar portadores de tumor W256, amostras de sangue e microdialisado de tecido do centro e periferia do tumor foram coletadas simultaneamente, até 7 h pós-dose, para determinar o fator de penetração no tumor. Um método analítico por CLAE-UV foi desenvolvido e validado para quantificação do etoposídeo nas amostras de plasma e dialisado. Os experimentos de farmacodinâmica foram conduzidos em ratos portadores de tumor W256 que receberam ETO 5 e 10 mg/kg i.v. bolus uma vez ao dia por 8 e 4 dias, respectivamente. O volume dos tumores foram monitorados diariamente durante 30 dias. Análise não-compartimental dos dados de PK foi realizada no WinNonlin®. A modelagem dos dados PK e PK/PD foi realizada no Monolix®, utilizando abordagem populacional. Os dados PK/PD foram analisados usando o modelo Simeoni TGI modificado através da introdução de uma função Emax para descrever a relação nãolinear entre a concentração plasmática e tumoral e o efeito. Resultados e Discussão: O método por CLAE-UV foi desenvolvido e validado para quantificar as amostras de ETO em plasma e tecido. A penetração do ETO no tumor foi maior na periferia (61 ± 15 % e 61 ± 29 %) do que no centro do tumor (34 ± 6 % e 28 ± 11 %) após administração das doses 10 e 20 mg/kg, respectivamente (ANOVA, α = 0.05). Um modelo de 4 compartimentos compreendendo uma distribuição saturável (cinética de Michaelis-Menten) nos compartimentos tumorais a partir do compartimento central modelou simultaneamente os perfis de concentração-tempo do ETO em plasma e em ambas regiões do tumor. O modelo populacional PK/PD Simeoni TGI–Emax foi capaz de descrever o efeito antitumoral dependente do regime de administração do ETO utilizando concentrações totais plasmáticas ou livres no tumor, resultando em um maior k2max (potência máxima) para as concentrações livres (25,8 mL.μg-1.dia-1 - intratumoral vs. 12,6 mL.μg-1.dia-1 - plasma total). Conclusões: Os resultados mostram que a utilização das concentrações livres do fármaco no tumor para a modelagem PK/PD pode fornecer um melhor entendimento da relação farmacocinética e farmacodinâmica e melhoram a capacidade de previsão do modelo, considerando que a eficácia dos fármacos antineoplásicos no tratamento de tumores sólidos é dependente da capacidade do fármaco em se distribuir no tecido tumoral. / Objective: The aim of this study was to describe the relationship between total plasma and free interstitial tumor etoposide (ETO) concentrations and the drug tumor growth inhibition observed in a Walker-256 (W256) tumor-bearing Wistar rat model using the pharmacokinetic/pharmacodynamic (PK/PD) modeling. Methods: The experiments with animals were approved by CEUA/UFRGS (protocol number 22302). Pharmacokinetic experiments were conducted to determine total plasma and free intratumoral concentrations in two regions of W256 solid tumor by microdialysis. After administration of ETO 10 or 20 mg/kg i.v. bolus to W256 tumorbearing Wistar rats, blood and tissue microdialysate samples from tumor center and periphery were simultaneously collected up to 7h to determine the tumor penetration factor. An analytical HPLC-UV method was developed and validated for quantification of ETO in plasma and microdialysate samples. The pharmacodynamic experiments were conducted in W256 tumor-bearing rats that received ETO 5 or 10 mg/kg i.v. bolus every day for 8 and 4 days, respectively. Tumor volumes were monitored daily for 30 days. Non-compartmental analysis of PK data was performed in WinNonlin®. The PK and PK/PD modeling by population approach were performed using Monolix®. PK/PD data were analyzed using a modification of Simeoni TGI model by introducing an Emax function to describe the nonlinear relationship between tumor and plasma concentrations and effect. Results and Discussion: The HLPCUV method was developed and validated to determine plasma and tissue samples of ETO. ETO tumor penetration was higher in the tumor periphery (61 ± 15 % and 61 ± 29 %) than center (34 ± 6 % and 28 ± 11 %) following 10 and 20 mg/kg doses, respectively (ANOVA, α = 0.05). A 4-compartment structural model comprising a saturable distribution (Michaelis-Menten kinetics) into the tumor compartments from the central compartment simultaneously described the ETO concentration–time profiles in plasma and both tumor regions. The PK/PD population Simeoni TGI–Emax model was capable of describing the schedule-dependent antitumor effects of ETO using total plasma or free tumor concentrations obtained in a W256-tumor bearing Wistar rat model, resulting in higher k2max (maximal potency) for free concentrations (25.8 mL.μg-1.day-1 - intratumoral vs. 12.6 mL.μg-1.day-1 total plasma). Conclusions: The results showed that the use of free intratumoral drug concentrations in the PK/PD modeling can provide a better understanding of the pharmacokinetics and pharmacodynamics relationship and improve the forecasting ability of the models considering that the efficacy of antineoplastic drugs in the treatment of solid tumors is dependent on the drug ability to distribute into the tumor.

Objetivo: O objetivo do presente estudo foi descrever a relação entre as concentrações plasmáticas totais e livres tumorais do etoposídeo (ETO) e a inibição do crescimento do tumor observada em ratos Wistar portadores de tumor Walker- 256 (W256) utilizando a modelagem farmacocinética/farmacodinâmica (PK/PD). Métodos: Os procedimentos com animais foram aprovados no CEUA/UFRGS sob o número 22302. Os experimentos de farmacocinética foram realizados para determinar concentrações plasmáticas e livres em duas regiões do tumor sólido W256 através de microdiálise. Após a administração do ETO nas doses de 10 ou 20 mg/kg i.v. bolus em ratos Wistar portadores de tumor W256, amostras de sangue e microdialisado de tecido do centro e periferia do tumor foram coletadas simultaneamente, até 7 h pós-dose, para determinar o fator de penetração no tumor. Um método analítico por CLAE-UV foi desenvolvido e validado para quantificação do etoposídeo nas amostras de plasma e dialisado. Os experimentos de farmacodinâmica foram conduzidos em ratos portadores de tumor W256 que receberam ETO 5 e 10 mg/kg i.v. bolus uma vez ao dia por 8 e 4 dias, respectivamente. O volume dos tumores foram monitorados diariamente durante 30 dias. Análise não-compartimental dos dados de PK foi realizada no WinNonlin®. A modelagem dos dados PK e PK/PD foi realizada no Monolix®, utilizando abordagem populacional. Os dados PK/PD foram analisados usando o modelo Simeoni TGI modificado através da introdução de uma função Emax para descrever a relação nãolinear entre a concentração plasmática e tumoral e o efeito. Resultados e Discussão: O método por CLAE-UV foi desenvolvido e validado para quantificar as amostras de ETO em plasma e tecido. A penetração do ETO no tumor foi maior na periferia (61 ± 15 % e 61 ± 29 %) do que no centro do tumor (34 ± 6 % e 28 ± 11 %) após administração das doses 10 e 20 mg/kg, respectivamente (ANOVA, α = 0.05). Um modelo de 4 compartimentos compreendendo uma distribuição saturável (cinética de Michaelis-Menten) nos compartimentos tumorais a partir do compartimento central modelou simultaneamente os perfis de concentração-tempo do ETO em plasma e em ambas regiões do tumor. O modelo populacional PK/PD Simeoni TGI–Emax foi capaz de descrever o efeito antitumoral dependente do regime de administração do ETO utilizando concentrações totais plasmáticas ou livres no tumor, resultando em um maior k2max (potência máxima) para as concentrações livres (25,8 mL.μg-1.dia-1 - intratumoral vs. 12,6 mL.μg-1.dia-1 - plasma total). Conclusões: Os resultados mostram que a utilização das concentrações livres do fármaco no tumor para a modelagem PK/PD pode fornecer um melhor entendimento da relação farmacocinética e farmacodinâmica e melhoram a capacidade de previsão do modelo, considerando que a eficácia dos fármacos antineoplásicos no tratamento de tumores sólidos é dependente da capacidade do fármaco em se distribuir no tecido tumoral. / Objective: The aim of this study was to describe the relationship between total plasma and free interstitial tumor etoposide (ETO) concentrations and the drug tumor growth inhibition observed in a Walker-256 (W256) tumor-bearing Wistar rat model using the pharmacokinetic/pharmacodynamic (PK/PD) modeling. Methods: The experiments with animals were approved by CEUA/UFRGS (protocol number 22302). Pharmacokinetic experiments were conducted to determine total plasma and free intratumoral concentrations in two regions of W256 solid tumor by microdialysis. After administration of ETO 10 or 20 mg/kg i.v. bolus to W256 tumorbearing Wistar rats, blood and tissue microdialysate samples from tumor center and periphery were simultaneously collected up to 7h to determine the tumor penetration factor. An analytical HPLC-UV method was developed and validated for quantification of ETO in plasma and microdialysate samples. The pharmacodynamic experiments were conducted in W256 tumor-bearing rats that received ETO 5 or 10 mg/kg i.v. bolus every day for 8 and 4 days, respectively. Tumor volumes were monitored daily for 30 days. Non-compartmental analysis of PK data was performed in WinNonlin®. The PK and PK/PD modeling by population approach were performed using Monolix®. PK/PD data were analyzed using a modification of Simeoni TGI model by introducing an Emax function to describe the nonlinear relationship between tumor and plasma concentrations and effect. Results and Discussion: The HLPCUV method was developed and validated to determine plasma and tissue samples of ETO. ETO tumor penetration was higher in the tumor periphery (61 ± 15 % and 61 ± 29 %) than center (34 ± 6 % and 28 ± 11 %) following 10 and 20 mg/kg doses, respectively (ANOVA, α = 0.05). A 4-compartment structural model comprising a saturable distribution (Michaelis-Menten kinetics) into the tumor compartments from the central compartment simultaneously described the ETO concentration–time profiles in plasma and both tumor regions. The PK/PD population Simeoni TGI–Emax model was capable of describing the schedule-dependent antitumor effects of ETO using total plasma or free tumor concentrations obtained in a W256-tumor bearing Wistar rat model, resulting in higher k2max (maximal potency) for free concentrations (25.8 mL.μg-1.day-1 - intratumoral vs. 12.6 mL.μg-1.day-1 total plasma). Conclusions: The results showed that the use of free intratumoral drug concentrations in the PK/PD modeling can provide a better understanding of the pharmacokinetics and pharmacodynamics relationship and improve the forecasting ability of the models considering that the efficacy of antineoplastic drugs in the treatment of solid tumors is dependent on the drug ability to distribute into the tumor.

Objetivo: O objetivo do presente estudo foi descrever a relação entre as concentrações plasmáticas totais e livres tumorais do etoposídeo (ETO) e a inibição do crescimento do tumor observada em ratos Wistar portadores de tumor Walker- 256 (W256) utilizando a modelagem farmacocinética/farmacodinâmica (PK/PD). Métodos: Os procedimentos com animais foram aprovados no CEUA/UFRGS sob o número 22302. Os experimentos de farmacocinética foram realizados para determinar concentrações plasmáticas e livres em duas regiões do tumor sólido W256 através de microdiálise. Após a administração do ETO nas doses de 10 ou 20 mg/kg i.v. bolus em ratos Wistar portadores de tumor W256, amostras de sangue e microdialisado de tecido do centro e periferia do tumor foram coletadas simultaneamente, até 7 h pós-dose, para determinar o fator de penetração no tumor. Um método analítico por CLAE-UV foi desenvolvido e validado para quantificação do etoposídeo nas amostras de plasma e dialisado. Os experimentos de farmacodinâmica foram conduzidos em ratos portadores de tumor W256 que receberam ETO 5 e 10 mg/kg i.v. bolus uma vez ao dia por 8 e 4 dias, respectivamente. O volume dos tumores foram monitorados diariamente durante 30 dias. Análise não-compartimental dos dados de PK foi realizada no WinNonlin®. A modelagem dos dados PK e PK/PD foi realizada no Monolix®, utilizando abordagem populacional. Os dados PK/PD foram analisados usando o modelo Simeoni TGI modificado através da introdução de uma função Emax para descrever a relação nãolinear entre a concentração plasmática e tumoral e o efeito. Resultados e Discussão: O método por CLAE-UV foi desenvolvido e validado para quantificar as amostras de ETO em plasma e tecido. A penetração do ETO no tumor foi maior na periferia (61 ± 15 % e 61 ± 29 %) do que no centro do tumor (34 ± 6 % e 28 ± 11 %) após administração das doses 10 e 20 mg/kg, respectivamente (ANOVA, α = 0.05). Um modelo de 4 compartimentos compreendendo uma distribuição saturável (cinética de Michaelis-Menten) nos compartimentos tumorais a partir do compartimento central modelou simultaneamente os perfis de concentração-tempo do ETO em plasma e em ambas regiões do tumor. O modelo populacional PK/PD Simeoni TGI–Emax foi capaz de descrever o efeito antitumoral dependente do regime de administração do ETO utilizando concentrações totais plasmáticas ou livres no tumor, resultando em um maior k2max (potência máxima) para as concentrações livres (25,8 mL.μg-1.dia-1 - intratumoral vs. 12,6 mL.μg-1.dia-1 - plasma total). Conclusões: Os resultados mostram que a utilização das concentrações livres do fármaco no tumor para a modelagem PK/PD pode fornecer um melhor entendimento da relação farmacocinética e farmacodinâmica e melhoram a capacidade de previsão do modelo, considerando que a eficácia dos fármacos antineoplásicos no tratamento de tumores sólidos é dependente da capacidade do fármaco em se distribuir no tecido tumoral. / Objective: The aim of this study was to describe the relationship between total plasma and free interstitial tumor etoposide (ETO) concentrations and the drug tumor growth inhibition observed in a Walker-256 (W256) tumor-bearing Wistar rat model using the pharmacokinetic/pharmacodynamic (PK/PD) modeling. Methods: The experiments with animals were approved by CEUA/UFRGS (protocol number 22302). Pharmacokinetic experiments were conducted to determine total plasma and free intratumoral concentrations in two regions of W256 solid tumor by microdialysis. After administration of ETO 10 or 20 mg/kg i.v. bolus to W256 tumorbearing Wistar rats, blood and tissue microdialysate samples from tumor center and periphery were simultaneously collected up to 7h to determine the tumor penetration factor. An analytical HPLC-UV method was developed and validated for quantification of ETO in plasma and microdialysate samples. The pharmacodynamic experiments were conducted in W256 tumor-bearing rats that received ETO 5 or 10 mg/kg i.v. bolus every day for 8 and 4 days, respectively. Tumor volumes were monitored daily for 30 days. Non-compartmental analysis of PK data was performed in WinNonlin®. The PK and PK/PD modeling by population approach were performed using Monolix®. PK/PD data were analyzed using a modification of Simeoni TGI model by introducing an Emax function to describe the nonlinear relationship between tumor and plasma concentrations and effect. Results and Discussion: The HLPCUV method was developed and validated to determine plasma and tissue samples of ETO. ETO tumor penetration was higher in the tumor periphery (61 ± 15 % and 61 ± 29 %) than center (34 ± 6 % and 28 ± 11 %) following 10 and 20 mg/kg doses, respectively (ANOVA, α = 0.05). A 4-compartment structural model comprising a saturable distribution (Michaelis-Menten kinetics) into the tumor compartments from the central compartment simultaneously described the ETO concentration–time profiles in plasma and both tumor regions. The PK/PD population Simeoni TGI–Emax model was capable of describing the schedule-dependent antitumor effects of ETO using total plasma or free tumor concentrations obtained in a W256-tumor bearing Wistar rat model, resulting in higher k2max (maximal potency) for free concentrations (25.8 mL.μg-1.day-1 - intratumoral vs. 12.6 mL.μg-1.day-1 total plasma). Conclusions: The results showed that the use of free intratumoral drug concentrations in the PK/PD modeling can provide a better understanding of the pharmacokinetics and pharmacodynamics relationship and improve the forecasting ability of the models considering that the efficacy of antineoplastic drugs in the treatment of solid tumors is dependent on the drug ability to distribute into the tumor.

Afin d'analyser les données in vitro de l'effet de l'erlotinib sur la croissance des cellules A431 suivie par impédance-métrie, nous avons développé un modèle pharmacocinétique - pharmacodynamique non linéaire décrivant simultanément la diminution de la concentration d'erlotinib et son effet sur la croissance cellulaire au cours du temps. La non-linéarité du modèle imposant le recours à des méthodes itératives pour l'estimation des paramètres, plusieurs étapes de la procédure d'identification du modèle ont été étudiées et des solutions proposées, avec des exemples d'application à des molécules utilisées en oncologie:Choix du critère d'optimisation à employer - supériorité de la relation fonctionnelle de la moyenne géométrique pour l'identification de modèles non linéaires. Application données réelles: courbe de calibration d'une expérience de dosage ELISA du bevacizumab. Choix de l'algorithme d'optimisation le plus approprié au problème d'identification du processus pharmacocinétique. Les algorithmes dérivatifs sont les plus performants. Application données réelles: estimation simultanée des paramètres du modèle pharmacocinétique du 5-fluorouracile et de son métabolite principal.Transformation de la forme différentielle initiale du modèle en temps continu vers un modèle récursif en temps discret. Par ce moyen le modèle devient linéaire en ses paramètres, ce qui permet d'estimer directement les paramètres sans utiliser d'algorithme d'optimisation. Il devient également possible de suivre les variations des paramètres au cours du temps. Application données réelles: pharmacocinétique de la fotemustine, de la mitoxantrone et du 5-fluorouracile. / A nonlinear pharmacokinetic-pharmacodynamic model has been devised do simultaneously describe the loss of erlotinib and its effect on the cell growth over time, in order to analyze impedance-based data of erlotinib effect on A431 cells growth in vitro over time. The model non-linearity requiring the use of iterative methods for parameter estimation, several steps of the model identification were studied, and solutions proposed, with application examples to cancer drugs :Choice of the optimization criterion - superiotity of the geometric mean functionnal relationship for non-linear model identification. Real data application : calibration curve of a bevacizumab ELISA quantification experiment.Choice of the most appropriate algorithm for the pharmacokinetic process identification problem. The derivative algorithms perform better. Real data application : simultaneous identification of the 5-fluorouracil and of its main metabolite pharmacokinetic system.Transform of the differential initial continuous-time model in a recursive discrete time model. The transformed model becomes linear with respect to its parameters, allowing straightforward parameter estimation without using any optimization algorithm. It is then also possible to track the parameter variations over time. Real data application : pharmacokinetic model parameter estimation of fotemustine, mitoxantrone and 5-fluorouracil.