Why Structural Checks Matter
An antibody PBPK model can produce a smooth plasma concentration-time profile and still be structurally wrong. A duplicated return flux can create antibody, a misplaced degradation term can remove it twice, and a concentration calculated with the wrong tissue volume can misstate the amount available for transport. Adjusting parameters may improve the fit without fixing any of these problems.
The useful starting point is an amount balance. A transport flux leaving one compartment must enter another compartment unless it represents a defined loss from the modeled system. This tutorial uses that principle to check four recurring errors in antibody disposition models. Amounts are denoted by \(A\), concentrations by \(C\), and fluxes by \(J\); amounts are in mg, concentrations in mg/L, and fluxes in mg/day unless stated otherwise.
The checks below assume that compartment amounts are the states, each compartment is well mixed at its defined concentration, and each internal transfer is represented once in each affected balance. These bookkeeping assumptions do not imply that a particular compartment structure is biologically complete. They do make it possible to detect whether the equations conserve antibody as intended.
Use the Whole-System Balance as a Reference
Consider a tissue with vascular amount \(A_V\), interstitial amount \(A_I\), lymph amount \(A_L\), and endosomal amount \(A_E\). Let \(J_{\mathrm{cap}}\) be net transcapillary flux, positive from vascular space to interstitium; \(J_{I,L}\) be interstitial-to-lymph flux; \(J_{L,V}\) be lymph-to-vascular return; and \(J_{\mathrm{up}}\) be uptake into endosomes. Let \(J_{\mathrm{salv}}\) be endosomal salvage returned to vascular space, and \(J_{\mathrm{deg}}\) be irreversible degradation. A simple balance is:
Adding these equations cancels every internal transfer: capillary exchange, lymphatic transport, uptake, and salvage. The remaining whole-system balance is:
Thus, without degradation or another explicitly modeled external loss, the sum of modeled amounts must remain constant after dosing stops. With degradation, the total falls at the degradation rate. This simple sum is a targeted structural check, not a claim that every tissue must use exactly these four compartments.
Check 1: Lymphatic Flux Must Have One Route
Lymphatic drainage can be represented with an explicit lymph compartment or as a direct return from interstitium to vascular space. Either choice can be internally consistent. The error is combining both representations for the same physical flow: for example, transferring antibody from interstitium into lymph and from lymph to vascular space, while also adding a direct interstitial-to-vascular return for that same lymph flow.
Worked Example: Detecting Double-Counted Lymph Return
Suppose the model has an explicit lymph compartment. At a particular time, \(J_{I,L}=0.08\) mg/day and \(J_{L,V}=0.08\) mg/day. The model also adds a direct return \(J_{I,V}=0.08\) mg/day to vascular space, without removing that direct amount from interstitium. The intended path is interstitium to lymph to vascular space; the direct term is a duplicate.
For the intended path, the contributions to the three relevant amount balances are \(\frac{dA_I}{dt}=-0.08\), \(\frac{dA_L}{dt}=0.08-0.08=0\), and \(\frac{dA_V}{dt}=0.08\), all in mg/day. Their sum is zero, as an internal transfer requires.
With the extra direct return added only to the vascular equation, the vascular contribution becomes \(0.08+0.08=0.16\) mg/day, while interstitium still loses 0.08 mg/day and lymph has no net change. The total derivative is therefore \(-0.08+0+0.16=0.08\) mg/day. The model creates 0.08 mg/day through an internal route. The correction is to remove the direct return or to use a direct-return formulation instead of the explicit lymph path—not to tune the lymph flow parameter.
A practical diagnostic is to draw the route and label each arrow once. For every arrow, identify the exact source and destination state variables and check that its source term is negative and its destination term is positive with the same magnitude. If a lymph compartment is explicit, a direct return term should not duplicate its return path.
Check 2: Concentrations Must Use the Matching Volume
For a well-mixed compartment, concentration is amount divided by that compartment’s volume: \(C_i=A_i/V_i\). A tissue’s total anatomical volume is not automatically the interstitial fluid volume. If an exchange or uptake flux depends on interstitial concentration, the amount must be divided by the interstitial volume used to define that state. Using a different volume changes the implied concentration and therefore changes the calculated flux.
Worked Example: An Inconsistent Interstitial Volume Basis
Suppose a tissue interstitial state contains \(A_I=0.80\) mg and its interstitial volume is \(V_I=0.20\) L. The appropriate concentration is \(C_I=A_I/V_I=0.80/0.20=4.00\) mg/L. If the model mistakenly divides by a 1.00 L total tissue volume, it obtains \(0.80\) mg/L instead.
Suppose an illustrative uptake expression is \(J_{\mathrm{up}}=PS\,C_I\), with \(PS=0.050\) L/day. Using the correct interstitial concentration gives \(J_{\mathrm{up}}=(0.050)(4.00)=0.200\) mg/day. Using the total tissue volume gives \(J_{\mathrm{up}}=(0.050)(0.80)=0.0400\) mg/day. The error understates the flux by a factor of five for these values.
The transfer may still appear dimensionally valid because both calculations yield mg/day. Dimensional analysis alone cannot verify that the correct physical volume was used. Check the definition of each concentration, the volume associated with its amount state, and whether a parameter such as \(PS\) was defined for that same concentration basis.
When a model uses concentration states rather than amount states, the volume basis also enters the differential equations. If \(A_i=V_iC_i\) and \(V_i\) is constant, then \(dA_i/dt=V_i\,dC_i/dt\). Consequently, equal-and-opposite amount fluxes do not generally appear as equal-and-opposite concentration derivatives when compartment volumes differ. The conservation check should be performed on amounts, or on concentration derivatives multiplied by their corresponding volumes.
Check 3: Endosomal Salvage Is Not the Same as Degradation
Endosomal handling separates antibody fates. Uptake moves antibody into the endosomal compartment; salvage returns a portion to the vascular or other specified extracellular space; degradation is an irreversible loss from the modeled antibody system. If degradation is intended but absent, the model retains antibody that should have been removed. Conversely, subtracting degradation from both endosomal and systemic amounts counts the same loss twice.
Worked Example: Finding a Missing Endosomal Degradation Term
Suppose endosomal uptake is \(J_{\mathrm{up}}=0.12\) mg/day and salvage is \(J_{\mathrm{salv}}=0.09\) mg/day. The unbound endosomal antibody amount is \(A_{E,u}=0.060\) mg, and the specified first-order degradation rate constant is \(k_{\mathrm{deg}}=0.50\) day\(^{-1}\). The degradation flux is \(J_{\mathrm{deg}}=k_{\mathrm{deg}}A_{E,u}=(0.50)(0.060)=0.030\) mg/day.
The correct endosomal balance at this instant is \(\frac{dA_E}{dt}=0.12-0.09-0.03=0\) mg/day. If degradation is omitted, the equation instead gives \(0.12-0.09=0.030\) mg/day, so endosomal antibody accumulates. Across the whole modeled system, uptake and salvage cancel as internal transfers; omitting degradation changes the total derivative from \(-0.030\) mg/day to zero. The equations therefore fail to represent the stated irreversible catabolic loss.
This check does not require all endosomal antibody to be degraded. It requires the degradation term to be included only for the modeled fraction and to be removed from the endosomal amount balance, not also subtracted a second time from the vascular or whole-system balance.
Check 4: Exchange Terms Must Cancel
For passive transcapillary exchange, a common illustrative form is \(J_{\mathrm{cap}}=PS(C_V-C_I)\), with positive flux defined from vascular space to interstitium. Here \(PS\) has units L/day and the concentrations have units mg/L, giving \(J_{\mathrm{cap}}\) in mg/day. The vascular balance receives \(-J_{\mathrm{cap}}\), and the interstitial balance receives \(+J_{\mathrm{cap}}\). Other models may use more detailed exchange terms, but the internal amount accounting still has to cancel.
Worked Example: Detecting a Sign Error in Exchange
Let \(PS=0.10\) L/day, \(C_V=5.0\) mg/L, and \(C_I=2.0\) mg/L. Then \(J_{\mathrm{cap}}=0.10(5.0-2.0)=0.30\) mg/day, positive from vascular to interstitium. The correct exchange contributions are \(-0.30\) mg/day to \(dA_V/dt\) and \(+0.30\) mg/day to \(dA_I/dt\), which sum to zero.
If the interstitial equation mistakenly uses \(-J_{\mathrm{cap}}\), both compartments lose 0.30 mg/day. Their combined exchange contribution is \(-0.30-0.30=-0.60\) mg/day. This is not a new elimination pathway: it is a sign error that destroys antibody during internal exchange. Reversing the sign in one of the two balances restores conservation.
A Short Derivation and a Practical Audit
The cancellation rule follows directly from the definition of an internal flux. For one transfer \(J_{i,j}\) from compartment \(i\) to compartment \(j\), its contributions are \(-J_{i,j}\) to \(dA_i/dt\) and \(+J_{i,j}\) to \(dA_j/dt\). Summing the two contributions gives:
Applying this to every internal transfer leaves only inputs and losses that cross the boundary of the full modeled system. This yields a direct audit procedure:
Write down every vascular, interstitial, lymphatic, and endosomal amount included in the model, including any explicitly separated endosomal species.
Record its source, destination, direction convention, units, and whether it is internal or crosses the system boundary. Follow lymphatic return as a route, not merely as an extra plasma input.
For each concentration used in a flux, verify \(C_i=A_i/V_i\) with the volume belonging to that compartment and confirm the flux parameters use compatible units and definitions.
Add all amount balances and simplify. Internal transfers should cancel; any remainder must match defined inputs, degradation, or other external losses.
Diagnostics, Alternatives, and Limits
A mass-balance audit is complementary to fitting plasma data. A good fit does not prove that lymph return is represented once, that tissue concentrations use the right volumes, or that endosomal degradation is present. Conversely, a mass-conserving model may still have incorrect parameters, an inappropriate compartment structure, or poor predictions. Conservation is a necessary bookkeeping check, not a sufficient validation criterion.
A useful numerical diagnostic is to evaluate the summed derivative at a chosen state and compare it with the specified external input and loss. A nonzero residual after internal terms are summed points to a structural or coding discrepancy. Unit checks can reveal incompatible flow and concentration units, while route diagrams help expose duplicate pathways that equations alone may obscure. A conservation residual near zero does not rule out two errors that happen to cancel; inspect each transfer individually as well as the total.
An explicit lymph compartment preserves information about lymph transit and its amount. A direct interstitial-to-vascular return can be a simpler approximation when that transit is not represented. These alternatives are not interchangeable within the same balance: the chosen structure must define one route and account for its source and destination. Similarly, amount-based states make conservation especially transparent, whereas concentration-based states can be used reliably if every volume conversion is applied consistently.
The checks do not establish that a chosen endosomal model captures all cellular handling, nor do they determine the biological value of a degradation rate. They identify whether the model equations implement the stated mechanisms without creating or losing antibody unintentionally. Separate evidence and model evaluation are needed to support the mechanisms and parameter values.
Common Mistakes
- “Both lymph terms are small, so keeping both is harmless.” Even a small duplicated flux creates a systematic amount-balance error. Remove the duplicate route rather than assuming it is negligible.
- “The concentration has the right units, so the volume is correct.” Units do not identify the physical compartment represented by a volume. Use the volume matching the amount state.
- “Salvage handles endosomal antibody, so a degradation term is unnecessary.” Salvage is return to another modeled compartment; degradation is irreversible loss. Include each specified fate once.
- “The plasma curve fits, so exchange must be balanced.” Parameter adjustment can conceal structural errors over a limited observation window. Sum the amount equations directly.
Check Your Understanding
For each question, identify the relevant balance check before consulting the solutions.
- An explicit lymph state receives 0.05 mg/day from interstitium and returns 0.05 mg/day to vascular space. The vascular equation also receives a direct 0.05 mg/day interstitial return, but interstitium has no corresponding direct outflow. What is the whole-system balance error?
- An interstitial amount is 1.2 mg in a 0.30 L interstitial volume. A model instead divides by a 1.5 L total tissue volume. Calculate both concentrations and state the consequence for a flux proportional to interstitial concentration.
- Endosomal uptake is 0.20 mg/day, salvage is 0.14 mg/day, and degradation is 0.02 mg/day. What should \(dA_E/dt\) be, and what is the whole-system amount derivative due to these processes?
- A vascular-to-interstitial exchange flux is 0.40 mg/day. The vascular equation includes \(-0.40\) mg/day and the interstitial equation includes \(+0.40\) mg/day. Does this exchange change the total modeled amount? Explain.
Solutions. (1) The model creates 0.05 mg/day: the direct return adds vascular amount without removing the same amount from interstitium, while the explicit lymph route already accounts for its own transfer. (2) The correct concentration is \(1.2/0.30=4.0\) mg/L; the incorrect value is \(1.2/1.5=0.80\) mg/L. A flux proportional to concentration is underestimated by a factor of five if its other terms are unchanged. (3) \(dA_E/dt=0.20-0.14-0.02=0.04\) mg/day. Uptake and salvage cancel across compartments, so the whole-system derivative is \(-0.02\) mg/day from degradation. (4) No. The two contributions sum to \(-0.40+0.40=0\) mg/day; exchange redistributes antibody internally without changing the total.