Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Testing Extended Clearance Models Against In Vivo Evidence

Hepatic Uptake and Efflux Transporters · Tutorial 431 of 799

Testing Extended Clearance Models Against In Vivo Evidence

Learn to compare an extended clearance model’s pathway-level predictions with in vivo measurements and distinguish parameter mismatch from evidence against its assumptions.

Advanced 11 min read

What You'll Learn

  • Calculate effective hepatic intrinsic clearance from uptake and competing intracellular pathways.
  • Relate effective intrinsic clearance to hepatic extraction under a well-stirred liver assumption.
  • Predict the relative amounts returning to blood, entering bile, and undergoing metabolism.
  • Check whether in vivo biliary recovery and intracellular exposure agree with model predictions.
  • Identify measurement and structural failures that can make apparent disagreement misleading.

Why Test the Pathways, Not Just the Total Clearance?

A model can reproduce plasma concentrations or total hepatic clearance while assigning drug to the wrong routes. For example, an overestimated uptake process might be offset by an overestimated return of drug from hepatocytes to blood. Agreement in one aggregate measure would not establish that the pathway assumptions are sound. Testing an extended clearance model therefore means comparing its predictions with several kinds of in vivo evidence: hepatic extraction, biliary recovery, and unbound intracellular exposure.

The central question is not whether one calculated number exactly matches one measurement. It is whether the model, with stated assumptions and uncertainty, can account for the observed evidence across pathways. A discrepancy can indicate an incorrect parameter, an unreliable measurement, or a structural omission such as an unrepresented elimination route. These possibilities need to be distinguished before changing the model.

Model and Assumptions

Consider a linear hepatic model with uptake from blood into hepatocytes and three competing intracellular fates: return to blood, biliary efflux, and metabolism. Define \(CL_{\mathrm{uptake}}\) as the uptake clearance, \(CL_{\mathrm{return}}\) as the clearance for return to blood, \(CL_{\mathrm{bile}}\) as biliary efflux clearance, and \(CL_{\mathrm{met}}\) as metabolic clearance. These clearances act on the relevant unbound concentrations. Let \(C_{\mathrm{u,b}}\) be unbound drug concentration in blood and \(C_{\mathrm{u,h}}\) unbound intracellular concentration in hepatocytes.

Definition: In the linear extended clearance model, uptake is driven by unbound blood concentration, and intracellular drug is removed by return to blood, biliary efflux, and metabolism. The model assumes that each pathway is proportional to its relevant unbound concentration over the exposure range being considered. Its predictions also depend on consistent concentration, flow, and clearance conventions.

For the basic calculations below, assume hepatic blood is well mixed, the liver is at steady state for the clearance calculation, and the intracellular pathways are linear and time-invariant. Assume the measured clearances apply to the same population and conditions as the in vivo evidence. The well-stirred liver relation further assumes that drug leaving the liver is represented by a single hepatic outflow concentration. These are testable simplifications, not universal properties of hepatic disposition.

Let \(Q_{\mathrm{h}}\) be hepatic blood flow and \(f_{\mathrm{u,b}}\) the unbound fraction in blood. All clearances and concentrations must use compatible blood or plasma conventions; the equations here use blood flow and blood-referenced concentrations. If the underlying data use plasma concentrations, the conversion between blood and plasma must be handled consistently.

From Intracellular Pathways to Hepatic Extraction

At steady state, the rate entering the hepatocyte equals the total rate leaving it. Uptake is \(CL_{\mathrm{uptake}}C_{\mathrm{u,b}}\), while total intracellular loss is \((CL_{\mathrm{return}}+CL_{\mathrm{bile}}+CL_{\mathrm{met}})C_{\mathrm{u,h}}\). Equating these rates gives the predicted intracellular concentration:

$$ C_{\mathrm{u,h}} = \frac{CL_{\mathrm{uptake}}} {CL_{\mathrm{return}}+CL_{\mathrm{bile}}+CL_{\mathrm{met}}} C_{\mathrm{u,b}} $$

Only biliary efflux and metabolism are irreversible removal from the body in this simplified liver model. Return to blood is not counted as hepatic elimination: the returned drug can circulate and reach the liver again. Thus the effective hepatic intrinsic clearance is uptake multiplied by the fraction of intracellular loss that follows either elimination route.

$$ CL_{\mathrm{int,h}} = CL_{\mathrm{uptake}} \frac{CL_{\mathrm{bile}}+CL_{\mathrm{met}}} {CL_{\mathrm{return}}+CL_{\mathrm{bile}}+CL_{\mathrm{met}}} $$

Under the well-stirred liver assumption, hepatic clearance and extraction ratio are:

$$ CL_{\mathrm{h}} = \frac{Q_{\mathrm{h}} f_{\mathrm{u,b}}CL_{\mathrm{int,h}}} {Q_{\mathrm{h}}+f_{\mathrm{u,b}}CL_{\mathrm{int,h}}}, \qquad E_{\mathrm{h}}=\frac{CL_{\mathrm{h}}}{Q_{\mathrm{h}}} $$

This formulation makes the model’s structure visible. Increasing uptake can raise effective intrinsic clearance, but that increase is reduced when a larger share of intracellular drug returns to blood. The extraction ratio is bounded between zero and one under positive finite clearances and flow.

Worked Example: Predicting Hepatic Extraction

Suppose hypothetical parameter estimates are \(Q_{\mathrm{h}}=90\) L/h, \(f_{\mathrm{u,b}}=0.10\), \(CL_{\mathrm{uptake}}=300\) L/h, \(CL_{\mathrm{return}}=120\) L/h, \(CL_{\mathrm{bile}}=20\) L/h, and \(CL_{\mathrm{met}}=60\) L/h. Calculate the model-predicted extraction ratio.

The total intracellular exit clearance is \(120+20+60=200\) L/h. The fraction of intracellular loss that is elimination is \((20+60)/200=0.400\). Therefore:

$$ CL_{\mathrm{int,h}} = 300\ \mathrm{L/h}\times\frac{80\ \mathrm{L/h}}{200\ \mathrm{L/h}} = 120\ \mathrm{L/h} $$

Applying the well-stirred relation:

$$ CL_{\mathrm{h}} = \frac{90\times0.10\times120}{90+0.10\times120} = \frac{1080}{102} = 10.59\ \mathrm{L/h} $$

Thus \(E_{\mathrm{h}}=10.59/90=0.1177\), or about 11.8%. If a compatible in vivo estimate were 0.12, this calculation would be consistent with that aggregate observation, subject to uncertainty. It would not, by itself, confirm the separate uptake, return, bile, or metabolism pathways.

Testing Biliary Recovery and Pathway Fractions

Under the linear competing-pathway assumption, the fraction of intracellular drug leaving by a particular route equals that route’s clearance divided by the sum of all intracellular exit clearances. For a known amount of drug entering hepatocytes, this gives the expected partition among blood return, bile, and metabolism. This comparison requires care: biliary recovery as a fraction of an administered dose is not automatically the same as biliary fraction of drug taken up by the liver. The administered dose may include drug that never entered hepatocytes, and recovery may be incomplete.

Worked Example: Comparing Predicted and Observed Biliary Recovery

Use the hypothetical intracellular clearances from the prior calculation. Suppose a tracer mass-balance experiment estimates that 100 μmol entered hepatocytes over the observation period. Calculate the expected route amounts, then compare with a hypothetical measured biliary recovery of 9 μmol.

The total intracellular exit clearance is 200 L/h. The pathway fractions are \(120/200=0.600\) for return to blood, \(20/200=0.100\) for bile, and \(60/200=0.300\) for metabolism. Therefore, the predicted amounts are:

$$ \begin{aligned} A_{\mathrm{return}} &= 100\times0.600=60\ \mathrm{\mu mol},\\ A_{\mathrm{bile}} &= 100\times0.100=10\ \mathrm{\mu mol},\\ A_{\mathrm{met}} &= 100\times0.300=30\ \mathrm{\mu mol}. \end{aligned} $$

The predicted amounts sum to 100 μmol, as required by the assumed complete accounting. The measured 9 μmol in bile is close to the predicted 10 μmol, but the comparison is meaningful only if the 100 μmol uptake estimate and the 9 μmol biliary recovery refer to compatible time windows and if bile collection captures the relevant drug or metabolites. If 9 μmol instead represented 9% of the administered dose, it could not be compared directly with the predicted 10% of hepatocyte uptake without accounting for the fraction of the dose that entered the liver.

Testing Intracellular Exposure

Intracellular exposure tests a different implication of the model. At steady state, the ratio of unbound intracellular to unbound blood concentration is determined by uptake relative to the sum of intracellular exit clearances. A high predicted intracellular concentration can coexist with modest hepatic extraction if uptake is high but return to blood is also substantial. Therefore, total clearance alone cannot determine intracellular exposure.

Worked Example: Using Intracellular Exposure to Challenge the Model

Suppose the hypothetical clearances above predict exposure at an unbound blood concentration of \(C_{\mathrm{u,b}}=0.20\) μmol/L. Calculate the predicted \(C_{\mathrm{u,h}}\). Then compare it with a hypothetical in vivo estimate of 0.14 μmol/L.

The total intracellular exit clearance is 200 L/h, so:

$$ C_{\mathrm{u,h}} = \frac{300\ \mathrm{L/h}}{200\ \mathrm{L/h}} \times0.20\ \mathrm{\mu mol/L} = 0.300\ \mathrm{\mu mol/L} $$

The predicted intracellular concentration is about 2.14 times the observed estimate, since \(0.300/0.14=2.14\). This is a material disagreement to investigate, not proof by itself that uptake is wrong. The intracellular estimate may have uncertainty, may represent total rather than unbound drug, or may not reflect steady state. If measurement and binding corrections are reliable and the model conditions match, the discrepancy challenges at least one of the assumed pathway clearances or the linear steady-state structure.

How to Judge Agreement and Locate Failure

A strong evaluation checks multiple predictions and their conditions rather than applying a single pass/fail threshold. Compare model predictions with observations on the same scale and at compatible times, and carry uncertainty in measured quantities and parameters into the comparison. If the evidence is imprecise, a numerical difference may be uninformative; if uncertainty is small and predictions repeatedly miss in distinct measures, structural assumptions deserve scrutiny.

1
Align the quantities.
Confirm whether observations refer to blood or plasma, unbound or total concentration, administered dose or hepatocyte uptake, and parent drug or total drug-related material.
2
Calculate pathway predictions.
Use the fitted or independently measured clearances to predict extraction, route fractions, and intracellular exposure. Do not tune a parameter to one observation and then treat that same observation as independent validation.
3
Check mass balance and timing.
Ask whether observed parent drug and metabolites account for the expected material, whether biliary collection is complete, and whether measurements cover the same disposition period.
4
Interpret joint discrepancies.
A mismatch in only one measurement may reflect that measurement’s limitations. Discordance in extraction, biliary recovery, and intracellular exposure together is stronger evidence that the model or its inputs are inadequate.

Several patterns are especially informative. Correct extraction but excessive predicted intracellular concentration may point to an incorrect uptake-to-exit balance. Correct intracellular exposure but too much predicted biliary recovery may implicate the bile-versus-metabolism partition. Extraction and biliary recovery that are both underpredicted could indicate an omitted elimination route, under-scaled uptake, or incompatible blood-flow assumptions. These patterns guide diagnosis but do not uniquely identify a parameter; identifiability must be considered.

Alternatives, Diagnostics, and Limitations

The extended model can be compared with a simpler well-stirred model using a single effective intrinsic clearance. The simpler model may be adequate for predicting hepatic clearance, but it cannot independently predict intracellular exposure or the split between bile and metabolism. Conversely, a compartmental hepatic PBPK model can represent changing concentrations, zonation, and additional compartments, but requires more parameters and more informative data. A detailed model is not automatically better: it is useful only when the evidence can constrain its added structure.

Look for practical warning signs: predicted extraction outside the physically possible range; route fractions that do not sum to one; a mismatch between calculated clearance units or concentration bases; sensitivity to a single poorly measured parameter; and disagreement that disappears only after inconsistent definitions of recovery are used. The steady-state formulas are not appropriate when pathways are saturable over the tested range, when concentrations change too rapidly for the assumed balance, or when transport direction and intracellular binding are misrepresented. Time-course modeling or nonlinear pathway descriptions may then be needed.

A common mistake is to treat agreement in hepatic extraction as proof of the transporter mechanism. Extraction is an aggregate outcome: different uptake and intracellular loss combinations can yield the same effective clearance. The correction is to test additional, mechanistically distinct observations—especially pathway-specific recovery and unbound intracellular exposure—and to ask whether the measurements can actually distinguish competing explanations.

Key takeaway: An extended clearance model is credible only to the extent that its pathway assumptions jointly explain compatible in vivo evidence. Matching total hepatic extraction is necessary for a good clearance prediction, but it is not sufficient to establish the proposed intracellular route structure.

Check Your Understanding

Use the linear, steady-state pathway model defined above. Answers are provided after each question so the calculations can be checked directly.

  1. If \(CL_{\mathrm{uptake}}=200\) L/h, \(CL_{\mathrm{return}}=100\) L/h, \(CL_{\mathrm{bile}}=30\) L/h, and \(CL_{\mathrm{met}}=70\) L/h, what fraction of intracellular exit is elimination? Solution: \((30+70)/(100+30+70)=100/200=0.500\), or 50%.
  2. For the same clearances, what is \(C_{\mathrm{u,h}}/C_{\mathrm{u,b}}\)? Solution: \(200/(100+30+70)=1.00\). The predicted unbound concentrations are equal under the stated steady-state model.
  3. Why can biliary recovery reported as a fraction of administered dose not automatically be compared with the model’s biliary pathway fraction? Solution: The pathway fraction is conditional on drug entering hepatocytes. The administered-dose denominator also includes drug that may not enter hepatocytes, and observed recovery may be incomplete.
  4. A model matches observed extraction but predicts twice the measured unbound intracellular concentration. Does the extraction match establish that the uptake and efflux pathways are correct? Solution: No. Different pathway combinations can produce the same aggregate extraction. The exposure discrepancy challenges the model or measurement assumptions and requires investigation.
  5. Name one condition under which the linear steady-state calculations may fail. Solution: Examples include saturable transport or metabolism over the concentration range, rapidly changing concentrations that do not satisfy the steady-state approximation, or incompatible concentration and clearance conventions.