Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Complex Internalization, Degradation, and Recycling

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

Complex Internalization, Degradation, and Recycling

Learn how to represent complex internalization with separate target-recycling, target-degradation, and intracellular drug-degradation processes, and how to check the resulting mass balances.

Advanced 11 min read

What You'll Learn

  • Write full TMDD equations with intracellular target and drug states
  • Distinguish target recycling from target degradation and intracellular drug degradation
  • Derive drug and target balance equations from the model
  • Calculate pathway fluxes and their effects on target recovery
  • Diagnose omitted terms, inconsistent units, and over-lumped internalization models

Why Internalization Needs More Than One Fate

When a drug binds its target, the resulting complex may be internalized. Internalization is not, by itself, a complete description of what happens next. The target may return to the cell surface, be degraded inside the cell, or remain temporarily in an intracellular pool; the drug may be degraded intracellularly or, in a different biological system, be returned to the extracellular space. These possibilities affect target abundance and drug disposition in different ways.

The model below makes those distinctions explicit. It extends the full TMDD mass-balance framework and target synthesis and turnover by adding intracellular target and drug states. Its central simplifying assumption is that an internalized complex promptly separates into target and drug, which then follow their own intracellular fates. This is a useful model structure, not a universal description of intracellular trafficking.

Definition: Let \(D\) be free drug concentration in the extracellular compartment, \(R\) be free surface-target concentration, and \(C\) be surface drug-target complex concentration. Let \(Q\) be intracellular target awaiting recycling or degradation, and \(I\) be intracellular drug awaiting degradation. The model includes target synthesis at rate \(k_{\mathrm{syn}}\), free-target degradation at rate \(k_R R\), reversible binding with association rate \(k_{\mathrm{on}}DR\) and dissociation rate \(k_{\mathrm{off}}C\), and complex internalization at rate \(k_{\mathrm{int}}C\). Internalized target recycles at rate \(k_{\mathrm{rec}}Q\) or degrades at rate \(k_{\mathrm{deg},T}Q\); intracellular drug degrades at rate \(k_{\mathrm{deg},D}I\).

All rate constants are nonnegative. Concentrations may be expressed in nM and time in hours, for example; then first-order rate constants have units h\(^{-1}\), \(k_{\mathrm{on}}\) has units nM\(^{-1}\)h\(^{-1}\), and \(k_{\mathrm{syn}}\) has units nM/h. Let \(u(t)\) denote drug input into the extracellular compartment in concentration per unit time, and let \(k_{\mathrm{el}}D\) denote nonspecific elimination of free drug.

Equations and Their Assumptions

Applying mass action to binding and first-order rates to the other processes gives the following equations:

$$ \begin{aligned} \frac{dD}{dt} &= u(t)-k_{\mathrm{el}}D-k_{\mathrm{on}}DR+k_{\mathrm{off}}C,\\ \frac{dR}{dt} &= k_{\mathrm{syn}}-k_RR-k_{\mathrm{on}}DR+k_{\mathrm{off}}C+k_{\mathrm{rec}}Q,\\ \frac{dC}{dt} &= k_{\mathrm{on}}DR-(k_{\mathrm{off}}+k_{\mathrm{int}})C,\\ \frac{dQ}{dt} &= k_{\mathrm{int}}C-(k_{\mathrm{rec}}+k_{\mathrm{deg},T})Q,\\ \frac{dI}{dt} &= k_{\mathrm{int}}C-k_{\mathrm{deg},D}I. \end{aligned} $$

The model assumes well-mixed compartments, one-to-one drug-target binding, and constant rate parameters over the modeled period. It treats target synthesis as input to the surface free-target pool and free-target degradation as a first-order process. Internalization removes surface complex and supplies one target equivalent to \(Q\) and one drug equivalent to \(I\). Thus it assumes prompt intracellular separation for purposes of the modeled kinetics. The model does not include intracellular complex accumulation, drug recycling, target transport delays beyond \(Q\), or drug elimination from any compartment other than free-drug elimination and degradation of \(I\).

The absence of a direct \(k_{\mathrm{int}}C\) term in the \(D\) equation is deliberate. Internalization removes complex, not free extracellular drug. It affects free drug over time by changing how much complex remains available to dissociate and return drug to the free pool. Similarly, recycling adds target back to \(R\), but does not return drug under these assumptions.

Mass-Balance Derivation

A useful check is to sum the state variables that represent each moiety. Define the modeled target amount, in concentration equivalents, as \(M_T=R+C+Q\). Add the equations for \(R\), \(C\), and \(Q\). The association terms cancel the matching dissociation terms, internalization removes \(C\) while adding \(Q\), and recycling removes \(Q\) while adding \(R\). The result is:

$$ \begin{aligned} \frac{dM_T}{dt} &=\frac{dR}{dt}+\frac{dC}{dt}+\frac{dQ}{dt}\\ &=k_{\mathrm{syn}}-k_RR-k_{\mathrm{deg},T}Q. \end{aligned} $$

This is a stoichiometric balance: target can enter through synthesis and leave through the two target-degradation pathways, but binding, internalization, and recycling only move target between modeled states. The target pool need not be constant; it changes whenever synthesis and degradation fluxes do not balance.

Likewise define \(M_D=D+C+I\), the drug amount in concentration equivalents across the modeled drug states. Summing the corresponding equations gives:

$$ \frac{dM_D}{dt}=u(t)-k_{\mathrm{el}}D-k_{\mathrm{deg},D}I. $$

Binding and dissociation cancel because they redistribute drug between \(D\) and \(C\). Internalization transfers drug from \(C\) to \(I\), so it also cancels in the sum. Drug leaves the modeled system through free-drug elimination and intracellular degradation. These balances are strong implementation checks: if the summed equations retain a binding, internalization, or recycling term, at least one state equation is inconsistent with the stated reactions.

Worked Example: Checking the Target Balance

Consider a hypothetical instant with \(R=2.00\) nM and \(Q=0.400\) nM. Suppose \(k_{\mathrm{syn}}=0.100\) nM/h, \(k_R=0.0200\) h\(^{-1}\), \(k_{\mathrm{rec}}=0.300\) h\(^{-1}\), and \(k_{\mathrm{deg},T}=0.200\) h\(^{-1}\). Let association, dissociation, and internalization fluxes at this instant be \(3.00\), \(0.400\), and \(0.600\) nM/h, respectively. Calculate the derivatives of \(R\), \(C\), and \(Q\), then check their sum.

For free target, include synthesis, free-target degradation, binding, dissociation, and recycling:

$$ \frac{dR}{dt} =0.100-(0.0200)(2.00)-3.00+0.400+(0.300)(0.400) =-2.42\ \text{nM/h}. $$

For complex, association adds complex while dissociation and internalization remove it:

$$ \frac{dC}{dt}=3.00-0.400-0.600=2.00\ \text{nM/h}. $$

For intracellular target, internalization supplies \(Q\), while recycling and degradation remove it:

$$ \frac{dQ}{dt}=0.600-(0.300+0.200)(0.400)=0.400\ \text{nM/h}. $$

Therefore \(dM_T/dt=-2.42+2.00+0.400=-0.0200\) nM/h. The balance equation gives the same result: \(0.100-(0.0200)(2.00)-(0.200)(0.400)=-0.0200\) nM/h. Although complex and intracellular target are increasing at this instant, total modeled target is decreasing slightly. The quantity \(k_{\mathrm{syn}}-k_RR\) alone is not the full-model derivative of free target or total target.

What Recycling and Degradation Rates Mean

If \(Q\) is held at a constant input flux for long enough to reach its local steady state, and its only exits are recycling and degradation, then setting \(dQ/dt=0\) gives:

$$ Q_{\mathrm{ss}}=\frac{k_{\mathrm{int}}C}{k_{\mathrm{rec}}+k_{\mathrm{deg},T}}. $$

At that condition, the fraction of incoming internalized target flux that recycles is \(k_{\mathrm{rec}}/(k_{\mathrm{rec}}+k_{\mathrm{deg},T})\), and the fraction that degrades is \(k_{\mathrm{deg},T}/(k_{\mathrm{rec}}+k_{\mathrm{deg},T})\). These fractions describe competing first-order exits from \(Q\); they are not the fractions of all surface target that recycle or degrade. Their interpretation depends on the specified pathway structure and the assumption that the rates are constant.

Worked Example: Quantifying Target Recycling and Degradation

Suppose complex internalization supplies \(Q\) at a constant rate of \(0.400\) nM/h. Let \(k_{\mathrm{rec}}=0.300\) h\(^{-1}\) and \(k_{\mathrm{deg},T}=0.200\) h\(^{-1}\). Find the steady intracellular target concentration and the recycling and degradation fluxes.

The total exit rate from \(Q\) is \(0.300+0.200=0.500\) h\(^{-1}\). Thus:

$$ Q_{\mathrm{ss}}=\frac{0.400\ \text{nM/h}}{0.500\ \text{h}^{-1}}=0.800\ \text{nM}. $$

The recycling flux is \((0.300\ \text{h}^{-1})(0.800\ \text{nM})=0.240\) nM/h. The target-degradation flux is \((0.200\ \text{h}^{-1})(0.800\ \text{nM})=0.160\) nM/h. They sum to the \(0.400\) nM/h entering flux. Equivalently, 60.0% of this incoming target flux recycles and 40.0% degrades. This calculation concerns the local \(Q\) compartment; it does not establish that the whole target system is at steady state.

Worked Example: Intracellular Drug Degradation and Free Drug

At a hypothetical instant, let \(D=5.00\) nM, \(C=0.400\) nM, \(k_{\mathrm{el}}=0.100\) h\(^{-1}\), and \(u(t)=0\). Suppose the binding flux is \(k_{\mathrm{on}}DR=3.00\) nM/h and dissociation flux is \(k_{\mathrm{off}}C=0.400\) nM/h. Calculate the instantaneous free-drug derivative. Then compare complex derivatives for \(k_{\mathrm{int}}=0.500\) h\(^{-1}\) and \(0.100\) h\(^{-1}\).

The free-drug equation gives:

$$ \frac{dD}{dt}=0-(0.100)(5.00)-3.00+0.400=-3.10\ \text{nM/h}. $$

At the stated \(D\), \(R\), and \(C\), this instantaneous derivative does not contain \(k_{\mathrm{int}}\). With \(k_{\mathrm{int}}=0.500\) h\(^{-1}\), the internalization flux is \((0.500)(0.400)=0.200\) nM/h and \(dC/dt=3.00-0.400-0.200=2.40\) nM/h. With \(k_{\mathrm{int}}=0.100\) h\(^{-1}\), the internalization flux is \(0.0400\) nM/h and \(dC/dt=3.00-0.400-0.0400=2.56\) nM/h. The immediate free-drug derivative is the same in both calculations, but complex changes at different rates. Subsequently, that difference can change dissociation back to free drug and the supply of intracellular drug for degradation. A claim that internalization instantaneously removes free drug directly would not match these equations.

Model Choice, Diagnostics, and Limitations

A simpler alternative is to remove complex through a single effective loss term, such as \(k_{\mathrm{int}}C\), without tracking \(Q\) or \(I\). This can be adequate when the scientific question concerns only the disappearance of surface complex, or when intracellular fates are not measured and need not be distinguished. It cannot, by itself, predict how much target returns to the surface, how much target is degraded after internalization, or how intracellular drug accumulates and degrades. The explicit-compartment model is more mechanistic, but it also adds parameters and states that may be difficult to estimate from limited data.

Useful checks include verifying that all rates have compatible units, that all states remain nonnegative under nonnegative inputs and initial conditions, and that the target and drug balance equations reduce to the expected synthesis and elimination fluxes. Compare model predictions with data over time: persistent errors in target recovery after complex declines may indicate that a single irreversible complex-loss term is inadequate. A delay in target recovery may require an intracellular target pool such as \(Q\), while systematic drug-persistence errors may indicate that drug recycling, intracellular complex, or another drug compartment is missing. Such patterns are clues, not unique diagnoses; assay error and other forms of model misspecification can produce similar discrepancies.

A common mistake is to count complex internalization as target degradation and then also include degradation from \(Q\). That removes target twice: once when \(C\) is internalized and again when \(Q\) degrades. In the model here, internalization transfers target from \(C\) to \(Q\); only \(k_{\mathrm{deg},T}Q\) and \(k_RR\) remove target. A related mistake is to omit the recycling term from \(dR/dt\) while retaining it as a loss from \(Q\). Recycling is a transfer, so it must appear with opposite signs in the two equations and cancel from the total-target balance.

The model also has important limits. The separate \(Q\) and \(I\) pools assume rapid dissociation of internalized complex and do not represent intact intracellular complex. A biological system in which the drug is recycled, target returns with drug bound, or intracellular processing is saturable needs a different reaction scheme. Finally, these equations describe disposition and target turnover; they do not specify how target engagement produces a pharmacodynamic response.

Key takeaway: Internalization is a transfer out of surface complex, not automatically a loss of both drug and target from the modeled system. Explicit intracellular states let recycling, target degradation, and drug degradation be represented as distinct fluxes, and the summed mass balances provide a direct check on the equations.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. In the stated model, which processes change total modeled target \(R+C+Q\), and which only redistribute it?
  2. If \(k_{\mathrm{rec}}=0.600\) h\(^{-1}\) and \(k_{\mathrm{deg},T}=0.400\) h\(^{-1}\), what fraction of incoming \(Q\) flux recycles?
  3. Why is there no direct internalization term in \(dD/dt\)?
  4. For \(k_{\mathrm{int}}C=0.300\) nM/h and \(k_{\mathrm{deg},D}=0.150\) h\(^{-1}\), what is \(I_{\mathrm{ss}}\) if intracellular drug is at steady state?
  5. What is the mass-balance error if internalization is included as a target loss from \(C\) but no corresponding target input to \(Q\) is included?

Solutions. (1) Synthesis adds total target; free-target degradation and degradation from \(Q\) remove it. Binding, dissociation, internalization, and recycling redistribute target among \(R\), \(C\), and \(Q\). (2) The recycling fraction is \(0.600/(0.600+0.400)=0.600\), or 60.0%. (3) Internalization removes complex, whereas \(D\) is free extracellular drug; the drug moves from \(C\) to \(I\) under the model assumptions. (4) At steady state, \(0=k_{\mathrm{int}}C-k_{\mathrm{deg},D}I\), so \(I_{\mathrm{ss}}=(0.300\ \text{nM/h})/(0.150\ \text{h}^{-1})=2.00\) nM. (5) The summed target balance would incorrectly retain a loss term \(-k_{\mathrm{int}}C\), treating internalization as destruction rather than transfer. Add \(+k_{\mathrm{int}}C\) to \(dQ/dt\) to represent the target entering the intracellular pool.

Clinical Biostats