Deriving the Parallel-Tube Liver Model

Why Model Drug Removal Along Flow Paths?

The liver receives drug in blood flowing through a network of sinusoids. A model of hepatic clearance must describe how drug concentration changes between the hepatic inlet and outlet. The parallel-tube model represents this passage as flow through many parallel paths: blood enters each path at the same concentration, drug is removed along the path, and the effluent leaves with a lower concentration.

The key distinction from the well-stirred liver model is the assumed pattern of mixing. In the parallel-tube model, blood moves through the liver with plug flow: it is not instantaneously mixed throughout the organ while it is being cleared. Instead, concentration declines along each path. The derivation below turns that first-order loss into a compact expression for hepatic extraction and clearance.

Definition: The hepatic extraction ratio, $E_H$, is the fraction of drug entering the liver that is removed in one passage through it: $$ E_H=\frac{C_{\mathrm{in}}-C_{\mathrm{out}}}{C_{\mathrm{in}}}. $$ Here $C_{\mathrm{in}}$ and $C_{\mathrm{out}}$ are hepatic inlet and outlet blood concentrations. Hepatic clearance, $CL_H$, is the hepatic elimination rate divided by $C_{\mathrm{in}}$. At steady flow, $CL_H=Q_H E_H$, where $Q_H$ is hepatic blood flow.

Model Assumptions and Notation

Let $CL_{\mathrm{int}}$ denote the whole-liver intrinsic clearance, expressed on the unbound-drug basis, and let $f_u$ be the unbound fraction in blood. Under the assumptions below, the elimination rate is proportional to total blood concentration, with effective proportionality $f_u CL_{\mathrm{int}}$. Thus, for a small amount of liver with volume $dV$, the local clearance is proportional to $dV$, and its drug-removal rate is proportional to the concentration passing through it.

The model relies on a steady hepatic blood flow, first-order unsaturated removal, and constant $f_u$ and $CL_{\mathrm{int}}$ during passage. It also assumes homogeneous perfusion and metabolic capacity across the parallel paths, so each path has the same ratio of intrinsic clearance to flow. There is no axial mixing between paths and no delay or storage that changes the concentration during transit except for elimination. These assumptions are idealizations; they make the model analytically tractable.

For the derivation, represent the liver as $N$ identical paths. Let total liver volume be $V_H$, so each path has volume $V_H/N$, and let each path carry flow $Q_H/N$. If intrinsic clearance is distributed uniformly, each path has intrinsic clearance $CL_{\mathrm{int}}/N$. The concentration in any path therefore follows the same first-order loss law.

Deriving the Outlet Concentration

Consider a small flowing element of blood moving through one path. Its concentration falls in proportion to its current concentration. Because the local clearance is distributed uniformly through liver volume, the first-order rate constant for the element is $f_u CL_{\mathrm{int}}/V_H$. This quantity has units of inverse time: clearance has units of volume per time, and liver volume has units of volume.

Writing $C(t)$ for the concentration in that element at time $t$, its mass-balance equation is

\frac{dC(t)}{dt}=-\frac{f_u CL_{\mathrm{int}}}{V_H}C(t).

Separate variables and integrate from entry, when $C(0)=C_{\mathrm{in}}$, to exit after transit time $\tau$:

\int_{C_{\mathrm{in}}}^{C(\tau)}\frac{dC}{C} = -\frac{f_u CL_{\mathrm{int}}}{V_H}\int_0^\tau dt, \qquad \ln\!\left(\frac{C(\tau)}{C_{\mathrm{in}}}\right) = -\frac{f_u CL_{\mathrm{int}}\tau}{V_H}.

In the homogeneous parallel-path representation, transit time is path volume divided by path flow. Both are the same fractions of their whole-liver totals, so this is $V_H/Q_H$. The concentration at exit is the common outlet concentration, giving

\tau=\frac{V_H}{Q_H}, \qquad \frac{C_{\mathrm{out}}}{C_{\mathrm{in}}} = \exp\!\left(-\frac{f_u CL_{\mathrm{int}}}{Q_H}\right).

The volume cancels: for a given flow and total intrinsic clearance, the predicted extraction does not depend separately on the liver volume in this formulation. This cancellation follows from the reciprocal effects of volume on the local removal rate and on transit time.

Result: Under the parallel-tube assumptions, the hepatic extraction ratio and clearance are $$ E_H=1-\exp\!\left(-\frac{f_u CL_{\mathrm{int}}}{Q_H}\right), \qquad CL_H=Q_H\left[1-\exp\!\left(-\frac{f_u CL_{\mathrm{int}}}{Q_H}\right)\right]. $$ The exponent is dimensionless because $f_u$ is dimensionless and both $CL_{\mathrm{int}}$ and $Q_H$ have units of volume per time.

To obtain the clearance expression, use the steady-flow mass balance. The amount removed per unit time is $Q_H(C_{\mathrm{in}}-C_{\mathrm{out}})$. Dividing by $C_{\mathrm{in}}$ gives $CL_H=Q_H E_H$. For positive finite $f_u CL_{\mathrm{int}}$ and $Q_H$, the exponential is between zero and one. Consequently, $0<E_H<1$ and $0<CL_H<Q_H$. If $f_u CL_{\mathrm{int}}=0$, both extraction and clearance are zero.

Worked Example: Calculate Extraction and Clearance

Suppose a hypothetical liver has $Q_H=90\ \mathrm{L/h}$, $f_u=0.080$, and $CL_{\mathrm{int}}=450\ \mathrm{L/h}$. Assume steady flow, first-order unsaturated removal, and homogeneous parallel paths with constant $f_u$ and intrinsic clearance. Calculate extraction and hepatic clearance.

First calculate the dimensionless exponent:

\frac{f_u CL_{\mathrm{int}}}{Q_H} = \frac{0.080(450\ \mathrm{L/h})}{90\ \mathrm{L/h}} = \frac{36\ \mathrm{L/h}}{90\ \mathrm{L/h}} = 0.400.

Then the outlet-to-inlet ratio and extraction are

\frac{C_{\mathrm{out}}}{C_{\mathrm{in}}}=e^{-0.400}=0.6703, \qquad E_H=1-0.6703=0.3297.

Thus about $32.97\%$ of drug entering the liver is extracted in one passage under these hypothetical assumptions. Clearance is

CL_H=(90\ \mathrm{L/h})(0.3297)=29.67\ \mathrm{L/h}.

The result is below $Q_H=90\ \mathrm{L/h}$, as required by the model’s physical bound.

Interpreting the Parallel Paths

The parallel paths matter because each path receives the same inlet concentration, but removal occurs during transit rather than after instantaneous mixing throughout the liver. When paths are identical, each has the same outlet-to-inlet ratio. Their combined outflow therefore has that same concentration ratio; summing the flows does not change the extraction ratio.

More generally, if paths do not have identical flows or enzyme distributions, each path can have a different concentration ratio. The total outlet concentration must then be calculated as a flow-weighted mixture of the path outlet concentrations. The single exponential formula above applies directly when the paths share the same ratio of effective intrinsic clearance to flow. It should not be assumed to hold unchanged for arbitrary heterogeneous paths.

Worked Example: Many Identical Tubes Give the Same Total Extraction

Consider a hypothetical liver represented by 100 identical tubes. Total hepatic flow is $80\ \mathrm{L/h}$, total intrinsic clearance is $240\ \mathrm{L/h}$, and $f_u=0.25$. Assume the flow and intrinsic clearance are divided equally among the tubes and the other parallel-tube assumptions hold. Calculate the extraction ratio using one tube and using whole-liver quantities.

Each tube carries $80/100=0.800\ \mathrm{L/h}$ and has intrinsic clearance $240/100=2.40\ \mathrm{L/h}$. Its exponent is

\frac{f_u(CL_{\mathrm{int}}/100)}{Q_H/100} = \frac{0.25(2.40\ \mathrm{L/h})}{0.800\ \mathrm{L/h}} = 0.750.

The tube extraction is $1-e^{-0.750}=1-0.4724=0.5276$. Using whole-liver quantities gives the same exponent:

\frac{f_u CL_{\mathrm{int}}}{Q_H} = \frac{0.25(240\ \mathrm{L/h})}{80\ \mathrm{L/h}} = 0.750.

Therefore the whole-liver extraction is also $E_H=0.5276$, and clearance is

CL_H=(80\ \mathrm{L/h})(0.5276)=42.21\ \mathrm{L/h}.

The number of tubes cancels because each tube receives the same proportional share of flow and intrinsic clearance. If those shares were not proportional, the path-specific outlet concentrations could differ.

Comparing the Model With the Well-Stirred Liver

The well-stirred model assumes rapid mixing within the liver, so the organ is represented by a single concentration rather than a concentration that declines along a flow path. Its standard clearance expression is $CL_H=Q_H f_u CL_{\mathrm{int}}/(Q_H+f_u CL_{\mathrm{int}})$. The parallel-tube expression instead results from integrating first-order removal during plug-flow transit. The two equations arise from different assumptions about how blood and drug are distributed during passage, not from different definitions of hepatic clearance.

Feature	Parallel-tube model	Well-stirred model
Mixing assumption	Plug flow along parallel paths; no axial mixing	Rapid mixing within the liver
Concentration during passage	Declines along each path as drug is removed	Represented by a single well-mixed concentration
Clearance expression	$Q_H[1-\exp(-f_u CL_{\mathrm{int}}/Q_H)]$	$Q_H f_u CL_{\mathrm{int}}/(Q_H+f_u CL_{\mathrm{int}})$
Key structural requirement	Homogeneous paths, or a common effective clearance-to-flow ratio	Rapid mixing and the well-stirred organ balance

Choosing between the models is a structural modeling decision. A fitted clearance value alone does not establish that the liver is literally well stirred or composed of identical tubes. Model choice should be guided by the biological question and by whether the assumptions are reasonable for the intended application.

Diagnostics, Pitfalls, and Limitations

Check units. The ratio $f_u CL_{\mathrm{int}}/Q_H$ must be dimensionless. A unit mismatch between intrinsic clearance and hepatic flow makes the exponential invalid.
Check physical bounds. For nonnegative $f_u CL_{\mathrm{int}}$ and positive $Q_H$, the model requires $0\leq E_H<1$ and $0\leq CL_H<Q_H$. A prediction outside these bounds signals an arithmetic, coding, or parameterization problem.
Do not confuse the two model equations. Using the well-stirred clearance formula while claiming a parallel-tube derivation applies the wrong mixing assumption. The correction is to use the exponential expression derived from first-order loss over transit time.
Do not treat heterogeneous paths as identical without justification. Large differences in path flow or metabolic capacity can produce different outlet concentrations. The simple equation may then be inadequate; a path-specific calculation or a richer liver representation may be needed.
Recognize first-order limits. Saturable metabolism, changing protein binding, or time-varying flow violates assumptions used in the derivation. In such settings, the local loss rate may not remain proportional to concentration with a constant coefficient.

The parallel-tube model is useful because it connects a mechanistic picture of removal during transit to an explicit extraction ratio and clearance. Its simplicity is also its limitation: it summarizes a spatially distributed organ with homogeneous paths and a single effective intrinsic clearance. The formula is appropriate only to the extent that those assumptions are acceptable.

Check Your Understanding

Use the parallel-tube assumptions for each question. State units where relevant.

For $Q_H=75\ \mathrm{L/h}$, $f_u=0.10$, and $CL_{\mathrm{int}}=300\ \mathrm{L/h}$, calculate the exponent, extraction ratio, and hepatic clearance.
If $E_H=0.40$ and $Q_H=100\ \mathrm{L/h}$, what is $CL_H$?
Why does liver volume cancel from the final extraction expression in this derivation?
What assumption allows the concentration ratio from a single tube to represent the whole liver?
A calculation gives $CL_H=112\ \mathrm{L/h}$ when $Q_H=100\ \mathrm{L/h}$. Is this compatible with the parallel-tube model for nonnegative parameters? Explain.

Solutions. (1) The exponent is $0.10(300)/75=0.400$. Thus $E_H=1-e^{-0.400}=0.3297$, and $CL_H=75(0.3297)=24.73\ \mathrm{L/h}$. (2) $CL_H=Q_HE_H=(100\ \mathrm{L/h})(0.40)=40.0\ \mathrm{L/h}$. (3) The local first-order rate constant is proportional to $1/V_H$, while transit time is $V_H/Q_H$; multiplying them leaves $f_u CL_{\mathrm{int}}/Q_H$. (4) The paths must have the same effective intrinsic-clearance-to-flow ratio, as in the homogeneous-path assumption. (5) No. The model gives $CL_H=Q_H E_H$ with $E_H<1$ for finite nonnegative parameters, so clearance must be less than $100\ \mathrm{L/h}$; check the inputs, units, and calculation.