Find Important Parameters with Sensitivity Analysis Using SimBiology Model Analyzer App

This example shows how to identify important model parameters for a tumor growth model [1]. In this example, you compute local time-dependent sensitivities of tumor growth with respect to model parameters for an anticancer drug.

Tumor Growth Model

The model used in this example is a SimBiology implementation of the pharmacokinetic/pharmacodynamic (PK/PD) model by Simeoni et al. It quantifies the effect of anticancer drugs on tumor growth kinetics from in vivo animal studies. The drug pharmacokinetics are described by a two-compartment model with IV bolus dosing and linear elimination (ke) from the Central compartment. Tumor growth is described as a biphasic process with an initial exponential growth followed by linear growth. The growth rate of the proliferating tumor cells is described by

L0*x1[1+(L0L1*w)ψ]1ψ

.

Here, L0, L1, and Ψ are tumor growth parameters, x1 is the weight of the proliferating tumor cells, and w is the total tumor weight. In the absence of any drugs, tumor is comprised only of proliferating cells, that is, w = x1. In the presence of an anticancer agent, it is assumed that a fraction of the proliferating cells is transformed into non-proliferating cells. The rate of this transformation is assumed to be a function of the drug concentration in the plasma and an efficacy factor k2. The nonproliferating cells, x2, go through a series of transit stages (x3 and x4) and are eventually cleared from the system. The flow through of the transit compartments is modeled as a first-order process with the rate constant k1.

Some adjustments were made to the pharmacodynamics of tumor growth.

  • Instead of defining the tumor weight as the sum of x1, x2, x3, and x4, the tumor weight is determined by the reaction named Increase, null → tumor_weight, with the reaction rate (2*L1*L0*x12L1+2*L0*x1)*tumor_weight.

    tumor_weight is the total tumor weight, x1 is the weight of the proliferating tumor cells, and L0, and L1 are tumor growth parameters.

  • Similarly, the decrease in tumor weight is defined by the reaction named Decay, tumor_weight → null, with the reaction rate k1*x4. k1 is the forward rate parameter and x4 is the last species in the series of transit reductions in tumor weight.

  • ke is defined as a function of the clearance and the volume of the central compartment: ke = Cl_Central/Central.

Sensitivity Analysis

Sensitivity analysis lets you determine which rate constants and concentrations in a model have significant influence on the overall behavior of the model. This example shows how to investigate which model parameters are sensitive to the tumor growth using the built-in sensitivity analysis program. Specifically, this example shows you how to calculate the sensitivity of tumor_weight species with respect to the model parameters as follows:

(tumor_weight)(L0),(tumor_weight)(L1),(tumor_weight)(k1),...

Calculate Sensitivities

Suppose that you have already calibrated tumor growth model parameters to experimental PK/PD data and saved parameter estimates and initial conditions as a variant called parameterEstimates in your model. As a follow-up analysis, you want to find out which model parameters are sensitive to tumor growth. The following steps show how to do so using the built-in sensitivity analysis program.

Load Tumor Growth Model

  1. Open the SimBiology Model Analyzer app by typing simBiologyModelAnalyzer at the MATLAB® command line or by clicking SimBiology Model Analyzer on the Apps tab.

  2. On the Home tab, click Open. Navigate to the folder matlabroot\examples\simbio, where matlabroot is the folder where MATLAB is installed. Open the SimBiology project file named tumor_growth_vpop_sa.sbproj.

    If you are using a Mac, press Command+Shift+G in the File Browser dialog box, and enter the full path to the folder.

Configure Calculate Sensitivities Program

  1. On the Home tab, select Program > Calculate Sensitivities. A new program opens.

  2. In the Variants section in the Model step of the program, click the option to view the variants to be applied. Then select parameterEstimates. This variant contains previously estimated parameter values and initial conditions for the drug.

  3. In the Doses section, select interval_dose. This repeat dose applies 30 mg of the drug every 4 days starting at day 7 for a total of five times.

  4. In the States To Log section, select [Tumor Growth Model].tumor_weight only.

    Tip

    The context menu of the table contains options to add species, parameters, or compartments, and lets you make selections. Access the context menu by right-clicking the table.

  5. In the Simulation step (an execution step), by default, the normalization method is set to Full (full dedimensionalization), meaning that the app fully normalize the sensitivities so that they can be compared to each other. For details, see Normalization.

    Tip

    Every program execution step can be run separately. An execution step is marked by the run icon at the top of the step. Running an individual step is helpful especially if the program contains multiple steps and you want to see the intermediate results from a particular step. By doing so, you can make adjustments as needed before running the next program step or the whole program. To run the whole program, click the Run button on the Home tab.

  6. In the Sensitivities to Compute section, right-click anywhere in the table and select Add All Constant Parameters. Clear ke because it is defined as a function of the Central volume and Cl_Central, which is already an input. Click the last empty row and type tumor. Then select [Tumor Growth Model].tumor_weight from the list. The tumor_weight species is then added as the output.

  7. On the Home tab, click Run.

Visualize Sensitivity Results

  1. After the sensitivity analysis, the program automatically generates two plots. Click the Plot1 tab to view a plot of the time course of the tumor weight and other sensitivity values d[tumor_weight]/d[parameter]. The time plots for d[tumor_weight]/d[k21] and d[tumor_weight]/d[k12] stay relatively flat compared to other sensitivities, indicating that the tumor growth is not sensitive to k21 and k12 parameters. The Property Editor lets you select which plots to display.

    Tip

    Plots are backed by data that are currently present in the app workspace. Plots are not snapshots. When the data (either experimental data or simulation results) is removed or changed, the plots are also updated according to the changes in the underlying data.

  2. You can also examine the same sensitivity data by checking the magnitude of the computed sensitivities integrated over time. The second generated plot (Plot2) shows a bar plot containing such information. The plot shows that the tumor weight is sensitive to L0, L1, w0, k2, k1, and Cl_Central, but not to k12 and k21.

  3. Save the sensitivity analysis results in a separate folder.

    1. Right-click the LastRun folder of the program in the Browser pane.

    2. Select Save Data.

    3. In the Save Data dialog, enter sa_data as the data name.

Investigate Variability of k12 and k21

The sensitivity analysis results indicate that tumor weight is not sensitive to the model values of k12 and k21. Explore whether varying these parameter values results in different sensitivity results.

Vary k12.  First, vary values for the k12 parameter.

  1. Click the Program1 tab. Click the (+) plus at the top of the program and select Generate Samples.

  2. The Generate Samples step appears. In the Parameter Set section of the step, double-click the empty cell in the Component Name, and enter k12. Set the following options:

    • Type — Range Of Values

    • Spacing — linear

    • Min — 1e-3

    • Max — 1

    • # Of Steps — 25

  3. Disable default plot generation by clicking the plot button at the top of the step.

  4. In the Simulation step, under Sensitivities to Compute, clear all the inputs except k12. Keep tumor_weight as the output. Also, disable the default plot generation by clicking the plot button at the top of the step.

  5. On the Home tab, click the Run button to run the whole program.

  6. Once the simulation finishes, the sensitivity results are stored in the LastRun folder. Expand the LastRun folder. Click results. Then select time from the Plot section on the Home tab.

    The Responses in the plot correspond to the tumor weight and the sensitivity of the tumor weight with respect to k12. Each response is plotted using a different line style. The Scenarios correspond to each parameter scan (simulation scenario). Each scenario is plotted using a different color.

  7. Customize the plot so that the plot uses a different color for each response instead. In the Slice Data table, clear the Style of Scenarios by selecting empty. Set Responses to Color.

  8. In the Responses table, clear tumor_weight and keep only the sensitivity results.

    The plot suggests that the tumor weight is sensitive to some k12 parameter values when model variants different than the estimated parameter values are explored.

  9. Plot the k12 sensitivities together with the previous sensitivity results to see how sensitive the tumor growth is to k12 relative to other model parameters. In the Browser pane, expand the folder sa_data that you saved previously. Then expand results.

  10. Multiselect (Ctrl + Click) the sensitivity data entries for L0, L1, w0, k2, k1, and Cl_Central. Then drag them on to the k12 sensitivity plot Plot3.

    The plot suggests that the tumor growth is still more sensitive to most of the other parameters than k12.

  11. Save the results in a separate folder. Right-click LastRun and select Save Data. Enter k12_data as the data name.

Vary k21.  You can perform a similar analysis by varying k21 and comparing its results with other parameter sensitivities.

  1. Return to Program1. At the bottom of the Generate Samples step, click Add parameter set to scan. A new Parameter Set section appears.

  2. Double-click the empty cell in the Component Name, and enter k21. Set the same options as k12:

    • Type — Range of Values

    • Spacing — linear

    • Min — 1e-3

    • Max — 1

    • # Of Steps — 25

  3. Disable the first parameter set (PS1). Point to the top corner of Parameter Set section of PS1 and click the action button. Then select Disable Parameter Set.

    The Parameter Set is then grayed out to indicate you have successfully disabled it.

  4. In the Simulation step, select k21 as the only input and tumor_weight as the output.

  5. Hit Run from the Home tab to run the program again.

  6. Once the simulation finishes, the sensitivity results are stored in the LastRun folder. Expand the LastRun folder. Click results. Then select time from the Plot section on the Home tab.

  7. In the Slice Data table, clear the Style of Scenarios by selecting empty. Set Responses to Color.

  8. In the Responses table, clear tumor_weight and keep only the sensitivity results.

    The plot suggests that the tumor weight is sensitive to some k21 parameter values when model variants different than the estimated parameter values are explored.

  9. Plot k21 sensitivities together with the previous sensitivity results to compare. Click Plot3. In the Browser pane, expand the LastRun folder. Then expand results.

  10. Drag the k21 sensitivity data d[tumor_weight]/d[k21] on to Plot3.

    The plot suggests that the tumor growth is relatively more sensitive to most other parameters than k21 and k12.

This example shows how to perform sensitivity analysis to find important model parameters to which tumor growth is sensitive. The initial analysis concludes that the tumor weight is not sensitive to some of the parameters. The example then explores the parameter space of these less sensitive parameters and compares the sensitivity results with those of other model parameters. You can use such sensitivity information for other analyses, such as to perform Monte Carlo simulations by varying sensitive parameters to explore the model dynamics and biological variability.

References

[1] Simeoni, M., P. Magni, C. Cammia, G. De Nicolao, V. Croci, E. Pesenti, M. Germani, I. Poggesi, and M. Rocchetti. 2004. Predictive pharmacokinetic-pharmacodynamic modeling of tumor growth kinetics in xenograft models after administration of anticancer agents. Cancer Research. 64:1094-1101.

See Also

|