# DDMODEL00000274: Terranova_2017_oncology_TGI

Short description:
PKPD model of tumor growth inhibition and toxicity outcome after administration of anticancer agents in xenograft mice
 Format: PharmML 0.8.x (0.8.1) Related Publication: .hiddenContent {display:none;} Evaluation of a PK/PD DEB-based model for tumor-in-host growth kinetics under anticancer treatment E.M.Tosca, E.Borella, N.Terranova, M.Rocchetti, P.Magni PAGE 25, 6/2015 Affiliation: Department of Electrical, Computer and Biomedical Engineering, University of Pavia, via Ferrata 5, Pavia, I-27100, Italy Abstract: Objectives: Mathematical models for describing the tumor growth in animals often neglect the relationship between tumor and host organism. To overcome this limitation, a more mechanistic model, based on energy balance between tumor and host, was developed. This PK/PD model, combining the Dynamic Energy Budget (DEB) theory with the Simeoni tumor growth inhibition (TGI) model, describes both the dynamics of the tumor-host interaction and the effect of anticancer treatments. Here a slightly revised model formulation and a new implementation are proposed. Moreover, a comparative study on the tumor growth in control groups between the DEB-TGI model and the widely used Simeoni TGI model is presented. Methods: Data used for model validation refer to xenograft experiments conducted on Harlan Sprague Dawley mice. Average data of tumor weight and mice net body weight were considered for the control and treated groups. The PKs were derived from separated studies. Monolix 4.3.3 was used for model identification, while Simulx was used to confirm the hypothesis emerged from a dynamic system analysis. Results: First of all, the model was identified on different experimental datasets with the following strategy: 1) physiological parameters of the tumor-free model were estimated on growth data of typical HSD mice; 2) estimated values were used to find the initial value for the energy reserve at the beginning of the experiment; 3) once fixed the tumor-free model parameters and energy initial value, the tumor-related and the drug-related parameters were simultaneously estimated. The mathematical analysis of the dynamic system showed that, as the Simeoni model, the DEB-TGI model predicts an exponential growth of the tumor in the early phases of its development. The exponential growth rate depends on several model parameters some of them related to the tumor cell lines and other to the host. We investigated also the relationship between the DEB-TGI model parameters and the decreasing of the tumor growth rate. Conclusions: The tumor-in-host DEB-based model confirmed its good capability in describing tumor growth and host body growth even when an anticancer drug is administered. Moreover, the affinities emerged from the comparative analysis with the Simenoni model provide a possible biological interpretation of the assumptions underlying the Simeoni model unperturbed (control) growth curve. Contributors: Elena Maria Tosca
 Context of model development: Candidate Comparison, Selection, Human Dose Prediction; Discrepancy between implemented model and original publication: Among the drugs considered in the paper, only PACLITAXEL has been used; Model compliance with original publication: Yes; Model implementation requiring submitter’s additional knowledge: No; Modelling context description: Host features, such as cell proliferation rates, caloric intake, metabolism and energetic conditions, significantly influence tumor growth; at the same time, tumor growth may have a dramatic impact on the host conditions. For example, in clinics, at certain stages of the tumor growth, cachexia (body weight reduction) may become so relevant to be considered as responsible for around 20% of cancer deaths. Unfortunately, anticancer therapies may also contribute to the development of cachexia due to reduced food intake (anorexia), commonly observed during the treatment periods. For this reason, cachexia is considered one of the major toxicity findings to be evaluated also in preclinical studies. However, although various pharmacokinetic-pharmacodynamic (PK-PD) tumor growth inhibition (TGI) models are currently available, the mathematical modelling of cachexia onset and TGI after an anticancer administration in preclinical experiments is still an open issue. To cope with this, a new PK-PD model, based on a set of tumor-host interaction rules taken from Dynamic Energy Budget (DEB) theory and a set of drug tumor inhibition equations taken from the well-known Simeoni TGI model, was developed. The model is able to describe the body weight reduction, splitting the cachexia directly induced by tumor and that caused by the drug treatment under study. It was tested in typical preclinical studies, essentially designed for efficacy evaluation and routinely performed as a part of the industrial drug development plans. For the first time, both the dynamics of tumor and host growth could be predicted in xenograft mice untreated or treated with different anticancer agents and following different schedules. ; Modelling task in scope: simulation; estimation; Nature of research: Preclinical development; In vivo; Therapeutic/disease area: Oncology;
 Validation Status: Annotations are correct. Certification Comment: This model is not certified.
• Model owner: Elena Maria Tosca
• Submitted: Dec 22, 2017 8:10:56 AM
##### Revisions
• Version: 5
• Submitted on: Dec 22, 2017 8:10:56 AM
• Submitted by: Elena Maria Tosca
• With comment: Model revised without commit message

### Name

Terranova_2017_oncology_TGI

### Description

PKPD model of tumor growth inhibition and toxicity outcome after administration of anticancer agents in xenograft mice

 T

### Function Definitions

 $\mathrm{additiveError}:\mathrm{real}\left(\mathrm{additive}:\mathrm{real}\right)=\mathrm{additive}$

### Parameter Model: $\mathrm{pm}$

#### Random Variables

${\mathrm{eps_RES_W}}_{\mathrm{vm_err.DV}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=1\right)$
${\mathrm{eps_RES_Wu}}_{\mathrm{vm_err.DV}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=1\right)$

#### Population Parameters

$\mathrm{K10_POP}$
$\mathrm{K12_POP}$
$\mathrm{K21_POP}$
$\mathrm{V1_POP}$
$\mathrm{En_initial_POP}$
$\mathrm{rho_b_POP}$
$\mathrm{xi_POP}$
$\mathrm{ni_POP}$
$\mathrm{gr_POP}$
$\mathrm{V1inf_POP}$
$\mathrm{mu_POP}$
$\mathrm{mu_u_POP}$
$\mathrm{gu_POP}$
$\mathrm{delta_Vmax_POP}$
$\mathrm{W_initial_POP}$
$\mathrm{Vu1_initial_POP}$
$\mathrm{IC50_POP}$
$\mathrm{k1_POP}$
$\mathrm{k2_POP}$
$\mathrm{b_W}$
$\mathrm{b_Wu}$

#### Individual Parameters

$\mathrm{K10}=\mathrm{pm.K10_POP}$
$\mathrm{K12}=\mathrm{pm.K12_POP}$
$\mathrm{K21}=\mathrm{pm.K21_POP}$
$\mathrm{V1}=\mathrm{pm.V1_POP}$
$\mathrm{En_initial}=\mathrm{pm.En_initial_POP}$
$\mathrm{rho_b}=\mathrm{pm.rho_b_POP}$
$\mathrm{xi}=\mathrm{pm.xi_POP}$
$\mathrm{ni}=\mathrm{pm.ni_POP}$
$\mathrm{gr}=\mathrm{pm.gr_POP}$
$\mathrm{V1inf}=\mathrm{pm.V1inf_POP}$
$\mathrm{mu}=\mathrm{pm.mu_POP}$
$\mathrm{mu_u}=\mathrm{pm.mu_u_POP}$
$\mathrm{gu}=\mathrm{pm.gu_POP}$
$\mathrm{delta_Vmax}=\mathrm{pm.delta_Vmax_POP}$
$\mathrm{W_initial}=\mathrm{pm.W_initial_POP}$
$\mathrm{Vu1_initial}=\mathrm{pm.Vu1_initial_POP}$
$\mathrm{IC50}=\mathrm{pm.IC50_POP}$
$\mathrm{k1}=\mathrm{pm.k1_POP}$
$\mathrm{k2}=\mathrm{pm.k2_POP}$
$\mathrm{density_V}=1$
$\mathrm{density_Vu}=1$
$\mathrm{omeg}=0.75$
$m=\frac{\mathrm{pm.ni}}{{\mathrm{pm.V1inf}}^{\frac{1}{3}}\cdot \mathrm{pm.gr}}$
$\mathrm{Z_initial}=\frac{\mathrm{pm.W_initial}}{\left(1+\mathrm{pm.En_initial}\cdot \mathrm{pm.xi}\right)}$

### Structural Model: $\mathrm{sm}$

#### Variables

$C=\frac{\mathrm{sm.Q1}}{\mathrm{pm.V1}}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Q1}=\mathrm{pm.K21}\cdot \mathrm{sm.Q2}-\left(\mathrm{pm.K10}+\mathrm{pm.K12}\right)\cdot \mathrm{sm.Q1}\\ \mathrm{Q1}\left(T=0\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Q2}=\mathrm{pm.K12}\cdot \mathrm{sm.Q1}-\mathrm{pm.K21}\cdot \mathrm{sm.Q2}\\ \mathrm{Q2}\left(T=0\right)=0\end{array}$
$\mathrm{rho}=\mathrm{pm.rho_b}\cdot \left(1-\frac{\mathrm{sm.C}}{\left(\mathrm{pm.IC50}+\mathrm{sm.C}\right)}\right)$
$\mathrm{ku}=\frac{\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}}{\left(\mathrm{sm.Z}+\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}\right)}$
$\mathrm{switch1}=\frac{\left(\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{pm.ni}\cdot \mathrm{sm.En}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}-\mathrm{pm.gr}\cdot \mathrm{pm.m}\cdot \mathrm{sm.Z}\right)}{\left(\mathrm{pm.gr}+\left(1-\frac{\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}}{\left(\mathrm{sm.Z}+\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}\right)}\right)\cdot \mathrm{sm.En}\right)}$
$\mathrm{switch2}=\frac{\left(\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{pm.ni}\cdot \mathrm{sm.En}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}-\mathrm{pm.gr}\cdot \mathrm{pm.m}\cdot \mathrm{sm.Z}\right)}{\left(1-\frac{\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}}{\left(\mathrm{sm.Z}+\mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}\right)}\right)\cdot \left(\mathrm{sm.En}+\mathrm{pm.omeg}\cdot \mathrm{pm.gr}\right)}$
$\mathrm{Wu}=\mathrm{pm.density_Vu}\cdot \left(\mathrm{sm.Vu1}+\mathrm{sm.Vu2}+\mathrm{sm.Vu3}+\mathrm{sm.Vu4}\right)$
$W=\mathrm{pm.density_V}\cdot \left(1+\mathrm{pm.xi}\cdot \mathrm{sm.En}\right)\cdot \mathrm{sm.Z}$
$\mathrm{W_err}=\mathrm{pm.b_W}\cdot \sqrt{\mathrm{sm.W}}$
$\mathrm{Wu_err}=\mathrm{pm.b_Wu}\cdot \sqrt{\mathrm{sm.Wu}}$
$\begin{array}{c}\frac{d}{dT}Z=\left\{\begin{array}{lll}\frac{\left(\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{pm.ni}\cdot \mathrm{sm.En}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}-\mathrm{pm.gr}\cdot \mathrm{pm.m}\cdot \mathrm{sm.Z}\right)}{\left(\mathrm{pm.gr}+\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{sm.En}\right)}& \text{if}& \mathrm{sm.switch1}\ge 0\\ \frac{\left(\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{pm.ni}\cdot \mathrm{sm.En}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}-\mathrm{pm.gr}\cdot \mathrm{pm.m}\cdot \mathrm{sm.Z}\right)}{\left(1-\mathrm{sm.ku}\right)\cdot \left(\mathrm{sm.En}+\mathrm{pm.omeg}\cdot \mathrm{pm.gr}\right)}& \text{if}& \mathrm{sm.switch1}<0\wedge \mathrm{sm.switch2}\le -\mathrm{pm.delta_Vmax}\\ -\mathrm{pm.delta_Vmax}& \text{otherwise}& \end{array}\\ Z\left(T=0\right)=\mathrm{pm.Z_initial}\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{En}=\frac{\mathrm{pm.ni}}{{\mathrm{sm.Z}}^{\frac{1}{3}}}\cdot \left(\mathrm{sm.rho}\cdot {\frac{\mathrm{pm.V1inf}}{\left(\mathrm{sm.Vu1}+\mathrm{sm.Z}\right)}}^{\frac{2}{3}}-\mathrm{sm.En}\right)\\ \mathrm{En}\left(T=0\right)=\mathrm{pm.En_initial}\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Vu1}=\left\{\begin{array}{lll}\frac{\left(\mathrm{pm.ni}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}+\mathrm{pm.m}\cdot \mathrm{sm.Z}\right)\cdot \mathrm{pm.gr}\cdot \mathrm{sm.ku}\cdot \mathrm{sm.En}}{\left(\mathrm{pm.gr}\cdot \mathrm{pm.gu}+\left(1-\mathrm{sm.ku}\right)\cdot \mathrm{pm.gu}\cdot \mathrm{sm.En}\right)}-\mathrm{pm.mu}\cdot \mathrm{sm.Vu1}-\mathrm{pm.k2}\cdot \mathrm{sm.Vu1}\cdot \mathrm{sm.C}& \text{if}& \mathrm{sm.switch1}\ge 0\\ \frac{\mathrm{pm.gr}\cdot \mathrm{pm.m}\cdot \mathrm{pm.mu_u}\cdot \mathrm{sm.Vu1}}{\mathrm{pm.gu}}-\mathrm{pm.mu}\cdot \mathrm{sm.Vu1}-\mathrm{pm.k2}\cdot \mathrm{sm.C}\cdot \mathrm{sm.Vu1}& \text{if}& \mathrm{sm.switch1}<0\wedge \mathrm{sm.switch2}<0\wedge \mathrm{sm.switch2}\ge -\mathrm{pm.delta_Vmax}\\ \frac{\mathrm{sm.ku}}{\mathrm{pm.gu}}\cdot \left(\mathrm{sm.En}\cdot \mathrm{pm.ni}\cdot {\mathrm{sm.Z}}^{\frac{2}{3}}+\mathrm{pm.delta_Vmax}\cdot \mathrm{sm.En}+\mathrm{pm.delta_Vmax}\cdot \mathrm{pm.omeg}\cdot \mathrm{pm.gr}\right)-\mathrm{pm.mu}\cdot \mathrm{sm.Vu1}-\mathrm{pm.k2}\cdot \mathrm{sm.C}\cdot \mathrm{sm.Vu1}& \text{otherwise}& \end{array}\\ \mathrm{Vu1}\left(T=0\right)=\mathrm{pm.Vu1_initial}\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Vu2}=\mathrm{pm.k2}\cdot \mathrm{sm.C}\cdot \mathrm{sm.Vu1}-\mathrm{pm.k1}\cdot \mathrm{sm.Vu2}\\ \mathrm{Vu2}\left(T=0\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Vu3}=\mathrm{pm.k1}\cdot \mathrm{sm.Vu2}-\mathrm{pm.k1}\cdot \mathrm{sm.Vu3}\\ \mathrm{Vu3}\left(T=0\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Vu4}=\mathrm{pm.k1}\cdot \mathrm{sm.Vu3}-\mathrm{pm.k1}\cdot \mathrm{sm.Vu4}\\ \mathrm{Vu4}\left(T=0\right)=0\end{array}$

### Observation Model: $\mathrm{om1}$

#### Continuous Observation

$\mathrm{Y1}=\mathrm{sm.W}+\mathrm{additiveError}\left(\mathrm{additive}=\mathrm{sm.W_err}\right)+\mathrm{pm.eps_RES_W}$

### Observation Model: $\mathrm{om2}$

#### Continuous Observation

$\mathrm{Y2}=\mathrm{sm.Wu}+\mathrm{additiveError}\left(\mathrm{additive}=\mathrm{sm.Wu_err}\right)+\mathrm{pm.eps_RES_Wu}$

## External Dataset

 OID $\mathrm{nm_ds}$ Tool Format NONMEM

### File Specification

 Format $\mathrm{csv}$ Delimiter comma File Location Simulated_DEB_TGI_data.csv

### Column Definitions

Column ID Position Column Type Value Type
$\mathrm{ID}$
$1$
$\mathrm{id}$
$\mathrm{int}$
$\mathrm{TIME}$
$2$
$\mathrm{idv}$
$\mathrm{real}$
$\mathrm{DV}$
$3$
$\mathrm{dv}$
$\mathrm{real}$
$\mathrm{DVID}$
$4$
$\mathrm{dvid}$
$\mathrm{int}$
$\mathrm{AMT}$
$5$
$\mathrm{dose}$
$\mathrm{real}$
$\mathrm{EVID}$
$6$
$\mathrm{evid}$
$\mathrm{real}$
$\mathrm{CMT}$
$7$
$\mathrm{cmt}$
$\mathrm{int}$

### Column Mappings

Column Ref Modelling Mapping
$TIME$
$T$
$DV$
$\left\{\begin{array}{lll}\mathrm{om1.Y1}& \text{if}& \mathrm{DVID}=1\\ \mathrm{om2.Y2}& \text{if}& \mathrm{DVID}=2\end{array}$
$AMT$
$\left\{\begin{array}{lll}\mathrm{sm.Q1}& \text{if}& \mathrm{AMT}>0\end{array}$

## Estimation Step

 OID $\mathrm{estimStep_1}$ Dataset Reference $\mathrm{nm_ds}$

### Parameters To Estimate

Parameter Initial Value Fixed? Limits
pm.K10_POP
$20.832$
true
$\left(,\right)$
pm.K12_POP
$0.144$
true
$\left(,\right)$
pm.K21_POP
$2.011$
true
$\left(,\right)$
pm.V1_POP
$813.1$
true
$\left(,\right)$
pm.En_initial_POP
$1.3$
true
$\left(,\right)$
pm.xi_POP
$0.184$
true
$\left(,\right)$
pm.ni_POP
$1.2242$
true
$\left(,\right)$
pm.gr_POP
$12.2$
true
$\left(,\right)$
pm.V1inf_POP
$22.6$
true
$\left(,\right)$
pm.rho_b_POP
$1$
true
$\left(,\right)$
pm.mu_POP
$0.0223$
false
$\left(0,\right)$
pm.mu_u_POP
$13.3$
false
$\left(0,\right)$
pm.gu_POP
$11.7$
false
$\left(0,\right)$
pm.delta_Vmax_POP
$0.185$
false
$\left(0,\right)$
pm.W_initial_POP
$21.2$
false
$\left(0,\right)$
pm.Vu1_initial_POP
$0.0023$
false
$\left(0,\right)$
pm.IC50_POP
$0.461$
false
$\left(0,\right)$
pm.k1_POP
$0.462$
false
$\left(0,\right)$
pm.k2_POP
$6.53E-4$
false
$\left(0,\right)$
pm.b_W
$0.101$
false
$\left(0,\right)$
pm.b_Wu
$0.134$
false
$\left(0,\right)$

### Operations

#### Operation: $1$

 Op Type generic
##### Operation Properties
Name Value
algo
$\text{foce}$

## Step Dependencies

Step OID Preceding Steps
$\mathrm{estimStep_1}$