Why Parameterization Errors Matter
A hepatic transporter model can be structurally plausible and still predict the wrong liver exposure if its parameters mix incompatible concentration bases or count the same process more than once. Four errors are especially consequential: correcting a clearance for binding in the wrong direction, treating uptake as net hepatic elimination, representing passive flux twice, and applying abundance scaling to a value that has already been scaled.
The earlier material on uptake assays, binding corrections, transporter activity scaling, and hepatic transport models provides the relevant measurement context. Here the emphasis is on checking the bookkeeping: what concentration a parameter uses, which flux it represents, and what biological amount it has already been scaled to. The examples use hypothetical values.
Total and Unbound Concentrations Are Not Interchangeable
In a hepatocyte incubation, some compound may bind to proteins, plastic, or other assay components. Let \(C_{\mathrm{tot,inc}}\) be the total incubation concentration and let \(f_{u,\mathrm{inc}}\) be the unbound fraction in that incubation. Under the assumption that binding is at equilibrium and the unbound fraction is constant over the concentration range considered,
If an uptake rate \(v\) is divided by the total concentration, the resulting total-based uptake clearance is \(CL_{\mathrm{uptake,tot}}=v/C_{\mathrm{tot,inc}}\). Dividing the same rate by unbound concentration gives \(CL_{\mathrm{uptake,u}}=v/C_{u,\mathrm{inc}}\). Therefore,
This conversion applies only when the original clearance used total concentration as its denominator and the rate is unchanged. If the rate was already divided by unbound concentration, applying the correction again is an error. If the unbound fraction varies with concentration, a single constant correction may not be adequate; measure or model the relevant unbound concentrations instead. For a saturable uptake model, the concentration basis of the fitted parameters must also be explicit rather than silently switching between total and unbound concentration.
Worked Example: Correcting an Uptake Clearance for Incubation Binding
A hypothetical assay has \(C_{\mathrm{tot,inc}}=2.00\ \mu\mathrm{M}\), \(f_{u,\mathrm{inc}}=0.200\), and a measured uptake rate of \(12.0\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\). Calculate the total-based and unbound-based uptake clearances.
Because \(1\ \mu\mathrm{M}=1\ \mathrm{pmol\,\mu L^{-1}}\), the total concentration is \(2.00\ \mathrm{pmol\,\mu L^{-1}}\). The total-based clearance is
The unbound concentration is \(0.200\times2.00=0.400\ \mu\mathrm{M}\), or \(0.400\ \mathrm{pmol\,\mu L^{-1}}\). Thus,
The larger unbound-based value is consistent with dividing the same rate by a smaller concentration. This correction changes the concentration basis of the clearance; it does not create additional uptake.
Uptake Is a Flux Into the Cell, Not Automatically Net Clearance
Uptake clearance describes entry from the extracellular or blood-facing space into hepatocytes, relative to the concentration used in its definition. It does not by itself state how much drug is ultimately eliminated from the body. Once inside a hepatocyte, compound may be metabolized, excreted into bile, or returned to blood. Return to blood can reduce net hepatic removal even when uptake is rapid.
A simple linear model makes this distinction precise. Let \(C_p\) be the relevant unbound plasma concentration, \(A\) the intracellular amount, and \(CL_{\mathrm{uptake}}\) the uptake clearance on a matching concentration basis. Let \(k_{\mathrm{met}}\), \(k_{\mathrm{bile}}\), and \(k_{\mathrm{return}}\) be first-order rate constants for intracellular metabolism, biliary loss, and return to blood, respectively. Their units are inverse time. Assume constant parameters, a well-mixed intracellular pool, and steady state for that pool. Then the uptake input rate is \(CL_{\mathrm{uptake}}C_p\), and the intracellular balance is
At steady state, set \(dA/dt=0\) and solve for \(A\). Net irreversible loss is the sum of metabolic and biliary output, not uptake alone. The resulting clearance from plasma into those loss pathways is
In this simplified model, the fraction of uptake that is ultimately lost through metabolism or bile is less than or equal to one. The result depends on the stated steady-state, linear-rate assumptions; more detailed models may include multiple intracellular spaces, saturable processes, or time-varying concentrations. The central bookkeeping point remains: an uptake parameter is one component in a pathway, not automatically the net clearance.
Worked Example: Distinguishing Uptake Clearance From Net Clearance
Suppose a hypothetical model has \(CL_{\mathrm{uptake}}=24.0\ \mathrm{L\,h^{-1}}\), \(k_{\mathrm{met}}=0.300\ \mathrm{h^{-1}}\), \(k_{\mathrm{bile}}=0.100\ \mathrm{h^{-1}}\), and \(k_{\mathrm{return}}=0.600\ \mathrm{h^{-1}}\). Assume the linear, steady-state model just defined.
The total intracellular exit rate constant is \(0.300+0.100+0.600=1.000\ \mathrm{h^{-1}}\). The fraction of uptake entering irreversible loss pathways is \((0.300+0.100)/1.000=0.400\). Therefore,
The uptake clearance is \(24.0\ \mathrm{L\,h^{-1}}\), while the net clearance through metabolism and bile in this model is \(9.60\ \mathrm{L\,h^{-1}}\). The difference reflects return to blood. Equating the two would overstate net loss under these assumptions.
Represent Passive Flux Once
A transporter assay may measure total accumulation that includes both transporter-mediated uptake and passive entry. At a specified concentration and time, write the measured control rate as \(J_{\mathrm{control}}=J_{\mathrm{active}}+J_{\mathrm{passive}}\). If an inhibitor fully and selectively blocks the transporter without changing passive permeability, the inhibited rate estimates \(J_{\mathrm{passive}}\). The transporter-dependent rate is then the difference \(J_{\mathrm{control}}-J_{\mathrm{inhibited}}\).
A model can represent these components either as an active term plus an independently parameterized passive term, or as a measured total term that already includes both. It must not use the total as the active term and then add passive flux again. The same principle applies when passive flux is written explicitly as a permeability term driven by the unbound concentration difference across a membrane: if another parameter already represents that same passive contribution, adding both counts one physical process twice.
The subtraction interpretation requires its assumptions. If the inhibitor is incomplete, the difference underestimates transporter-mediated flux. If it alters passive permeability, the inhibited rate is not a clean passive estimate. Other active pathways may also remain in the inhibited condition. Those limitations should be handled in the assay interpretation, rather than hidden by labeling the residual as pure passive flux.
Worked Example: Avoiding Double-Counting of Passive Uptake
In a hypothetical assay at the same unbound concentration, the control uptake rate is \(16.0\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\), and the rate with a fully effective, selective transporter inhibitor is \(6.00\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\). Assume the inhibitor leaves passive flux unchanged and there are no other inhibitor-sensitive pathways.
The transporter-dependent rate is \(16.0-6.00=10.0\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\), and the passive rate represented by the inhibited condition is \(6.00\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\). If the concentration is \(2.00\ \mathrm{pmol\,\mu L^{-1}}\), the corresponding clearances are \(5.00\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\) for the transporter-dependent component and \(3.00\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\) for passive flux. Their sum is the observed total clearance, \(8.00\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\).
A model using the transporter-dependent value of \(5.00\) plus passive clearance of \(3.00\) reconstructs the total of \(8.00\). In contrast, using the total value of \(8.00\) as if it were active uptake and adding \(3.00\) again predicts \(11.0\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\). That excess is the passive component counted twice.
Scale Transporter Activity Only Once
In vitro transporter activity may be reported per unit assay protein, per unit transporter abundance, or for a reference expression system. The appropriate scaling depends on which of these bases the reported value already uses. For example, if a measured activity is expressed per milligram of membrane protein, scaling to liver may involve the membrane protein per gram of liver and liver mass. If activity is expressed per amount of transporter, it may instead be multiplied by the target liver transporter amount. Do not apply both routes to the same activity unless the quantities are defined to avoid overlap.
Abundance-ratio scaling follows the same rule. If a reference-system activity has already been multiplied by the target-to-reference abundance ratio, the result is already adjusted for that abundance difference. Multiplying the adjusted result by the ratio again applies the same factor twice. Units provide a practical check: multiplying a whole-liver activity by an abundance ratio should not be done merely because the ratio is available; first determine whether the activity already represents the target abundance.
Worked Example: Checking an Abundance Scaling Factor
A hypothetical reference system has transporter activity of \(0.0200\ \mathrm{nmol\,min^{-1}}\) per picomole of transporter and contains \(200\ \mathrm{pmol}\) of transporter. The target liver is estimated to contain \(600\ \mathrm{pmol}\), so its target-to-reference abundance ratio is \(600/200=3.00\). Assume the activity per picomole is applicable to both systems.
The reference-system activity is \(0.0200\times200=4.00\ \mathrm{nmol\,min^{-1}}\). Scaling that activity once by the abundance ratio gives
The same result follows directly from the per-transporter activity and target abundance: \(0.0200\times600=12.0\ \mathrm{nmol\,min^{-1}}\). Multiplying the already scaled value by 3.00 again would give \(36.0\ \mathrm{nmol\,min^{-1}}\), which is inconsistent with the target abundance of \(600\ \mathrm{pmol}\) under the stated activity assumption.
Practical Checks and Common Failure Modes
Before using a transporter parameter in a PBPK model, trace its measurement and transformation from the assay to the model. The following checks help localize errors without mistaking one problem for another.
| Check | Warning sign | Correction |
|---|---|---|
| Concentration denominator | A clearance is labeled unbound-based, but its original calculation used total concentration without a documented conversion. | Confirm the rate denominator and apply the binding correction once, if its assumptions hold. |
| Direction of binding correction | With \(f_{u,\mathrm{inc}}<1\), a corrected unbound clearance is larger than its total-based counterpart. | Check for a reversed correction or another inconsistency: under the stated conversion, \(CL_{\mathrm{uptake,u}}=CL_{\mathrm{uptake,tot}}/f_{u,\mathrm{inc}}\), so the unbound-based value is larger. |
| Repeated binding correction | The unbound-based clearance is implausibly large relative to the assay rate and concentrations. | Check whether division by \(f_{u,\mathrm{inc}}\) was applied more than once. Repeating that correction makes the value larger still. |
| Meaning of uptake | Uptake clearance is entered as net elimination even though intracellular return is represented. | Keep uptake distinct from downstream metabolism, biliary loss, and return flux. |
| Passive component | A total measured uptake parameter is combined with a separate passive term. | Use a transporter-dependent component plus passive flux, or use the total without adding its passive component again. |
| Scaling basis | An activity already adjusted to target abundance is multiplied by the same abundance ratio again. | Write down the units and biological level of the activity at each scaling step. |
A smaller corrected unbound clearance than the total-based value is a warning of a reversed correction or another inconsistency; it is not evidence of repeated division by \(f_{u,\mathrm{inc}}\). Repeating that division produces an excessively large value. Keeping these diagnoses separate is important because they point to different corrections.
There are reasonable alternatives to each calculation route. Direct measurement of unbound incubation concentration avoids relying on a binding correction, while a measured \(f_{u,\mathrm{inc}}\) can be useful when direct measurement is impractical and binding is adequately characterized. In uptake assays, inhibitor subtraction can isolate a transporter-dependent component under restrictive assumptions; explicit passive-flux modeling can be preferable when permeability data and intracellular concentrations are available. For scaling, one can use protein-per-liver scaling or transporter-abundance scaling, provided the activity units and abundance basis match. None is automatically superior: the data and assumptions determine which is defensible.
These checks do not establish that a chosen transporter mechanism is correct or that the resulting model is supported by in vivo evidence. They address parameter consistency. A well-scaled, correctly corrected parameter can still be uncertain, poorly identifiable, or unsuitable outside the assay conditions; those limitations should remain visible in model evaluation.
Check Your Understanding
Work through each question, then compare your reasoning with the solutions.
- A total-based uptake clearance is \(8.00\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\) and \(f_{u,\mathrm{inc}}=0.250\). What is the unbound-based clearance, assuming the correction applies?
- In the linear steady-state model, what fraction of uptake enters irreversible loss pathways if \(k_{\mathrm{met}}=0.200\ \mathrm{h^{-1}}\), \(k_{\mathrm{bile}}=0.100\ \mathrm{h^{-1}}\), and \(k_{\mathrm{return}}=0.700\ \mathrm{h^{-1}}\)?
- A control assay rate is \(20.0\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\) and the inhibited rate is \(8.00\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\). Under the assumptions of selective, complete inhibition with unchanged passive flux, what are the transporter-dependent and passive rates?
- A reference activity of \(5.00\ \mathrm{nmol\,min^{-1}}\) has already been scaled by a target-to-reference abundance ratio of 2.00. What value results if the ratio is incorrectly applied a second time, and what is the correctly scaled value?
Solutions. (1) \(8.00/0.250=32.0\ \mathrm{\mu L\,min^{-1}\,mg^{-1}}\). The corrected value is larger because it divides the same rate by the smaller unbound concentration. (2) The irreversible fraction is \((0.200+0.100)/(0.200+0.100+0.700)=0.300\), or 30.0%. (3) The transporter-dependent rate is \(20.0-8.00=12.0\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\); the passive rate is \(8.00\ \mathrm{pmol\,min^{-1}\,mg^{-1}}\). (4) The correctly scaled value is \(5.00\ \mathrm{nmol\,min^{-1}}\), because scaling has already been applied. Incorrectly applying the ratio again gives \(5.00\times2.00=10.0\ \mathrm{nmol\,min^{-1}}\).