Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Plausibility Audits for PBPK Parameter Sets

Biological Plausibility and Parameter Uncertainty · Tutorial 155 of 799

Plausibility Audits for PBPK Parameter Sets

Learn to screen PBPK parameter sets for implausible values, unit errors, double counting, and mismatches between a value and its source.

Advanced 12 min read

What You'll Learn

  • Separate hard domain checks from context-dependent physiological screening.
  • Check organ masses and tissue fractions against body mass without double counting.
  • Reconcile blood-flow values with their units, network structure, and cardiac output.
  • Audit drug-input units, ranges, and source metadata before model use.
  • Record a flag as a question to investigate, not automatic proof that a parameter is wrong.

Why Audit a Parameter Set?

A PBPK model can be coded correctly and still give misleading predictions if its inputs are implausible, inconsistently defined, or transcribed with the wrong units. A plausibility audit is a structured review of the parameter record before model evaluation. It asks whether each value has a clear meaning and source, falls within a reasonable domain, and fits with the other inputs that describe the same physiology or compound.

The audit is a screening step, not a method for estimating the “true” value. Physiological quantities vary with age, body size, health, posture, and measurement conditions. A value outside a broad reference range should prompt investigation: it might reflect a special population, a different tissue definition, or an error. It should not be rejected solely because it differs from a familiar reference value.

Definition: A plausibility audit is a documented review of PBPK inputs for valid domains, appropriate units and definitions, credible source context, and consistency with basic physiological or chemical constraints. An audit flag identifies an input that needs explanation or correction; it is not, by itself, evidence that the value is false.

The audit assumes that the model’s intended population and parameter definitions are specified. For example, “liver volume” might refer to an anatomical volume, a tissue mass converted using a density, or an effective model compartment. Those are not automatically interchangeable. Likewise, blood flow may refer to total hepatic inflow or only hepatic arterial flow. The audit must preserve the definitions used by the model rather than compare superficially similar labels.

A Structured Audit

A useful audit proceeds from simple, objective checks toward more contextual judgments. Begin with data integrity: confirm parameter names, units, species, population, and source. Then check hard mathematical domains, such as positive organ mass and fractions between zero and one. Next, compare magnitudes with broad physiological expectations and examine whether quantities with shared definitions can be reconciled. Finally, record what is flagged, why, and what evidence would resolve it.

1
Confirm meaning and provenance.
Record the parameter definition, units, population, measurement conditions, source, and any transformation used to obtain the model value.
2
Check domains and units.
Identify impossible signs or bounds, and convert units explicitly before comparing values or adding quantities.
3
Check physiological scale and closure.
Compare with an appropriate reference context and test basic totals, while ensuring that compartments or flows are not counted twice.
4
Classify and document each flag.
Distinguish a definite data error from an unexplained difference or a population-specific value; retain the source and rationale for any change.

Organ Volumes and Tissue Composition

Organ volume and tissue composition checks begin with units and definitions. If the model stores tissue mass but a source reports volume, conversion requires a stated density and compatible definitions. A mass in grams must not be compared directly with a volume in liters. For an adult reference subject, organ masses such as liver around 1.5 kg, brain around 1.4 kg, and total kidneys around 0.3 kg are useful only as rough orientation; values depend on the individual and on what the measurement includes.

For a set of mutually exclusive tissue compartments, their masses should sum to the modeled body mass, apart from explicitly represented residual compartments and rounding. If the model uses fractions of body mass, each fraction is dimensionless and lies between zero and one; fractions for a complete, non-overlapping partition sum to one. These checks are valid only when the listed compartments are truly disjoint. Adding an organ mass to a tissue category that already includes that organ creates an apparent excess even if neither input is individually wrong.

Audit principle: A total is meaningful only when its components share units, refer to the same subject and definition, and do not overlap. Before interpreting a sum as inconsistent, verify the partition represented by the model.

Worked Example: Checking a Tissue-Mass Partition

Consider a hypothetical 70 kg adult model. Its mutually exclusive compartments are entered as liver 1.50 kg, kidneys 0.30 kg, heart 0.32 kg, lungs 0.50 kg, brain 1.40 kg, skeletal muscle 28 kg, adipose tissue 12 kg, and other tissues 26 kg. Check the total and the corresponding fractions.

The summed mass is

$$ 1.50+0.30+0.32+0.50+1.40+28+12+26=70.02\ \mathrm{kg}. $$

The difference from the stated body mass is \(70.02-70.00=0.02\) kg, or 20 g. This is small relative to the rounded inputs and is not, on its own, a meaningful discrepancy. The corresponding total fraction is

$$ \frac{70.02\ \mathrm{kg}}{70.00\ \mathrm{kg}}=1.0003. $$

If the skeletal-muscle entry were instead 38 kg, the total would be 80.02 kg and the total fraction would be \(80.02/70.00=1.143\). That is a strong audit flag for a supposedly exhaustive, non-overlapping partition. It suggests a possible transcription error, a mismatch in the population or definition, or double counting. The arithmetic does not identify which explanation is correct, so the source and compartment definitions must be checked.

Blood Flows: Check Units and Network Meaning

Blood-flow inputs should be compared in compatible units and with attention to how blood moves through the modeled network. Resting cardiac output in a typical adult is often on the order of 4–7 L/min, but body size, activity, clinical state, and other conditions matter. Approximate values such as renal flow near 1.0–1.3 L/min, cerebral flow around 0.7–0.9 L/min, and total hepatic inflow around 1.2–1.8 L/min are broad screening references, not universal acceptance limits.

A common structural problem is treating serial and parallel flows as though they were all independent branches. Blood passing from the gut through the portal vein to the liver is part of hepatic inflow. If a parameter list includes total hepatic flow and also adds portal-drained gut flow as a separate independent systemic outflow, the same flow may be counted twice. Check the diagram and definitions before using a total as a consistency test.

Worked Example: Investigating an Apparent Flow Excess

A hypothetical adult model lists cardiac output as 5.0 L/min and gives flows of liver 1.5 L/min, kidneys 1.1 L/min, brain 0.75 L/min, muscle 0.90 L/min, skin 0.30 L/min, gut 1.0 L/min, and other tissues 0.45 L/min. The values appear to be independent organ outflows. Add them and compare with cardiac output.

The listed sum is

$$ 1.5+1.1+0.75+0.90+0.30+1.0+0.45=6.00\ \mathrm{L/min}. $$

This is 1.00 L/min above the stated cardiac output, and \(6.00/5.00=1.20\), or 120% of cardiac output. That is not a plausible sum of simultaneous, independent systemic outflows under the stated definitions. However, the excess does not prove that an individual flow value is wrong. If the 1.5 L/min liver value is total hepatic inflow and the 1.0 L/min gut value represents portal flow entering the liver, then those entries overlap. The audit flag is therefore “check flow topology and definitions,” not “reduce a flow by 1.0 L/min.” If instead liver flow means hepatic arterial flow only, the interpretation and reconciliation change. The parameter metadata and model connections determine the appropriate conclusion.

Drug Inputs: Domains, Units, and Source Context

Drug input audits begin with the same basic questions: what is the measured or predicted quantity, what units are expected, and does the value fall within its mathematical domain? For example, molecular weight must be positive and needs units such as g/mol; a fraction unbound is dimensionless and lies between zero and one; and an acid dissociation constant is often represented as pKa, a logarithmic quantity. A numerical value cannot be judged from its magnitude alone if its unit or encoding is unclear.

Do not infer a universal “plausible” range for every chemical property. A wide range can be legitimate across compounds, and values from different methods may not represent the same conditions. Instead, verify the compound identity, assay or prediction method, experimental conditions, and any conversions. The earlier material on binding and blood-to-plasma partitioning provides definitions for those inputs; in an audit, preserve those definitions and check the recorded value and representation before attempting a mechanistic interpretation.

Worked Example: Finding a Fraction and Unit-Scale Error

A hypothetical compound record gives molecular weight as “350 mg/mol,” pKa as 7.4, logP as 2.1, plasma fraction unbound as 8.0 (fraction), and tissue fraction unbound as 0.02. Identify entries that need review.

The molecular-weight number may have been copied with the wrong unit. If the intended value is 350 g/mol, the equivalent is \(350{,}000\) mg/mol, not 350 mg/mol. Conversely, 350 mg/mol is 0.350 g/mol. The record’s units therefore imply a thousand-fold difference, so the compound source must be checked before the parameter is used.

A fraction unbound of 8.0 is outside the required domain \(0\leq f_u\leq1\), so it cannot be used as a fraction. If “8.0” was reported as a percentage, the corresponding fraction is \(8.0/100=0.080\). The tissue value 0.02 is within the mathematical domain, but domain validity alone does not verify its source or conditions. The pKa and logP values are not rejected merely because of their numerical size; their compound identity, definition, and source still need confirmation.

Why the Checks Work—and Where They Stop

Mass closure follows from the definition of a complete partition. If \(M_i\) is the mass of compartment \(i\), and the compartments are mutually exclusive and exhaustive for a subject of mass \(M\), then their masses must add to that subject’s total mass. Defining \(w_i=M_i/M\), division by \(M\) gives the corresponding fraction constraint.

$$ \sum_i M_i=M \quad\Longrightarrow\quad \sum_i w_i=\sum_i\frac{M_i}{M}=\frac{\sum_i M_i}{M}=1. $$

The argument depends on the assumptions stated: the compartments refer to the same subject, use compatible mass definitions, and neither overlap nor omit mass. A numerical total near one does not establish that each compartment is correct; opposing errors can cancel. Similarly, a flow sum can diagnose a mismatch only after the relevant flows have been defined as independent branches over the same interval.

A plausibility audit is complementary to, rather than a substitute for, other approaches. A reference-range comparison is quick but may overlook unusual populations and definitional differences. A formal mass-balance or network reconciliation is more informative when the model structure is explicit, but it cannot repair inaccurate source data automatically. Expert review can interpret context that a simple rule misses, though it should be documented so the decision is reproducible. For uncertain inputs, retain the evidence and uncertainty rather than changing values merely to make the model look tidy.

Diagnostics and Common Pitfalls

Several patterns signal that an audit is not yet reliable: the same parameter appears with different units in different files; fractions exceed one or are entered as percentages without conversion; organ totals are far above or below the subject’s body mass; flows appear to exceed cardiac output; or a value has no traceable source or population definition. These are reasons to pause, inspect the raw record, and ask whether the discrepancy is a unit, definition, topology, or transcription issue.

  • Rejecting a value just because it is outside a familiar range. A reference value may not fit the modeled age, body size, or clinical state. Confirm context before rejecting it.
  • Adding overlapping compartments. A tissue category may already include an organ or its blood content. Establish a non-overlapping partition before interpreting its total.
  • Comparing unlike flow definitions. Total hepatic inflow and hepatic arterial flow are not interchangeable. Check the network before summing.
  • Confusing a percentage with a fraction. Eight percent is 0.08 as a fraction, not 8.0.
  • Treating an audit flag as a correction. A flag identifies a question. Any change needs a defensible source or explicit rationale.
Key takeaway: A credible PBPK parameter set needs more than plausible-looking numbers. Every input should have a defined meaning, compatible units, a traceable source, and a documented check against relevant physiological or mathematical constraints.

Check Your Understanding

Use the definitions and audit principles above. Answers follow the questions.

  1. A complete set of non-overlapping tissue fractions sums to 1.08. What should you check before changing any entries?
  2. A model records cardiac output as 5.5 L/min and one organ flow as 5,500 mL/min. Are the values necessarily inconsistent? Explain.
  3. A plasma fraction unbound is reported as 12%, but the model field expects a fraction. What value should be entered?
  4. Why does an organ-mass total of 103% of body mass not, by itself, identify which parameter is incorrect?
  5. A total hepatic inflow and portal gut flow are both included in a sum of independent outflows. What is the first audit question?

Solutions. (1) Check whether all fractions refer to the same subject and denominator, whether the categories are mutually exclusive, and whether rounding or a transcription error explains the excess. (2) Not necessarily: \(5{,}500\ \mathrm{mL/min}=5.5\ \mathrm{L/min}\). The units must be converted before comparison, and the organ-flow definition still matters. (3) Enter \(12/100=0.12\). (4) The excess identifies a partition problem, but several entries could be wrong, compartments could overlap, or definitions could differ; the total alone cannot locate the cause. (5) Determine whether portal gut flow is already included in total hepatic inflow, and inspect the model’s flow topology before summing them as independent outflows.