From Clinical Observations to PBPK Inputs
A hepatic PBPK model needs quantities such as functional liver mass, hepatic blood flow, enzyme capacity, plasma protein binding, and biliary transport capacity. Clinical records more often contain aminotransferases, bilirubin, albumin, international normalized ratio (INR), imaging findings, and a severity score. The central task is to decide what each observation can support, then express the resulting model input with a clear definition, unit, reference, and uncertainty.
A laboratory result is not automatically a measurement of a PBPK parameter. Alanine aminotransferase (ALT) and aspartate aminotransferase (AST) primarily signal hepatocellular injury; bilirubin and alkaline phosphatase can inform assessment of cholestasis or excretory dysfunction; albumin and INR can inform assessment of synthetic function but have important nonhepatic influences. None of these is, by itself, a direct measurement of hepatic enzyme abundance, hepatic flow, or functional liver mass. The earlier material on Child-Pugh and MELD scores provides clinical severity context; here the emphasis is on translating clinical evidence into model covariates.
Define the Covariate Before Estimating It
For each proposed covariate, state what parameter it modifies, its reference condition, and its direction. A ratio is often convenient: it is dimensionless, has a defined reference of 1, and makes the intended comparison explicit. For a physiological quantity \(\theta\), define the impairment scaling factor as
The numerator and denominator must refer to comparable quantities and compatible units. For example, if \(\theta\) is hepatic blood flow in litres per hour, both values must represent the same flow definition. A ratio of 0.70 means the specified quantity is estimated as 70% of its reference value; it does not mean that every hepatic process is reduced by 30%.
The clinical-to-covariate mapping should distinguish at least four evidence types:
- Direct or near-direct physiological measurements: for example, imaging-based liver volume or a suitably measured flow quantity. These may still require assumptions about functional tissue, measurement error, and comparability to the model parameter.
- Functional readouts: for example, a probe-substrate phenotype or a clearance measurement. These can reflect several processes together, so they should not be assigned to a single mechanism without justification.
- Clinical biomarkers: laboratory tests that reflect injury, excretion, or synthesis indirectly and are influenced by factors beyond liver physiology.
- Clinical scores and categories: summaries of multiple observations that can stratify patients, but generally do not uniquely identify individual physiological parameters.
This distinction matters because the same score can arise from different combinations of findings. If one patient has low albumin and another has ascites, the same aggregate severity category does not establish that their enzyme capacity, flow, or functional mass is equal. A score may be useful for defining a population stratum or a conditional distribution of covariates, but it is not itself a measurement of each covariate.
Worked Example: Reading Routine Tests Without Overassigning Enzyme Capacity
Suppose a hypothetical patient has ALT 156 U/L, AST 120 U/L, alkaline phosphatase 108 U/L, total bilirubin 0.8 mg/dL, albumin 4.0 g/dL, and INR 1.0. The laboratory upper limits of normal are 40 U/L for ALT and AST and 120 U/L for alkaline phosphatase. What can be assigned to a hepatic PBPK model from these results alone?
First express the enzyme results relative to their stated reference limits:
Alkaline phosphatase is \(108/120=0.900\) times its upper limit. These calculations describe laboratory results; they are not scaling factors for enzyme capacity. The elevated ALT and AST support an injury signal, while the other listed results do not show an elevation relative to their stated references. These results alone do not quantify hepatic blood flow, functional liver mass, or the capacity of a particular drug-metabolizing enzyme.
A defensible model record would preserve the values, units, reference limits, and clinical interpretation, but leave unsupported mechanistic scaling factors unassigned. If an enzyme-capacity factor is required, it needs an additional evidence source or an explicitly uncertain assumption. Assigning \(S_{\mathrm{enzyme}}=0.39\) merely because ALT is 3.90 times its limit would reverse the meaning of the ratio and invent a relationship that these measurements do not establish.
Estimating Covariates From Available Physiology
When a measurement corresponds closely to a model quantity, estimate a ratio only after checking how the measurement was obtained and what the model parameter represents. Imaging-based liver volume, for example, can inform an anatomical mass estimate. It does not automatically measure the fraction of that tissue that is functional, nor does it directly determine the abundance of a specific enzyme. If functional fraction is estimated separately, the assumptions and sources should be recorded rather than hidden inside a single number.
Worked Example: Calculating a Liver-Mass Scaling Factor
In a hypothetical assessment, imaging gives a patient liver volume of 1.20 L. A suitable reference liver volume for the intended comparison is 1.50 L. For illustration, an independently justified estimate assigns 0.80 of the patient's measured volume to functional tissue. Assume the reference value represents fully functional tissue. Calculate the functional-mass scaling factor and show how uncertainty in functional fraction changes it.
Under these assumptions, estimated functional volume is \(1.20\ \mathrm{L}\times0.80=0.960\ \mathrm{L}\). The scaling factor is therefore
The model input is 0.640 relative to this defined reference. The units cancel, and the value does not imply that enzyme capacity or hepatic flow is also 64% of reference. If the functional fraction is uncertain between 0.70 and 0.90, the corresponding range is
Thus, the estimated scaling factor ranges from 0.560 to 0.720 under the stated functional-fraction bounds. These bounds are conditional on the chosen reference and assumptions; they are not a confidence interval unless a statistical procedure and sampling model justify that interpretation.
Other physiological inputs need their own mappings. A flow measurement can inform a hepatic-flow ratio if it measures the relevant flow and is comparable with the reference. Evidence of portosystemic shunting should be represented separately when the model distinguishes shunted from liver-perfusing flow; it should not be silently folded into a single flow reduction. The earlier material on hepatic blood flow and shunting establishes why those are distinct features. Similarly, a probe clearance can reflect flow, shunting, binding, and intrinsic metabolic capacity together, so using it as a pure enzyme-capacity measurement requires additional evidence.
For binding, a measured albumin concentration may support a disease description or a binding model, but concentration alone does not establish the unbound fraction for every drug. Drug-specific binding information and relevant competing factors are needed. For biliary disposition, bilirubin, alkaline phosphatase, and related clinical findings can inform the assessment of cholestasis, but they do not directly quantify a particular transporter’s capacity. The earlier material on hepatic enzyme capacity, binding, and cholestasis provides the mechanistic distinctions needed to avoid treating these as interchangeable effects.
Missing Measurements and Uncertainty
A blank laboratory result is not a zero, a normal value, or evidence that the corresponding process is absent. First determine whether the measurement was not ordered, was unavailable, was technically unreliable, or was omitted for another reason. Missingness can be associated with observed patient characteristics or with the unobserved value itself. Consequently, the reason and pattern of missingness affect how defensible an imputation is.
If a covariate is missing, reasonable options include leaving it unestimated and reporting that limitation, imputing from a model that conditions on relevant observed information, or evaluating a set of plausible scenarios. A single substituted value conceals uncertainty. When imputation is justified, the imputation model should preserve relevant relationships among variables and should be compatible with the patient population and measurement process. A clinical score can be an input to such a model, but it cannot replace a missing physiological measurement by definition.
where \(M\) denotes additional information used in the imputation or conditional model. This is the law of total variance: average the variance conditional on both \(Z\) and \(M\), then add the variance of the conditional mean across plausible \(M\). The first term represents residual uncertainty even after the additional information is known; the second represents uncertainty because that information is not fully known. The practical implication is that a model input derived from incomplete records should generally carry uncertainty forward rather than be treated as exact.
Worked Example: Handling a Missing Hepatic-Flow Measurement
A hypothetical patient has no hepatic-flow measurement. For illustration only, suppose an appropriate external reference for patients with comparable observed clinical features supports three equally weighted plausible flow values: 54, 65, and 77 L/h. The reference flow is 90 L/h. Calculate the corresponding flow ratios and summarize the uncertainty without claiming that any one ratio is observed.
For each plausible value, divide by the same reference flow:
Under the stated equal weights, the mean ratio is
The population standard deviation of these three illustrative equally weighted possibilities is
The ratio is therefore represented here by three plausible values, with an illustrative mean of 0.726 and standard deviation of about 0.104. The patient did not have a measured flow of \(0.726\times90=65.3\) L/h; that is only the value corresponding to the mean ratio. In an actual analysis, the conditional distribution and weights must be supported by suitable data, and the uncertainty should be propagated through the PBPK analysis. If no defensible conditional reference exists, report the flow as missing and examine justified scenarios rather than manufacturing a point estimate.
Check the Mapping and Its Failure Modes
A practical covariate record should make the translation auditable. For each input, document the target model parameter, observed source, units and timing, reference definition, transformation, assumptions, and uncertainty. This also helps identify whether two observations are independent evidence or different manifestations of the same underlying finding. Counting the same information twice—for example, using a score that already includes a laboratory result and then applying that result again as an independent severity adjustment—can exaggerate the implied disease effect.
Diagnostics should focus on both the data and the resulting covariate set. Check that units, reference ranges, dates, and patient identifiers are correct; compare values with clinically plausible ranges; and investigate contradictions rather than averaging them away. A large aminotransferase elevation paired with stable synthetic markers may be consistent with a different clinical pattern than low albumin and prolonged INR, but neither pattern uniquely determines the PBPK parameters. If assigned covariates imply impossible values, reverse the expected direction, or vary implausibly with small changes in a laboratory result, revisit the mapping and the model’s parameter definitions.
The main alternatives are not competing formulas so much as different evidence strategies. Direct physiological measurement is preferable when it measures the intended quantity reliably. A validated functional readout can be useful when its relationship to the target parameter is established. A severity category can support stratification or conditional distributions when direct measurements are unavailable. Expert assumptions or scenario ranges may be necessary when evidence is sparse, but their status should remain explicit. These strategies are not interchangeable: a biomarker that predicts a clinical category may still be a poor estimator of a mechanistic parameter.
Specify the PBPK quantity, its units, and the reference condition before looking for a proxy.
Separate direct measurement, functional readout, clinical biomarker, and severity score; record timing and limitations.
Calculate a ratio only when numerator and denominator represent comparable quantities. Otherwise use a validated conditional model, a clearly labelled scenario, or leave the covariate unassigned.
Represent plausible values or distributions, examine contradictions and dependence among inputs, and document what the data do not identify.
The central limitation is identifiability: routine clinical data often support a description of disease severity more strongly than a unique numerical estimate of each mechanism. Several different combinations of functional mass, flow, binding, and enzyme capacity can produce similar clinical findings. Additional measurements, prior physiological knowledge, or explicit assumptions are needed to distinguish those combinations. A covariate estimate should therefore communicate both what the data support and what remains uncertain.
- Common mistake: Treating an elevated ALT as a proportional decrease in every hepatic enzyme capacity. Correction: ALT is an injury biomarker; a quantitative enzyme-capacity factor requires evidence linking the available data to that particular capacity.
- Common mistake: Setting a missing flow or albumin value to zero, or silently replacing it with a typical value. Correction: Record it as missing, investigate why it is missing, and use justified conditional imputation or sensitivity scenarios that preserve uncertainty.
- Common mistake: Applying both a severity-score adjustment and all of its component findings as if they were independent. Correction: identify overlap and avoid double counting the same clinical information.
Check Your Understanding
For each item, give a concise answer and compare it with the solutions below.
- A patient’s ALT is four times the upper limit of normal. Does this establish an enzyme-capacity scaling factor of 0.25? Why or why not?
- A measured liver volume is 1.10 L, the reference volume is 1.40 L, and an independently estimated functional fraction is 0.75. Calculate the functional-mass scaling factor.
- A hepatic-flow result is absent from the record. What is wrong with entering a flow ratio of zero?
- Why might two patients in the same clinical severity category need different hepatic PBPK covariates?
Solutions. (1) No. ALT indicates an injury signal but does not quantify the capacity of a particular drug-metabolizing enzyme; the proposed value has no justified mapping. (2) \(1.10\times0.75/1.40=0.589\) to three significant figures, relative to the specified reference. (3) Missing means unobserved, not zero flow; zero would assert an extreme physiological state and could create a false model input. Use justified conditional information or scenarios, or leave the covariate unassigned. (4) A category can summarize different combinations of clinical findings, and those findings do not uniquely determine functional mass, flow, binding, or enzyme capacity.