Tutorials › Quantitative Systems Pharmacology (QSP) › Resistance From Altered Target Sensitivity

Oncology: Tumor Growth and Therapy Response · Tutorial 263 of 799

Resistance From Altered Target Sensitivity

Compare potency shifts and reduced maximal drug effect in tumor-growth models, and see why these mechanisms need not predict the same response as a distinct resistant clone.

Advanced 12 min read

What You'll Learn

  • Represent altered target sensitivity as a change in EC50 or maximal drug effect.
  • Calculate drug effect and net tumor growth at a specified concentration.
  • Derive the concentration required to achieve a target effect and identify unattainable targets.
  • Compare a single altered-sensitivity population with a mixture of sensitive and resistant clones.
  • Recognize assumptions, diagnostic clues, and limitations of these model representations.

Representing Altered Target Sensitivity

A tumor may respond less to a drug even when the exposure is unchanged. One compact representation is to keep a single tumor population but alter the parameters linking concentration to effect. A shift in potency means that a higher concentration is needed to produce a given effect. A reduction in maximal effect means that even very high concentrations cannot produce the same maximum effect. These representations describe different changes in the concentration–effect relationship; neither, by itself, specifies the biological cause of resistance.

To make the distinction concrete, let \(C\) be the effective drug concentration at the target, measured in nM, and let \(E(C)\) be the drug-induced tumor-cell killing rate, measured in day\(^{-1}\). Consider a saturable Emax function:

$$ E(C)=\frac{E_{\max}C}{EC_{50}+C}. $$

Here \(E_{\max}\) is the maximum killing rate and \(EC_{50}\) is the concentration producing half of that maximum. A potency shift changes \(EC_{50}\); a reduction in maximal effect changes \(E_{\max}\). For a less-sensitive target, the potency shift is typically upward, so \(EC_{50,R}>EC_{50,S}\), where \(S\) and \(R\) denote sensitive and less-sensitive parameter values. Alternatively, a reduced maximum has \(E_{\max,R}<E_{\max,S}\).

Definition: In a single-population altered-sensitivity model, the tumor is represented by one burden \(N\), and resistance is represented by a change in the parameters of its drug-effect function. A potency shift changes \(EC_{50}\); an efficacy reduction changes \(E_{\max}\). These changes can be modeled separately or together, but their implications differ.

For the tumor dynamics, suppose untreated tumor growth follows a Gompertz law and drug-induced killing is proportional to tumor burden:

$$ \frac{dN}{dt} = rN\ln\!\left(\frac{K}{N}\right)-E(C)N. $$

In this expression, \(N\) is tumor burden, \(r\) is a growth-rate parameter in day\(^{-1}\), and \(K\) is the Gompertz asymptotic size. The per-capita growth rate at the current burden is \(r\ln(K/N)\), while \(E(C)\) is the per-capita killing rate. The net per-capita rate is therefore \(r\ln(K/N)-E(C)\): a negative value predicts shrinkage at that instant, and a positive value predicts growth. This is a local prediction; it does not guarantee the same sign at later times if burden, concentration, or model parameters change.

The model assumes that concentration is known and adequately represents exposure at the target, the effect is immediate, and the same tumor growth law applies before and during treatment. It also assumes a single effective population, with killing proportional to burden. These are modeling choices, not universal biological facts. The earlier material on exposure-effect links and tumor-growth laws supplies the relevant building blocks; here the focus is how altered sensitivity changes their predictions.

How Potency and Maximum Effect Change Predictions

For a potency shift, write the sensitive and less-sensitive effects as:

$$ E_S(C)=\frac{E_{\max}C}{EC_{50,S}+C}, \qquad E_R(C)=\frac{E_{\max}C}{EC_{50,R}+C}, \qquad EC_{50,R}>EC_{50,S}. $$

Both curves have the same upper limit, \(E_{\max}\), but the less-sensitive curve produces less effect at any finite positive concentration. The difference becomes small at concentrations far above both \(EC_{50}\) values, where both effects approach the common maximum. By contrast, if sensitivity changes through a reduced maximum, then increasing concentration cannot restore the lost effect.

Worked Example: A Potency Shift Changes the Net Growth Rate

Suppose a hypothetical tumor has a current Gompertz per-capita growth rate of \(0.2993\) day\(^{-1}\). At a constant concentration of \(2\) nM, let the sensitive parameters be \(E_{\max}=0.80\) day\(^{-1}\) and \(EC_{50,S}=1\) nM. Let the less-sensitive target have \(EC_{50,R}=4\) nM, with the same \(E_{\max}\). Calculate the drug effects and the resulting instantaneous net per-capita rates.

For the sensitive case:

$$ E_S(2)=\frac{(0.80\ \text{day}^{-1})(2\ \text{nM})}{1\ \text{nM}+2\ \text{nM}} =\frac{1.60}{3} =0.5333\ \text{day}^{-1}. $$

The net per-capita rate is \(0.2993-0.5333=-0.2340\) day\(^{-1}\), so the model predicts instantaneous tumor shrinkage.

For the less-sensitive target:

$$ E_R(2)=\frac{(0.80\ \text{day}^{-1})(2\ \text{nM})}{4\ \text{nM}+2\ \text{nM}} =\frac{1.60}{6} =0.2667\ \text{day}^{-1}. $$

The net rate is \(0.2993-0.2667=0.0326\) day\(^{-1}\), so the model predicts instantaneous growth. This difference follows from the higher \(EC_{50}\), not a change in the assumed growth rate or maximum killing rate. It does not mean that every tumor with this concentration will grow: the prediction depends on the specified parameters and current burden.

Deriving the Concentration Needed for a Target Effect

The Emax equation also makes the distinction useful for dose or exposure reasoning. If a desired killing rate \(E^*\) is below \(E_{\max}\), solve the Emax equation for the required concentration:

$$ E^*=\frac{E_{\max}C}{EC_{50}+C} \quad\Longrightarrow\quad E^*(EC_{50}+C)=E_{\max}C \quad\Longrightarrow\quad C=\frac{EC_{50}E^*}{E_{\max}-E^*}. $$

The algebra requires \(E_{\max}-E^*>0\). Thus, a finite positive concentration can achieve a positive target effect only when \(0<E^*<E_{\max}\). If \(E^*\geq E_{\max}\), the target is unattainable at finite concentration. Under a potency shift with unchanged \(E_{\max}\), the concentration required for any attainable target increases in direct proportion to \(EC_{50}\). Under a lower \(E_{\max}\), some targets that were previously attainable may become impossible.

Worked Example: Separating a Potency Shift From a Lower Maximum

Use a hypothetical sensitive profile with \(EC_{50}=1\) nM and \(E_{\max}=0.80\) day\(^{-1}\). Compare it with (a) a potency-shifted profile having \(EC_{50}=4\) nM and the same maximum, and (b) a reduced-maximum profile having \(EC_{50}=1\) nM and \(E_{\max}=0.40\) day\(^{-1}\). Find the concentration needed for \(E^*=0.30\) day\(^{-1}\).

For the sensitive profile:

$$ C_S=\frac{(1\ \text{nM})(0.30\ \text{day}^{-1})} {0.80\ \text{day}^{-1}-0.30\ \text{day}^{-1}} =\frac{0.30}{0.50}\ \text{nM} =0.600\ \text{nM}. $$

For the potency-shifted profile:

$$ C_R=\frac{(4\ \text{nM})(0.30\ \text{day}^{-1})} {0.80\ \text{day}^{-1}-0.30\ \text{day}^{-1}} =\frac{1.20}{0.50}\ \text{nM} =2.40\ \text{nM}. $$

The fourfold increase in \(EC_{50}\) produces a fourfold increase in the concentration needed for this same effect. For the reduced-maximum profile:

$$ C=\frac{(1\ \text{nM})(0.30\ \text{day}^{-1})} {0.40\ \text{day}^{-1}-0.30\ \text{day}^{-1}} =\frac{0.30}{0.10}\ \text{nM} =3.00\ \text{nM}. $$

This profile can still reach \(0.30\) day\(^{-1}\), but it requires a higher concentration because the target effect is closer to its lower maximum. A target of \(0.45\) day\(^{-1}\), however, exceeds its \(0.40\) day\(^{-1}\) maximum and cannot be achieved by increasing concentration in this model.

These calculations assume the concentration can be increased as specified and that the fitted effect function remains valid over that range. They do not establish that the required concentration is safe, clinically achievable, or appropriate. A concentration–effect model describes a response relationship; it is not, by itself, a dosing recommendation.

Why a Resistant Clone Is a Different Representation

A distinct resistant clone represents two populations with separate burdens, rather than one population with a shifted effect curve. To isolate the difference, suppose both clones share the same per-capita growth rate \(g(N)\), but have different killing rates \(E_S(C)\) and \(E_R(C)\). Their equations are:

$$ \frac{dN_S}{dt}=g(N)N_S-E_S(C)N_S, \qquad \frac{dN_R}{dt}=g(N)N_R-E_R(C)N_R, \qquad N=N_S+N_R. $$

Adding the equations gives the total-burden rate. If \(f=N_R/N\) is the resistant fraction, then \(1-f=N_S/N\), and the instantaneous total per-capita killing rate is:

$$ \frac{E_S(C)N_S+E_R(C)N_R}{N} =(1-f)E_S(C)+fE_R(C). $$

Thus, a clone mixture has a weighted effect at a particular concentration and composition. That weighted value can match a single altered-sensitivity model at one instant. But it need not match the same single Emax curve across concentrations, and the fraction \(f\) can change as the clones respond differently to treatment. A fixed shifted \(EC_{50}\) does not explicitly track that changing composition.

Worked Example: A Mixture Need Not Behave Like One Potency-Shifted Population

Consider a hypothetical fixed snapshot with equal sensitive and resistant fractions, \(f=0.50\). Both have \(E_{\max}=0.80\) day\(^{-1}\), but \(EC_{50,S}=1\) nM and \(EC_{50,R}=9\) nM. Calculate the mixture effect at \(C=1\) nM and \(C=9\) nM, then test whether one single-population curve with the same \(E_{\max}\) and a fixed \(EC_{50}\) can match both.

At \(C=1\) nM, the clone-specific effects are \(0.80(1)/(1+1)=0.400\) day\(^{-1}\) and \(0.80(1)/(9+1)=0.0800\) day\(^{-1}\). The mixture effect is:

$$ E_{\mathrm{mix}}(1) =(0.50)(0.400)+(0.50)(0.0800) =0.240\ \text{day}^{-1}. $$

A single curve with \(E_{\max}=0.80\) day\(^{-1}\) matching this point must have \(EC_{50}=0.80/0.240-1=2.333\) nM. At \(C=9\) nM, the mixture effect is:

$$ E_{\mathrm{mix}}(9) =(0.50)\frac{(0.80)(9)}{1+9} +(0.50)\frac{(0.80)(9)}{9+9} =0.560\ \text{day}^{-1}. $$

The single curve calibrated at \(1\) nM instead predicts \(0.80(9)/(2.333+9)=0.6353\) day\(^{-1}\). It does not match the mixture at \(9\) nM. This example does not prove that every mixture is distinguishable from every flexible single-population model: allowing other parameters to vary can improve a fit. It shows that a fixed-\(E_{\max}\), shifted-\(EC_{50}\) curve is not generally an exact substitute for explicitly representing two populations.

Model Checks, Comparisons, and Limitations

The choice between altered sensitivity and a distinct clone depends on the question and the available data. A single altered-sensitivity model is compact when the goal is to represent a changed concentration–effect relationship without tracking subpopulations. An explicit clone model is useful when clone fractions, clone-specific responses, or treatment-driven changes in composition are central to the prediction. The two representations can yield similar predictions over a limited exposure range, so a good fit alone does not necessarily identify the underlying mechanism.

FeatureAltered sensitivity in one populationDistinct resistant clone
Population structureOne total burden \(N\)Separate burdens \(N_S\) and \(N_R\)
How reduced response is representedChanged \(EC_{50}\), \(E_{\max}\), or bothDifferent clone-specific effects and abundances
Effect of changing clone fractionsNot represented explicitlyCan change the weighted total response
Data and parameters neededConcentration–response information and growth dataThose data plus information or assumptions about clone composition and dynamics

Useful checks include plotting observed or simulated effect against concentration, checking whether the response approaches a plausible plateau, and verifying that units are consistent. If an estimated \(EC_{50}\) is far outside the observed concentration range, potency may be poorly identified. If data show a plateau below the expected maximum, a reduced \(E_{\max}\) may be a candidate explanation, but incomplete exposure coverage or other model misspecification can produce similar patterns. If the apparent concentration–effect relationship changes over time, a constant single-population parameter may fail to describe the data; a changing clone composition is one possible explanation, not the only one.

A further diagnostic is to compare predictions at multiple concentrations, times, and tumor burdens rather than judging the model from one endpoint. A model can match one measured tumor change while giving a different prediction under another exposure or after treatment continues. Inspect residual patterns and parameter uncertainty, and check whether the estimated maximum is supported by observations near saturation. Where the data do not cover a broad concentration range or do not distinguish populations, report that limitation rather than assigning a unique mechanism.

  • Mistake: treating potency and maximum effect as interchangeable. Both can lower effect at a particular concentration, but a potency shift preserves the upper limit whereas a reduced maximum changes it.
  • Mistake: inferring shrinkage from drug effect alone. Compare \(E(C)\) with the current per-capita growth rate; tumor burden shrinks instantaneously only if \(E(C)>r\ln(K/N)\) in the stated model.
  • Mistake: interpreting a fitted shifted curve as proof of a resistant clone or a target alteration. A concentration–effect fit describes model behavior; biological mechanism requires supporting evidence.
  • Mistake: assuming greater concentration can always overcome resistance. That may follow from a potency shift within the model, but it does not follow when the target effect is at or above a reduced \(E_{\max}\), and model calculations do not address safety.
Key takeaway: Altered target sensitivity changes a concentration–effect relationship: a higher \(EC_{50}\) shifts potency, while a lower \(E_{\max}\) limits attainable effect. A resistant clone instead adds a distinct population whose contribution to total response can change with its fraction.

Check Your Understanding

Use the definitions and equations above. Each solution is included after its question.

  1. A hypothetical drug has \(E_{\max}=0.60\) day\(^{-1}\), \(EC_{50}=2\) nM, and \(C=2\) nM. What is \(E(C)\)?
    Solution: \(E(C)=0.60(2)/(2+2)=0.300\) day\(^{-1}\).
  2. If \(E_{\max}=0.50\) day\(^{-1}\), can this Emax model reach a target effect of \(0.50\) day\(^{-1}\) at a finite concentration? Explain.
    Solution: No. The equation approaches \(E_{\max}=0.50\) day\(^{-1}\) only as concentration tends to infinity; no finite concentration reaches or exceeds the maximum.
  3. At a current per-capita growth rate of \(0.25\) day\(^{-1}\), an effect of \(0.20\) day\(^{-1}\) predicts instantaneous shrinkage or growth? What about an effect of \(0.30\) day\(^{-1}\)?
    Solution: The net rates are \(0.25-0.20=0.05\) day\(^{-1}\), predicting growth, and \(0.25-0.30=-0.05\) day\(^{-1}\), predicting shrinkage.
  4. For a mixture with \(f=0.20\), \(E_S=0.40\) day\(^{-1}\), and \(E_R=0.10\) day\(^{-1}\), calculate the instantaneous weighted killing rate.
    Solution: \((1-f)E_S+fE_R=(0.80)(0.40)+(0.20)(0.10)=0.34\) day\(^{-1}\).
  5. State one reason a single shifted-\(EC_{50}\) model may fail to predict a resistant-clone mixture across time.
    Solution: The resistant fraction can change under treatment, changing the mixture’s weighted effect, while a fixed single-population \(EC_{50}\) does not explicitly track that changing composition.
Clinical Biostats