Modeling Hepatic Uptake-Efflux Interactions

Why Uptake and Efflux Must Be Considered Together

A hepatic uptake transporter moves drug from plasma into hepatocytes. Once inside, drug may return to plasma through basolateral efflux, leave toward bile through canalicular efflux, or be metabolized. These routes compete: changing one pathway can alter both hepatocyte concentration and the fraction of drug that ultimately leaves the plasma compartment. Therefore, the effect of inhibiting a transporter cannot generally be inferred from its direction of transport alone.

The model in this tutorial isolates those competing processes in a single hepatocyte compartment. It is useful for interpreting the direction and approximate size of an interaction, and for checking whether a PBPK model’s pathway-level behavior is internally consistent. It is not a substitute for a full liver model when blood-flow limitation, multiple cell types, nonlinear transport, or other disposition routes are important.

A Steady-State Uptake-Efflux Model

The amount of drug entering hepatocytes per unit time is \(CL_u C_{p,u}\). The three exit fluxes are \(CL_b C_{h,u}\), \(CL_c C_{h,u}\), and \(CL_{met} C_{h,u}\). At steady state, entry equals total exit. Define the total hepatocyte exit clearance as \(D=CL_b+CL_c+CL_{met}\). The steady-state formulas below require \(D>0\).

The ratio \(C_{h,u}/C_{p,u}\) describes hepatocyte exposure relative to unbound plasma exposure in this model. It is not, by itself, a total-tissue-to-total-plasma ratio; binding in plasma and hepatocytes can make total concentrations behave differently.

Uptake into hepatocytes is not necessarily irreversible removal from plasma. Drug returning through the basolateral route re-enters plasma. The net disappearance flux from plasma due to the liver is uptake minus this return flux. Dividing that net flux by \(C_{p,u}\) gives the hepatic clearance referenced to unbound plasma concentration.

This expression makes the competing roles explicit. Canalicular efflux and metabolism are irreversible exits from the modeled hepatocyte compartment, whereas basolateral efflux is return to plasma. Inhibiting uptake reduces entry. Inhibiting canalicular efflux can increase hepatocyte concentration while reducing the fraction of uptake that exits irreversibly; the net result depends on the other exit pathways.

Deriving the Clearance Expression

The result follows directly from mass balance, but the denominator condition matters. At a positive plasma concentration, steady state requires \(CL_u C_{p,u}=D C_{h,u}\). If \(D>0\), division by \(D C_{p,u}\) yields the concentration ratio. Substituting that ratio into the net plasma disappearance flux gives the clearance equation.

Two limiting cases illustrate why the condition \(D>0\) cannot be omitted. If \(CL_c=CL_{met}=0\) and \(CL_b>0\), then \(D=CL_b>0\): all drug entering hepatocytes eventually returns to plasma, so \(CL_h=0\). If \(CL_b=0\) and \(CL_c+CL_{met}>0\), then \(D>0\) and \(CL_h=CL_u\): all uptake is followed by an irreversible exit in this model. But if all three exit clearances are zero, \(D=0\); with positive uptake and plasma concentration, drug accumulates in the hepatocyte compartment and no finite steady-state ratio or clearance formula applies.

Representing Transporter Inhibition

For a specified inhibitor exposure, represent the remaining activity of each affected pathway by a fraction between zero and one. For example, if uptake is inhibited by 50%, set \(CL_{u,I}=0.50CL_u\). If canalicular efflux is inhibited by 75%, set \(CL_{c,I}=0.25CL_c\). The other pathway clearances remain unchanged unless there is evidence that they are also affected. These adjusted clearances are then inserted into the same mass-balance equations; inhibition does not justify deleting the unaffected pathways.

This calculation assumes that the inhibition fractions are appropriate at the relevant transporter site and that the clearances remain linear at the concentrations considered. It also assumes the inhibitor does not independently alter metabolism, protein binding, or another clearance route. Those assumptions should be checked against the interaction mechanism and the PBPK model inputs.

Assumptions, Diagnostics, and Alternatives

The equations are most directly applicable when disposition is linear over the concentration range, the hepatocyte compartment can be treated as well mixed, and the clearances can be regarded as constant. They assume a quasi-steady hepatocyte concentration for the comparison being made. They do not explicitly represent hepatic blood-flow limitation, concentration gradients, transporter saturation, time-varying inhibition, intracellular binding, or changes in enzyme activity. In a full PBPK model, those features may require explicit compartments, concentration-dependent processes, and dynamic simulation.

A useful diagnostic is to inspect all pathway clearances before evaluating the ratios. Confirm that units agree, that the inhibited pathway has been changed in the intended direction, and that \(D>0\). Then check whether the predicted hepatocyte ratio responds consistently: reducing \(CL_u\) alone lowers it; reducing an exit clearance with uptake held fixed raises it, provided the remaining total exit clearance stays positive. The hepatic clearance should not exceed \(CL_u\) under this model, because \(0\leq(CL_c+CL_{met})/D\leq1\).

The compartment calculation also differs from treating uptake as if it were irreversible hepatic clearance. That simpler approximation uses \(CL_u\) as net removal and can overpredict hepatic clearance when basolateral return is substantial. Conversely, adding separate percentage changes for uptake and efflux ignores the shared denominator \(D\) and generally does not reproduce the mass balance. A more detailed PBPK representation is preferable when flow, multiple intracellular regions, or nonlinear kinetics determine exposure; this compact model is best used when its assumptions are defensible and its parameters are interpretable.

Common Mistake: Ignoring the Steady-State Condition

A tempting but incorrect shortcut is to set \(CL_c=CL_{met}=0\), conclude that all uptake returns through basolateral efflux, and immediately state that \(CL_h=0\). This conclusion is valid only when \(CL_b>0\). If \(CL_b=0\) as well, then \(D=0\), so dividing by \(D\) is invalid and there is no finite steady state with positive uptake. Similarly, the limiting conclusion \(CL_h=CL_u\) when basolateral return is zero requires \(CL_c+CL_{met}>0\). Always verify a positive total exit clearance before using the steady-state ratio or clearance formula.