Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Matching Model Detail to the Prediction Target

Compartmental Versus Whole-Body Modeling · Tutorial 37 of 799

Matching Model Detail to the Prediction Target

Choose a lumped or organ-resolved PBPK structure by asking which predictions must be reliable, not simply which model has more compartments.

Beginner 10 min read

What You'll Learn

  • Define a prediction target before selecting model structure.
  • Distinguish plasma exposure from average tissue burden and organ-specific concentration.
  • Calculate a volume-weighted average from organ concentrations.
  • Identify when a pooled tissue value cannot answer an organ-specific question.
  • Compare model predictions using target-specific error and tolerances.
  • Diagnose common mismatches between model detail and intended use.

Choose Detail for the Prediction You Need

A PBPK model can represent tissues as individual organs or combine several tissues into a smaller number of pooled compartments. Neither choice is automatically better. The relevant question is whether the chosen structure can provide the prediction needed for a particular decision. A model intended to predict plasma exposure may not need to report every organ concentration. A model intended to predict concentration in a named organ must, however, preserve enough information to distinguish that organ from the tissues pooled with it.

The target should be stated before comparing structures: specify the quantity, its location or aggregation, the time period, and any acceptable prediction error. “Predict tissue distribution” is too vague. “Predict the volume-weighted mean concentration across liver, muscle, and adipose tissue at 24 hours” is a defined target, as is “predict kidney concentration at 24 hours.”

Definition: A prediction target is the precisely specified model output to be used or evaluated, including its quantity, anatomical scope, and time or time interval. A lumped model represents multiple tissues with one pooled state and generally provides a combined prediction for them. An organ-resolved model represents the relevant organs as separate states, allowing separate organ predictions. Model detail is adequate when the structure can represent the target at the required accuracy under the intended conditions.

Three common targets illustrate why the choice matters. Plasma concentration or plasma area under the concentration–time curve (AUC) is a systemic target. Average tissue burden combines information across a defined set of tissues. Concentration in a particular organ is a local target. These outputs are related, but one cannot generally be substituted for another: a correct plasma profile does not by itself establish correct tissue concentrations.

Define the Aggregation Before Pooling

Suppose a target averages concentrations across \(m\) tissues. Let \(V_i\) be the volume of tissue \(i\), and let \(C_i(t)\) be its concentration at time \(t\). If the target is the mean concentration over the combined tissue volume, its value is

$$ \bar{C}_{\mathrm{tissue}}(t) = \frac{\sum_{i=1}^{m} V_i C_i(t)} {\sum_{i=1}^{m} V_i}. $$

The weights are tissue volumes because this definition averages concentration across the combined volume. The total amount in those tissues is a different quantity, \(\sum_i V_i C_i(t)\), with units of amount rather than concentration. An unweighted mean, \(\sum_i C_i(t)/m\), answers yet another question: it gives every tissue the same weight regardless of size. The target definition, not convenience, determines which summary is appropriate.

The same weighting applies to AUC when the set of tissues and their volumes are fixed. Integrating the volume-weighted concentration over time gives

$$ \mathrm{AUC}_{\mathrm{tissue}} = \int_0^T \bar{C}_{\mathrm{tissue}}(t)\,dt = \frac{\sum_{i=1}^{m} V_i\mathrm{AUC}_i} {\sum_{i=1}^{m} V_i}, \qquad \mathrm{AUC}_i=\int_0^T C_i(t)\,dt. $$

This follows because the volumes are constant: the finite weighted sum can be taken outside the integral. It does not imply that a single pooled compartment will reproduce each tissue’s trajectory. It establishes how to calculate the specified aggregate from organ-resolved predictions.

Worked Example: Calculating Average Tissue Burden

Hypothetical data. At a specified time, an organ-resolved model predicts concentrations of 4.0 mg/L in liver, 1.0 mg/L in muscle, and 2.0 mg/L in adipose tissue. Their volumes are 1.5 L, 20 L, and 15 L, respectively. The target is the mean concentration over the combined volume of these three tissues.

First calculate each tissue’s amount: liver, \(1.5\ \mathrm{L}\times4.0\ \mathrm{mg/L}=6.0\ \mathrm{mg}\); muscle, \(20\ \mathrm{L}\times1.0\ \mathrm{mg/L}=20\ \mathrm{mg}\); and adipose tissue, \(15\ \mathrm{L}\times2.0\ \mathrm{mg/L}=30\ \mathrm{mg}\). The total amount is \(56\ \mathrm{mg}\), and the combined volume is \(1.5+20+15=36.5\ \mathrm{L}\).

Therefore,

$$ \bar{C}_{\mathrm{tissue}} = \frac{(1.5)(4.0)+(20)(1.0)+(15)(2.0)} {1.5+20+15} = \frac{56}{36.5} = 1.534\ \mathrm{mg/L}. $$

The unweighted mean would be \((4.0+1.0+2.0)/3=2.333\ \mathrm{mg/L}\), but that is not the defined target. It gives the small liver the same weight as the much larger muscle and adipose tissue. If a lumped model predicts \(1.60\ \mathrm{mg/L}\), its difference from this target is \(1.60-1.534=0.066\ \mathrm{mg/L}\), or approximately \(4.30\%\) of the organ-resolved value. Whether that difference is acceptable depends on the prespecified use and tolerance.

Match Structure to Three Common Targets

For a plasma target, a lumped tissue representation may be sufficient if it reproduces plasma predictions accurately across the relevant dose, time range, and conditions. Plasma data alone do not establish that tissue predictions are accurate: different tissue distributions can be compatible with similar plasma profiles. If the application only uses plasma exposure, that limitation may be acceptable. If tissue concentrations affect the decision, plasma agreement is not enough.

For an average tissue-burden target, pooling can be appropriate when the required output is specifically the defined average, and the pooled representation preserves that average over the relevant conditions. It may be unnecessary to estimate separate organ concentrations when no decision depends on them. But a pooled state is not automatically equal to a volume-weighted mean: the model’s pooled quantity and the target’s weights must actually correspond.

For a named-organ target, the model must distinguish that organ or provide an independently justified mapping to its concentration. A pooled prediction can sometimes approximate one organ under restricted conditions, but it does not generally identify that organ’s value. When organs differ in flow, tissue partitioning, or kinetic behavior, a shared state can conceal differences that matter to the target.

Worked Example: Selecting Structure for a Plasma Target

Hypothetical data. The target is plasma concentration at four specified times after the same dose and under the same conditions. A prespecified benchmark provides values of 10.0, 6.0, 3.0, and 1.0 mg/L. A lumped model predicts 10.1, 5.9, 3.2, and 0.8 mg/L. An organ-resolved model predicts 10.05, 6.05, 3.1, and 0.9 mg/L. Assume the predictions are evaluated on the same benchmark, and the practical tolerance is an RMSE no greater than 0.20 mg/L.

For \(n\) paired predictions, define root mean squared error as the square root of the mean squared prediction error:

$$ \mathrm{RMSE} = \sqrt{\frac{1}{n}\sum_{j=1}^{n} \left(\widehat{C}_{j}-C_{j,\mathrm{benchmark}}\right)^2}. $$

For the lumped model, the errors are \(0.1,-0.1,0.2,-0.2\ \mathrm{mg/L}\). Thus,

$$ \mathrm{RMSE}_{\mathrm{lumped}} = \sqrt{\frac{0.1^2+(-0.1)^2+0.2^2+(-0.2)^2}{4}} = \sqrt{0.025} = 0.1581\ \mathrm{mg/L}. $$

For the organ-resolved model, the errors are \(0.05,0.05,0.1,-0.1\ \mathrm{mg/L}\), giving

$$ \mathrm{RMSE}_{\mathrm{organ}} = \sqrt{\frac{0.05^2+0.05^2+0.1^2+(-0.1)^2}{4}} = \sqrt{0.00625} = 0.07906\ \mathrm{mg/L}. $$

Both models meet the stated plasma-target tolerance; the organ-resolved model has the lower RMSE, but that improvement is not by itself a reason to require its extra detail. If the intended use is limited to these plasma predictions and the simpler model is otherwise adequate, the lumped model may be a defensible choice. This comparison says nothing about which model better predicts an individual organ concentration.

Worked Example: When a Named-Organ Target Requires Resolution

Hypothetical data. The target is kidney concentration at a specified time. Two equal-volume tissues, kidney and another tissue, have the same volume-weighted mean concentration of \(5.0\ \mathrm{mg/L}\) in each of two possible distributions: in distribution A, their concentrations are 8.0 and 2.0 mg/L; in distribution B, they are 2.0 and 8.0 mg/L. An independent hypothetical benchmark for kidney concentration is \(7.5\ \mathrm{mg/L}\).

For distribution A, the pooled concentration is \(((1)(8.0)+(1)(2.0))/(1+1)=5.0\ \mathrm{mg/L}\). For distribution B, it is \(((1)(2.0)+(1)(8.0))/(1+1)=5.0\ \mathrm{mg/L}\) as well. Thus the pooled value cannot distinguish the two kidney concentrations, even though they are \(8.0\) and \(2.0\ \mathrm{mg/L}\).

If a lumped model reports \(5.0\ \mathrm{mg/L}\) as a proxy for kidney concentration, its absolute error against the benchmark is \(|5.0-7.5|=2.5\ \mathrm{mg/L}\). An organ-resolved model predicting \(7.8\ \mathrm{mg/L}\) has absolute error \(|7.8-7.5|=0.3\ \mathrm{mg/L}\). For this target and these hypothetical results, the organ-resolved prediction is much closer. The central point is structural: the pooled average alone does not contain enough information to recover a particular organ’s concentration.

Justify the Choice and Check for Failure

A useful selection procedure is to make the target explicit, test whether the simpler structure preserves it, and add resolution only when that structure cannot meet the target’s requirements. This is consistent with the earlier material on model structure as a prediction choice, pooling tissues, and structural error: the comparison must focus on the output needed, rather than on compartment count alone.

1
Specify the output.
Name the quantity, anatomical scope, time or interval, conditions, and practical tolerance. Define any average with its weights.
2
Check what the structure can identify.
Determine whether the model has a separate state or justified output mapping for the target. A pooled state ordinarily identifies a pooled quantity, not each component value.
3
Evaluate predictions for that target.
Compare structures against relevant observations or a suitable benchmark using target-specific errors over the intended time and conditions. Do not infer organ accuracy from plasma agreement alone.
4
Choose the least detailed adequate structure.
Prefer a simpler model if it meets the target’s accuracy needs and the extra organ outputs are not required. Retain additional resolution when it changes a decision-relevant prediction.

The comparison is only informative when the predictions are evaluated on the same target, conditions, and benchmark. A model that is calibrated and assessed on the same observations may appear accurate without demonstrating performance for new predictions. Check the time course, not just one summary value, when the intended target is an exposure profile or AUC. For average burden, verify that the tissue set, volumes, and weighting in the calculation match the target definition.

Warning signs include a pooled concentration being labelled as if it were an organ-specific measurement; a model matching plasma while missing tissue observations; or an average that changes materially when the weights are corrected. Large, systematic discrepancies for one organ suggest that the pooled structure may be hiding relevant heterogeneity. Conversely, if organ-resolved outputs are not supported by relevant data or do not improve the target prediction, their apparent detail should not be mistaken for established accuracy.

Common Mistake and Practical Limits

Wrong reasoning: “The organ-resolved model has more compartments, so it must be the better model for every prediction.” More detail can provide outputs unavailable from a pooled model, but it does not guarantee that those outputs are accurate. Correction: judge both structures against the stated prediction target and its required accuracy. If only plasma AUC is needed, organ-specific states may add complexity without improving the relevant prediction. If kidney concentration is needed, plasma fit or a whole-body average cannot establish that kidney concentration is correct.

This approach has limits. A target-specific comparison depends on having suitable observations or a defensible benchmark, and a tolerance appropriate to the intended application. Results may not carry over to different doses or physiological conditions. A structure adequate for one target may fail for another: a model that predicts plasma exposure adequately may not predict an organ concentration, while a model that captures a named organ may not automatically provide a correct average over a differently defined tissue set. State the scope of the conclusion rather than calling one structure universally adequate.

Key takeaway: Match model detail to the output that will be used. Plasma exposure, a defined average tissue burden, and a particular organ’s concentration are distinct targets. Pool tissues only when the pooled structure can preserve the required prediction to an acceptable degree.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A decision uses only plasma AUC over a specified interval. What should determine whether a lumped model is adequate?
  2. Two tissues have volumes of 2 L and 8 L, with concentrations of 6 mg/L and 1 mg/L. Calculate the volume-weighted mean concentration across them.
  3. Why can a pooled concentration fail to answer a question about concentration in one named organ?
  4. A model matches plasma concentration but differs substantially from observed liver concentrations. Is the plasma agreement sufficient to establish suitability for a liver-concentration target?

Solutions. (1) Evaluate its plasma AUC predictions over the relevant interval, conditions, and intended use against an appropriate benchmark and tolerance; organ detail is not required solely because it exists. (2) The weighted mean is \(((2)(6)+(8)(1))/(2+8)=20/10=2.0\ \mathrm{mg/L}\). (3) Pooling combines tissues into an aggregate and can discard the differences needed to identify an individual organ’s value. (4) No. Suitability must be assessed against the liver target; plasma agreement alone does not establish accurate liver concentrations.