Deriving an Endosomal FcRn Recycling Model

Why Derive the Endosomal Model Explicitly?

FcRn can protect IgG from degradation by binding it in the acidic endosome and returning it to the circulation. To represent that mechanism, a model must track not just circulating IgG, but also endosomal free IgG, available FcRn, and their complex. The equations then make clear which IgG is taken up, which is salvaged, and which is lost.

The introductory material on pH-dependent FcRn binding establishes the key distinction: binding is favored in the acidic endosome, whereas IgG is released as the complex reaches the more neutral cell surface. Here that mechanism is translated into compartmental ordinary differential equations (ODEs). The derivation uses amount-based states so that transfer between plasma and the endosomal compartment can be checked directly by mass balance.

Building the Fluxes Before the ODEs

Writing the fluxes first makes the signs in the ODEs easier to verify. Let \(u(t)\) be an external IgG input into plasma, in amount per unit time. Plasma IgG is taken up into the endosomal pool at first-order rate constant \(k_u\), giving an uptake flux \(k_uA\). This is an internal transfer: it leaves \(A\) and enters \(I\).

In the acidic endosome, free IgG and free FcRn bind with association rate constant \(k_{\mathrm{on}}^{a}\), and the complex dissociates with rate constant \(k_{\mathrm{off}}^{a}\). The superscript \(a\) denotes the acidic condition. With amounts as states, the association flux is \(k_{\mathrm{on}}^{a}IR/V_e\); dividing by \(V_e\) converts the amounts to concentrations for the mass-action reaction. The dissociation flux is \(k_{\mathrm{off}}^{a}C\).

Free endosomal IgG is degraded at rate \(k_{\mathrm{deg},I}I\). The complex can recycle at rate \(k_{\mathrm{rec}}C\), returning its IgG to plasma after neutral-pH release, or undergo complex-associated IgG degradation at rate \(k_{\mathrm{deg},C}C\). The latter represents loss of the IgG cargo, not loss of FcRn. Free FcRn has a source \(q_R\) and a first-order loss \(k_{\mathrm{deg},R}R\). All rate constants have units of inverse time, except \(k_{\mathrm{on}}^{a}\), which has units of volume per amount per time; \(q_R\) and \(u(t)\) have units of amount per time.

Deriving the Coupled ODEs

Each derivative is the sum of the fluxes entering that state minus the fluxes leaving it. The plasma equation receives the external input and recycled IgG, and loses IgG through systemic clearance and endosomal uptake. Here \(k_{\mathrm{cl}}\) denotes first-order systemic clearance other than the modeled endosomal uptake.

Endosomal free IgG gains the uptake flux and the products of complex dissociation and loses IgG through binding and free-IgG degradation. The complex loses IgG through recycling and complex-associated degradation.

Free FcRn is consumed by complex formation and replenished by complex dissociation. It also receives FcRn returned after either recycling or complex-associated IgG degradation, and loses FcRn through its own turnover. Thus:

The last two terms in the FcRn equation follow from the stated assumption that FcRn is intact and returned to the free pool after either complex fate. An alternative biological assumption—for example, loss of FcRn with complex degradation—would change these terms and the FcRn balance. That choice should be explicit rather than hidden in a rate constant.

Mass-Balance Checks and a Short Derivation

A useful structural check is to add the three IgG-state equations. Uptake cancels between plasma and endosome. Binding and dissociation also cancel because they only transfer IgG between free and complexed forms. Recycling cancels because it removes complexed IgG from the endosome and returns the same IgG to plasma. What remains is:

This derivation shows that the modeled internal transfers do not create or destroy IgG. The total modeled IgG amount changes only through external input and the specified clearance or degradation processes. It also gives a practical diagnostic: if a coded model does not reduce to this expression, a transfer term may have the wrong sign, be missing, or be counted twice.

Similarly, add the FcRn and complex equations. Association, dissociation, recycling, and complex-associated degradation cancel because, under the stated assumption, these processes move FcRn between \(R\) and \(C\) or return it intact. The result is:

This second check confirms that complex processing does not destroy FcRn in this model. It does not say that total FcRn is constant: receptor synthesis and free-receptor turnover can change it.

How pH Enters the Model

The pH dependence enters through the binding rate constants and through the fate assigned to a complex that reaches the cell surface. For the acidic endosomal reaction, the dissociation constant is \(K_d^{a}=k_{\mathrm{off}}^{a}/k_{\mathrm{on}}^{a}\), in concentration units. A lower acidic \(K_d^{a}\) corresponds to stronger equilibrium binding. Neutral-pH release is represented here by the recycling flux, rather than by adding a separate neutral-pH complex state.

Assumptions, Diagnostics, and Alternatives

This model is useful when the aim is to represent the main endosomal salvage mechanism while retaining explicit states for IgG, FcRn, and complex. Its assumptions should be checked against the intended application:

• Well-mixed pools: each compartment is treated as spatially uniform. Strong spatial gradients or distinct trafficking subpopulations may require additional compartments.
• Mass-action binding: association depends on free IgG and free FcRn, with acidic rate constants. The model does not represent multiple binding sites or heterogeneous affinities.
• Lumped recycling and release: surface transit and neutral-pH dissociation are represented by \(k_{\mathrm{rec}}C\). If surface residence or incomplete release matters, a separate surface state is needed.
• Explicit receptor-fate assumption: FcRn is returned intact after recycling or IgG degradation. If receptor is lost with cargo, the receptor equation must change.

During implementation, check that \(A\), \(I\), \(R\), and \(C\) remain nonnegative for nonnegative initial amounts and inputs, and verify the two summed balances numerically. A negative state often indicates a sign error or an unsuitable numerical solver step. Also check units: the association term must have units of amount per time, which it does only when \(V_e\) converts the amount states to concentrations. If predicted complex exceeds available receptor or IgG, inspect the states and equations; the physical constraints \(C\leq R+C\) and \(C\leq I+C\) should hold for nonnegative states initialized consistently.

There are simpler and more detailed alternatives. A single effective clearance term can summarize net IgG loss but cannot separately predict free IgG, FcRn occupancy, and recycling flux. A fixed-receptor approximation removes the receptor ODE and assumes an effectively constant receptor pool; it is appropriate only when that assumption is defensible. At the other extreme, explicit surface and trafficking compartments can represent delayed release and multiple routes, at the cost of more states and parameters. The present model is a mechanistic middle ground, not a universal description of all intracellular processing.