Worked Construction of an IgG Tissue-Distribution Model

From Mechanisms to a Checkable Model

A mechanistic antibody model is useful only if its fluxes connect consistently: antibody that leaves one compartment must enter another, and the whole-body amount should change only through explicitly represented elimination. This tutorial constructs a hypothetical one-organ IgG model that combines perfusion, transcapillary exchange, interstitial convection, lymph return, and FcRn-mediated recycling. The aim is not to represent every tissue in detail, but to make the accounting and numerical checks transparent.

The model follows the compartment and flux ideas introduced in the earlier material on antibody amount balance, transcapillary exchange, interstitial convection, and FcRn salvage. Here those pieces are assembled into one set of differential equations and evaluated with numerical values. All quantities and parameter values below are hypothetical and illustrative, not estimates for a particular antibody or organ.

Model Structure and Assumptions

The central plasma and organ vascular spaces exchange blood at flow \(Q\). Perfusion moves IgG from plasma to organ vasculature at rate \(Q C_p\), and returns it at rate \(Q C_v\); the net transfer is therefore \(Q(C_p-C_v)\). Across the organ capillary, passive exchange is represented by \(PS(C_v-C_i)\), where \(PS\) is a permeability-surface area product. Positive values indicate net movement from organ vasculature to interstitium.

A filtration flow \(Q_f\) carries a fraction \(1-\sigma\) of the vascular concentration into the interstitium. Here \(\sigma\) is the vascular reflection coefficient: \(\sigma=1\) gives no convective IgG transfer in this representation, whereas a lower value permits more solvent-drag transfer. Interstitial fluid and its IgG leave at lymph flow \(Q_L\), with IgG flux \(Q_L C_i\), and lymph returns to central plasma. For the numerical example, filtration and lymph flows are equal, \(Q_f=Q_L\), so the assumed fluid balance is steady.

Interstitial IgG enters the endosomal compartment at first-order rate \(k_u A_i\). Endosomal IgG returns to central plasma at rate \(k_r A_e\), or is degraded at rate \(k_d A_e\). The model also includes systemic clearance from plasma, at rate \(CL C_p\). These are simplifying assumptions: concentrations are well mixed within each compartment, the listed parameters are constant, and no binding, spatial gradients within a compartment, or additional elimination pathways are represented.

The terms are amount rates. For example, \(PS(C_v-C_i)\) has units \((\mathrm{L/h})(\mathrm{mg/L})=\mathrm{mg/h}\), as does \(k_u A_i\). This dimensional check helps catch errors such as multiplying a flow by an amount instead of a concentration.

Why the Whole-Body Balance Is a Strong Check

Add the four compartment equations. The perfusion terms cancel because the amount lost from plasma is gained by organ vasculature, and vice versa. Capillary transfer \(J\) cancels between vascular and interstitial compartments. Lymph loss from the interstitium cancels against lymph return to plasma. Endosomal uptake and recycling also cancel as internal transfers. What remains is

Thus total modeled IgG decreases only through the two explicit elimination routes: systemic clearance and endosomal degradation. This result is a bookkeeping identity, not a claim that those are the only biological elimination mechanisms. If a simulation shows an additional change in total amount, a flux has been omitted, duplicated, or assigned an inconsistent sign.

For a simple numerical solution, the explicit Euler update over a time step \(\Delta t\) is \(A_j(t+\Delta t)=A_j(t)+\Delta t\,(dA_j/dt)\), for each compartment \(j\). It approximates the amount change using the rates at the beginning of the step. It is most appropriate here for illustrating short steps; it is not an exact solution over a long interval.

Checking the Coupled Fluxes at the Next Time Point

After the first step, the compartments no longer have equal concentrations. The next calculation checks that the model uses the changed concentrations, rather than reusing initial fluxes. This distinction matters in dynamic models: gradients and resulting fluxes evolve as amounts move between spaces.

What FcRn Salvage Means in This Model

In this formulation, FcRn-mediated salvage is represented by an endosomal return flux competing with endosomal degradation. If an amount enters the endosome and remains there until one of these two first-order exits occurs, the probability of recycling is \(k_r/(k_r+k_d)\), and the probability of degradation is \(k_d/(k_r+k_d)\). These fractions follow from competing exponential exit rates; they apply to the simplified endosomal compartment, not necessarily to every biological route taken by IgG.

Diagnostics, Alternatives, and Limitations

A useful numerical check is to calculate each transfer once, name its direction, and use that same value with opposite signs in the two affected compartments. Then sum all compartment rates and compare with the elimination terms. Also inspect concentrations and amounts at every step: negative amounts or unexpectedly large oscillations indicate that the numerical step may be too large, the signs may be wrong, or units may be inconsistent. Here the organ vascular space has a relatively rapid perfusion turnover, \(Q/V_v=4.00\) h\(^{-1}\), so the step size should be small enough to resolve its changes. A smaller Euler step or a suitable ODE solver can improve numerical integration, but neither can repair an incorrect model or flux definition.

This lumped model is one modeling choice, not a universal representation of antibody distribution. A model that omits explicit organ vascular space could combine vascular exchange differently; a model with multiple tissues would allocate perfusion and lymph return across organs. More detailed representations may separate endothelial and other cellular uptake, include binding or spatial gradients, or use concentration-dependent processes. Those alternatives need additional assumptions and data. The compact version here is valuable because its fluxes and mass balance are easy to inspect, but its predictions should not be treated as organ-specific evidence without appropriate parameterization and evaluation.

• Sign error: Giving capillary transfer the same sign in vascular and interstitial balances creates or destroys antibody internally. The correction is to subtract \(J\) from \(A_v\) and add \(J\) to \(A_i\).
• Wrong concentration for a flow: Lymph flow is multiplied by \(C_i\), not by \(A_i\). The correction follows from flow times concentration having units of amount per time.
• Misreading salvage: Recycling returns an amount to plasma; it is not an additional elimination term. Only the degraded fraction contributes to \(k_dA_e\) in this model.
• Overinterpreting a coarse step: Euler updates use beginning-of-step rates. A large step can misrepresent evolving gradients, even when the continuous equations are correctly written.