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.”
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
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
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,
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:
For the lumped model, the errors are \(0.1,-0.1,0.2,-0.2\ \mathrm{mg/L}\). Thus,
For the organ-resolved model, the errors are \(0.05,0.05,0.1,-0.1\ \mathrm{mg/L}\), giving
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.
Name the quantity, anatomical scope, time or interval, conditions, and practical tolerance. Define any average with its weights.
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.
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.
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.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A decision uses only plasma AUC over a specified interval. What should determine whether a lumped model is adequate?
- 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.
- Why can a pooled concentration fail to answer a question about concentration in one named organ?
- 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.