Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Integrating Age-Dependent Hepatic Intrinsic Clearance

Pediatric Ontogeny and Scaling · Tutorial 587 of 799

Integrating Age-Dependent Hepatic Intrinsic Clearance

Learn to combine pediatric enzyme and liver physiology inputs into hepatic intrinsic clearance, then account for blood binding and hepatic blood flow to estimate hepatic clearance.

Advanced 12 min read

What You'll Learn

  • Distinguish intrinsic clearance per unit liver mass from total hepatic intrinsic clearance.
  • Combine enzyme abundance and catalytic activity across parallel metabolic pathways.
  • Convert plasma unbound fraction to an unbound fraction in blood using the blood-to-plasma ratio.
  • Apply the well-stirred liver model to estimate age-specific hepatic clearance.
  • Diagnose double counting, mismatched units, and inappropriate clearance approximations.

From Ontogeny Inputs to Hepatic Clearance

Age-specific enzyme abundance is not itself a clearance. To predict hepatic drug removal, it must be combined with the activity of the enzyme, the amount of liver tissue, binding in blood, and hepatic blood flow. These quantities affect different parts of the calculation: abundance and catalytic activity determine intrinsic clearance, while binding and blood flow determine how that intrinsic capacity translates into hepatic clearance.

The earlier material on age-dependent hepatic enzyme abundance and capacity, developmental enzyme profiles, and plasma binding supplies relevant age-specific inputs. Here the aim is to integrate those inputs consistently. Maturation values should be evaluated at the age being modeled, and liver mass and hepatic blood flow should also correspond to that age.

Definition: Hepatic intrinsic clearance, \(CL_{int,u}\), is the liver’s intrinsic capacity to eliminate unbound drug, before accounting for hepatic blood-flow limitation. In this tutorial, \(CL_{int,u}\) is the total clearance capacity of the whole liver, with units of volume per time. A per-unit-liver-mass intrinsic clearance has units of volume per time per mass of liver and must be multiplied by liver mass to obtain the whole-liver value.

For parallel metabolic pathways, the whole-liver intrinsic clearances are added. For enzyme \(j\), let \(c_{int,u,j,ref}\) be its reference intrinsic clearance per unit liver mass. Let \(R_{E,j}(a)\) be the ratio of enzyme abundance at age \(a\) to reference abundance, and \(R_{A,j}(a)\) the ratio of catalytic activity per unit enzyme to its reference value. Then:

$$ CL_{int,u}(a) = M_{liver}(a) \sum_j c_{int,u,j,ref}\,R_{E,j}(a)\,R_{A,j}(a) $$

Here, \(M_{liver}(a)\) is liver mass at age \(a\). Each pathway-specific product has units of clearance per unit liver mass; multiplying by liver mass gives clearance for the whole liver. The abundance and activity ratios are dimensionless. When several enzymes contribute to one pathway, or a drug has multiple metabolic routes, include the relevant contributions once each, provided their intrinsic clearances are defined on a compatible basis.

Assumptions Behind the Integration

This scaling expression treats the reference clearance as proportional to enzyme abundance and to activity per unit enzyme. It assumes that the measured or specified abundance is relevant to the drug’s metabolic pathway, that the activity adjustment captures age-related changes not already included in abundance, and that the pathway contributions can be added. It also assumes that other important determinants—such as substrate availability, cofactors, inhibition, and induction—are either represented in the inputs or do not change materially over the modeled comparison.

The distinction between abundance and activity matters. If an age-specific clearance measurement already incorporates both enzyme amount and catalytic function, multiplying it by separate abundance and activity ratios would count those effects twice. Conversely, an abundance ratio alone does not capture a change in activity per enzyme molecule. The inputs must be defined before they are combined.

Worked Example: Scaling Two Pathways to a Three-Year-Old Child

Hypothetical inputs. Suppose reference intrinsic clearances per kilogram of liver are 60 L/h/kg for a CYP pathway and 15 L/h/kg for a UGT pathway. For a three-year-old child, assume CYP abundance and activity ratios of 0.50 and 0.90, and UGT ratios of 0.70 and 1.10. Let liver mass be 0.75 kg.

For CYP, the age-adjusted intrinsic clearance per kilogram of liver is:

$$ 60\ \mathrm{L/h/kg}\times 0.50\times 0.90 =27.00\ \mathrm{L/h/kg} $$

For UGT, it is:

$$ 15\ \mathrm{L/h/kg}\times 0.70\times 1.10 =11.55\ \mathrm{L/h/kg} $$

The summed intrinsic clearance per kilogram of liver is \(27.00+11.55=38.55\) L/h/kg. Scaling by the child’s liver mass gives:

$$ CL_{int,u} =0.75\ \mathrm{kg}\times 38.55\ \mathrm{L/h/kg} =28.91\ \mathrm{L/h} $$

Thus the hypothetical whole-liver intrinsic clearance is 28.91 L/h. This is not yet the predicted hepatic clearance from blood: binding and hepatic blood flow have not been applied.

Binding and the Well-Stirred Liver Model

The well-stirred model represents the liver as a uniformly mixed compartment. Drug enters the liver in blood, can be eliminated according to its unbound intrinsic capacity, and leaves in blood. Let \(Q_h\) be hepatic blood flow, \(C_{in,b}\) the incoming whole-blood drug concentration, \(C_{out,b}\) the outgoing whole-blood concentration, and \(f_{u,b}\) the unbound fraction in blood. Under the model, the unbound concentration in the liver is taken to equal the unbound concentration in exiting blood. The elimination rate is therefore \(CL_{int,u} f_{u,b} C_{out,b}\), while the net rate delivered by blood is \(Q_h(C_{in,b}-C_{out,b})\). Equating these rates gives:

$$ Q_h(C_{in,b}-C_{out,b}) = CL_{int,u}f_{u,b}C_{out,b} $$

Rearranging, \(C_{out,b}/C_{in,b}=Q_h/(Q_h+f_{u,b}CL_{int,u})\). Hepatic clearance from blood is blood flow multiplied by the fraction extracted, \(1-C_{out,b}/C_{in,b}\). Thus:

$$ CL_h = Q_h\frac{f_{u,b}CL_{int,u}}{Q_h+f_{u,b}CL_{int,u}} $$
Model assumptions: The well-stirred equation assumes uniform mixing within the liver, a steady flow of blood, linear elimination over the concentration range considered, and no separate flow bypassing the modeled liver. It treats the unbound drug concentration as the concentration available for elimination. It is a model approximation; it does not represent every spatial or transport process in the liver.

If plasma unbound fraction \(f_{u,p}\) and the blood-to-plasma concentration ratio \(R_b=C_b/C_p\) are available, and unbound drug concentrations in plasma and blood are assumed equal at binding equilibrium, then:

$$ f_{u,b}=\frac{f_{u,p}}{R_b} $$

This conversion is important because \(CL_h\) in the well-stirred equation is clearance from blood and \(Q_h\) is blood flow. Using \(f_{u,p}\) directly in place of \(f_{u,b}\) is justified only when the relevant blood-to-plasma ratio makes them equal, for example when \(R_b=1\).

Worked Example: Converting Intrinsic to Hepatic Clearance

Continue the hypothetical three-year-old example, with \(CL_{int,u}=28.91\) L/h. Assume \(Q_h=22.0\) L/h, \(f_{u,p}=0.080\), and \(R_b=0.80\). These values are hypothetical and are assumed to apply to this child and drug.

First calculate the unbound fraction in blood:

$$ f_{u,b}=\frac{0.080}{0.80}=0.100 $$

The product representing unbound intrinsic capacity is \(f_{u,b}CL_{int,u}=0.100\times28.91=2.891\) L/h. Applying the well-stirred equation:

$$ CL_h = \frac{22.0\ \mathrm{L/h}\times2.891\ \mathrm{L/h}} {22.0\ \mathrm{L/h}+2.891\ \mathrm{L/h}} = \frac{63.60\ \mathrm{L^2/h^2}}{24.89\ \mathrm{L/h}} =2.555\ \mathrm{L/h} $$

The predicted hepatic clearance is 2.555 L/h. The corresponding extraction ratio is \(CL_h/Q_h=2.555/22.0=0.1161\), or about 11.61%. Because \(f_{u,b}CL_{int,u}\) is small relative to \(Q_h\), clearance is below blood flow and is relatively sensitive to intrinsic capacity and binding.

Why the Integration Is Nonlinear

The equation shows why a proportional change in intrinsic clearance does not always produce the same proportional change in hepatic clearance. If \(f_{u,b}CL_{int,u}\) is much smaller than \(Q_h\), the denominator is approximately \(Q_h\), so \(CL_h\) is approximately \(f_{u,b}CL_{int,u}\). In that low-extraction regime, hepatic clearance is sensitive to both unbound fraction and intrinsic capacity. If \(f_{u,b}CL_{int,u}\) is much larger than \(Q_h\), clearance approaches \(Q_h\), so increased intrinsic capacity has a smaller effect on hepatic clearance.

A short sensitivity justification follows directly by differentiating with respect to intrinsic clearance while holding other quantities fixed:

$$ \frac{\partial CL_h}{\partial CL_{int,u}} = \frac{Q_h^2 f_{u,b}}{(Q_h+f_{u,b}CL_{int,u})^2} $$

This derivative is positive, so increasing intrinsic clearance increases hepatic clearance, but the increase becomes less pronounced as \(f_{u,b}CL_{int,u}\) grows relative to flow. The same expression also makes clear that hepatic clearance cannot exceed \(Q_h\) under this model.

Worked Example: A Low-Extraction Newborn Prediction

Consider a second hypothetical child, a newborn with liver mass 0.30 kg and hepatic blood flow 10.0 L/h. Use the same reference per-kilogram clearances as before. Assume newborn CYP abundance and activity ratios of 0.12 and 0.80, and UGT ratios of 0.25 and 0.90. Also assume \(f_{u,p}=0.15\) and \(R_b=0.75\).

The CYP contribution per kilogram of liver is \(60\times0.12\times0.80=5.760\) L/h/kg. The UGT contribution is \(15\times0.25\times0.90=3.375\) L/h/kg. Their sum is 9.135 L/h/kg, so whole-liver intrinsic clearance is:

$$ CL_{int,u}=0.30\ \mathrm{kg}\times9.135\ \mathrm{L/h/kg} =2.741\ \mathrm{L/h} $$

Convert binding to blood: \(f_{u,b}=0.15/0.75=0.200\). Then \(f_{u,b}CL_{int,u}=0.200\times2.741=0.5482\) L/h, and:

$$ CL_h = \frac{10.0\times0.5482}{10.0+0.5482}\ \mathrm{L/h} = 0.5197\ \mathrm{L/h} $$

The extraction ratio is \(0.5197/10.0=0.05197\), or 5.197%. This predicted clearance reflects the specified enzyme scaling, liver mass, binding, and flow together. The difference from the three-year-old prediction cannot be attributed to age alone: several age-specific inputs differ in these hypothetical calculations.

Diagnostics, Alternatives, and Limitations

Before interpreting an age-specific prediction, check that each input has a clear definition and compatible units. In particular, verify whether intrinsic clearance is per unit liver mass or for the whole liver, whether the abundance ratio is already embedded in the activity measurement, and whether binding is expressed for plasma or whole blood. Confirm that \(Q_h\) and \(CL_h\) use the same volume-per-time units. A computed \(CL_h\) greater than \(Q_h\) signals a calculation or unit error under the well-stirred model.

The well-stirred model is commonly contrasted with the parallel-tube model, which represents blood flowing past hepatocytes rather than being uniformly mixed. The models make different assumptions about intrahepatic concentration profiles and can yield different predictions, particularly when extraction is high. A more spatially detailed dispersion model represents mixing between these idealizations but requires additional information. These alternatives do not remove the need for age-appropriate intrinsic clearance, binding, and flow inputs; they change how liver passage is represented.

The integrated calculation is only as reliable as its components. Enzyme abundance may not capture changes in cellular localization or functional activity. Transporter-mediated uptake or efflux can affect access to metabolic enzymes and is not automatically represented by an enzyme-only intrinsic clearance. Changes in blood binding, nonlinear metabolism, extrahepatic elimination, or disease-related physiology may also require model structure beyond the simple calculation shown here. Treat a prediction as conditional on its assumptions, not as an observed clinical clearance.

Common pitfall: A tempting but incorrect calculation is to multiply adult whole-liver clearance by an enzyme maturation ratio and then multiply by liver mass again. This mixes whole-liver and per-unit-mass quantities, potentially counting organ size twice. Instead, scale a per-unit-liver-mass clearance by the relevant abundance and activity ratios, sum pathways, and multiply by liver mass once. Then apply blood binding and flow through the chosen hepatic model.
1
Define the reference basis.
Check whether each pathway clearance is per unit liver mass or already scaled to a whole liver.
2
Scale pathway capacity.
Apply age-specific abundance and activity ratios to each pathway, then sum the pathway contributions.
3
Scale organ amount.
Multiply the summed per-mass capacity by age-specific liver mass to obtain \(CL_{int,u}\).
4
Apply binding and flow.
Convert \(f_{u,p}\) to \(f_{u,b}\) if needed, then calculate \(CL_h\) using a stated liver model.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A reference pathway clearance is 40 L/h/kg liver. At a specified age, abundance and activity ratios are 0.60 and 0.80, and liver mass is 0.50 kg. What is the whole-liver intrinsic clearance for this pathway?
  2. If \(f_{u,p}=0.12\) and \(R_b=0.60\), what is \(f_{u,b}\), assuming equal unbound concentrations in plasma and blood?
  3. Why does the well-stirred model predict \(CL_h\) to be less than or equal to \(Q_h\)?
  4. A prediction gives \(CL_{int,u}\) for the whole liver, but the analyst multiplies it by liver mass again. What is wrong?

Solutions. (1) Per-kilogram clearance is \(40\times0.60\times0.80=19.2\) L/h/kg; whole-liver clearance is \(19.2\times0.50=9.600\) L/h. (2) \(f_{u,b}=0.12/0.60=0.200\). (3) From \(CL_h=Q_h[f_{u,b}CL_{int,u}/(Q_h+f_{u,b}CL_{int,u})]\), the bracketed fraction is nonnegative and less than or equal to 1, so \(CL_h\leq Q_h\). (4) Liver mass has already been incorporated into whole-liver intrinsic clearance; multiplying by mass again double-counts organ size and gives incompatible scaling.