Why a Tissue Can Have a Binding-Site Barrier
An antibody entering tissue from the vasculature does not necessarily distribute uniformly through the interstitium. If it encounters abundant targets near the tissue surface, binding followed by internalization can remove antibody before much of it reaches deeper regions. The result is a steep spatial concentration gradient: high concentration near the entry boundary and much lower concentration farther in. This is called a binding-site barrier.
The key balance is between movement through tissue and local removal. Diffusion carries free antibody deeper; target-mediated uptake removes it along the way. A characteristic penetration length summarizes this competition. It is not a fixed property of the antibody alone: it depends on tissue diffusivity, target density, and the rates of binding and internalization.
From Target Binding to a Tissue Sink
Consider a uniform tissue represented as a one-dimensional semi-infinite space. Let \(c(x)\) be free-antibody concentration at depth \(x\), and let \(b(x)\) be target-bound antibody concentration. The total target concentration is \(R_T\). Antibody binds at association rate \(k_{\mathrm{on}}\), dissociates at rate \(k_{\mathrm{off}}\), and target-bound antibody is internalized at rate \(k_{\mathrm{int}}\). Assume a constant effective diffusion coefficient \(D\), a fixed free concentration \(C_0\) at the tissue boundary, and a steady spatial profile. Ignore convection, spatial variation in target density, and return of internalized antibody to the extracellular space.
Under local quasi-steady binding, the bound-target balance is association minus dissociation and internalization. Solving that balance gives the bound concentration as a function of local free concentration. Internalization of this bound pool is the local removal flux.
Here \(J(c)\) is the amount of free antibody removed per tissue volume per unit time. If free antibody is low enough that \(k_{\mathrm{on}}c \ll k_{\mathrm{off}}+k_{\mathrm{int}}\), target occupancy is low and the removal flux is approximately proportional to \(c\). The proportionality constant \(k_s\) is an effective first-order sink rate:
The dimensions are consistent: \(k_{\mathrm{on}}\) has units of concentration\(^{-1}\) time\(^{-1}\), \(R_T\) has units of concentration, and the rate constants have units of inverse time. Thus \(k_s\) has units of inverse time. In the linear approximation, the steady diffusion-and-removal balance is \(D\,d^2c/dx^2-k_sc=0\). This approximation describes a sustained sink only when internalization (or another process that removes bound antibody) operates on the timescale of interest. Reversible binding without removal can delay distribution, but cannot by itself maintain this steady-state sink indefinitely.
Deriving the Exponential Profile
For constant \(D\) and \(k_s\), the steady equation has solutions formed from exponentials. The boundary conditions are \(c(0)=C_0\) and \(c(x)\to0\) as \(x\to\infty\). The solution that decays into the tissue is the exponential below; substituting it into the differential equation also identifies the length scale.
Indeed, differentiating \(C_0e^{-x/\lambda}\) twice gives \(d^2c/dx^2=c/\lambda^2\). Substitution yields \((D/\lambda^2-k_s)c=0\), so the nonzero profile requires \(\lambda^2=D/k_s\). Larger \(D\) increases penetration length, whereas a larger effective sink decreases it. In this linear model, multiplying \(R_T\) by four multiplies \(k_s\) by four and divides \(\lambda\) by two. This square-root dependence is important: a fourfold increase in target density does not make the penetration length four times shorter.
Worked Example: Concentration at a Given Tissue Depth
Hypothetical values. Let \(D=10\ \mu\mathrm{m}^2/\mathrm{s}\), \(R_T=1.00\ \mathrm{nM}\), \(k_{\mathrm{on}}=0.100\ \mathrm{nM}^{-1}\mathrm{s}^{-1}\), \(k_{\mathrm{off}}=0.100\ \mathrm{s}^{-1}\), and \(k_{\mathrm{int}}=0.400\ \mathrm{s}^{-1}\). Set the boundary concentration to \(C_0=0.0100\ \mathrm{nM}\). The low-occupancy approximation is appropriate because \(k_{\mathrm{on}}C_0=0.00100\ \mathrm{s}^{-1}\), much less than \(k_{\mathrm{off}}+k_{\mathrm{int}}=0.500\ \mathrm{s}^{-1}\).
First calculate the effective sink rate:
The penetration length is \(\lambda=\sqrt{10/0.0800}=\sqrt{125}=11.1803\ \mu\mathrm{m}\), or \(11.18\ \mu\mathrm{m}\) to four significant figures. At \(x=20.0\ \mu\mathrm{m}\), the fraction of boundary concentration remaining is
Thus \(c(20.0)\approx(0.0100)(0.1672)=0.001672\ \mathrm{nM}\). The tissue concentration is about 16.72% of the boundary concentration at this depth, under the stated steady, linear model.
How Target Density Changes Penetration
Because the effective sink rate is proportional to \(R_T\), target-rich tissue can capture more antibody per unit time near the entry region. The next calculation changes only target density; all other assumptions and parameter values remain as in the preceding example. Holding \(D\), \(k_{\mathrm{on}}\), \(k_{\mathrm{off}}\), and \(k_{\mathrm{int}}\) fixed makes the density effect directly interpretable.
Worked Example: Fourfold Higher Target Density
Hypothetical comparison. Increase target concentration from \(1.00\) to \(4.00\ \mathrm{nM}\), retaining the same kinetic constants. The effective sink rate becomes \(k_s=(0.400)(0.100)(4.00)/(0.500)=0.320\ \mathrm{s}^{-1}\). Therefore,
The penetration length is half the baseline value, as predicted by its inverse square-root dependence on target density. At \(20.0\ \mu\mathrm{m}\), the concentration fraction is \(\exp[-20.0/\sqrt{31.25}]=\exp(-3.577709\ldots)\approx0.02793\). The predicted free concentration is \(0.0100\times0.02793=0.0002793\ \mathrm{nM}\). This is lower than in the baseline case, but it is not zero: the exponential profile declines continuously.
Binding Kinetics Matter Even at the Same Affinity
The equilibrium dissociation constant is \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\), but equal \(K_D\) does not guarantee equal tissue penetration. In the sink expression, both rates affect how much target-bound antibody is available for internalization. For example, with a fixed \(k_{\mathrm{int}}\), a slower association rate can reduce capture before antibody diffuses onward, even if dissociation slows proportionally and \(K_D\) stays the same. This comparison assumes low occupancy and the same target density and diffusivity.
Worked Example: Equal Affinity, Different Kinetic Rates
Hypothetical comparison. Use \(D=10\ \mu\mathrm{m}^2/\mathrm{s}\), \(R_T=1.00\ \mathrm{nM}\), and \(k_{\mathrm{int}}=0.400\ \mathrm{s}^{-1}\) in both cases. For antibody A, let \(k_{\mathrm{on}}=0.100\ \mathrm{nM}^{-1}\mathrm{s}^{-1}\) and \(k_{\mathrm{off}}=0.100\ \mathrm{s}^{-1}\); for antibody B, use \(k_{\mathrm{on}}=0.0100\ \mathrm{nM}^{-1}\mathrm{s}^{-1}\) and \(k_{\mathrm{off}}=0.0100\ \mathrm{s}^{-1}\). Both have \(K_D=1.00\ \mathrm{nM}\). At \(C_0=0.0100\ \mathrm{nM}\), the low-occupancy condition holds in both cases.
For A, \(k_s=(0.400)(0.100)(1.00)/(0.100+0.400)=0.0800\ \mathrm{s}^{-1}\), so \(\lambda_A=\sqrt{10/0.0800}=\sqrt{125}=11.1803\ \mu\mathrm{m}\). For B,
At \(20.0\ \mu\mathrm{m}\), the fractions are \(\exp[-20.0/\sqrt{125}]\approx0.1672\) for A and \(\exp[-20.0/\sqrt{1025}]\approx0.5354\) for B. Under these hypothetical assumptions, B penetrates farther despite having the same equilibrium affinity. The distinction arises because the slower association of B produces a weaker internalization sink when \(k_{\mathrm{int}}\) is held fixed.
Assumptions, Diagnostics, and Alternatives
The exponential profile is useful because it gives a transparent scale, but its assumptions should be checked before interpreting a fitted or calculated \(\lambda\). In particular, target density and diffusivity are treated as spatially uniform; the boundary concentration is fixed; binding is locally quasi-steady; and uptake is approximately linear in free concentration. If the antibody concentration is high enough to occupy a substantial fraction of targets, the sink saturates. From the full expression for \(J(c)\), increasing \(c\) eventually produces diminishing increases in uptake, so a single concentration-independent \(k_s\) no longer describes the whole profile. The low-concentration tail may still be approximately exponential even when the region near the boundary is not.
Useful warning signs include a predicted boundary concentration for which \(k_{\mathrm{on}}C_0\) is not small relative to \(k_{\mathrm{off}}+k_{\mathrm{int}}\), large differences in measured concentrations at nearby tissue depths, or systematic disagreement between the exponential profile and spatial observations. These indicate that the linear, uniform-sink approximation may be inadequate; they do not by themselves identify which assumption failed. Target heterogeneity, finite tissue dimensions, spatially varying diffusion, convection, reversible retention, and time-dependent exposure can all change the profile. A steady-state penetration length should not be interpreted as a transient penetration distance without checking that the steady profile is reached on the timescale studied.
A well-mixed tissue compartment is a useful alternative when the question concerns average tissue amount or exposure and spatial gradients are unimportant. It cannot represent a near-vessel-to-deep-tissue gradient. An equilibrium partition approach can summarize an eventual distribution under appropriate reversible-binding assumptions, but it does not represent a sustained internalization sink or the distance-dependent profile. The reaction–diffusion approximation addresses that spatial competition, at the cost of stronger assumptions and the need for spatial information or justified parameter values. The choice of spatial discretization is a separate modeling decision.
| Approach | What it represents | Main limitation for a binding-site barrier |
|---|---|---|
| Well-mixed tissue compartment | Average tissue concentration or amount | Does not resolve depth-dependent gradients |
| Equilibrium partitioning | Distribution under equilibrium assumptions | Does not by itself represent sustained spatial removal |
| Linear reaction–diffusion model | Steady spatial gradient from diffusion and first-order uptake | Requires low occupancy and uniform tissue properties |
| Saturable or time-dependent spatial model | Concentration-dependent uptake or evolving profiles | Requires more kinetic detail and data to constrain it |
Check Your Understanding
Use the linear steady model unless a question asks about its assumptions. Each solution follows the question.
- If \(R_T\) increases ninefold while all other parameters stay fixed, how does \(\lambda\) change? Solution: Since \(\lambda\propto1/\sqrt{R_T}\), the new length is the old length divided by \(\sqrt{9}=3\).
- If \(D\) increases fourfold with the sink unchanged, how does \(\lambda\) change? Solution: Since \(\lambda\propto\sqrt{D}\), the new penetration length is twice the old length.
- Why can reversible binding without internalization not maintain the same steady sink indefinitely? Solution: Reversible binding only transfers antibody between free and bound pools. Without a removal process, it does not continually eliminate antibody from the tissue; it can retard or redistribute antibody transiently but does not provide the sustained sink assumed in the steady model.
- For \(D=8.00\ \mu\mathrm{m}^2/\mathrm{s}\) and \(k_s=0.0800\ \mathrm{s}^{-1}\), calculate \(\lambda\). Solution: \(\lambda=\sqrt{D/k_s}=\sqrt{8.00/0.0800}=\sqrt{100}=10.0\ \mu\mathrm{m}\).