Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Diagnosing Implausible Obesity-Related PBPK Inputs

Obesity, Aging and Disease States · Tutorial 675 of 799

Diagnosing Implausible Obesity-Related PBPK Inputs

Learn a practical sequence for detecting and tracing implausible obesity-related PBPK inputs while distinguishing bookkeeping errors from genuine physiological variation.

Advanced 11 min read

What You'll Learn

  • Check whether modeled organ and tissue volumes form a coherent fraction of the body-volume convention used.
  • Verify that flows to mutually exclusive destinations sum to cardiac output.
  • Distinguish a tissue’s share of total flow from its flow per unit tissue volume.
  • Detect inappropriate linear scaling of every organ and flow with body weight.
  • Trace failed checks to possible duplicate, omitted, stale, or mismatched inputs.
  • Use independent physiological constraints without forcing values to close a balance artificially.

Why Input Audits Matter

A PBPK simulation can be numerically stable and still represent an implausible person. In obesity-related models, this can happen when inputs assembled from different sources use incompatible body-size descriptors, when a flow is counted twice, or when an organ volume is scaled as if every tissue grew in direct proportion to body weight. These errors affect the model before any drug-specific parameter or exposure prediction is considered.

A useful audit asks several distinct questions. Do the volumes fit the model’s body-volume bookkeeping? Do blood flows leaving the arterial supply sum to cardiac output? Are tissue-specific flows plausible relative to the sizes of those tissues? Finally, do the changes with body weight agree with the physiological scaling assumptions used to construct the individual? A failed check is a signal to investigate, not an automatic diagnosis of which input is wrong.

Definition: For a modeled tissue \(i\), its body-volume fraction is \(f_i=V_i/V_{\mathrm{body}}\), where \(V_i\) is the tissue volume and \(V_{\mathrm{body}}\) is the volume denominator specified by the model. Its flow fraction is \(q_i=Q_i/Q_{\mathrm{CO}}\), where \(Q_i\) is the flow assigned to that tissue and \(Q_{\mathrm{CO}}\) is cardiac output. Tissue-specific perfusion is \(Q_i/V_i\), with units of inverse time when both quantities use volume units.

The denominator and compartment definitions matter. A sum of organ fractions can be expected to equal one only when the listed compartments are an exhaustive, mutually exclusive partition of the model’s stated volume. Some implementations include additional tissues, blood pools, or residual compartments; others use a different volume convention. Likewise, flows should be summed over destinations that partition cardiac output. A total liver flow and one of its component inflows cannot both be counted as separate destinations.

Deriving the Two Balance Checks

Suppose a model represents the body as \(k\) mutually exclusive tissue compartments that together exhaust its defined body-volume denominator. Then \(V_{\mathrm{body}}=\sum_{i=1}^{k}V_i\). Dividing both sides by \(V_{\mathrm{body}}\) gives the volume-fraction identity:

$$ \sum_{i=1}^{k} f_i = \sum_{i=1}^{k}\frac{V_i}{V_{\mathrm{body}}} = \frac{\sum_{i=1}^{k}V_i}{V_{\mathrm{body}}} =1. $$

This is a bookkeeping identity under the stated partition assumption, not a universal rule that every list of organ volumes must sum to an anatomically exact whole-body volume. If the modeled list is incomplete, the remainder must be represented explicitly or the expected sum must reflect the model’s convention.

For a steady circulation in which the listed, mutually exclusive destinations receive all cardiac output, conservation of flow gives \(Q_{\mathrm{CO}}=\sum_{i=1}^{k}Q_i\). Dividing by \(Q_{\mathrm{CO}}\) yields \(\sum_i q_i=1\). If a row is duplicated with flow \(\delta\), the apparent sum becomes \(1+\delta/Q_{\mathrm{CO}}\). Thus the size of a flow-balance discrepancy can help identify a duplicated or misclassified component, although it does not prove that this is the cause.

Worked Example: Auditing a Body-Volume Partition

Worked Example: Auditing a Body-Volume Partition

Hypothetical model inputs. Consider an adult with body weight 120 kg. For this audit only, the model defines its body-volume budget as 120 L, using a simplifying one-kilogram-per-litre convention. Its exhaustive compartment list is liver 2.2 L, kidneys 0.45 L, brain 1.4 L, heart 0.45 L, muscle 32 L, adipose tissue 72 L, and other tissue 17.5 L. An independent body-composition input used by the same model gives adipose tissue as 40 L.

Check the total. Summing the listed compartment volumes gives \(2.2+0.45+1.4+0.45+32+72+17.5=126\) L. Relative to the stated 120 L budget, the total fraction is \(126/120=1.05\), or 105%. The partition exceeds its defined budget by \(126-120=6\) L.

Inspect the fractions. The listed adipose fraction is \(72/120=0.600\), or 60.0%. The independently supplied adipose volume corresponds to \(40/120=0.333\), or 33.3%. These values are inconsistent with one another; the listed adipose input is a strong place to investigate, but the total alone does not establish that adipose is the sole source of the error.

Replacing 72 L with 40 L gives a new listed total of \(126-72+40=94\) L, leaving \(120-94=26\) L unallocated. That replacement does not produce a complete partition. The appropriate finding is that the assembled inputs are inconsistent: the adipose values disagree, and the remaining compartments also fail to account for the stated budget. The modeler should check compartment definitions, source records, and whether a residual compartment has been omitted rather than altering a value merely to force the sum to 120 L.

Checking Cardiac Output and Tissue Perfusion

The flow audit begins with a clear destination list. Confirm whether each row represents a complete organ flow or a component of another row’s flow. Then sum the destination flows and compare them with the cardiac output assigned to that individual. Record the discrepancy in both absolute units and as a fraction of cardiac output. A total near the target can still conceal offsetting errors, so row definitions and units need review as well.

A flow fraction and tissue-specific perfusion answer different questions. \(q_i=Q_i/Q_{\mathrm{CO}}\) describes how much of total cardiac output is assigned to tissue \(i\). \(Q_i/V_i\) describes flow relative to that tissue’s volume. A large tissue can receive a substantial share of total flow without having an unusually high flow per unit volume; conversely, a small tissue can have a small flow fraction but high tissue-specific perfusion. Compare the latter with the model’s intended physiological values and units, not with an arbitrary universal cutoff.

Worked Example: Finding a Duplicate Flow

Hypothetical inputs. A model assigns cardiac output \(Q_{\mathrm{CO}}=7.5\) L/min. Its flow table lists total liver flow 1.5 L/min, hepatic arterial flow 0.3 L/min, kidney flow 1.0 L/min, brain flow 0.75 L/min, heart flow 0.25 L/min, muscle flow 1.25 L/min, adipose flow 0.90 L/min, skin flow 0.65 L/min, and other flow 1.2 L/min. The table’s intended destinations are mutually exclusive.

Sum the recorded rows. The total is \(1.5+0.3+1.0+0.75+0.25+1.25+0.90+0.65+1.2=7.8\) L/min. The excess is \(7.8-7.5=0.3\) L/min, and the ratio of listed flow to cardiac output is \(7.8/7.5=1.04\). The apparent flow fractions therefore sum to 1.04 rather than 1.

Inspect the definitions. The hepatic arterial flow is listed separately even though total liver flow is already included. If the 0.3 L/min hepatic arterial flow is a component of the 1.5 L/min total liver flow, counting both duplicates 0.3 L/min. Removing the component row from the destination sum gives \(7.8-0.3=7.5\) L/min, matching cardiac output. This is a plausible bookkeeping explanation; it should be confirmed against the model’s flow definitions before the table is changed.

For example, the adipose flow fraction is \(0.90/7.5=0.120\), or 12.0%. If adipose volume is 38 L, its tissue-specific perfusion is \(0.90/38=0.0237\) min\(^{-1}\). Those two summaries are not interchangeable: the first is a share of cardiac output, while the second is flow per unit tissue volume.

Testing Scaling With Body Weight

A common input error is to multiply every reference volume and flow by the same ratio of body weights. That operation assumes direct proportionality to total body weight for every quantity. Obesity-related physiology does not generally justify that blanket assumption: adipose tissue, lean tissues, individual organs, and cardiac output need not change in the same proportion. The earlier material on selecting body-size descriptors and constructing consistent physiological scaling provides the modeling context; here the aim is to test whether the assembled values obey their stated scaling rules.

For each quantity, compare the obese individual with the reference individual using the ratio \(R_i=V_{i,\mathrm{obese}}/V_{i,\mathrm{reference}}\) for volumes or \(R_Q=Q_{\mathrm{CO,obese}}/Q_{\mathrm{CO,reference}}\) for cardiac output. Compare these ratios with the specified physiological model, not automatically with the body-weight ratio. Also check that related inputs come from compatible individuals and definitions: an organ volume from one body-size method paired with a cardiac output from another can create a plausible-looking but incoherent set.

Worked Example: Detecting Blanket Body-Weight Scaling

Hypothetical inputs. A 70 kg reference adult has liver volume 1.5 L, kidney volume 0.30 L, brain volume 1.4 L, muscle volume 28 L, adipose volume 2 L, and cardiac output 5.0 L/min. A separately specified physiological model for a 140 kg adult gives liver 2.0 L, kidneys 0.36 L, brain 1.4 L, muscle 34 L, adipose tissue 38 L, and cardiac output 6.5 L/min. These are illustrative model inputs, not universal physiological reference values.

The body-weight ratio is \(140/70=2.00\). If every input were scaled linearly by that ratio, the predicted values would be liver \(1.5(2)=3.0\) L, kidneys \(0.30(2)=0.60\) L, brain \(1.4(2)=2.8\) L, muscle \(28(2)=56\) L, adipose tissue \(2(2)=4.0\) L, and cardiac output \(5.0(2)=10.0\) L/min.

Compare the blanket-scaled values with the specified model values. The liver ratio is \(2.0/1.5=1.33\); the kidney ratio is \(0.36/0.30=1.20\); the brain ratio is \(1.4/1.4=1.00\); the muscle ratio is \(34/28=1.21\); the adipose ratio is \(38/2=19.0\); and the cardiac-output ratio is \(6.5/5.0=1.30\). These differ markedly from the blanket factor of 2.00. In particular, doubling adipose tissue from a small reference value gives 4.0 L, far below the specified 38 L, whereas doubling brain volume gives 2.8 L, twice the specified 1.4 L.

The diagnosis is not that one ratio is physiologically correct for every model. It is that blanket scaling conflicts with the specified individual-level inputs and should not be used as a substitute for the model’s tissue-specific rules. The inputs must be traced to the intended scaling method before simulation.

A Practical Diagnostic Sequence

When an audit fails, avoid changing values until the source of the discrepancy is understood. A useful sequence separates arithmetic checks from physiological judgment:

1
Confirm definitions and units.
Check whether volumes are tissue volumes, body-volume fractions, or mass-derived estimates, and whether flows are in consistent units. Identify what the model includes in its volume denominator and flow destinations.
2
Recompute balances.
Sum the applicable volume fractions and organ flows from the raw inputs. Retain the unrounded values for arithmetic; rounding can obscure small discrepancies but cannot explain a material imbalance.
3
Inspect individual ratios.
Compare tissue volumes and flows with their model-specific expectations and with the corresponding reference individual. Check both flow fraction and flow per unit tissue volume.
4
Trace discrepancies to their source.
Look for duplicated component flows, omitted compartments, stale reference values, mismatched body-size methods, and unintended unit conversions. Recalculate after correcting a confirmed source error.

Several approaches can support this sequence. A strict balance check is effective for detecting bookkeeping inconsistencies when the model’s compartments are exhaustive. A comparison with tissue-specific physiological inputs is more informative when the model intentionally includes a residual compartment or uses a nonstandard denominator. A reference-ratio comparison helps reveal inappropriate scaling, but it does not establish that an individual value is wrong unless the model’s scaling assumptions are known. These checks complement one another; none replaces review of definitions and provenance.

Diagnostics also have limits. Biological variation means that an unusual value is not automatically an error. A balance can close even when two inputs are both wrong in offsetting directions. Conversely, a sum below one may reflect an omitted modeled compartment rather than a faulty recorded value. These audits assess internal consistency and plausibility, not whether the PBPK model predicts clinical exposure accurately.

Common Mistakes and Limitations

  • Forcing a sum to one. Adjusting one volume or flow simply to close a balance hides the source of the inconsistency. First establish that the list is exhaustive and that the quantities have compatible definitions.
  • Counting components and totals together. A vascular component such as hepatic arterial flow must not be added as a separate destination when the listed organ flow already includes it.
  • Treating body weight as a universal scaling factor. The same weight ratio need not apply to adipose tissue, organs, muscle, and cardiac output. Use the specified tissue-specific assumptions.
  • Confusing flow fraction with perfusion. A tissue’s percentage of cardiac output does not describe flow per unit tissue volume. Calculate each quantity with its own denominator and units.
  • Using a balance as proof of physiological truth. A perfectly closed set of inputs can still be implausible. Independent physiological comparisons and input provenance remain necessary.
Key takeaway: Treat volume closure, flow conservation, tissue-specific perfusion, and body-weight scaling as complementary audits. A failed check identifies an inconsistency to trace; it does not, by itself, identify the correct replacement value.

Check Your Understanding

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

  1. A model’s exhaustive compartment list totals 108 L against a defined body-volume denominator of 120 L. What fraction is accounted for, and what should be checked before concluding that a compartment is missing?
  2. Cardiac output is 6.0 L/min, while mutually exclusive destination flows sum to 6.6 L/min. What is the relative excess, and name one bookkeeping error to investigate.
  3. A tissue receives 0.4 L/min of flow and has volume 20 L. Calculate its tissue-specific perfusion and state its units.
  4. Body weight doubles, but a model’s specified brain volume is unchanged. Why does this not necessarily indicate an error?
  5. Why is it unsafe to change one organ flow until the total matches cardiac output?

Solutions. (1) The accounted fraction is \(108/120=0.900\), or 90.0%. First confirm the denominator, units, and whether the list is truly exhaustive under the model’s compartment definitions. (2) The excess is \(6.6-6.0=0.6\) L/min, or \(0.6/6.0=0.100\), a 10.0% excess. A duplicated component flow is one possibility to check. (3) \(0.4/20=0.020\) min\(^{-1}\). (4) Different tissues need not scale in direct proportion to total body weight; the value should be assessed against the model’s specified tissue-specific rule. (5) The mismatch may arise from a duplicated row, an omitted or misclassified destination, or inconsistent definitions. Arbitrarily changing a flow can hide the actual cause.