Why Plasma Binding Can Depend on Concentration
A fraction unbound measured at one drug concentration is not necessarily a fixed property of a compound. If binding sites have finite capacity, increasing the drug concentration fills a larger share of those sites. The bound amount then approaches a maximum, and a progressively larger share of additional drug remains unbound. For PBPK parameterization, this means that a single fraction-unbound value may not describe the whole concentration range of interest.
The model below describes reversible binding to one class of identical sites at equilibrium. It separates the free drug concentration from the drug bound to protein and uses mass balance to connect them to the measured total plasma concentration. As in the earlier material on recovery and artifacts in plasma-binding assays, the concentration values used for interpretation need to represent the intended equilibrium system; apparent concentration dependence caused by assay problems is not evidence of saturation.
Deriving the Saturable-Binding Relationship
The model assumes that free drug binds reversibly to a fixed population of sites, that binding has reached equilibrium, and that all sites in the population have the same affinity. The plasma conditions and binding-protein abundance are treated as constant over the concentration range considered. There is no competition from other ligands, and there are no additional binding processes in the model. The bound concentration is given by the one-site equilibrium expression:
This expression is approximately proportional to \(C_u\) when \(C_u\) is small relative to \(K_d\), but it approaches \(B_{max}\) as \(C_u\) becomes large. Substituting it into the mass balance gives the total concentration as a function of the unbound concentration:
To calculate \(C_u\) from an observed or specified \(C_{tot}\), multiply both sides by \(K_d+C_u\) and collect terms. This produces a quadratic equation:
For positive \(C_{tot}\), the constant term is negative, so the two roots have opposite signs. Only the positive root is physically meaningful. Thus the unbound concentration and fraction unbound can be calculated as:
This solution gives \(C_u\) in the same concentration units as \(C_{tot}\), \(K_d\), and \(B_{max}\). If concentrations are expressed in mg/L, for example, all three model parameters must also use mg/L. \(B_{max}\) represents the binding capacity in drug-equivalent concentration units, not necessarily the molar concentration of the protein itself.
Worked Example: Calculating Fraction Unbound at One Total Concentration
Hypothetical data: Suppose a plasma-binding system has \(K_d=2.00\) mg/L and \(B_{max}=18.0\) mg/L. Calculate \(C_u\) and \(f_u\) when \(C_{tot}=10.0\) mg/L, using the one-site model.
Substitute the values into the quadratic equation:
so \(C_u^2+10.0C_u-20.0=0\). The positive root is:
The resulting fraction unbound is \(f_u=1.708/10.0=0.1708\). As a mass-balance check, \(C_b=10.0-1.708=8.292\) mg/L. The binding equation also gives \(18.0(1.708)/(2.00+1.708)=8.292\) mg/L, agreeing with the mass balance. Under this model and at this total concentration, about 17.1% of the drug is unbound.
How Fraction Unbound Changes as Sites Fill
The model has two useful limits. At low unbound concentrations, the denominator \(K_d+C_u\) is approximately \(K_d\). The bound concentration is then nearly proportional to the free concentration, and the fraction unbound approaches a constant:
At high unbound concentrations, the binding capacity is nearly filled: \(C_b\) approaches \(B_{max}\). In that limit, \(C_{tot}\) is approximately \(C_u+B_{max}\), and \(f_u\) approaches 1 as total concentration increases. The fraction unbound therefore generally rises with concentration in this model. The change need not be large over every observed range; it depends on the relation between the concentrations and the binding parameters.
A short derivative makes the changing response explicit. Differentiating the mass-balance expression with respect to \(C_u\) gives \(dC_{tot}/dC_u=1+B_{max}K_d/(K_d+C_u)^2\). Reversing this derivative gives the local change in free concentration per unit change in total concentration:
For positive parameter values this slope is greater than 0 and less than 1. At low concentrations, binding absorbs a relatively large part of an increase in total concentration, so the free concentration rises more slowly. As sites approach saturation, the slope moves toward 1: less of an additional increase in total concentration can be accommodated by additional binding. The fraction unbound can rise even though the local absolute increase in \(C_u\) per unit increase in \(C_{tot}\) remains at most 1.
Worked Example: Comparing Free Concentration Before and Near Saturation
Hypothetical data: Keep \(K_d=2.00\) mg/L and \(B_{max}=18.0\) mg/L. Compare the free concentrations at total concentrations of 10.0 and 30.0 mg/L.
At 10.0 mg/L, the previous calculation gave \(C_u=1.708\) mg/L and \(f_u=0.1708\). At 30.0 mg/L, substitution gives:
or \(C_u^2-10.0C_u-60.0=0\). Taking the positive root:
Therefore \(f_u=14.22/30.0=0.4740\). The bound concentration is \(30.0-14.22=15.78\) mg/L, or \(15.78/18.0=87.7\%\) of the available capacity. Tripling total concentration from 10.0 to 30.0 mg/L increased the modeled free concentration from 1.708 to 14.22 mg/L. This is a nonlinear response: the binding sites are much closer to full occupancy at the higher concentration.
Worked Example: Using the Local Slope Near Half-Occupancy
Hypothetical data: Use the same binding parameters, and start at \(C_u=2.00\) mg/L. Calculate the corresponding total concentration and estimate the effect of a small increase in total concentration. Then compare the estimate with an exact calculation for an increase of 1.00 mg/L.
At \(C_u=2.00\) mg/L, the bound concentration is \(18.0(2.00)/(2.00+2.00)=9.00\) mg/L, so \(C_{tot}=2.00+9.00=11.0\) mg/L. The local slope is:
Near this starting point, a small increase of 1.00 mg/L in total concentration is therefore estimated to increase \(C_u\) by about \(0.3077\) mg/L. For the exact result at \(C_{tot}=12.0\) mg/L, solve:
The positive root is \(C_u=(-8.00+\sqrt{160})/2=2.325\) mg/L. The exact increase is \(2.325-2.00=0.325\) mg/L, close to the local estimate. The derivative is a local approximation; over a larger concentration change, the slope itself changes as binding becomes more saturated.
Assumptions, Alternatives, and Failure Modes
The model is useful when equilibrium binding is reversible and a single effective site class adequately describes the data. Its assumptions should be considered before using the quadratic solution. Changes in protein abundance, pH, competing ligands, or other binding mechanisms can alter the observed relationship. A measurement made before equilibrium, or one affected by the assay artifacts addressed in the earlier material, can also produce apparent concentration dependence that this equilibrium model does not explain.
A constant \(f_u\) is a simpler alternative. It may be a reasonable approximation over a concentration range where binding is effectively linear, such as when \(C_u\) is small relative to \(K_d\). It is not a general replacement for the saturable model across a range that approaches site capacity. Conversely, a one-site model may be too simple if different proteins or site classes contribute substantially. In that case, a model with multiple binding terms can represent more than one affinity or capacity, but it requires enough informative data to estimate the additional parameters.
Useful checks include verifying that predicted concentrations are nonnegative, that \(C_u\leq C_{tot}\), and that \(C_b=C_{tot}-C_u\) agrees with the binding equation. The bound concentration should approach, but not exceed, \(B_{max}\). A systematic change in measured \(f_u\) with concentration may be consistent with saturation, but by itself does not establish that the one-site mechanism is correct. Large inconsistencies with the model, or parameter estimates that change markedly when the sampled concentration range changes, call for examining the data range, measurement quality, and model adequacy.
- Pitfall: “The fraction unbound was 0.17 at one concentration, so multiplying every total concentration by 0.17 gives the free concentration.” This assumes that \(f_u\) stays constant. Under saturable binding, that fraction generally increases as total concentration increases.
- Correction: Use the one-site mass-balance equation and solve for \(C_u\) at each relevant \(C_{tot}\), or justify a constant-\(f_u\) approximation over the concentration range of interest.
- Pitfall: Treating \(B_{max}\) as an unlimited proportional binding coefficient. It is a capacity; the model predicts \(C_b\) cannot exceed it.
Check Your Understanding
Use the one-site equilibrium model unless a question says otherwise. Concentrations are in consistent units.
- For low \(C_u\), what value does \(f_u\) approach in terms of \(K_d\) and \(B_{max}\)?
- With \(K_d=2.00\) mg/L and \(B_{max}=18.0\) mg/L, the model gives \(C_u=2.00\) mg/L. Calculate \(C_{tot}\) and \(f_u\).
- What happens to \(C_b\) and \(f_u\) as \(C_u\) becomes very large?
- Why is a constant fraction unbound measured at low concentration not necessarily suitable at a much higher concentration?
Solutions. (1) At low \(C_u\), \(C_b\approx (B_{max}/K_d)C_u\), so \(C_{tot}\approx C_u(1+B_{max}/K_d)\) and \(f_u\approx K_d/(K_d+B_{max})\). (2) \(C_b=18.0(2.00)/(2.00+2.00)=9.00\) mg/L, so \(C_{tot}=2.00+9.00=11.0\) mg/L and \(f_u=2.00/11.0=0.1818\). (3) \(C_b\) approaches \(B_{max}\), while \(f_u\) approaches 1 as total concentration increases. (4) At higher concentrations, binding sites may be more occupied, so the bound amount does not keep rising in proportion to free concentration; \(f_u\) can therefore increase.