Noncompartmental Analysis
Noncompartmental analysis (NCA) lets you compute pharmacokinetic (PK) parameters of a drug from the time course of measured drug concentrations. This approach does not require the assumption of a specific compartmental model. NCA is often used to determine the degree of exposure following administration of a drug, such as AUC, and other PK parameters, such as the clearance and the terminal halflife.
Data
SimBiology^{®} lets you calculate NCA parameters from concentration–time data. The data must contain a time column, a concentration column, and a dose column that defines dose amounts. Three types of drug administration routes are supported: IV bolus, IV infusion, and Extravascular. You can have a column for each type. For infusion doses, an infusion rate column is also needed.
If you have data containing multiple groups of observations, you can define a group column. If needed, you can use two levels of hierarchy to specify grouping. Specify the outer level of grouping using the group column, and specify the inner level (subgroups) in the ID column. Consider data that contains three groups, where each group contains four patients. The group column labels the three groups, and the ID column labels each patient.
Dosing
Singledosing data contains a single dose amount for each individual. Multipledosing data has several doses at different times for each individual. There are common parameters calculated for either type of dosing data, and parameters that are specific to single or multiple dosing.
Common Parameters for Single and Multiple Dosing
SimBiology computes some common parameters for single or multipledosing data. This figure represents the concentrationtime profile after a single dose. For multiple dosing, the same principles apply, except that SimBiology uses a steady state dosing period.
Figure A shows concentration–time data in a
linear scale and illustrates how the AUC from time 0 to infinity is calculated.
Figure B shows the same data in a
semilogarithmic scale. To compute the terminal rate constant
(Lambda_z), SimBiology performs a set of linear regressions of the logtransformed data
using each of the last n points (n =
3, 4, 5, ...
) from the terminal portion of the curve.
Lambda_z is chosen from the regression that uses the most
points and has the maximum
adjusted_R^{2}.
This table describes the common parameters for single and multiple dosing.
Parameter  Description 

AUC_0_last  Area under the concentration–time curve from time = 0 to the last time point. $$AUC\_0\_last={\displaystyle \underset{0}{\overset{Tlast}{\int}}C(t)dt}$$, where C(t) is the plasma concentration at time t. SimBiology provides two ways to calculate the AUC.

AUC_infinity  Total area under the concentration–time curve extrapolating to Inf using the terminal rate constant Lambda_z. $$AUC\_infinity=AUC\_0\_last+\frac{C\_last}{Lambda\_z}$$, where C_last is the last observed concentration and Lambda_z is the terminal rate constant. 
AUC_infinity_dose  $$AUC\_infinity\_dose=\frac{AUC\_infinity}{DM}$$. 
AUCx_y  Partial AUC computed for a custom time range, where x and y are the start and end times, respectively. 
AUC_extrap_percent  Fraction of total AUC_infinity obtained from extrapolation. $$AUC\_extrap\_percent=\frac{AUC\_infinityAUC\_0\_last}{AUC\_infinity}*100$$. 
Lambda_z  To calculate the terminal rate constant
(Lambda_z), SimBiology performs a set of linear regressions of the
log(concentration)–time data using each of the last
n points
( Lambda_z is chosen from the regression that uses the most points and has the maximum adjusted_R^{2} among all regressions. $$adjusted\_{R}^{2}=1\frac{(1{R}^{2})*(n1)}{n2}$$ 
R2  Coefficient of determination for the linear regressions (Statistics and Machine Learning Toolbox)
used in the Lambda_z calculation. 
Num_points  Number of data points from the profile used in the determination of Lambda_z. 
C_0  Extrapolated concentration at time = 0, computed using a regression of the first two data points in a profile. This parameter is for IV Bolus doses only. 
C_max  Maximum observed concentration. 
C_max_Dose  $$C\_\mathrm{max}\_Dose=\frac{C\_max}{DM}$$. 
C_max_x_y  Maximum observed concentration within a given time range,
specified by the start time x and the end
time y. This parameter is computed when you
specify a custom time range in the Cmax Time
Range box in the SimBiology Model
Analyzer app or set the
C_max_ranges property of the options
object created by sbioncaoptions . 
T_{last}  Time of the last observed concentration value above the lower limit of quantization (LOQ). 
T_half  Terminal halflife of the drug. $$T\_half=\frac{\mathrm{ln}(2)}{Lambda\_z}$$. 
T_max  T_max is the time point at which the maximum concentration (C_max) is observed. 
T_max_x_y  Time point at which maximum concentration is observed within
a given time range, specified by the start time
x and the end time y.
This parameter is computed when you specify a custom time range
in the Cmax Time Range box in the
SimBiology Model Analyzer or set the
C_max_ranges property of the options
object created by sbioncaoptions . 
V_ss  Apparent volume of distribution at equilibrium. This parameter is for IV Bolus doses only. $$V\_ss=MRT*CL$$. 
V_z  Volume of distribution during the terminal phase. $$V\_z=\frac{DM}{AUC\_infinity*Lambda\_z}$$. 
DM  Administered dose amount. For multiple dosing, the last administered dose is reported. 
doseSchedule  Single or multipledosing data. 
administrationRoute  Dose administration route. Supported routes are
IVBolus , IVInfusion ,
ExtraVascular . 
Parameters for Single Dosing
In addition to the common parameters, SimBiology reports parameters for singledosing data.
Parameter  Description 

AUMC_0_last  Area under the first moment of the concentration–time curve from time 0 to the last time point Tlast. $$AUMC\_0\_last={\displaystyle \underset{0}{\overset{Tlast}{\int}}t*C(t)dt}$$. 
AUMC  Total area under the first moment of the concentration–time curve extrapolating to Inf using Lambda_z. $$AUMC=AUMC\_0\_last+\frac{C\_last}{Lambda\_{z}^{2}}+\frac{Tlast*C\_last}{Lambda\_z}$$. 
AUMC_extrap_percent  Fraction of total AUMC obtained from extrapolation. $$AUMC\_extrap\_percent=\frac{AUMCAUMC\_0\_last}{AUMC}*100$$. 
CL  Total drug clearance. $$CL=\frac{DM}{AUC\_infinity}$$, where DM is the dose amount. 
MRT  Mean residence time. $$MRT=\frac{AUMC}{AUC\_infinity}$$. 
Parameters for Multiple Dosing
This figure shows the concentrationtime profile after multiple doses. SimBiology uses a steady state dosing period to compute the following NCA parameters for multipledosing data, in addition to the common parameters listed previously. In the following figure, the last dosing period is used for illustration purposes.
Parameter  Description 

AUC_Tau  Area under the concentration–time curve during a dosing period of length Tau. SimBiology uses a steadystate dosing period (SS_period). $$AUC\_Tau={\displaystyle {\int}_{{T}_{SS\_period}}^{{T}_{SS\_period}+Tau}C(t)dt}$$. 
Tau  Dosing interval. 
AUMC_Tau  Area under the first moment of the concentration–time curve during a steadystate dosing period of length Tau. $$AUMC\_Tau={\displaystyle {\int}_{{T}_{SS\_period}}^{{T}_{SS\_period}+Tau}t*C(t)dt}$$. 
C_avg  Average concentration over one period. $$C\_avg=\frac{AUC\_Tau}{Tau}$$. 
C_min  Minimum observed concentration during the first period,
that is, 
PTF_percent  Peak trough fluctuation over one dosing interval at steady state. $$PTF\_Percent=\frac{C\_\mathrm{max}C\_\mathrm{min}}{C\_Avg}*100$$. 
Accumulation_Index  Theoretical accumulation ratio. $$Accumulation\_Index=\frac{1}{1{e}^{Lambda\_z*Tau}}$$. 
T_min  Time at which the minimum concentration is reached in a dosing interval. 
MRT  Mean residence time. $$MRT=\frac{AUMC\_Tau+Tau*\left(AUC\_infinityAUC\_0\_last\right)}{AUC\_Tau}$$. Note that for drugs with
prolonged halflives, the extrapolation necessary to compute
the term

CL  Total drug clearance. $$CL=\frac{DM}{AUC\_Tau}$$ Here, DM is the dose amount. 
Sparse Sampling. To calculate PK parameters, measured concentrations at multiple time points for each individual is needed after the drug administration. Under certain circumstances, it is not feasible or not practical to obtain such longitudinal data on a single subject. In such cases, concentration data is collected from multiple individuals at each time point and then averaged to calculate NCA parameters for each group instead. SimBiology performs such sparse sampling by taking the average of the dependent variable for all individuals at the same time point. It then returns the values of NCA parameters for each group. Time values for each measurement across individuals (IDs) within a group must be identical.
Calculate NCA Parameters
You can calculate NCA parameters using the sbionca
function in the command line or using the SimBiology
Model Analyzer app.
Use sbionca
sbionca
provides command line functionality to compute
NCA parameters. Define the data classification options and parameter calculation
options using an option object created by sbioncaoptions
. For an example, see Compute NCA Parameters from ConcentrationTime Data.
Use SimBiology Model Analyzer
After you import the data, select Program > NonCompartmental Analysis on the Home tab. You can classify your data column in the NCA step of the program. If your data has a grouping column, specify it using Group. Use ID to specify the inner level of grouping. Specify the dosing data column (IV Bolus Dose or Extravascular Dose). Lower limit of quantization (LOQ) is a threshold below which the values of dependent variables are truncated to zero.
Lambda Time Range lets you specify a custom time range to compute the terminal rate constant (Lambda_z). The time range applies to all groups; you cannot specify a different time range for each group.
Cmax Time Range lets you specify a custom time range to report the maximum observed concentration within the range (C_max) and the time (T_max) when it is observed. You can specify a different time range for each group.
Partial AUC lets you specify a custom time range to compute the partial AUC bounded by the start and end times. You can specify a different time range for each group.
You can export the NCA results to MATLAB^{®} workspace. By default, the data is exported as a table. To convert
it to a dataset, use table2dataset
(Statistics and Machine Learning Toolbox).
For a workflow example, see Calculate NCA Parameters and Fit Model to PK/PD Data Using SimBiology Model Analyzer.
Note
In the NCA program, you can use simulation results as the input data only
if the underlying simulation used some doses to generate such results. If
there are multiple simulation results (such as nby1
SimData
for simulating n number of
groups), each simulation result (or each group) must use one dose of type
SimBiology.RepeatDose
only.