Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Provenance for Derived and Model-Dependent Inputs

Parameter Sourcing, Provenance and Documentation · Tutorial 194 of 799

Provenance for Derived and Model-Dependent Inputs

You will learn to trace a PBPK input from its source measurements through its calculations and assumptions, while clearly labeling scaled and model-dependent values.

Intermediate 12 min read

What You'll Learn

  • Distinguish directly measured values from calculated, scaled, and model-estimated inputs.
  • Record the source, units, equation, assumptions, and applicability of each derivation.
  • Scale intrinsic clearance using enzyme abundance and physiological scaling factors.
  • Derive a pH-dependent distribution coefficient and document its chemical assumptions.
  • Recognize when a parameter is conditional on a model structure or fitted quantity.
  • Identify provenance gaps that can make an apparently precise input difficult to reproduce.

Why Derived Inputs Need a Traceable Lineage

A PBPK input may look like a single number in a parameter table, but it can represent a chain of evidence and decisions. A value reported by an assay may be transformed into a tissue-level quantity, scaled to an organ, or estimated conditional on a model equation. If that chain is hidden, a reviewer cannot tell what the number means, whether its assumptions apply, or which part should change when new evidence becomes available.

The key distinction is not simply “published” versus “unpublished.” A published value can still be a model-derived estimate, while a quantity calculated from an assay may be reproducible if its inputs, equation, and assumptions are recorded. The provenance of an input is its traceable account of where the evidence came from and how it became the value used in the model.

Definition: A measured input is a quantity observed or reported from an experiment or other measurement process. A derived input is calculated from one or more reported quantities using a stated equation. A scaled input is derived by applying a scaling factor to translate an assay-level quantity to a different biological amount or system. A model-dependent input is estimated or defined conditional on a specified model structure, equations, and assumptions. These categories can overlap: a scaled value is derived, and a fitted estimate is model-dependent.

A useful provenance record therefore identifies the value actually used, its unit and status, the primary evidence for every source quantity, the calculation or model that connects those quantities to the input, and the assumptions that make that connection defensible. Earlier material on Documenting Transformations From Measurement to Input addresses the general transformation record; here the emphasis is on making the derivation’s evidentiary and model dependence explicit.

What to Record for a Derived or Model-Dependent Value

A record should let another analyst reconstruct the input without guessing. For each quantity in the calculation, preserve its reported value and qualifier, unit, experimental system, and source. Then give the equation and substitution, including unit conversions. State the assumptions that connect the source measurement to the target quantity, and say whether the input is measured, calculated, scaled, fitted, or some combination. Finally, record applicability limits: for example, the enzyme, species, tissue, pH, or concentration range to which the source evidence applies.

The source of an equation is not necessarily the source of its inputs. A publication might provide a scaling equation, while separate assay reports provide the activity and abundance values inserted into it. Both links belong in the lineage. If a value is chosen from among plausible alternatives, record that decision and its rationale rather than presenting the choice as though it were measured.

Provenance elementWhat it answersExample entry
Input and unitWhat quantity enters the model?Hepatic intrinsic clearance, mL/min
Evidence and sourceWhere did each starting value come from?Enzyme abundance reported in a specified microsomal preparation
Status and derivationWas it observed, calculated, scaled, or fitted?Scaled from an assay value by an abundance ratio
Assumptions and limitsWhen is the transformation defensible?Activity assumed proportional to abundance for the enzyme and substrate

Scaling an Enzyme-Dependent Assay Result

A common derivation translates an intrinsic clearance measured per amount of microsomal protein to a target tissue abundance. If the assay clearance is attributable to a particular enzyme, and activity per unit of that enzyme is assumed unchanged between the assay preparation and target tissue, an enzyme-abundance ratio can scale the assay value. Further physiological factors can then scale the result to an organ.

$$ CL_{\mathrm{int,target\ per\ mg}} = CL_{\mathrm{int,assay\ per\ mg}} \left(\frac{A_{\mathrm{target}}}{A_{\mathrm{assay}}}\right) $$

Here \(A\) is enzyme abundance per milligram of microsomal protein in the stated preparation. The ratio is dimensionless only when both abundances use compatible units and refer to the same enzyme. Multiplication by microsomal protein per gram of tissue and tissue mass then yields an organ-level clearance. This scaling is not itself a new measurement of organ clearance; it is a calculation conditional on the abundance and transfer assumptions.

Worked Example: Scaling Microsomal Clearance to a Liver

Suppose a hypothetical assay reports \(8.00\ \mu\text{L/min/mg}\) microsomal protein for a drug’s enzyme-specific intrinsic clearance. Enzyme abundance is \(0.25\ \text{pmol/mg}\) in the assay preparation and \(0.50\ \text{pmol/mg}\) in the target liver microsomes. For illustration, assume \(40.0\ \text{mg}\) microsomal protein per gram of liver and a liver mass of \(1{,}500\ \text{g}\). These are hypothetical values, not measurements from a real study.

First calculate the abundance ratio and the target clearance per milligram of microsomal protein:

$$ \frac{A_{\mathrm{target}}}{A_{\mathrm{assay}}} = \frac{0.50\ \text{pmol/mg}}{0.25\ \text{pmol/mg}} =2.00, \qquad CL_{\mathrm{int,target\ per\ mg}} = 8.00\times 2.00 =16.0\ \mu\text{L/min/mg} $$

Now scale by microsomal protein per gram and liver mass:

$$ 16.0\ \frac{\mu\text{L}}{\text{min}\cdot\text{mg}} \times 40.0\ \frac{\text{mg}}{\text{g}} \times 1{,}500\ \text{g} = 960{,}000\ \mu\text{L/min} = 960\ \text{mL/min} $$

The derived organ-level value is \(960\ \text{mL/min}\). Its record should retain the assay clearance, both abundance values, the microsomal protein scaling factor, liver mass, units, and their sources. It should also state the assumptions: the clearance is attributable to the enzyme whose abundance is compared; activity scales linearly with abundance; the preparations are comparable; and the physiological scaling factors represent the modeled population. If any of those conditions is not supported, the calculated value may not represent the target liver.

Documenting Equation-Based Physicochemical Inputs

A physicochemical input can also be derived from measured or estimated molecular properties using an equation. The resulting value depends on both the source quantities and the chemical assumptions built into that equation. A record should not label the output as a direct measurement merely because its inputs came from experiments.

Worked Example: Deriving a pH-Specific Distribution Coefficient

Consider a hypothetical monoprotic weak acid with a neutral-species partition coefficient \(P=100\), \(pK_a=4.40\), and a specified aqueous pH of \(6.40\). Assume that only the neutral form partitions appreciably into the organic phase and that the phases are at equilibrium. Under this simplified assumption, the distribution coefficient is:

$$ D(\mathrm{pH}) = \frac{P}{1+10^{\mathrm{pH}-pK_a}} $$

Substituting the hypothetical values gives:

$$ D(6.40) = \frac{100}{1+10^{6.40-4.40}} = \frac{100}{1+100} = \frac{100}{101} \approx 0.9901 $$

Therefore, \(\log_{10}D(6.40)\approx -0.00432\). This is a calculated, pH-specific value, not a pH-independent measurement of partitioning. The provenance record should identify the source and status of \(P\) and \(pK_a\), the specified pH, the equation, and the assumption that ionized species make negligible contribution to organic-phase partitioning. If that assumption is inappropriate, this equation does not support the reported result; a more suitable chemical model or direct measurement would be needed.

When an Input Depends on the PBPK Model

Some quantities used as inputs are inferred through model equations rather than observed directly. In that case, a numerical estimate is conditional on the model structure and on the values of other parameters. This conditionality should be visible in the record, especially if the estimated quantity is later treated as fixed in a simulation.

For example, under the well-stirred liver model, hepatic availability \(F_h\) and intrinsic clearance \(CL_{\mathrm{int}}\) are related by:

$$ F_h = \frac{Q_h}{Q_h+f_{u,b}CL_{\mathrm{int}}} $$

Here \(Q_h\) is hepatic blood flow and \(f_{u,b}\) is the unbound fraction in blood. Rearranging gives an estimate of intrinsic clearance from an estimated \(F_h\):

$$ CL_{\mathrm{int}} = \frac{Q_h}{f_{u,b}} \left(\frac{1}{F_h}-1\right) $$

Worked Example: Recovering a Model-Dependent Clearance

Suppose a hypothetical analysis estimates hepatic availability as \(F_h=0.750\), using a specified model and data. Let hepatic blood flow be \(90.0\ \text{L/h}\) and the unbound fraction in blood be \(0.100\). Treat these as fixed for this calculation. The estimated \(F_h\) is not a direct measurement of intrinsic clearance; it is an estimate conditional on the analysis.

Substitute into the rearranged well-stirred equation:

$$ CL_{\mathrm{int}} = \frac{90.0\ \text{L/h}}{0.100} \left(\frac{1}{0.750}-1\right) = 900\ \text{L/h}\times 0.3333 \approx 300\ \text{L/h} $$

A check by substitution gives \(F_h=90.0/(90.0+0.100\times300)=90.0/120=0.750\). Thus \(300\ \text{L/h}\) is consistent with the stated model and inputs. Its provenance must include the estimate of \(F_h\), the data and model used to obtain it, \(Q_h\), \(f_{u,b}\), the equation, and the assumptions that the well-stirred representation and stated conditions apply. It should not be described as an experimentally measured \(CL_{\mathrm{int}}\).

Justifying the Scaling and Recognizing Failure

The abundance scaling equation follows from a proportionality assumption. If enzyme-specific intrinsic clearance per milligram of microsomal protein is proportional to enzyme abundance, write \(CL_{\mathrm{int}}=kA\), where \(k\) is activity per unit abundance under the relevant assay conditions. For assay and target preparations, this gives \(CL_{\mathrm{int,assay}}=kA_{\mathrm{assay}}\) and \(CL_{\mathrm{int,target}}=kA_{\mathrm{target}}\). Dividing the second relation by the first cancels \(k\), provided it is the same in both preparations, yielding the abundance-ratio formula. If substrate conditions, enzyme activity, or contributions from other enzymes differ, that cancellation is not justified and the ratio alone does not establish the target clearance.

Several checks can expose a weak or broken provenance chain:

  • Check dimensions. Confirm that units cancel as claimed and that the final quantity has the unit required by the model. An abundance ratio should be dimensionless only after compatible units are verified.
  • Check definitions and scope. A value for one enzyme, species, tissue, or assay condition cannot silently become a value for another. Preserve qualifiers such as “total,” “unbound,” or “per milligram.”
  • Check arithmetic independently. Recalculate each step from the recorded values, including conversions between microliters and milliliters or between minutes and hours.
  • Check sensitivity to assumptions. If modest changes in abundance, binding, or a model estimate cause large changes in the derived input, report that dependence rather than presenting excess precision.
  • Check for circularity. If a quantity was estimated by fitting the same model that later uses it, identify the data and model connection. Do not imply that the estimate is independent experimental confirmation.

A direct measurement may avoid some scaling assumptions but may not exist for the target tissue or population. An empirical scaled value may be more feasible, but depends on transferability assumptions. A mechanistic or fitted estimate can use several evidence sources coherently, yet is conditional on model structure and other inputs. These are alternatives with different evidence and limitations, not interchangeable labels for the same certainty.

Common Provenance Pitfalls

Common mistake: Reporting a calculated value as though it were directly measured because every number in the calculation came from a publication. Correction: Label the output according to how it was obtained, show the equation and substitutions, and preserve the source and status of each input.

Another frequent error is recording only the final value and a citation. That does not reveal which reported quantity was selected, how it was transformed, whether the units were compatible, or what assumptions were added by the analyst. A citation is a path to evidence, not a substitute for documenting the derivation. Similarly, a value labeled “scaled” is not fully interpretable unless the reference system, target system, scaling factor, and transfer assumptions are stated.

The practical standard is reproducibility with honest qualification: another modeler should be able to reconstruct the calculation and understand what would invalidate it. This does not require claiming that every assumption is certain. It requires distinguishing observed evidence from analyst choices and model conditions, and stating where uncertainty or applicability limits enter the lineage.

Check Your Understanding

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

  1. A clearance value is calculated by multiplying an assay result by a tissue-to-assay enzyme-abundance ratio. Is the final value measured, derived, or both? What assumption is central to this scaling?
  2. A parameter was estimated by fitting a PBPK model to concentration data. Name two items that should be recorded so its model dependence is clear.
  3. If a distribution coefficient calculation uses a pH that differs from the intended model condition, can the reported coefficient be treated as the same input? Explain.
  4. An assay clearance is \(5.00\ \mu\text{L/min/mg}\), and the target-to-assay abundance ratio is \(3.00\). What is the scaled clearance per milligram, and what units does it have?

Solutions. (1) It is derived and scaled, not directly measured in the target tissue. The key assumption is that enzyme-specific activity is proportional to abundance and transfers between the assay and target preparations. (2) Record the data and analysis that produced the estimate, the model structure and equations, and the other fixed or estimated parameters on which it depends. (3) No. The coefficient is condition-specific; the pH used in the derivation must correspond to the intended condition, or the value must be recalculated under an appropriate equation. (4) \(5.00\times3.00=15.0\ \mu\text{L/min/mg}\). The units remain clearance per milligram because the abundance ratio is dimensionless.