Simulations
Last updated
Last updated
Simulations can easily be performed in PK-Sim® after having defined at least one building block for Individuals, Compounds and Administration Protocols in the Building Block panel. Optionally, populations, formulations, events and observed data may be defined. Details on these building blocks can be found in the respective sections of this manual:
Even if not all required building blocks are defined you may still set up a simulation. You will then be systematically led through each missing building block where you can specify data or parameters. For details on each building block please refer to the respective chapters of this manual.
Simulations are performed in 3 steps: Creating a simulation (see Creating new simulations in PK-Sim®), running a simulation (see Running a simulation in an individual) and displaying the results in the chart window (see Analyzing results for a simulation in an individual). Experimental data can also be imported in the chart window to enable comparison to the simulated results and to facilitate model refinement.
To create a new simulation:
Click Simulation in the Create Group of the Modeling Tab or
Use the short cut Ctrl+Alt+S
In the Create Simulation dialog window the simulation is set up by specifying Individual or Population simulation, selecting a compound from a list and choosing the model settings:
In order to create a new simulation, you are guided through five tabs of the Create Simulation window:
The Model Structure tab to define model properties
The Compounds tab to select the compounds used for the simulation
The Processes tab to assign the relevant biological processes
The Administration tab to select between different predefined applications
The Events tab to choose a specific event defined in the Event Building Block
In the Model Structure tab of the Create Simulation window you need to:
Enter a name for the simulation
Specify an individual or population by doing one of the following:
Select a previously defined individual or population from the drop-down menu
Allow aging: If a human individual or population is selected the growth of the human individual(s) during the simulation time will be taken into account when5 choosing this option. Based on the human growth and maturation functions available for most parameters in PK-Sim® (e.g. organ volumes, blood flow rates, organ composition, etc.) the parameters are updated along the time scale of the simulation. This is important for multiple drug administration to
e.g. preterm and term neonates, for which the rapid changes in anatomical and physiological properties can influence the pharmacokinetics during the simulated study circle. Please note that the use of the growth function will need additional simulation time so that it is recommended to use this option only if really needed.
Specify compound(s) by doing one of the following:
Select a compound from the list of previously created compounds in the Building Block window.
Select the model settings, i.e. either select the standard model for small molecules or the model for proteins and large molecules
| Model Settings | Description |
|---|---|
Standard model for small molecules | The PK-Sim® standard distribution model assumes 4 subcompartments per organ, i.e. compartments for blood cells, plasma, interstitial space, and cellular space. This model type considers a permeation barrier between plasma and interstitial space and takes into account that the interstitial space has a lower protein content than the plasma. It is especially suitable for small molecules. The plasma-interstitial partition coefficients result from the lower protein concentration in the interstitial space compared to plasma and the unbound fraction in plasma. It is assumed, that the drug has the same affinity to plasma as to interstitial proteins. Thus, effects from the partitioning between plasma and interstitial space may (depending on other physico- chemical data) become important for compounds with a low fraction unbound. The rate of permeation through the endothelial barrier between plasma and interstitial space is determined by the product of endothelial permeability and surface area. The drug dependent permeability can be defined within the Distribution Tab and the individual dependent surface area can be defined for the individual. In the present version of PK-Sim®, the endothelial permeability parameters of the plasma-interstitial barriers are not calculated from physico-chemical compound data. The default value for the plasma-interstitial permeabilities is very large, i.e. the permeabilities have to be adjusted manually, if the permeation across the plasma-interstitial barriers of the organs is expected to be restricted. Using the large default value for plasma-interstitial permeabilities, the exchange rate between plasma and interstitial in 4 subcompartments model are almost instantaneous. |
Model for proteins and large molecules | This model type was developed to simulate the pharmacokinetics of macromolecules such as therapeutic proteins. It is an extension of the standard model for small molecules and contains a description of the transcapillary drug exchange, lymph flows and an additional sub-compartment per organ that represents the endosomal space. The endosomal compartment contains FcRn receptors that are able to protect the macromolecules from catabolism via binding to these receptors (important for e.g. IgG antibodies). |
Please note that currently only for mouse, monkey and human species specific values for the concentration of the FcRn receptor, the concentration of the endogenous IgG and the affinity of the endogenous IgG to the FcRn receptor (needed for the Model for proteins and large molecules) are available. In case of all other animal species just reasonable values are used as physiological default values.
In the Compounds tab of the Create Simulation window you can review the Calculation methods selected for the partition coefficients and the cellular permeabilities. Also, you can review and override values for the Parameter Alternatives.
Click Next in order to assign the relevant biological processes. Note that you can switch between the Tabbed view and the Accordion View in the Utilities/ Options menu if you selected more than one compound.
In the Processes tab of the Create Simulation window, processes like e.g. metabolic, transport, excretion and specific binding processes defined in the selected individual (or population) can be linked to the ADME properties of the selected compounds. A detailed description of how to set up interaction of a drug with enzymes/transporters (e.g. via a drug-drug interaction or induction) can be found in PK-Sim® Compounds: Defining Inhibition/Induction Processes.
If, and only if, processes in individuals are linked to properties of the compounds, these processes are modeled in the simulation.
A green check mark indicates that the process is modeled when running the simulation whereas the yellow exclamation mark indicates that there is no match between the expression data defined in the individual and active processes defined in the compounds, as shown in the following screenshot.
Example 2:
Renal excretion of the compounds should be incorporated into the model. The following settings are required for the selected individual (the capability of renal excretion is a default setting):
Add a renal clearance process for the selected compound in Compound → Biological processes → Transport & Excretion → Renal clearances, as exemplified in the following screenshot:
Link the process in the individual with the biological process of the compound (if not yet automatically done):
If you wish to simulate a renally impaired individual, decrease the GFR properties in the selected individual, not in the compounds!
Click on Next in order to assign an application to your simulation
In the Administration tab of the Create Simulation window, the administration protocol can be loaded from the project database or from a template. Further, a novel administration protocol can be defined and used. In addition, in the case of oral and user defined administration types, a formulation is required for the simulation.
To insert an administration do one of the following:
Select a previously defined protocol from the Administration Protocol drop- down menu
In case of oral and user defined administration, you will have to additionally map the formulation(s) defined in the Formulation building block. In case of sophisticated dosing schedules, different formulations can be chosen for the various administration times. For an intravenous administration (Intravenous Bolus and Intravenous Infusion), the definition of a formulation name is not necessary, because the drug is always assumed to be given in solution. For these administration types, therefore, no match with formulations has to be made. For further information please see PK-Sim® - Administration Protocols.
None or exactly one administration protocol can be assigned to every compound used in the simulation. Each administration protocol can be assigned to maximally one compound.
In the Events tab of the Create Simulation window, events such as the administration of meals and/or discrete gallbladder emptying can be defined at various points of the simulation.
Events are optional for creating a simulation. Possible predefined events are currently the application of a meal to simulate fed conditions for an oral administration and gallbladder emptying independent from meal administration. To insert an event do one of the following:
If you have successfully created a simulation, it will appear in the simulation window with its name and its components:
The green check marks indicate that the building blocks used in the simulation have the same settings as the original building blocks saved in the building block explorer window. Changing the settings in the simulation is done locally without affecting the original building blocks in the building block explorer window.
PK-Sim® offers two alternatives to define drug metabolites. First, metabolites can be a "sink" which means that not actively or passively transported. They possess no physico-chemical or ADME properties and cannot be used as compounds in a simulation. Second, one of the compounds in a simulation can be assigned to be a metabolite or another compound. The metabolite possesses physico-chemical and ADME properties and be transported. In addition, the metabolite can be used in further metabolizations and thus a metabolic network can be built.
If you want the compound not be treated as sink, it needs to be included into the list of compounds in the building blocks.
If you click on Show Diagram, a reaction network of all reactions in the simulation is shown (not just the metabolic network).
In the Observed Data tab observed data sets can be mapped to the corresponding simulation outputs. This mapping is required e.g. for displaying goodness of fit plots and is used as default when adding the simulation to a parameter identification.
All observed data sets belonging to the opened simulation are listed in the mapping table. When adding or removing an observed data set to/from the simulation, the table is updated. Newly added observed data sets are automatically mapped to simulation outputs according to Organ, Compartment and Molecule meta data of the data set and path elements of the outputs. If no matching output can be found, the mapped simulation output is set to None. This means that the specified observed data set is not mapped. The user can also clear an output mapping by selecting the None entry from the Output dropdown. By clicking the "x" on the right side of the grid, the user can delete the observed data from the simulation.
Because meta data of observed data can be incomplete or wrong, you should check whether the right output is mapped to each observed data set. In case of different outputs with the same meta data (this can happen at least in MoBi), you should also check whether the automatically chosen output is correct.
In case of incomplete or missing meta data, it is recommended to correct the meta data first to enable automatic mapping.
For each mapping, the scaling can be defined as Linear or Log which determines the residual calculation.
Scaling
Linear | Residuals are calculated as: |
Log | Residuals are calculated as: |
To reflect the quality or importance of the Observed Data set you can edit the weights of each mapping.
If there is no simulated time value corresponding to an observed time value, linear interpolation between simulated points is used to calculate the missing simulated value. This happens in the calculation of the Residuals in the Simulation, opposed to the calculation of Residuals in Parameter Identification (where every observed time point is added to simulation output points). This might result in slightly different total error values calculated in "Simulation" view and in Parameter Identification.
Per default, the peripheral venous blood plasma is already selected. If desired, select further outputs and then press OK. The result of the individual simulation will be automatically displayed.
In order to select another or additional outputs for the current simulation, click on Define Settings and Run simulation button.
By clicking Save Settings into... at the bottom left corner of the Curve Selection window the settings can be saved as a default for the project or as a user-specific default.
Two venous blood outputs can be selected: “Venous Blood” and “Peripheral Venous Blood”. “Venous Blood” refers to the compartment “Venous Blood” representing the large veins. In clinical practice it is common to sample blood at patients superficial veins, e.g. the antecubital vein. Therefore, PK-Sim® offers the opportunity to also display the pharmacokinetics of the drug in the peripheral venous plasma in order to allow a more accurate description of clinical data. Per default “Peripheral Venous Blood” is a weighted mean of skin and muscle tissue blood (about 70% contribution from skin and about 30% contribution from muscle for all species). You can change the default contribution to “Peripheral Venous Blood” by adjusting the parameters “Fraction of peripheral blood flow in organ” at “Physiology” -> “Flow Rates” -> “Peripheral Blood Flow Fraction” (select “Advanced” view for parameters). The contributions can be defined for arterial blood, bone, fat, muscle, and skin, i.e. all compartments which could possibly contribute to “Peripheral Venous Blood”. The arterial plasma is also considered because of the arteriovenous anastomoses in e.g. the skin of the hand (shunts between arteries and veins involved in the regulation of body temperature). A similar approach to describe peripheral venous plasma concentrations can be found in literature [41].
As the observer for peripheral venous blood sampling represents a balanced mixture of plasma concentrations of the drug in arterial blood, bone, fat, muscle, and skin please do not use the peripheral venous blood observer in the case of e.g. subcutaneous or intramuscular drug administration. This will lead to an overestimation of the concentration in plasma.
For general information about the chart component, i.e. chart settings and options etc., please see Chart Component.
If you would like to cancel the running process, press the Stop button next to the Run button.
If you wish to change the settings of the simulation, click on the Parameters tab. If you have simulated a population, there are two more tabs User Defined Variability and Distribution in which the settings also might be changed. For both, individual and population simulations, there are three views on the parameter settings to select at the bottom of the window:
Simple view: in this view, not all parameters are displayed
Advanced view: all parameters are displayed
Hierarchy view: all parameters are displayed as located in the spatial model structure
You may change the parameter settings for your simulation in order to e.g. achieve a better fit to the observed data. The settings in the building blocks will be unaffected.
It is recommended to select all parameters under consideration as Favorites and to document the source of all parameter values changed from the default in the column Value Description. Then you have a comprehensive overview about the essential input of your simulation, which you can document by copying just the Favorites table.
Update from building block: The simulation settings (local) will be updated with the (global) settings of the building block. This is useful if you want to discard the settings of your simulation and get back to the original settings defined in the building block.\
Commit to building block: The (local) settings of the simulation will be committed to the building block (global settings). This is useful if you want to make these settings available in other simulations.\
Configure: Instead of updating a single building block within your simulation you can also Configure the entire simulation and update or exchange several building blocks at a time. To do so, right mouse click on the simulation and select Configure. The Create Simulation dialog will open where you can exchange the parameters and building blocks of your simulation. In the case of building blocks for which changes were made in the simulation, the name will be supplemented by the warning This is not the template building block!. To update the settings of the simulation select the appropriate building block from the drop-down menu.
As described in the previous chapter Run simulation, clicking on Run in the Modeling & Simulation ribbon or pressing F5 starts the calculation of the simulation. The results will be automatically displayed after finishing the calculation. The calculated results can the be exported to various file formats.
For more information on displaying and editing the charts, see Chart Component
For more information on uploading observed data, see Import and Edit of Observed Data
Generally, two different views are available and switching between these view can be done by clicking on Show PK-Analysis and Show Chart button on the bottom of the chart window:
The simulated curves can be displayed. If simulated curves selected in the curve selection window are missing after (re-)running a simulation, most probably these are not selected in the chart. The graphs can be selected by clicking on Chart Editor on the right hand side of the chart. Curve options, axis options and chart options can be selected here and the general layout of the chart can be configured.
For more information on displaying and editing the chart Display, go to Chart Component
For comparison, observed data in MS Excel® format (*.xls and *.xlsx) can be loaded and displayed in the chart as well. If observed data were added to the Observed Data Building Block, they can be added to the chart by dragging and dropping them onto the chart.
For more information on how to load observed data and compare them to your simulated data, go to Import and Edit of Observed Data
Multiple results windows
Special features for population simulations
If a simulation has been successfully created, click Show PK-Analysis on the bottom of the chart window. PK parameters are then calculated and displayed next to the chart.
The calculated pharmacokinetic parameters are:
In all simulations:
| PK parameter | Description |
|---|---|
AUC_tEnd | |
AUC_inf | |
%AUC(tlast-inf) | |
AUC_tEnd_norm | |
AUC_inf_norm | |
AUC Ratio (AUCR) | AUC_I/AUC (Area under the plasma concentration-time curve of the substrate in the presence of the inhibitor / Area under the plasma concentration-time curve of the substrate in the absence of the inhibitor) |
C_max | Maximum concentration |
C_max_norm | Maximum concentration normalized to dose dose in drug mass per body weight |
C_max Ratio (Cmax_R) | C_max_I/C_max (Maximum concentration of the plasma concentration-time curve of the substrate in the presence of the inhibitor / Maximum concentration of the plasma concentration-time curve of the substrate in the absence of the inhibitor) |
C_tEnd | concentration at the end of simulation |
t_max | Time at which the maximum concentration is assumed |
Half-Life | |
MRT |
In simulations with intravenous administration:
| PK parameter | Description |
|---|---|
VSS(plasma) | |
Vd(plasma) | |
Vss(phys-chem) | Predicted volume of distribution = Volume-weighted mean organ to plasma partition coefficient calculated from physico- chemical compound data. |
Total plasma clearance CL | |
Total body clearance | D/AUC - Total body clearance of drug or apparent clearance (CL/F for extravascular application), dimension is ml/min/kg |
In simulations with oral administration
| PK parameter | Description |
|---|---|
Vss(plasma)/F | Volume of distribution at steady state calculated from plasma curve (see above) divided by bioavailability |
Vd(plasma)/F | |
Total plasma clearance/F | |
Fraction absorbed | Absorbed fraction of applied oral dose. Please note that, e.g. in the case of enterohepatic circulation, this fraction may exceed 1 |
Bioavailability | The bioavailability is only calculated on request. After pressing the button Bioavailability a second simulation with an intravenous (i.v.) short infusion is carried out (internally without being displayed) using identical parameters to the last simulation with oral (p.o.) administration. The bioavailability is then calculated from AUCinf (p.o.)/AUCinf (i.v.) in the venous blood compartment. For a proper estimate of the AUCinf (p.o.) it is recommended to simulate as long as total gastrointestinal transit takes. After changing any parameter the results of previously run simulations are no longer valid. In such a situation, the Bioavailability button is deactivated until the p.o. simulation has been executed with the current parameters. |
In simulations with multiple administrations
| PK parameter | Description |
|---|---|
AUC_inf_tD1 | Area under the concentration vs. time curve from the first data point extrapolated to infinity (further administrations are not considered!) |
AUC_inf_tD1_n | oArmrea under the concentration vs. time curve from the first data point extrapolated to infinity (further administrations are not considered!) normalized to dose, with dose in drug mass per body weight |
....tDi-tDj | Respective PK parameter from the administration time of the first dose until the administration time of the second dose |
....tDlast-tDEnd | Respective PK parameter following the last application |
....tDlast-1- tDlast | Respective PK parameter in the interval between the (last -1) application and the last application |
C_trough_dDi | Trough concentration just before the i-th dose is administered |
C_trough_dlast | Trough concentration just before the last dose is administered |
All values are calculated using the standard equations for PK-values (see e.g. M. Rowland, T. N. Tozer, “Clinical Pharmacokinetic Concepts and Applications”, (1994) Lippincott Williams & Wilkins, Philadelphia). For extrapolation to infinity an exponential function is used on the basis of the last 10% of the calculated time steps. AUC is calculated by extrapolating the first time steps to t = 0. Depending on the curve shape the result of this extrapolation may be sensitive to time resolution. This may lead to some variability in AUC and thus, in clearance and distribution volumes.
PK-parameters for selected outputs of population simulations are shown in two ways:
Individual PK Values are calculated for all individuals within the range being analyzed and the median is presented
Aggregated PK Values are calculated from aggregated curve being analyzed
Globally calculated PK-parameters are always calculated for all individuals with the median being presented
By clicking on Export to Excel® the calculated PK-parameters (including the simulated concentration-time profiles) can be exported to MS Excel® format.
Population simulation analysis function offers a variety of analyses and graphical displays: Time Profile, Box Whisker, Scatter and Range Analysis that will be explained in detail in the following sections of the manual.
Unlike for simulations in single individuals, additional variability for parameters in population simulations may be defined. In the User Defined Variability tab, parameters from all building blocks that need to be varied may be selected and specific distributions/variabilities can be set. In the Distribution tab, these settings are graphically visualized.
For more information on adding variability to a defined population, go to PK-Sim® - Creating Populations.
Choose the organ and the compartment by expanding the respective drop-down menue and select by ticking the box in the right column of the table. Press OK to start the simulation run.
By clicking into the Save Settings at the bottom left corner of the Curve Selection window the settings can be saved as user-specific default.
After the simulation has been calculated, the Create Time Profile Analysis window opens up next in which the graphical display of simulated time profiles can be specified. If you do not want to plot time profiles, close this window and choose the analysis type you require from the Analysis button.
For plotting time profiles, use the Output tab to specify in which of the previously selected organs/compartment you want to analyze and plot time profiles. Highlight the respective organ/compartment and add it to the right panel. In the panel below, the Outputs: Select distribution statistics for display, you may specify which statistical quantity of your output will be graphically displayed.
You may change the name, the unit and the scaling of selected output in all plot types. Changing the name may be useful when graphically displaying several parameters and the parameter name is printed along the y-axis.
The calculated PK-parameters of all individual curves can be exported using the Export PK-Analyses to CSV... item from the context menu of the simulation.
The underlying data of population analysis can be exported to Excel® using the Export to Excel®**...** item from the context menu of the chart view. Two sheets are created containing the data as original and transposed matrix.
Output may be grouped using population parameters and PK-parameters. The Population Parameters tab offers grouping by population parameters e.g. Characteristics of individual (weight, age, etc.), anatomy or physiology. A comprehensive list itemizes each PK-parameter in the organs/compartments specified in the Outputs tab. In order to select a parameter as a grouping criterion, expand the tree underlying each group of criteria and use the Add or Removebutton to add or remove a specific grouping PK-Parameter to the right panel. Highlight the respective grouping and click on the Create Grouping button to specify the grouping intervals (binning).
For PK-parameters with numerical values like C_max or AUC grouping can be done into equally sized intervals using the User defined limits. Alternatively, grouping by numerical PK-parameters can be done into equally populated intervals of varying size using the User defined (equally populated) number of bins. For PK-parameters that assume string values like gender or race grouping is defined by Value mapping. Groupings have to be named and the name of the grouping appears later on as an item in the list of grouping parameters.
User defined limits
Grouping by customized intervals is available in both the Population Parameter and the PK-Parameter tabs. If you decide on defining the limits of each bin, the minimum and maximum values of your grouping criterion are automatically entered in the table listing the grouping intervals.
Continue as follows:
Define the maximum value of the first bin.
Add additional bins by clicking on the + button in the right column of the table that lists the binning intervals.
In your last bin, the upper limit of you last bin is the maximum value of the parameter.
For each interval, select a label and change color and symbol (used in scatter plots) by changing default settings.
User defined (equally populated) number of bins
If you decide on using equally populated bins for grouping of your output, labels for the individual bins have to be designed for the display.
In the Label Generation panel, do the following:
Decide on the number of bins used for grouping.
Choose from the drop-down menu a symbol that will be used in a scatter plot.
Choose a start and end color for plotting the output in the first and the last interval. The color gradient in between is set automatically. Colors in individual intervals can be changed in the Generated Label panel below.
Individual intervals are named using the selected Naming pattern that can be customized. Per default it consists of the selected Template e.g. alphabetic or numerical labelling together with the lower and upper limit of the interval. However, any of these elements may be omitted. More information on customizing the label generation are available directly within the software.
When you finish defining the grouping intervals, the grouping appears in the list of selected parameters. If you highlight a specific grouping, the lower panel of the Population Parameters and the PK-Parameters shows the distribution of the individuals of the population when stratified by this grouping along with the label of each bin.
Value mapping for string parameters
Population parameters such as gender, race and population name can be grouped by their string values. This grouping can be used to customize labels e.g. F or M instead of Female or Male and/or the relative order of these labels in the analysis (i.e. which label should come first) menu serves the design of customized labels.
In the Time Profile tab, concentration-time profiles or fraction-time profiles can be plotted. First, in the Selected Outputs window, select from list of previously specified organs/compartments the concentration-/fraction-time profile(s) that you want to plot. Please note that concentration and fraction cannot be used at the same time for a given analysis. Second, if you want to display by grouping, select from the Available Fields panel a grouping criterion specified in the Population Parameters or the PK-Parameters.
There are several ways to visualize stratified results:
If you drag a grouping criterion into the Panes panel, time-profiles in each subgroup will be displayed in a separate panel on the right. More than one grouping criterion may be selected and the respective number of panels will be generated.
If you drag a grouping criterion into the Colors panel, time-profiles in each subgroup will be displayed in one color. Colors were previously selected in the Outputs window. Only one grouping criterion may be selected here, to be able to differentiate the subgroups.
If you want to display results grouped by two criteria, use the Panes and the Colors panels simultaneously. In this case, it may wise to select only one output because curves from different compartments can only be differentiated by their progression.
For each curve PK parameters are calculated and can be shown by clicking on the Show PK-Analysis button.
If you hover over the simulated time-profiles, details of the curve are displayed, e.g. the organ/compartment in which the time-profile was simulated, the statistical quantity that is plotted and the x- and y-values of the underlying data points.
Details of your grouping can be changed:
Go back to the Population Parameters or the PK-Parameters tab.
Highlight the name of the grouping in the panel listing the selected parameters.
Right mouse click on the name of the grouping and select Edit from the menu.
Use the magnifying glass to zoom into an area that you wish to enlarge.
Use the magnifying glass to zoom into an area that you wish to enlarge.
In order to display Observed Data in your Time Profile Analysis do and note the following:
Drag and drop the observed data from the Building Block panel into the Analysis panel.
Either display specific files or folders of observed data in the building block.
Use the Chart Options to specify the settings for display of the observed data.
When observed data is displayed for the first time, it appears in all analysis panels by default.
In the Chart Options, observed data can be grouped by meta data and is displayed in the respective panel. To apply the grouping also to the observed data the meta data field of the observed data must be named like the field in the analysis and the value must match.
Once observed data is dragged into the analysis panel, it is associated to that simulation and can be removed by removing the observed data from the simulation tree.
For comparison charts observed data can not be removed from the analysis, but only made invisible by unticking the Visible box in the Chart Options.
In contrast to the Time Profile Analysis, population parameters and PK- parameters serve not only as grouping criteria in the Box Whisker Analysis, but may also be plotted. Any parameter selected in the Population Parameters tab and PK-Parameters tab is listed in the Selected Outputs panel of the Box Whisker Plot tab. Select a parameter for plotting by ticking the respective box. Grouping criteria are defined as described above and can be selected from the Available Fields. For grouping along the x-axis, several criteria may be selected by dragging to the X-Grouping panel. Grouping by color may be used for one criterion only and a corresponding legend is automatically created. The box comprises 25% (lower rim) to 75% (upper rim) of the values. The whiskers extend from 2.5% to 97.5% of the values in the population. The distance between the lower and upper box rim is called the inner quartile range (IQR). Outliers are values which lie outside the range from lower whisker limit - 1.5 time IQR to upper whisker limit + 1.5 time IQR.
A separate panel is created for each of the selected output parameters in the graphical display. The name and unit of the output parameters is printed along the Y-axis.
For each observed value a point is plotted with observed value as x-Value and corresponding simulated value as y-Value.
Adding Deviation Lines to the plot
In a Predicted vs. Observed plot the user can right click on the chart and add deviation lines:
This opens a dialog where the user can specify the fold value of the deviation curves.
This will create two deviation lines according to the given x-fold value which has to be greater than 1. An x-fold deviation range includes simulated values within x-fold and 1/x-fold of observed values.
In the Chart Editor the deviation lines are grouped under the Category Identity.
This chart is similar to the Time Profile chart, but on the y-axis the (absolute) residuals are plotted. The chart includes scaling and weights.
At the top of the chart, the total residual error is displayed.
In the Scatter Plot Analysis population parameters and PK-parameters can also be used as grouping criteria and be plotted. Grouping criteria are defined as described above. In addition to grouping in different panes or by different colors, you can group your output by symbols. Do so by dragging a grouping criterion into the Symbols panel. The symbols need to be defined along with the other grouping criteria in the Population Parameter or PK-Parameter tab.
Similar to the Scatter Plot Analysis, the Range Plot Analysis offers continuous plotting of parameters on both axis. Grouping by Population Parameters and/or PK-Parameters is done as described above. The range plot displays the median of the parameter as a thick line and the parameter range that comprises 5-95% of individuals of the population.
Details of your analysis can be changed by right mouse click into the plot and selecting Edit.
In certain cases, it can be helpful to clone a simulation, e.g. in order to keep all parameters defined in the simulation, such as the partition coefficient method, and exchange only one particular building block, e.g. the administration protocol. This can easily be done by cloning an existing simulation and subsequently exchanging the respective building block using the Configure functionality (see above).
To clone a simulation in the project:
Right mouse click on the respective simulation in the Simulations Explorer
The Cloning Simulation dialog, which is identical to the Create Simulation dialog for new simulations, will open. You will have to enter an alternative name for the simulation clone.
You can now go through the Cloning Simulation dialog and exchange parameters and/or entire building blocks of your simulation clone. To update the settings of the simulation select the appropriate building block from the drop-down menus.
Please note that a simulation can only be cloned if the local simulation settings are in agreement with the global settings of the building blocks, as indicated by green check marks.
The comparison chart function allows for comparison of results of different individual or population simulations in one plot.
To create a comparison plot for simulation in individuals proceed in the following way:
Open the comparison chart (if not already open). An empty chart will be displayed.
Drag and drop the simulations you want to compare from the Simulation window into the empty comparison chart.
Select the simulated curves you would like to display in the comparison chart.
If desired, add observed data from the building block by drag and drop to the chart.
To rename a comparison chart, right mouse click on the respective summary chart in the Simulation window, select Rename and enter the new name of the chart.
PK-Sim® offers a comparison chart also for population simulations. To create a comparison plot for simulations in populations proceed in the following way:
When you have approved your selection of the curves by clicking OK, the Comparison chart chart offers you to look at the distribution of population parameters in the set union of the populations.
You can then select the analysis type for your comparison from the Analysis button and then proceed with selection output and grouping criteria as described above.
For a comparison of populations, only the intersecting set of output appears, e.g. if you selected venous plasma concentrations as output in all populations, it appears in the list of Outputs in the time profile analysis tab.
If you define a population as a reference population, it will not be grouped into subpopulations. This might be useful when you want to compare variability in subgroups of one population with the variability of a reference population. This is illustrated below.
Define a new individual or population by clicking Create and follow the instructions of the dialog windows
Load a new individual or population from template by clicking Load and follow the instructions of the dialog windows
Create a new compound by clicking Create and follow the instructions in the dialog windows.
Load a new compound from template by clicking Load and follow the instructions in the dialog windows.
Example 1:
Define a new administration by clicking Create and follow the instructions of the dialog windows
Load an administration from the template database by clicking Load and follow the instructions of the dialog windows
Select a previously defined event from the Event drop-down menu
Define a new event by clicking Create and follow the instructions of the appearing dialog windows
Load an event from the template database by clicking and follow the instructions of the appearing dialog windows
Click OK to finish the Create Simulation process and to close the window.
If a simulation was successfully created, press the Run simulation button in the Modeling & Simulation ribbon or press the F5 key. If the simulation is run for the first time, ths following window will appear in which the simulation curves that will be generated in the simulation can be selected:
User Defined parameter node shows an overview of all parameters changed by the user in the simulation.
If you change the parameter settings in the Parameter tab, the green check marks (traffic lights) on the Results tab will turn red indicating that the displayed simulation results were not performed with the current settings: Press the Run simulation button in the Modeling & Simulation ribbon or press the F5 key again to re-run the simulation with the current settings and display the results.
If parameter settings were changed in the simulation, the red traffic lights in the Simulation window indicate that the local settings in the simulation are different from the settings in the Building Block, i.e. the global settings: A right click on the red traffic lights in the simulation window allows for two possibilities: \
Chart view (default): The simulated curves and the observed data are displayed
PK-Analysis view : The calculated PK parameters for the selected simulated curves are displayed.
If you wish to display different curves or use different settings for the same simulation, you can add new results charts windows by clicking on the Results button in the Ribbon group Simulation.
If you wish to display other percentiles or other curves without simulating again, click on the Results button in the ribbon group Simulation. You will be directed to the Curve selection for chart window (see above).
Area under curve from tstart to tend of the simulation.
Area under curve extrapolated to infinity (using the terminal 10% of data points).
, where λ is the terminal elimination rate (calculated from the terminal 10% of simulated points).
Percentage of AUCinf after end of simulation time.
Area under curve extrapolated to infinity and normalized to the dose dose in drug mass per body weight.
Area under curve extrapolated to infinity and normalized to the dose [dose in drug mass per body weight].
Terminal half life time (calculated from the terminal 10% of data points).
, where λ is the terminal elimination rate (calculated from the terminal 10% of simulated points).
Mean residence time is calculated by:
The infusion time is set to 0 for non intravenous administrations. The AUMCs are calculated according to:
, where λ is the terminal elimination rate (calculated from the terminal 10% of simulated points).
Volume of distribution at steady state calculated from the plasma curve according to:
Apparent volume of distribution is calculated from the plasma curve according to: , where λ is the terminal elimination rate (calculated from the terminal 10% of simulated points). If simulation time is short, the distribution volume may be incorrect. Run a longer simulation (to approach steady state) and reevaluate the distribution volume.
Total clearance calculated from plasma curve according to:
Apparent volume of distribution calculated from the plasma curve according to (see above) divided by bioavailability
Total clearance calculated from plasma curve according to divided by bioavailability
Once all parameters are set, click the Run simulation button in the Modeling & Simulation group or press the F5 key. The Population Simulation Settings window: Curve selection window will appear in which organs and compartments are selected in which time profiles and PK-parameters will be calculated.
Chart view (default): The calculated curves and the observed data are displayed
PK-Analysis view : The calculated PK parameters for the selected calculated curves are displayed. See “PK- Analysis view”.
Confirm and close the window by clicking OK
Select Clone...
Once all parameters and building blocks are appropriate confirm and close the window by clicking OK
Click on the Comparison Chart button in the Ribbon group Simulation and select For Individual Simulations. A new entry in the Simulation window will appear, Comparison Chart 1.
Click on the Comparison Chart button in the Ribbon group Simulation and select For Population Simulations. A new window Select simulations to use in comparison opens up. Select the simulations you want to compare by ticking the respective boxes. In addition, you can select a population as a reference population and specify labeling, symbol and color for plotting.
