Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Co-Maturation of Enzymes and Transporters in Hepatic Clearance

Pediatric Ontogeny and Scaling · Tutorial 580 of 799

Co-Maturation of Enzymes and Transporters in Hepatic Clearance

Learn to combine maturing hepatic uptake, efflux, and metabolic capacity in a clearance model, and interpret how their competing changes affect intracellular exposure and systemic clearance.

Advanced 10 min read

What You'll Learn

  • Derive effective hepatic intrinsic clearance from uptake, efflux, and metabolism.
  • Distinguish intracellular concentration from net removal of drug from blood.
  • Combine separate developmental profiles without assuming they mature together.
  • Calculate hepatic clearance under the well-stirred liver model.
  • Recognize boundary cases where changing a transporter or enzyme has no effect on net clearance.
  • Check assumptions and identify when a linear, steady-state hepatic model is inadequate.

Why Co-Maturation Matters

Hepatic uptake, efflux, and metabolism act in sequence but need not mature at the same rate. Uptake moves drug from blood into hepatocytes; efflux returns some of that drug to blood; metabolism removes drug from the intracellular pool. A developmental increase in one process can therefore be offset, amplified, or rendered irrelevant by the levels of the others. Predicting hepatic clearance from the maturation of a single enzyme or transporter can miss these interactions.

The earlier material on scaling in vitro transporter abundance to pediatric tissues provides a basis for age-specific transporter inputs. Here the focus is how those inputs combine with enzyme activity inside a simplified hepatic model. The model separates two outputs that should not be conflated: intracellular concentration relative to blood, and net removal from blood by the liver.

A Linear Hepatocyte Model

Consider one well-mixed hepatocyte compartment at a fixed age. Let \(C_b\) be the drug concentration in blood entering the liver and \(C_h\) the concentration in the hepatocyte compartment. Let \(CL_U\) denote hepatic uptake clearance, \(CL_E\) efflux clearance, and \(CL_M\) metabolic clearance. For this model, all three clearances use compatible concentration bases and units, such as mL/min/kg liver. The clearances are nonnegative.

Assume uptake flux is \(CL_U C_b\), while the efflux and metabolic fluxes are \(CL_E C_h\) and \(CL_M C_h\). At steady state, the intracellular amount is not changing, so uptake balances efflux plus metabolism. If \(CL_E+CL_M>0\), this gives:

$$ CL_U C_b=(CL_E+CL_M)C_h, \qquad \frac{C_h}{C_b}=\frac{CL_U}{CL_E+CL_M}. $$

Net removal from blood is uptake minus the drug returned by efflux. Substituting the steady-state intracellular concentration gives an effective intrinsic clearance:

$$ CL_{\mathrm{int,eff}} =\frac{CL_U C_b-CL_E C_h}{C_b} =CL_U\left(1-\frac{CL_E}{CL_E+CL_M}\right) =\frac{CL_U CL_M}{CL_E+CL_M}. $$
Definition: In this linear, steady-state model, \(CL_{\mathrm{int,eff}}\) is the net hepatic removal clearance resulting from uptake, efflux, and metabolism together. It is not the same quantity as \(CL_M\) alone, nor does it by itself determine \(C_h/C_b\).

The result has a useful interpretation. Uptake supplies drug to the hepatocyte, but only the fraction not returned through efflux is available for net removal through metabolism. If metabolism is absent, uptake may still occur, but at steady state all entering drug must return through efflux in this simplified model. If efflux is absent, all drug entering the cell is metabolized at steady state.

How Developmental Profiles Enter

For age \(a\), represent the three clearances as age-specific quantities, for example \(CL_U(a)\), \(CL_E(a)\), and \(CL_M(a)\). Each may be obtained from its own developmental profile and tissue-scaling inputs. The combined result is evaluated at that age:

$$ CL_{\mathrm{int,eff}}(a) =\frac{CL_U(a)CL_M(a)}{CL_E(a)+CL_M(a)}, \qquad \frac{C_h(a)}{C_b(a)} =\frac{CL_U(a)}{CL_E(a)+CL_M(a)}. $$

This formulation does not require uptake, efflux, and metabolism to share one maturation function. Indeed, assigning a common age multiplier to all three can conceal important differences in developmental timing. Keep body-size scaling conceptually separate: the clearances above can be expressed per unit liver mass, with liver mass and hepatic blood flow handled separately when calculating whole-organ or systemic quantities. In the examples below, values are normalized per kg liver and blood flow and binding are held fixed unless stated otherwise.

Worked Example: Combining Three Hypothetical Maturation Profiles

Suppose three hypothetical age-specific profiles give the following clearances, all in mL/min/kg liver: Profile A has \(CL_U=60\), \(CL_E=30\), and \(CL_M=30\); Profile B has \(CL_U=30\), \(CL_E=8\), and \(CL_M=48\). A hypothetical adult reference has \(CL_U=90\), \(CL_E=30\), and \(CL_M=60\). Calculate the effective intrinsic clearance and intracellular-to-blood concentration ratio for each.

For Profile A:

$$ CL_{\mathrm{int,eff}}=\frac{60(30)}{30+30} =\frac{1800}{60}=30.00\ \mathrm{mL/min/kg\ liver}, \qquad \frac{C_h}{C_b}=\frac{60}{60}=1.000. $$

For Profile B:

$$ CL_{\mathrm{int,eff}}=\frac{30(48)}{8+48} =\frac{1440}{56}=25.714\ \mathrm{mL/min/kg\ liver}, \qquad \frac{C_h}{C_b}=\frac{30}{56}=0.5357. $$

For the adult reference:

$$ CL_{\mathrm{int,eff}}=\frac{90(60)}{30+60} =\frac{5400}{90}=60.00\ \mathrm{mL/min/kg\ liver}, \qquad \frac{C_h}{C_b}=\frac{90}{90}=1.000. $$

Profile B has lower effective intrinsic clearance and a lower intracellular-to-blood concentration ratio than Profile A. These two outputs are related through the model, but they are not interchangeable: \(C_h/C_b\) describes the intracellular level at a given blood concentration, whereas \(CL_{\mathrm{int,eff}}\) describes net removal from blood.

From Intrinsic Clearance to Hepatic Clearance

For a well-stirred liver, hepatic blood flow and the unbound fraction in blood also affect systemic clearance. Let \(Q_H\) be hepatic blood flow and \(f_{u,b}\) the unbound fraction in blood. Under the well-stirred assumption, the hepatic clearance is:

$$ CL_H=\frac{Q_H f_{u,b}CL_{\mathrm{int,eff}}} {Q_H+f_{u,b}CL_{\mathrm{int,eff}}}. $$

This model treats the liver as a single, well-mixed compartment exposed to blood flow \(Q_H\). It assumes linear clearance over the concentration range considered and uses an intrinsic clearance defined on a compatible basis. The expression implies \(CL_H\) cannot exceed \(Q_H\). In the next calculation, blood flow, binding, and liver-mass normalization are held constant so that differences in \(CL_H\) arise from the combined uptake, efflux, and metabolic clearances alone.

Worked Example: Translating Co-Maturation into Hepatic Clearance

Use the three profiles above with hypothetical \(Q_H=90\) mL/min/kg liver and \(f_{u,b}=0.10\) for each. For Profile A, \(f_{u,b}CL_{\mathrm{int,eff}}=0.10(30)=3\) mL/min/kg liver, so:

$$ CL_H=\frac{90(3)}{90+3} =\frac{270}{93} =2.903\ \mathrm{mL/min/kg\ liver}. $$

For Profile B, the unrounded effective intrinsic clearance is \(1440/56=180/7\). Thus \(f_{u,b}CL_{\mathrm{int,eff}}=18/7\), and:

$$ CL_H=\frac{90(18/7)}{90+18/7} =\frac{1620/7}{648/7} =2.500\ \mathrm{mL/min/kg\ liver}. $$

For the adult reference, \(f_{u,b}CL_{\mathrm{int,eff}}=0.10(60)=6\), giving:

$$ CL_H=\frac{90(6)}{90+6} =\frac{540}{96} =5.625\ \mathrm{mL/min/kg\ liver}. $$

Under these fixed-flow and fixed-binding assumptions, Profile A predicts slightly greater hepatic clearance than Profile B, while the adult reference predicts greater clearance than either. These are model-based comparisons, not claims about real children or a particular drug. If hepatic blood flow, binding, liver mass, or other clearance pathways vary across the comparison, they must also be represented before interpreting systemic differences.

Why the Direction of Change Needs Conditions

The formula for effective intrinsic clearance makes it possible to state precisely how changing one pathway affects net removal while holding the other two fixed. For \(CL_E+CL_M>0\), its sensitivities are:

$$ \frac{\partial CL_{\mathrm{int,eff}}}{\partial CL_M} =\frac{CL_U CL_E}{(CL_E+CL_M)^2}, \qquad \frac{\partial CL_{\mathrm{int,eff}}}{\partial CL_E} =-\frac{CL_U CL_M}{(CL_E+CL_M)^2}. $$

Because the clearances are nonnegative, increasing \(CL_M\) weakly increases \(CL_{\mathrm{int,eff}}\), and increasing \(CL_E\) weakly decreases it. The change from increasing \(CL_M\) is strict when \(CL_U>0\) and \(CL_E>0\); the change from increasing \(CL_E\) is strict when \(CL_U>0\) and \(CL_M>0\). Similarly, increasing \(CL_U\) weakly increases effective intrinsic clearance, strictly so when \(CL_M>0\). The word “weakly” matters: at boundary values, a change in one pathway can leave net removal unchanged.

Worked Example: Checking Boundary Cases

Let \(CL_U=1\) and \(CL_E=0\), in arbitrary but compatible clearance units. If \(CL_M\) increases from 1 to 2, effective intrinsic clearance is:

$$ \frac{1(1)}{0+1}=1 \quad\longrightarrow\quad \frac{1(2)}{0+2}=1. $$

Metabolism increases, but net clearance does not: with no efflux, all drug entering the cell was already metabolized at steady state. Now let \(CL_U=1\) and \(CL_M=0\). Increasing \(CL_E\) from 1 to 2 gives:

$$ \frac{1(0)}{1+0}=0 \quad\longrightarrow\quad \frac{1(0)}{2+0}=0. $$

With no metabolism, efflux changes the intracellular concentration but cannot create net metabolic removal. These examples show why unqualified statements that more enzyme always raises net clearance or more efflux always lowers it are too strong.

Diagnostics, Alternatives, and Limitations

A useful diagnostic is to inspect both outputs at each age: \(CL_{\mathrm{int,eff}}(a)\) and \(C_h(a)/C_b(a)\). A profile that produces a large change in intracellular concentration but little change in net clearance may be plausible when metabolism or uptake is limiting. Conversely, a predicted clearance change driven by a pathway whose clearance is zero in the model deserves scrutiny. Check that each age-specific input uses the intended tissue and concentration basis, that all clearances share compatible units, and that the assumed steady state is meaningful. A nonzero uptake clearance with both efflux and metabolism equal to zero gives no finite steady-state intracellular concentration in this model.

The model is deliberately compact. It assumes one intracellular compartment, linear and concentration-independent clearances, immediate steady state, and no additional routes such as biliary loss or intracellular binding. If transport or metabolism is saturable, if intracellular concentration gradients matter, or if multiple hepatic cell types have distinct roles, a richer mechanistic model may be needed. The well-stirred liver expression is also not the only liver model: a parallel-tube model represents progressive extraction along the sinusoid, while a dispersion model represents flow heterogeneity. These alternatives change the relationship between intrinsic clearance and organ clearance; they do not remove the need to define how uptake, efflux, and metabolism combine.

A common mistake is to apply an enzyme maturation factor directly to systemic hepatic clearance while holding all transporter contributions implicit. The correction is to update each relevant age-dependent process on its own basis, combine them in the hepatocyte model, and then apply the chosen organ-clearance model with its blood-flow and binding assumptions. Another mistake is to interpret a lower \(C_h/C_b\) as necessarily meaning lower systemic clearance. The equations show that intracellular exposure and net removal are distinct outcomes and should be reported separately when both are relevant.

Key takeaway: Co-maturation is a joint-model problem. Age-specific uptake, efflux, and metabolism determine intracellular exposure and effective intrinsic clearance through different expressions; hepatic blood flow and binding then connect effective intrinsic clearance to hepatic clearance.

Check Your Understanding

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

  1. For \(CL_U=20\), \(CL_E=10\), and \(CL_M=30\) mL/min/kg liver, calculate \(CL_{\mathrm{int,eff}}\) and \(C_h/C_b\).
  2. With \(CL_U>0\) and \(CL_M>0\), what is the direction of change in \(CL_{\mathrm{int,eff}}\) when \(CL_E\) increases? Is the change strict?
  3. Why can increasing \(CL_M\) leave effective intrinsic clearance unchanged when \(CL_E=0\)?
  4. If \(Q_H=100\) mL/min/kg liver, \(f_{u,b}=0.20\), and \(CL_{\mathrm{int,eff}}=25\) mL/min/kg liver, calculate \(CL_H\).

Solutions. (1) \(CL_{\mathrm{int,eff}}=20(30)/(10+30)=600/40=15.00\) mL/min/kg liver, and \(C_h/C_b=20/40=0.500\). (2) It decreases strictly: the derivative is \(-CL_U CL_M/(CL_E+CL_M)^2<0\) under the stated positive conditions. (3) With no efflux, all uptake is ultimately metabolized at steady state, so \(CL_{\mathrm{int,eff}}=CL_U CL_M/CL_M=CL_U\) when \(CL_M>0\). (4) \(f_{u,b}CL_{\mathrm{int,eff}}=0.20(25)=5\), so \(CL_H=100(5)/(100+5)=500/105=4.762\) mL/min/kg liver.