Why Recognize Unnecessary Structure?
A PBPK model can distinguish organs or tissues for sound physiological reasons. But a compartment is not justified merely because an organ exists, or because a parameter can be assigned to it. The practical question is whether the available evidence and the prediction task support that distinction. If observations cannot inform tissue-specific parameters, a more detailed structure may give an appearance of precision without better-supported predictions.
This is a question of model adequacy for a specified use, not a general claim that fewer compartments are always preferable. A pooled model may be appropriate for predicting a plasma profile yet inadequate for a question about concentrations in a particular tissue. The introductory material on model structure as a prediction choice established why the target matters; here the focus is on recognizing signs that a proposed level of detail is not earning its place.
A useful assessment starts with a fair comparison. The candidate models should represent the same dosing scenario, prediction target, and observation process. Their parameters should be estimated using the same data and compatible assumptions about residual error. If one model is given different data, a different target, or a more favorable error treatment, the comparison cannot isolate the consequence of its structure.
What Evidence Would Support the Detail?
A tissue-resolved structure is more defensible when relevant observations distinguish the tissues or when physiological information constrains the parameters that govern their behavior. Tissue concentrations, for example, may directly inform tissue-specific distribution. Plasma observations can also inform tissue parameters indirectly when those parameters have distinct, detectable effects on the plasma predictions. Whether they do so depends on the model, the sampling times, and the variability of the measurements.
By contrast, when only a sparse plasma profile is observed, multiple combinations of tissue-specific parameters may explain the data similarly. A model can still produce numerical estimates, but an optimizer returning values is not proof that the values are well informed. Look for signs such as large changes in a tissue parameter producing negligible changes in observed predictions, or coordinated changes in several parameters leaving those predictions nearly unchanged. Estimates that repeatedly hit imposed bounds are another warning: the data may not locate the parameter within its allowed range.
A direct diagnostic is to perturb a parameter and recalculate the predictions for the observations and target of interest. If a substantial parameter change barely changes those predictions, the data or target provide little leverage on that parameter. Conversely, a clear response does not by itself establish that the parameter is identified: other parameters may produce a similar response. This diagnostic is a way to find questions for further assessment, not a substitute for considering measurement error, physiology, and the intended use.
Worked Example: A Tissue Parameter With Little Observable Influence
Hypothetical illustration. Suppose a simplified local approximation to a model's observed plasma summary is \(q=0.4K_A+0.6K_B\), where \(K_A\) and \(K_B\) are dimensionless tissue-specific parameters. The only available observation constrains \(q\) to 1.20. The approximation is used here to illustrate parameter trade-offs; it is not a universal PBPK equation.
One parameter pair, \(K_A=1.50\) and \(K_B=1.00\), gives \(q=0.4(1.50)+0.6(1.00)=0.60+0.60=1.20\). A different pair, \(K_A=0.75\) and \(K_B=1.50\), gives \(q=0.4(0.75)+0.6(1.50)=0.30+0.90=1.20\). Thus this observation alone cannot distinguish these pairs.
More generally, changing the first pair by \(\Delta K_A=0.30\) and \(\Delta K_B=-0.20\) changes the summary by \(\Delta q=0.4(0.30)+0.6(-0.20)=0.12-0.12=0\). The parameters have moved, but the observed summary has not. If the target is this plasma summary, the separate tissue values are weakly constrained by this evidence. If the target is a tissue-specific concentration, that target could differ across the parameter pairs, making tissue observations or other informative constraints important.
Compare Structure Against the Prediction Target
The relevant comparison is not “Which model has more compartments?” but “Does the added distinction matter for the prediction I need, and is that distinction supported?” The prediction target should be stated precisely: for example, plasma concentration at specified times, a tissue concentration, or an exposure summary over a stated interval. Then compare candidate structures on that target using observations appropriate to it.
For observed concentration data, one simple descriptive measure is the root mean squared error (RMSE), when the residuals are on a common scale and the observations are being compared consistently. Define residual \(e_i\) as observed concentration minus predicted concentration at observation \(i\). For \(n\) observations, the RMSE is the square root of the average squared residual. It summarizes discrepancy in the same units as concentration; it does not, by itself, establish that a model is physiologically correct or that a difference is important.
Worked Example: A Small Fit Improvement May Not Justify Extra Tissue Detail
Hypothetical data. Plasma concentrations are measured at 1, 2, 4, and 8 hours, in milligrams per litre. The observed values are 10.0, 7.0, 5.0, and 3.5. A pooled model predicts 9.8, 7.2, 4.8, and 3.7; a model with additional tissue detail predicts 9.9, 7.1, 4.9, and 3.6. Assume the measurements share a comparable error scale and both models use the same observations and prediction times.
For the pooled model, observed-minus-predicted residuals are \(0.2,-0.2,0.2,-0.2\) milligrams per litre. Therefore, \(\mathrm{RMSE}=\sqrt{(0.2^2+(-0.2)^2+0.2^2+(-0.2)^2)/4} =\sqrt{0.16/4}=\sqrt{0.04}=0.200\) milligrams per litre.
For the detailed model, the residuals are \(0.1,-0.1,0.1,-0.1\) milligrams per litre. In particular, at 8 hours the residual is \(3.5-3.6=-0.1\), not \(+0.1\). Its RMSE is \(\sqrt{(0.1^2+(-0.1)^2+0.1^2+(-0.1)^2)/4} =\sqrt{0.04/4}=\sqrt{0.01}=0.100\) milligrams per litre.
The detailed model halves this descriptive RMSE on these four observations. That is a real numerical difference, but it does not alone show that separate tissue parameters are well determined or that the extra detail meaningfully improves the intended prediction. If the target is this plasma profile and typical measurement variation is larger than 0.100 milligrams per litre, the difference may be practically modest. If the target is a tissue-specific concentration, these plasma residuals do not assess that target at all. The conclusion requires target-specific evidence, not just the residual calculation.
For a target that is not directly observed, compare the models' predictions over the relevant scenario and ask whether the difference is large enough to matter for the scientific or decision context. Where possible, use an independent data set or measurements that directly address that target. Do not treat agreement with the data used to fit the model as equivalent to independent validation.
Worked Example: An Added Compartment Does Not Change the Stated Target
Hypothetical scenario. A team wants the predicted plasma AUC from 0 to 24 hours, in milligram-hours per litre. The pooled model predicts 48.0 and the tissue-resolved model predicts 48.5. Assume the target has a prespecified practical tolerance of 2.0 milligram-hours per litre, based on the use of the prediction. No tissue concentration observations are available.
The absolute difference is \(|48.5-48.0|=0.5\) milligram-hours per litre. Relative to the stated tolerance, this is \(0.5/2.0=0.25\), or one quarter of that tolerance. On the specified AUC target, the extra structure has not produced a practically meaningful change under this criterion.
This result does not show that tissue detail is useless for every purpose. The detailed model may predict tissue concentrations that differ substantially from the pooled model, but without tissue observations those predictions are not supported by the evidence described here. If tissue exposure is a future target, suitable tissue data or defensible external constraints could change the assessment. For the present AUC target, the small difference offers little practical reason to retain extra detail by itself.
A Short Justification for Sensitivity Checking
The reason parameter perturbation can reveal weak constraint follows from a first-order approximation. Let \(\boldsymbol{\theta}\) be a vector of model parameters, and let \(\mathbf{f}(\boldsymbol{\theta})\) be the vector of predictions at the observed times. For a small parameter change \(\Delta\boldsymbol{\theta}\), differentiability gives \(\mathbf{f}(\boldsymbol{\theta}+\Delta\boldsymbol{\theta})\approx \mathbf{f}(\boldsymbol{\theta})+J\Delta\boldsymbol{\theta}\), where \(J\) is the matrix of prediction sensitivities, with entries \(J_{ij}=\partial f_i/\partial\theta_j\).
If a parameter change has little effect on the predictions, the relevant sensitivity is small. If a nonzero change satisfies \(J\Delta\boldsymbol{\theta}\approx\mathbf{0}\), different parameter combinations produce nearly indistinguishable predictions near the fitted values. This explains why a detailed model can have weakly constrained parameters even when its predictions fit the observed data. The approximation is local: it diagnoses behavior near the current parameter values and does not prove that predictions are identical everywhere.
Name the output, time range, scenario, and units that matter for the decision or scientific question.
Identify which observations and physiological constraints directly or indirectly inform each proposed tissue distinction.
Perturb tissue-specific parameters and check changes in fitted observations and target predictions; look for compensating parameter changes.
Compare candidate structures on the same target and decide whether any change is meaningful relative to measurement variability and the intended use.
Alternatives, Diagnostics, and Limits
Pooling similar tissues is one alternative when the target and evidence do not require separate representation; earlier material on pooling tissues with similar kinetics explains conditions under which predictions can be preserved. Keeping the resolved structure is more appropriate when distinct tissue behavior matters and can be supported. A third option is to collect informative observations, such as tissue concentrations or measurements at times when candidate structures predict different behavior, rather than forcing a choice from uninformative data.
Useful warning signs include parameters that change greatly across plausible fits while target predictions remain stable, parameter pairs that compensate for one another, estimates repeatedly reaching bounds, and added compartments that leave the target essentially unchanged. These signs do not automatically prove a compartment should be removed. They prompt a closer look at data quality, parameter constraints, numerical fitting, and whether the chosen target can reveal the proposed distinction.
A common mistake is to reason, “The detailed model fits the observed plasma data better, so its tissue parameters must be correct.” The correction is that improved fit on those observations supports only a narrower claim: the detailed model can reproduce those data more closely under the fitting setup. It does not establish that its added parameters are individually constrained, that the improvement matters for the target, or that tissue predictions are reliable. Likewise, a lack of improvement on one target is not a reason to discard detail needed for another target.
Any judgment is conditional on the evidence, physiological assumptions, measurement error, and intended prediction. A simple model can omit a meaningful process; a detailed one can be hard to constrain. The aim is not minimum complexity at any cost, but a structure whose distinctions are supported and useful for the stated prediction.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A tissue-specific parameter changes substantially during fitting, but the observed plasma predictions barely change. What does this suggest, and what does it not establish?
- A detailed model reduces RMSE from 0.30 to 0.20 milligrams per litre. What additional information is needed before concluding that its extra compartments are justified?
- Two models predict a target AUC of 30.0 and 30.4 milligram-hours per litre. The stated practical tolerance is 1.0 milligram-hours per litre. Calculate their difference as a fraction of the tolerance and interpret it.
- Why might a model with weakly supported tissue parameters still be appropriate for a tissue-specific prediction question?
Solutions. (1) The parameter may be weakly constrained by the plasma observations, or it may compensate with other parameters. This does not prove it is physiologically irrelevant or that it cannot affect an unobserved tissue target. (2) Check that both models use the same data, error assumptions, and target; assess whether the RMSE difference is meaningful relative to measurement variability; examine whether tissue parameters are constrained; and consider independent or tissue-specific observations. (3) The difference is \(|30.4-30.0|=0.4\) milligram-hours per litre, and \(0.4/1.0=0.40\) of the tolerance. Under the stated criterion, the difference is less than the practical tolerance, though other evidence may matter. (4) The target may require tissue-level predictions even if existing plasma data do not constrain the parameters. The appropriate response is to recognize the uncertainty and seek relevant tissue data or physiological constraints, not to assume that plasma fit alone validates those predictions.