Why pH-Dependent Binding Matters
Immunoglobulin G (IgG) can bind the neonatal Fc receptor, FcRn, in acidic endosomal conditions, while binding is much weaker near neutral pH. This pH dependence provides a molecular basis for capturing IgG in an acidic environment and releasing it when the environment becomes less acidic. A PBPK model needs to represent both sides of that behavior: binding cannot be treated as one fixed affinity if the relevant pH changes between compartments or over time.
The purpose here is to formulate the binding and dissociation process, not to derive the subsequent tissue-scale recycling fluxes. The essential distinction is between a kinetic description, which specifies association and dissociation rates, and an equilibrium description, which summarizes their ratio at a fixed pH. A model that uses only an equilibrium affinity can describe equilibrium occupancy, but it cannot by itself specify how quickly binding adjusts after pH changes.
Kinetic and Equilibrium Descriptions
For this model, the complex concentration changes according to the difference between the association and dissociation rates:
Here concentrations may be expressed, for example, in mol/L. The units of \(k_{\mathrm{on}}\) are concentration\(^{-1}\) time\(^{-1}\), whereas \(k_{\mathrm{off}}\) has units time\(^{-1}\). At equilibrium at a fixed pH, the two rates are equal. The equilibrium dissociation constant is therefore
A smaller \(K_D\) means tighter binding at equilibrium; a larger \(K_D\) means weaker binding. For the capture-and-release behavior of interest, the model should represent stronger binding at acidic pH than at near-neutral pH, so \(K_D\) is lower in the acidic condition and higher in the near-neutral condition. Both rate constants can change with pH. Specifying a pH-dependent \(K_D\) alone does not determine either rate separately: many pairs of rates can have the same ratio but imply different binding and release times.
A common simplifying assumption is a single reversible binding site with one effective complex, well-mixed concentrations within the modeled compartment, and rates that respond to the local pH. This is an effective representation, not a claim that molecular binding involves only one ionizable group. Protonation of ionizable residues on IgG and FcRn can alter their interaction; a fitted pH-response curve summarizes the combined behavior over the conditions used to estimate it.
Worked Example: Equilibrium Occupancy at Acidic and Near-Neutral pH
Suppose a hypothetical system has a free IgG concentration of \(0.100\ \mu\mathrm{M}\). At pH 6.0, take \(K_D=0.0258\ \mu\mathrm{M}\); at pH 7.4, take \(K_D=9.984\ \mu\mathrm{M}\). Assume free IgG is maintained at the stated concentration and is not appreciably depleted by binding.
Let \(\theta\) be the fraction of FcRn sites occupied by IgG. With free IgG held constant, equilibrium gives \([C]/[R]=[G]/K_D\). Since the total receptor concentration is \([R]+[C]\), the occupied fraction is
At pH 6.0,
Thus about \(79.5\%\) of the modeled FcRn sites are occupied at equilibrium under these hypothetical conditions. At pH 7.4,
That is about \(0.992\%\) occupancy. The contrast follows from the stipulated pH-dependent affinities; it does not establish that these particular values apply to a biological system. The calculation also relies on free IgG remaining at \(0.100\ \mu\mathrm{M}\). If receptor binding materially depletes free IgG, the simple occupancy expression is not sufficient.
Representing the pH Response
A PBPK model can represent pH dependence using measured or estimated values of \(k_{\mathrm{on}}(\mathrm{pH})\) and \(k_{\mathrm{off}}(\mathrm{pH})\), with \(K_D(\mathrm{pH})\) calculated from their ratio. If kinetic data are unavailable, it may instead use a fitted \(K_D(\mathrm{pH})\) curve for equilibrium calculations, while recognizing that this does not identify the kinetic time course.
For example, a smooth empirical curve can interpolate between a low acidic-pH dissociation constant and a high near-neutral-pH dissociation constant. One possible form uses a transition function \(q(\mathrm{pH})\), where \(q\) declines from near one to near zero as pH increases:
Here \(h\) controls the steepness and \(\mathrm{pH}_{50}\) locates the transition in this chosen parameterization. This is an empirical interpolation, not a direct count of protons involved in binding. Its parameters should be fitted to relevant data or treated as explicitly hypothetical. Extrapolating beyond the pH range used to support the curve is generally not justified without additional evidence.
Worked Example: A Smooth Interpolation Between pH Conditions
For illustration only, set \(K_{D,\min}=0.0100\ \mu\mathrm{M}\), \(K_{D,\max}=10.0\ \mu\mathrm{M}\), \(h=4\), and \(\mathrm{pH}_{50}=6.70\). At pH 6.00,
Therefore,
At pH 7.40,
so
These values match the hypothetical affinities used in the equilibrium occupancy example. The steepness here is a selected curve parameter, not evidence that the molecular response has a Hill coefficient of four or that four protons are involved.
Finite Ligand Concentration and Mass Balance
When both total receptor and total IgG are specified and neither is effectively in excess, binding changes the free concentrations of both. Let \(R_T=[R]+[C]\) and \(G_T=[G]+[C]\). Substituting \([R]=R_T-[C]\) and \([G]=G_T-[C]\) into the equilibrium definition gives
Rearranging produces a quadratic equation in \([C]\). The physically admissible solution is the smaller root, since a complex concentration cannot exceed either total concentration:
This expression assumes one-to-one binding, equilibrium at a fixed pH, and no additional binding partners or compartments in the mass balance. It is useful when the free-ligand assumption in the simpler occupancy formula does not hold.
Worked Example: Binding When IgG Is Not in Large Excess
Suppose a hypothetical endosomal binding pool has \(R_T=0.0800\ \mu\mathrm{M}\), \(G_T=0.100\ \mu\mathrm{M}\), and \(K_D=0.0258\ \mu\mathrm{M}\) at pH 6.0. Apply the finite-concentration expression:
The sum is \(0.2058\ \mu\mathrm{M}\), and the quantity under the square root is \(0.2058^2-0.0320=0.01035\ (\mu\mathrm{M})^2\). Thus,
The fraction of total receptor bound is \(0.0520/0.0800=0.650\), or about \(65.0\%\). This is lower than the \(79.5\%\) calculated when free IgG was fixed at \(0.100\ \mu\mathrm{M}\), because here \(0.100\ \mu\mathrm{M}\) is the total IgG concentration and binding depletes the free amount.
Why Binding Rates Matter During a pH Shift
An endosome changes its environment over time, so equilibrium occupancy is not necessarily reached instantly. If free IgG is held constant and total receptor is conserved, define \(\theta=[C]/R_T\). Then \([R]=R_T-[C]\), and the kinetic equation becomes a linear equation in \(\theta\). At constant pH its solution approaches the equilibrium value exponentially:
This follows by collecting the terms in the differential equation: \(d\theta/dt=k_{\mathrm{on}}[G]-(k_{\mathrm{on}}[G]+k_{\mathrm{off}})\theta\). Its steady value is \(k_{\mathrm{on}}[G]/\lambda\), which equals \([G]/([G]+K_D)\) because \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\). The characteristic adjustment time is \(1/\lambda\). This solution applies only while pH, free IgG, and the rate constants remain fixed.
Worked Example: Occupancy After a Shift to Near-Neutral pH
Use hypothetical kinetic constants consistent with the affinities above: at pH 6.0, \(k_{\mathrm{on}}=1.00\times10^5\ \mathrm{M}^{-1}\mathrm{s}^{-1}\) and \(k_{\mathrm{off}}=0.00258\ \mathrm{s}^{-1}\), giving \(K_D=2.58\times10^{-8}\ \mathrm{M}=0.0258\ \mu\mathrm{M}\). At pH 7.4, let \(k_{\mathrm{on}}=5.00\times10^4\ \mathrm{M}^{-1}\mathrm{s}^{-1}\) and \(k_{\mathrm{off}}=0.4992\ \mathrm{s}^{-1}\), giving \(K_D=9.984\times10^{-6}\ \mathrm{M}=9.984\ \mu\mathrm{M}\). Suppose free IgG remains \(0.100\ \mu\mathrm{M}\), and occupancy is \(0.800\) just before the pH shift.
At pH 7.4, \(k_{\mathrm{on}}=0.0500\ \mu\mathrm{M}^{-1}\mathrm{s}^{-1}\). Thus \(\lambda=(0.0500)(0.100)+0.4992=0.5042\ \mathrm{s}^{-1}\), and \(\theta_{\mathrm{eq}}=(0.0500)(0.100)/0.5042=0.009917\). Five seconds after the shift,
The modeled occupancy has fallen to about \(7.34\%\), but has not yet reached its near-neutral equilibrium value of about \(0.992\%\). This calculation assumes the pH shift changes the rates immediately and that free IgG remains fixed. If free IgG also changes, or if pH changes gradually, the fixed-condition exponential solution does not apply directly.
Assumptions, Diagnostics, and Alternatives
A pH-dependent binding model should be checked against the observations and the conditions it represents. If acidic and near-neutral data imply similar \(K_D\) values when the biology or experimental design supports a pH effect, inspect pH control, assay conditions, and the fitted rate constants. Check that each fitted \(K_D\) equals \(k_{\mathrm{off}}/k_{\mathrm{on}}\), that rates are positive, and that units are consistent. A fitted curve that changes sharply outside the measured range may be an interpolation artifact rather than a defensible prediction.
The one-site model can also fail if binding is heterogeneous, if multiple receptor or ligand states matter, if complexes have additional fates, or if the system is far from equilibrium while only equilibrium data are used. Time-course data can help distinguish a change in affinity from a change in the speed of reaching equilibrium. A model with separate protonation states may be appropriate when the data support those states and their transitions; it is more mechanistic but requires more parameters and evidence. A table of pH-specific rate constants is less smooth but can be preferable when measurements exist only at a few pH values. A single pH-independent \(K_D\) is simpler, but cannot represent different acidic capture and near-neutral release behavior.
The introductory material on antibody amount balance and tissue exchange provides the broader compartment context. Here the binding model supplies a pH-sensitive interaction; a whole-body PBPK model must still specify where that interaction occurs and how the bound and unbound forms are represented. The key limitation is that a good fit to equilibrium affinity alone does not validate the kinetic rates or the assumed molecular mechanism.
Check Your Understanding
Use the definitions and assumptions in this tutorial to answer each question. Solutions follow the questions.
- If \(k_{\mathrm{on}}=2.0\times10^5\ \mathrm{M}^{-1}\mathrm{s}^{-1}\) and \(k_{\mathrm{off}}=0.020\ \mathrm{s}^{-1}\) at a specified pH, what is \(K_D\) in \(\mu\mathrm{M}\)?
- With free IgG \(0.050\ \mu\mathrm{M}\) and \(K_D=0.050\ \mu\mathrm{M}\), what fraction of receptor sites is occupied at equilibrium under the fixed-free-IgG assumption?
- Why can two models with the same \(K_D\) predict different occupancy time courses?
- When is the finite-concentration binding equation preferable to \(\theta=[G]/([G]+K_D)\)?
- Does an empirical pH-response curve with steepness parameter \(h=4\) establish that four protons participate in binding? Explain.
Solutions. (1) \(K_D=0.020/(2.0\times10^5)=1.0\times10^{-7}\ \mathrm{M}=0.100\ \mu\mathrm{M}\). (2) \(\theta=0.050/(0.050+0.050)=0.500\), or \(50.0\%\). (3) The ratio \(k_{\mathrm{off}}/k_{\mathrm{on}}\) can be equal for different pairs of rates, while the adjustment rate \(\lambda=k_{\mathrm{on}}[G]+k_{\mathrm{off}}\) differs. (4) Use it when total receptor and total IgG are specified and binding appreciably depletes free IgG; the simple expression assumes free IgG is fixed. (5) No. In this formulation \(h\) controls the empirical curve’s steepness; it is not a direct count of participating protons.