Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Soluble and Membrane Target Pools in Disposition Models

Target-Mediated Drug Disposition and PK/PD · Tutorial 259 of 598

Soluble and Membrane Target Pools in Disposition Models

Learn to represent soluble and membrane target binding in separate compartments and determine how accessibility and complex internalization contribute to drug disposition.

Advanced 11 min read

What You'll Learn

  • Define separate free-target and drug–target complex states for soluble and membrane-associated targets
  • Write mass-balance equations that distinguish plasma and tissue drug concentrations
  • Interpret accessibility factors without double-counting limited tissue exposure
  • Derive how binding and complex internalization produce target-mediated drug loss
  • Compare location-specific models with a single lumped target pool and identify their limitations
  • Diagnose when the model's assumed binding, exchange, or target-recovery pathways are not supported

Why Target Location Belongs in the Disposition Model

A drug can encounter a soluble target in plasma and a membrane-associated target in tissue, but those targets need not be equally accessible. The free drug concentration at a tissue surface may differ from the plasma concentration, and binding can lead to different complex fates in the two locations. Treating both pools as one target may therefore obscure which binding events remove drug from circulation and which chiefly alter local exposure.

The full TMDD framework distinguishes free drug, free target, and drug–target complex states, with binding, dissociation, target turnover, and complex loss represented as rate processes. Here the same bookkeeping is applied to two locations: a plasma-accessible soluble pool and a tissue-accessible membrane pool. The model is mechanistic only to the extent that its locations and pathways correspond to the system being studied.

Definition: A soluble target pool is free target present in a fluid compartment, represented here in plasma. A membrane target pool is free target associated with a cell surface, represented here in a tissue compartment. Each pool has its own free-target state and drug–target complex state. Accessibility describes whether, and at what local free-drug concentration, drug can reach a target pool and bind it.

A Two-Location State Model

Let \(C_p\) be free drug concentration in plasma and \(C_t\) free drug concentration in the tissue region containing the membrane target. Let \(S\) and \(R\) be, respectively, free soluble-target concentration in plasma and free membrane-target concentration in tissue. Let \(C_s\) and \(C_m\) denote soluble and membrane drug–target complex concentrations. Concentrations of \(S\) and \(C_s\) use the plasma volume \(V_p\); concentrations of \(R\) and \(C_m\) use the tissue volume \(V_t\). Drug input \(R_{\mathrm{in}}\) is an amount per time, \(CL_p\) and \(CL_t\) are free-drug clearances, and \(Q\) is an intercompartmental exchange clearance.

For location \(j\), \(k_{\mathrm{on},j}^{\mathrm{eff}}\) is the effective association-rate constant, \(k_{\mathrm{off},j}\) is the dissociation-rate constant, and \(k_{\mathrm{int},j}\) is the rate constant for complex internalization or another complex-loss process that removes drug from the modeled free and bound states. Target synthesis and free-target turnover are represented by \(k_{\mathrm{syn},j}\) and \(k_{\mathrm{deg},j}\). A fraction \(\rho_j\) of target removed with complex loss may return promptly to the free-target pool; use \(\rho_j=0\) when no such recycling is represented. The effective association rate can include an accessibility factor, \(k_{\mathrm{on},j}^{\mathrm{eff}}=a_j k_{\mathrm{on},j}\), where \(0\leq a_j\leq1\). If local drug concentration already captures limited access, set \(a_j=1\) rather than accounting for the same restriction twice.

$$ \begin{aligned} V_p\frac{dC_p}{dt} &=R_{\mathrm{in}}-CL_pC_p-Q(C_p-C_t) -V_p k_{\mathrm{on},s}^{\mathrm{eff}}C_pS+V_p k_{\mathrm{off},s}C_s,\\ V_t\frac{dC_t}{dt} &=Q(C_p-C_t)-CL_tC_t -V_t k_{\mathrm{on},m}^{\mathrm{eff}}C_tR+V_t k_{\mathrm{off},m}C_m,\\ \frac{dC_s}{dt} &=k_{\mathrm{on},s}^{\mathrm{eff}}C_pS -(k_{\mathrm{off},s}+k_{\mathrm{int},s})C_s,\\ \frac{dC_m}{dt} &=k_{\mathrm{on},m}^{\mathrm{eff}}C_tR -(k_{\mathrm{off},m}+k_{\mathrm{int},m})C_m,\\ \frac{dS}{dt} &=k_{\mathrm{syn},s}-k_{\mathrm{deg},s}S -k_{\mathrm{on},s}^{\mathrm{eff}}C_pS +k_{\mathrm{off},s}C_s+\rho_s k_{\mathrm{int},s}C_s,\\ \frac{dR}{dt} &=k_{\mathrm{syn},m}-k_{\mathrm{deg},m}R -k_{\mathrm{on},m}^{\mathrm{eff}}C_tR +k_{\mathrm{off},m}C_m+\rho_m k_{\mathrm{int},m}C_m. \end{aligned} $$

Each association term removes free drug and free target while forming complex in the same location. Dissociation restores both. Complex loss removes the complex; only the fraction \(\rho_j\) of the target is returned directly to the free-target state. The equations do not imply that all complex loss is literally internalization: the chosen mechanism must match the available biological evidence.

The drug exchange term \(Q(C_p-C_t)\) moves free drug between locations and cancels when the two drug amounts are summed. This is a simplified representation of distribution. A spatially heterogeneous tissue, intracellular transport, or explicit surface-binding region may require additional states rather than one effective tissue concentration.

Worked Example: Comparing Initial Binding Rates by Location

Consider hypothetical initial conditions \(C_p=2.0\) nM, \(C_t=0.50\) nM, \(S=10\) nM, and \(R=40\) nM, with no complexes initially. Suppose \(k_{\mathrm{on},s}^{\mathrm{eff}}=0.10\) nM\(^{-1}\)day\(^{-1}\) and \(k_{\mathrm{on},m}^{\mathrm{eff}}=0.050\) nM\(^{-1}\)day\(^{-1}\). Calculate the initial complex-formation rates per unit volume.

For the soluble pool, the initial association rate is \(k_{\mathrm{on},s}^{\mathrm{eff}}C_pS\). Substitution gives \(0.10\times2.0\times10=2.0\) nM/day.

For the membrane pool, the rate is \(k_{\mathrm{on},m}^{\mathrm{eff}}C_tR\). Substitution gives \(0.050\times0.50\times40=1.0\) nM/day.

The soluble pool initially forms complex at twice the rate per unit volume in this hypothetical example, despite the larger membrane-target concentration. The relevant drug concentration for membrane binding is \(C_t\), not \(C_p\); the target concentration alone does not determine the binding rate. Total complex formation in amount per time would also depend on the corresponding compartment volumes.

How Complex Formation Becomes Drug Loss

Binding by itself is not necessarily irreversible clearance. A drug molecule in a complex can dissociate and return to the free pool. If complex loss competes with dissociation, only some binding events result in drug removal through that route. This distinction can be made precise using the complex equation.

$$ \frac{dC_j}{dt} =k_{\mathrm{on},j}^{\mathrm{eff}}D_jT_j -(k_{\mathrm{off},j}+k_{\mathrm{int},j})C_j $$

Here \(D_j\) is the local free-drug concentration and \(T_j\) the corresponding free-target concentration: \((D_s,T_s)=(C_p,S)\) and \((D_m,T_m)=(C_t,R)\). If \(D_j\) and \(T_j\) are approximately constant over the interval of interest and the complex rapidly reaches a quasi-steady value, set \(dC_j/dt\approx0\). Rearranging gives \(C_j\approx k_{\mathrm{on},j}^{\mathrm{eff}}D_jT_j/(k_{\mathrm{off},j}+k_{\mathrm{int},j})\). The drug-loss rate per unit volume is then \(k_{\mathrm{int},j}C_j\), or equivalently \(k_{\mathrm{on},j}^{\mathrm{eff}}D_jT_j [k_{\mathrm{int},j}/(k_{\mathrm{off},j}+k_{\mathrm{int},j})]\).

The ratio \(k_{\mathrm{int},j}/(k_{\mathrm{off},j}+k_{\mathrm{int},j})\) is the probability that complex loss occurs before dissociation when those two competing events have constant rates. The quasi-steady expression is not a general steady-state solution for the entire TMDD system: it holds only when the complex adjusts quickly enough and the local free-drug and target concentrations change slowly on that timescale.

Worked Example: Estimating Complex-Mediated Drug Loss

Use hypothetical local conditions for a soluble pool: \(C_p=1.0\) nM, \(S=5.0\) nM, \(k_{\mathrm{on},s}^{\mathrm{eff}}=0.20\) nM\(^{-1}\)day\(^{-1}\), \(k_{\mathrm{off},s}=0.30\) day\(^{-1}\), and \(k_{\mathrm{int},s}=0.20\) day\(^{-1}\). Assume the complex is at quasi-steady state and \(V_p=3.0\) L.

The association rate is \(0.20\times1.0\times5.0=1.0\) nM/day. The quasi-steady complex concentration is \(C_s=1.0/(0.30+0.20)=2.0\) nM. The fraction of complexes lost rather than dissociating is \(0.20/0.50=0.40\). Thus the drug-loss rate per unit volume is \(0.20\times2.0=0.40\) nM/day. Multiplying by volume gives \(0.40\) nmol/L/day \(\times3.0\) L \(=1.20\) nmol/day.

For a hypothetical membrane pool, suppose \(C_t=0.20\) nM, \(R=20\) nM, \(k_{\mathrm{on},m}^{\mathrm{eff}}=0.10\) nM\(^{-1}\)day\(^{-1}\), \(k_{\mathrm{off},m}=0.10\) day\(^{-1}\), \(k_{\mathrm{int},m}=0.40\) day\(^{-1}\), and \(V_t=10\) L. Association is \(0.10\times0.20\times20=0.40\) nM/day. The quasi-steady complex concentration is \(0.40/(0.10+0.40)=0.80\) nM. Drug loss is \(0.40\times0.80=0.32\) nM/day, or \(0.32\times10=3.20\) nmol/day.

Under these stated approximations, the membrane route removes more drug per day in total, although its local free-drug concentration is lower. The comparison depends on target abundance, binding and complex-loss rates, and volume, not simply on whether the target is soluble or membrane-associated.

Mass-balance check: Define modeled drug amount as \(A=V_p(C_p+C_s)+V_t(C_t+C_m)\). Adding the four drug and complex equations cancels association, dissociation, and intercompartmental exchange terms, giving \(\frac{dA}{dt}=R_{\mathrm{in}}-CL_pC_p-CL_tC_t-V_pk_{\mathrm{int},s}C_s-V_tk_{\mathrm{int},m}C_m\). Thus, in this model, complex-mediated drug loss enters through complex loss, not through the binding step itself.

Accessibility Can Change the Apparent Binding Contribution

A membrane target may be abundant but poorly exposed to drug. One way to represent this is an explicit tissue drug concentration determined by distribution. Another is an effective association rate that summarizes limited access not otherwise represented. These approaches describe different levels of detail: if tissue distribution already lowers \(C_t\), applying a second accessibility factor for that same limitation exaggerates the reduction in binding.

Worked Example: Tissue Exposure Versus Plasma Exposure

Suppose tissue free drug is governed, in a simplified no-binding steady state, by exchange from plasma and tissue free-drug clearance: \(0=Q(C_p-C_t)-CL_tC_t\). Let \(Q=2.0\) L/day and \(CL_t=6.0\) L/day. Solving gives \(C_t=Q C_p/(Q+CL_t)=2.0C_p/8.0=0.25C_p\).

For a hypothetical \(C_p=4.0\) nM, this predicts \(C_t=1.0\) nM under the stated no-binding approximation. If \(R=10\) nM and \(k_{\mathrm{on},m}^{\mathrm{eff}}=0.20\) nM\(^{-1}\)day\(^{-1}\), membrane complex formation is \(0.20\times1.0\times10=2.0\) nM/day. Using plasma concentration in place of tissue concentration would predict \(0.20\times4.0\times10=8.0\) nM/day, four times as large.

This is a distribution illustration, not a full TMDD prediction: adding target binding can change \(C_t\), and the tissue concentration may not be at steady state. The calculation shows why the model should use the drug concentration at the target location.

Model Choice, Diagnostics, and Limitations

A single lumped target pool can be a useful simpler alternative when the two locations have similar exposure and kinetics, or when the data cannot distinguish them. It uses fewer parameters and states, but cannot separately attribute binding and complex loss to soluble and membrane targets. A two-location model can express different local concentrations, target turnover, and complex fates, but requires enough information to estimate or constrain those additional quantities.

RepresentationUseful whenMain limitation
Single lumped target poolLocations behave similarly or only aggregate target-mediated effects are supported by dataCannot assign disposition effects to a specific target location or pathway
Separate soluble and membrane poolsLocal exposure, target abundance, or complex fate plausibly differs by locationMore parameters and states can be weakly identifiable from plasma drug data alone
Explicit tissue distribution or additional intracellular statesTransport delays, gradients, or internalized material matter to the questionRequires additional measurements and assumptions about transport or intracellular fate

Useful checks follow directly from the model structure. Confirm that each concentration is paired with the correct volume and that \(k_{\mathrm{on}}DT\) has units of concentration per time. Check that all rate constants and clearances are nonnegative and that the chosen initial conditions do not imply negative states. Compare predicted tissue exposure and soluble-target measurements with relevant observations where available. If the model requires an extreme accessibility factor, a large unexplained target pool, or compensating changes in association and internalization rates, the assumed location structure may be unsupported or insufficiently informed.

A common mistake is to count complex formation itself as irreversible clearance. That reasoning treats every binding event as drug loss, even though dissociation can restore free drug. The correction is to track complex as a state and count drug removal only through the modeled complex-loss pathway (or another explicitly defined removal process). Another frequent error is to apply both a low tissue concentration and an accessibility multiplier for the same barrier. Use one representation for each limitation unless the factors describe distinct mechanisms.

Separate pools do not automatically establish that a target is physically accessible, that internalization destroys the drug, or that a recycling fraction is known. Those are mechanistic assumptions, not consequences of fitting equations. If the data do not distinguish soluble from membrane pathways, the more detailed model may produce uncertain parameter estimates even when its aggregate predictions appear adequate.

Check Your Understanding

Answer each question using the stated model assumptions.

  1. If \(C_t\) is one quarter of \(C_p\), what should be used in the membrane-target association term?
  2. For a complex with \(k_{\mathrm{off}}=0.30\) day\(^{-1}\) and \(k_{\mathrm{int}}=0.20\) day\(^{-1}\), what fraction of complexes are expected to undergo complex loss before dissociation?
  3. At quasi-steady complex concentration, explain why binding alone does not imply net drug elimination.
  4. Give one reason to prefer a single lumped target pool over separate soluble and membrane pools.

Solutions. (1) Use \(C_t\), because it is the local free-drug concentration at the membrane target; substituting \(C_p\) would ignore the stated tissue exposure difference. (2) The fraction is \(k_{\mathrm{int}}/(k_{\mathrm{off}}+k_{\mathrm{int}})=0.20/(0.30+0.20)=0.40\), or 40%. (3) Association transfers drug from the free state to the complex state, and dissociation returns it; the two terms cancel in total drug amount. Net loss occurs through complex loss or another explicitly modeled removal route. (4) A lumped pool may be preferable when the locations have similar kinetics or the data cannot distinguish their contributions; it also requires fewer parameters.)

Clinical Biostats