Recovering Intrinsic Clearance From Hepatic Clearance

Why Invert Hepatic Clearance?

A hepatic clearance measurement summarizes the net removal of drug from blood passing through the liver. Intrinsic clearance, by contrast, represents the liver’s metabolic capacity before the model accounts for blood flow and the fraction of drug available to the eliminating enzymes. Recovering intrinsic clearance from observed hepatic clearance is therefore a model-based calculation: the same measured clearance can imply different intrinsic clearances under different liver models.

This tutorial gives the inverse equations for the well-stirred and parallel-tube models. The inversions require hepatic blood flow and the unbound fraction in blood as inputs, as well as assumptions that make the forward model applicable. All numerical values below are hypothetical.

Inverting the Well-Stirred Model

The well-stirred model treats the liver as a single, uniformly mixed compartment. At steady state, the rate at which drug enters in hepatic blood minus the rate leaving must equal the rate eliminated. With linear elimination and unbound drug available for metabolism, the model is

Here \(C_{\mathrm{in}}\) and \(C_{\mathrm{out}}\) are the incoming and outgoing blood concentrations. Clearance is defined by the rate removed divided by the incoming concentration, so \(CL_H=Q_H(C_{\mathrm{in}}-C_{\mathrm{out}})/C_{\mathrm{in}}\). Solving the steady-state balance gives the forward equation and its inverse:

The inverse requires \(f_{u,B}>0\) and \(0\leq CL_H<Q_H\). For a positive clearance below flow, it returns a finite, positive intrinsic clearance. As \(CL_H\) approaches \(Q_H\), the denominator \(Q_H-CL_H\) approaches zero and the inferred intrinsic clearance grows without bound. A value \(CL_H=Q_H\) is a limiting prediction, not a finite intrinsic-clearance estimate.

Inverting the Parallel-Tube Model

The parallel-tube model represents the liver as many parallel channels with drug concentration declining along each channel. It assumes plug-like flow through the channels rather than complete mixing throughout the organ. With uniform enzyme capacity and flow along the channels, the model predicts an exponential decline in concentration from inlet to outlet:

Because \(E_H=1-C_{\mathrm{out}}/C_{\mathrm{in}}\) and \(CL_H=Q_HE_H\), the forward and inverse equations are

For a finite, nonnegative intrinsic clearance, the model predicts \(0\leq CL_H<Q_H\). Thus, the inverse requires \(f_{u,B}>0\) and \(0\leq CL_H<Q_H\). As clearance approaches flow from below, the logarithm’s argument approaches zero and the inferred intrinsic clearance diverges. A clearance at or above flow cannot yield a finite nonnegative estimate under this model.

What the Inversions Assume

An algebraic inverse is only as useful as the model and inputs behind it. Both equations assume hepatic clearance can be represented by a steady-state, linear relationship, that the hepatic blood flow and unbound fraction correspond to the same physiological and experimental condition as the clearance, and that the clearance is expressed on a blood basis. The intrinsic clearance is an effective model parameter, not a direct measurement of enzyme abundance.

• Well-stirred model: The liver is treated as uniformly mixed, so the drug concentration available to the eliminating process is represented by the mixed hepatic concentration. The model also assumes a consistent effective intrinsic clearance and no unrepresented process that changes the input-output balance.
• Parallel-tube model: Drug passes through parallel channels with concentration declining along their length. The model assumes an appropriate uniform distribution of flow and intrinsic capacity, with no axial back-mixing between positions along a channel.
• Both models: The equations use \(f_{u,B}CL_{\mathrm{int,H}}\) as the unbound-capacity term. A plasma clearance cannot be inserted without first expressing clearance and related inputs on a compatible blood basis. Likewise, a value of \(f_u\) measured under different conditions may not represent \(f_{u,B}\) for the calculation.

These assumptions are not verified just because an inverse returns a numerical answer. A result should be checked against known physiology and the measurement context. If the observed \(CL_H\) is negative or at least as large as \(Q_H\), the stated models do not provide a finite, nonnegative inverse. Such a result calls for checking measurement uncertainty, units, blood-versus-plasma scaling, and model applicability rather than forcing the calculation.

Comparing the Estimates and Checking Their Behavior

The models have different forward relationships between capacity and clearance, so their inverse estimates generally differ even when \(Q_H\), \(f_{u,B}\), and \(CL_H\) are identical. At low extraction, both models have a similar first-order limit: when \(f_{u,B}CL_{\mathrm{int,H}}/Q_H\) is small, clearance is approximately \(f_{u,B}CL_{\mathrm{int,H}}\). At greater extraction, the predicted relationship diverges and model choice matters more.

For a concrete comparison, use the hypothetical inputs from the first example: \(Q_H=90\) L/h, \(f_{u,B}=0.10\), and \(CL_H=30\) L/h. The well-stirred estimate was 450 L/h. Under the parallel-tube model,

Both estimates reproduce the same clearance under their own equations, but they are not interchangeable. Neither estimate is a model-free recovery of a uniquely measured biological quantity; each is conditional on its model assumptions.

Diagnostics, Limitations, and a Common Pitfall

Before reporting an inverted estimate, check the calculation as well as the assumptions. Confirm that \(CL_H\) and \(Q_H\) use the same units and blood basis, that \(0<f_{u,B}\leq1\), and that the measured clearance lies in the model’s admissible range. Then substitute the inferred \(CL_{\mathrm{int,H}}\) into the corresponding forward equation. The reconstructed \(CL_H\) should agree with the input up to rounding. A disagreement usually indicates an algebra, unit, or transcription error.

An estimate may be mathematically valid but poorly constrained by the inputs. In particular, when \(CL_H\) is near \(Q_H\), the inverse depends strongly on the difference \(Q_H-CL_H\) in the well-stirred model and on \(1-CL_H/Q_H\) inside the logarithm in the parallel-tube model. Measurement errors in clearance or flow can therefore have a large effect on the result. An exact inverse does not remove uncertainty in measured inputs or uncertainty about the model.

A common mistake is to calculate \(CL_{\mathrm{int,H}}=CL_H/f_{u,B}\) for every extraction level. That expression is only a low-extraction approximation, not the general inverse: it omits the model’s flow-dependent correction. For example, in the first worked example it gives \(30/0.10=300\) L/h, rather than the well-stirred estimate of 450 L/h or the parallel-tube estimate of 364.8 L/h. Use the full inverse for the selected model, and label the resulting estimate with that model and its assumptions.

The central limitation is that inversion does not decide which liver model is appropriate. The well-stirred and parallel-tube equations are alternatives with distinct structural assumptions; agreement with one observed clearance alone does not establish that its model is correct. The practical result is therefore best reported as a model-conditional estimate, together with the blood flow, unbound fraction, clearance basis, and model used.