From Hepatic Mass Balance to Clearance
The well-stirred liver model converts three physiological and drug-specific quantities—hepatic blood flow, the unbound fraction in blood, and intrinsic clearance—into a prediction of hepatic clearance. Its central idea is a mass balance: the drug removed from blood as it passes through the liver must equal the amount eliminated by the liver. The resulting clearance is not generally equal to intrinsic clearance, because delivery by blood flow and binding in blood also constrain elimination.
The derivation uses the well-stirred assumption: blood entering the liver is treated as mixing with the hepatic compartment, and the concentration leaving the liver is the concentration to which the eliminating processes are exposed. This assumption and its physiological implications are the focus of the earlier material on assumptions behind the well-stirred liver. Here, the aim is to derive the equations and understand what each term represents.
Assumptions and the Elimination Rate
The standard derivation assumes steady flow and a steady mass balance across the liver; no net accumulation of drug within the liver during the balance interval; and the well-stirred mixing condition described above. It also assumes that elimination is linear in the unbound concentration over the concentrations considered. In this linear range, the elimination rate is intrinsic clearance multiplied by unbound concentration. The intrinsic clearance is treated as an aggregate hepatic property, not as the clearance of the whole body.
In the well-stirred model, the concentration exposed to elimination is \(C_{\mathrm{out}}\). The corresponding unbound concentration is \(f_{u,B}C_{\mathrm{out}}\). Thus the hepatic elimination rate can be written as:
This expression has units of amount per time: \(CL_{\mathrm{int}}\) has units of volume per time and \(f_{u,B}C_{\mathrm{out}}\) has units of amount per volume. It presumes that the unbound fraction is the relevant fraction for elimination and that \(CL_{\mathrm{int}}\) is defined with respect to unbound drug.
Deriving Extraction and Clearance
The mass entering the liver per unit time is \(Q_H C_{\mathrm{in}}\); the mass leaving per unit time is \(Q_H C_{\mathrm{out}}\). Their difference is the net amount removed from blood per unit time. At steady state, this difference equals the elimination rate. Combining the flow mass balance with the intrinsic-clearance expression gives:
Rearrange to collect the terms involving \(C_{\mathrm{out}}\). This step makes explicit how the outlet concentration depends on both delivery and elimination:
By definition, the extraction ratio is one minus the fraction of inlet concentration remaining in the outlet blood. Substituting the concentration ratio yields the well-stirred extraction equation:
Clearance is the elimination rate divided by the inlet concentration. The flow mass balance therefore gives \(CL_H=Q_H(C_{\mathrm{in}}-C_{\mathrm{out}})/C_{\mathrm{in}}\), or \(CL_H=Q_HE_H\). Substitution produces the well-stirred hepatic clearance equation:
The same result follows by dividing the elimination-rate expression by \(C_{\mathrm{in}}\) and substituting the derived outlet-to-inlet ratio. The units of \(CL_H\) are volume per time. It is a whole-organ clearance referenced to inlet blood concentration; it is not the same quantity as \(CL_{\mathrm{int}}\), which describes elimination capacity per unit unbound concentration.
Worked Calculations
Worked Example: Calculating Clearance and Extraction
Hypothetical data. Suppose \(Q_H=90\) L/h, \(f_{u,B}=0.10\), and \(CL_{\mathrm{int}}=600\) L/h. Assume that the linear, well-stirred model applies. First calculate the unbound intrinsic-clearance term: \(f_{u,B}CL_{\mathrm{int}}=0.10(600)=60\) L/h. Then:
Therefore \(CL_H=Q_HE_H=90(0.400)=36.0\) L/h. Equivalently, \(CL_H=90(60)/(90+60)=36.0\) L/h. The predicted extraction is 40.0% per pass, and the predicted hepatic clearance is 36.0 L/h. It is below the hepatic blood flow of 90 L/h, as required by the model.
Worked Example: Checking the Mass Balance With Concentrations
Hypothetical data. Use \(Q_H=90\) L/h, \(f_{u,B}=0.10\), and \(CL_{\mathrm{int}}=600\) L/h, and suppose the inlet concentration is \(C_{\mathrm{in}}=2.00\) mg/L. From the concentration-ratio equation, \(C_{\mathrm{out}}/C_{\mathrm{in}}=90/(90+60)=0.600\). Thus \(C_{\mathrm{out}}=0.600(2.00)=1.20\) mg/L.
The mass removed from blood per hour is:
The elimination-rate expression gives the same answer:
This agreement is the mass-balance check. The extraction ratio is \((2.00-1.20)/2.00=0.400\), consistent with the previous calculation. The clearance is the elimination rate divided by inlet concentration: \(72.0\ \mathrm{mg/h}/(2.00\ \mathrm{mg/L})=36.0\) L/h.
Worked Example: A Small Unbound Intrinsic-Clearance Term
Hypothetical data. Let \(Q_H=90\) L/h, \(f_{u,B}=0.050\), and \(CL_{\mathrm{int}}=18\) L/h. The product governing the equations is \(f_{u,B}CL_{\mathrm{int}}=0.050(18)=0.900\) L/h. Therefore:
The predicted extraction is about 0.9901%, and clearance is about 0.8911 L/h. As a check, \(Q_H f_{u,B}CL_{\mathrm{int}}/(Q_H+f_{u,B}CL_{\mathrm{int}})=90(0.900)/90.9=0.8911\) L/h. The small extraction arises because the unbound intrinsic-clearance term is small relative to hepatic blood flow.
Interpreting the Equation and Checking a Calculation
The equation has a useful structure: clearance is the product of flow and the fraction extracted, while the fraction extracted depends on the competition between \(Q_H\) and \(f_{u,B}CL_{\mathrm{int}}\). A reliable calculation therefore begins by checking that both terms in the denominator have the same units. \(Q_H\) and \(CL_{\mathrm{int}}\) must both be expressed as volume per time, and \(f_{u,B}\) is dimensionless.
Two limiting checks also help catch algebra or coding errors. When \(f_{u,B}CL_{\mathrm{int}}\) is much smaller than \(Q_H\), the denominator is close to \(Q_H\), so \(CL_H\) is close to \(f_{u,B}CL_{\mathrm{int}}\). When \(f_{u,B}CL_{\mathrm{int}}\) is much larger than \(Q_H\), extraction approaches one and clearance approaches \(Q_H\). These are limiting interpretations of this equation, not permission to substitute an approximation when the full formula is needed.
For positive inputs, the denominator exceeds \(f_{u,B}CL_{\mathrm{int}}\), so \(E_H<1\). Also, \(CL_H=Q_HE_H\), which means \(CL_H<Q_H\). A computed clearance greater than hepatic blood flow, or an extraction ratio greater than one, signals a likely error in units, parameter definitions, or implementation—or a mismatch between the data and the model assumptions.
Confirm that flow and intrinsic clearance use compatible volume-per-time units, that \(f_{u,B}\) is a fraction between zero and one, and that concentrations are on a consistent blood basis.
Compute \(f_{u,B}CL_{\mathrm{int}}\). This product, not \(CL_{\mathrm{int}}\) alone, enters the well-stirred equations.
Use the extraction equation, then multiply \(E_H\) by \(Q_H\) to obtain \(CL_H\). Alternatively, substitute directly into the clearance equation.
Verify that \(0\leq E_H<1\), \(0\leq CL_H<Q_H\), and that any calculated outlet concentration is no greater than the inlet concentration under the model.
Common Pitfalls and Model Limits
- Using intrinsic clearance as hepatic clearance. The quantities have different definitions. The model predicts \(CL_H=Q_H f_{u,B}CL_{\mathrm{int}}/(Q_H+f_{u,B}CL_{\mathrm{int}})\), not \(CL_H=CL_{\mathrm{int}}\).
- Omitting the unbound fraction. If intrinsic clearance is referenced to unbound concentration, the elimination-rate term is \(f_{u,B}CL_{\mathrm{int}}C_{\mathrm{out}}\). Omitting \(f_{u,B}\) changes the model and can greatly alter a prediction.
- Mixing units or concentration bases. A clearance in mL/min cannot be combined directly with a flow in L/h without conversion. Likewise, inlet, outlet, and binding quantities must refer to a consistent blood basis.
- Assuming linear elimination outside its range. If elimination saturates or other concentration-dependent processes matter, a constant \(CL_{\mathrm{int}}\) may not describe the elimination rate. The derivation then does not establish the same equations for that situation.
- Reading a model prediction as a measured fact. Calculated \(CL_H\) depends on the model assumptions and parameter inputs. If concentration profiles or independent flow information conflict with its predictions, investigate the assumptions and measurements rather than forcing the calculation to fit.
The well-stirred result is one way to translate inlet–outlet mass balance and an intrinsic elimination process into an organ-level clearance. A direct inlet–outlet calculation can estimate extraction from measured concentrations and flow without first estimating intrinsic clearance. The well-stirred model instead predicts that extraction from flow, binding, and intrinsic clearance. These approaches answer related but different questions: one summarizes observed removal, while the other supplies a mechanistic relationship under specified assumptions. The model does not describe every possible pattern of within-liver concentration or transport; those distinctions require other liver models.
Check Your Understanding
Answer each question, then compare your reasoning with the short solutions.
- If \(Q_H=100\) L/h and \(f_{u,B}CL_{\mathrm{int}}=25\) L/h, what are \(E_H\) and \(CL_H\)?
- For \(Q_H=80\) L/h, \(f_{u,B}=0.20\), and \(CL_{\mathrm{int}}=200\) L/h, calculate the unbound intrinsic-clearance term, extraction ratio, and hepatic clearance.
- Why does the model predict clearance below hepatic blood flow for positive finite parameters?
- A calculation gives \(E_H=1.15\) and \(CL_H=115\) L/h when \(Q_H=100\) L/h. Name a model-based plausibility problem and one item to check.
- What mass-balance equality connects hepatic flow, inlet and outlet concentrations, and intrinsic clearance?
Solutions. (1) \(E_H=25/(100+25)=0.200\), so \(CL_H=100(0.200)=20.0\) L/h. (2) \(f_{u,B}CL_{\mathrm{int}}=0.20(200)=40.0\) L/h; \(E_H=40/(80+40)=0.3333\); \(CL_H=80(0.3333)=26.67\) L/h. (3) Since \(E_H=f_{u,B}CL_{\mathrm{int}}/(Q_H+f_{u,B}CL_{\mathrm{int}})<1\), \(CL_H=Q_HE_H<Q_H\). (4) Both values exceed the model’s physical bounds. Check units and parameter definitions, and check whether the well-stirred assumptions are appropriate. (5) \(Q_H(C_{\mathrm{in}}-C_{\mathrm{out}})=f_{u,B}CL_{\mathrm{int}}C_{\mathrm{out}}\).