Tutorials › Pharmacometrics › PK/PD Models with Disease Progression
Pharmacokinetics · PK/PD Modeling

PK/PD Models with Disease Progression

Learn how PK/PD models can separate the natural course of disease from drug effects, represent disease progression over time, and connect drug exposure to changing biomarkers, symptoms, tumor burden, or other clinical outcomes.

Intermediate PK/PD Modeling Disease Progression Pharmacometrics
01 · The big picture

1. What Is a Disease-Progression PK/PD Model?

A conventional PK/PD model often describes how drug concentration changes over time and how that concentration produces a pharmacologic effect. In many diseases, however, the outcome being modeled is not static. The underlying disease itself may change with time, even in the absence of treatment.

A disease-progression PK/PD model explicitly represents this changing baseline and then describes how drug exposure modifies the disease trajectory.

Dose PK model C(t) Drug effect treatment signal Disease Progression natural history Outcome biomarker / symptom The model separates disease dynamics from treatment-related changes

Disease-progression PK/PD models combine pharmacokinetics, drug effects, and an explicit representation of the underlying disease trajectory.

Core idea: the observed clinical trajectory can reflect both what the disease would have done without treatment and how treatment changes that trajectory. Disease-progression models attempt to represent both components.
02 · Why progression matters

2. Why Is Disease Progression Important in PK/PD Modeling?

Many PD models assume that the baseline response is fixed or that the drug effect acts on an otherwise stable system. That assumption may be inappropriate when the disease changes continuously over the study period.

Examples include tumor growth, decline in pulmonary function, progression of neurodegenerative disease, worsening of a biomarker, or gradual changes in disease severity.

SituationPotential progression processWhy it matters
OncologyTumor growth or tumor burden changeA reduction in tumor size may reflect treatment opposing an underlying growth process.
Chronic diseaseGradual biomarker deteriorationTreatment may slow decline rather than produce an immediate improvement.
NeurodegenerationProgressive loss of functionA treatment effect may appear as a change in the slope of decline.
Respiratory diseaseTime-varying lung-function trajectoryTreatment may modify the rate of deterioration.
Metabolic diseaseChanging disease-state markerDrug effects may be superimposed on natural temporal changes.

Without a progression component, a model can incorrectly attribute changes caused by natural disease history to the drug—or overlook a treatment effect because the outcome is changing for reasons unrelated to treatment.

03 · Model structure

3. The Basic Structure of a Disease-Progression Model

A useful conceptual decomposition is:

\[ \text{Observed outcome} = \text{disease progression} + \text{drug effect} + \text{residual variability} \]

The exact mathematical form depends on the disease and endpoint. A simple linear progression model might begin with:

\[ D(t)=D_0+kt \]

Here, \(D(t)\) represents the disease state, \(D_0\) is the baseline disease state, and \(k\) is the natural progression rate.

If the drug changes the progression rate, one possible model is:

\[ \frac{dD}{dt}=k_{\text{prog}}-E(C(t)) \]

where \(E(C(t))\) is a concentration-dependent treatment effect.

Important distinction: the drug does not necessarily have to produce an immediate change in the disease marker. It may instead change the rate of progression, the rate of recovery, or another component of the disease dynamics.
04 · Natural history

4. Representing Natural Disease Progression

The first modeling task is often to describe what the disease would do over time in the absence of treatment.

Several functional forms are possible. A simple linear model is:

\[ D(t)=D_0+\text{Slope}\times t \]

For processes that accelerate or decelerate, nonlinear functions may be more appropriate. For example, exponential progression can be represented as:

\[ D(t)=D_0e^{kt} \]

Other disease-progression models may use turnover, logistic growth, indirect-response models, or mechanistic state-transition systems.

Progression modelGeneral formInterpretation
Linear\(D(t)=D_0+kt\)Constant absolute rate of progression.
Exponential\(D(t)=D_0e^{kt}\)Rate increases with the disease state.
Turnover\(\frac{dD}{dt}=k_{\text{in}}-k_{\text{out}}D\)Balance between production and loss.
LogisticGrowth slows as a limiting level is approached.Useful for bounded growth processes.
State-transitionMultiple disease states connected by transition rates.Represents movement between clinically meaningful states.

The progression model should be chosen according to the scientific behavior of the endpoint and the information contained in the data.

05 · Drug effect

5. Adding the Drug Effect

Once the underlying disease trajectory is represented, a drug effect can be incorporated into the model.

A common concentration-effect relationship is the \(E_{\max}\) model:

\[ E(C)=\frac{E_{\max}C}{EC_{50}+C} \]

Suppose disease progression is naturally worsening at rate \(k_{\text{prog}}\), while drug exposure reduces that rate. One possible model is:

\[ \frac{dD}{dt} = k_{\text{prog}} - \frac{E_{\max}C(t)}{EC_{50}+C(t)} \]

When concentration is high, the treatment effect approaches \(E_{\max}\). When concentration is low, the effect approaches zero.

This formulation illustrates an important feature of disease-progression PK/PD models: the pharmacologic effect can act on a disease process rather than simply shifting the observed response by a fixed amount.

06 · Slope effects

6. Drug Effects on the Rate of Disease Progression

For many chronic diseases, an especially useful question is whether treatment changes the slope of disease progression.

Consider a disease marker that naturally increases over time:

\[ D(t)=D_0+k_{\text{prog}}t \]

A treatment that slows progression could be represented as:

\[ D(t)=D_0+ \left( k_{\text{prog}}-E(C) \right)t \]

If \(E(C)>0\), the treatment reduces the effective progression rate.

Natural progression Treatment Time Disease measure

A treatment may reduce the rate of worsening without necessarily returning the disease marker to baseline.

This distinction is clinically important. A therapy that slows deterioration may have a meaningful effect even when the measured disease marker continues to worsen.

07 · Turnover systems

7. Disease Progression as a Turnover Process

Some disease biomarkers are naturally modeled as turnover systems, where the observed quantity reflects production and loss.

A basic turnover model is:

\[ \frac{dD}{dt}=k_{\text{in}}-k_{\text{out}}D \]

At steady state:

\[ D_{\text{ss}}=\frac{k_{\text{in}}}{k_{\text{out}}} \]

A drug can act on either process. For example, treatment might inhibit production:

\[ \frac{dD}{dt} = k_{\text{in}}\left(1-I(C)\right)-k_{\text{out}}D \]

Alternatively, the drug might increase loss or elimination of the disease marker.

Why this matters: a turnover model provides a mechanistic interpretation of delayed PD responses. Drug concentration can change rapidly while the disease marker changes more slowly because the marker itself has production and loss dynamics.
08 · Delayed effects

8. Disease Progression Can Produce Delayed Pharmacodynamic Effects

A drug concentration may reach its maximum quickly, while the clinical endpoint changes much more slowly. Disease progression can be one reason for this temporal separation.

One approach is to introduce an indirect-response or turnover component:

\[ \frac{dR}{dt} = k_{\text{in}}(C(t)) - k_{\text{out}}R \]

Another approach is to model a latent disease state that responds to drug exposure and then drives the observed endpoint.

\[ C(t) \rightarrow \text{drug effect} \rightarrow D(t) \rightarrow Y(t) \]

Here, \(D(t)\) represents a latent or mechanistic disease state and \(Y(t)\) is the measured clinical endpoint.

These models help distinguish PK delay, transduction delay, and disease-system delay, which can otherwise appear similar in concentration-effect plots.

09 · Biomarkers

9. Biomarkers as Disease-Progression Variables

Many disease-progression models use a biomarker as a quantitative representation of disease state.

For example, a biomarker might increase as disease worsens:

\[ B(t)=B_0+\text{progression component}-\text{treatment component} \]

The biomarker can then serve as an intermediate variable between drug exposure and a clinical endpoint.

Model componentExample role
PKDetermines the time-varying drug concentration.
Drug-effect modelConverts concentration into a pharmacologic signal.
Disease-state modelDescribes natural progression and treatment modification.
Biomarker modelLinks latent disease state to a measurable biomarker.
Clinical endpoint modelLinks disease state to symptoms, function, event risk, or another outcome.

This structure can make a PK/PD model more mechanistic because it separates drug exposure from the downstream biological process that ultimately determines the observed endpoint.

10 · Common structures

10. Common Disease-Progression Model Structures

Several broad structures appear repeatedly in pharmacometric applications.

StructureCore ideaTypical use
Baseline + slopeDisease changes at an approximately constant rate.Longitudinal biomarkers and functional measures.
Drug effect on slopeTreatment changes the rate of progression.Chronic progressive disease.
Drug effect on levelTreatment shifts the disease state.Endpoints with relatively rapid response.
Turnover modelProduction and loss determine the disease marker.Biomarkers with delayed response.
Transit compartmentsDrug effect propagates through intermediate states.Delayed pharmacodynamic responses.
Latent disease-state modelAn unobserved disease state drives measured outcomes.Mechanistic disease modeling.
State-transition modelPatients move between disease states over time.Discrete disease stages or severity categories.
11 · Worked example

11. Worked Example: A Drug That Slows Disease Progression

Consider a hypothetical disease biomarker for which higher values indicate greater disease burden. At baseline, the biomarker is 100 units. Without treatment, the disease progresses at approximately 4 units per month.

Step 1: Natural progression

The untreated trajectory can be represented as:

\[ D(t)=100+4t \]

After 6 months:

\[ D(6)=100+4(6)=124 \]

Thus, the model predicts a disease burden of 124 units under natural progression.

Step 2: Add a treatment effect

Suppose a sustained drug concentration produces an estimated treatment effect of 2 units/month on the progression slope.

\[ \frac{dD}{dt}=4-2=2 \]

The treated trajectory becomes:

\[ D_{\text{treated}}(t)=100+2t \]

Step 3: Predicted disease state after 6 months

\[ D_{\text{treated}}(6)=100+2(6)=112 \]

Step 4: Compare trajectories

TimeUntreated progressionTreated progressionDifference
0 months1001000
2 months1081044
4 months1161088
6 months12411212

The treatment does not cause the disease marker to decline below baseline. Instead, it reduces the rate of worsening. By month 6, the predicted disease burden is 12 units lower than it would have been under the untreated trajectory.

Interpretation: the important model quantity is not necessarily the absolute treatment response at one time point. In a progression model, the drug may exert its benefit by changing the trajectory of the disease over time.
12 · Connecting PK

12. Where Does Pharmacokinetics Enter the Model?

The progression model does not usually receive dose directly. Instead, the PK model converts dosing information into a time-varying concentration profile.

\[ \text{Dose} \rightarrow C(t) \rightarrow E(C(t)) \rightarrow \frac{dD}{dt} \rightarrow D(t) \]

For example, a simple one-compartment IV bolus model might provide:

\[ C(t)=\frac{D}{V}e^{-CLt/V} \]

The concentration then drives the drug-effect component of the disease-progression model:

\[ \frac{dD}{dt} = k_{\text{prog}} - \frac{E_{\max}C(t)}{EC_{50}+C(t)} \]

This creates a complete chain from dose to exposure, pharmacologic effect, disease trajectory, and ultimately the observed endpoint.

13 · Baseline and placebo

13. Why Placebo and Natural-History Data Are Valuable

Disease-progression models are particularly sensitive to information about what happens without active treatment.

Placebo or control data can help estimate the natural progression component, while treated data provide information about how treatment modifies that trajectory.

Modeling principle: when possible, separating natural disease progression from drug effect requires information that identifies both processes. A treatment-only dataset may contain limited information about what the disease would have done without treatment.

Longitudinal control data can therefore be valuable for estimating progression slopes, turnover parameters, baseline disease states, or other natural-history components.

14 · Population modeling

14. Disease Progression in Population PK/PD Models

In population pharmacometrics, disease progression parameters can vary between individuals just as PK parameters can.

For example, an individual progression rate might be represented as:

\[ k_{\text{prog},i} = k_{\text{prog,pop}} e^{\eta_{\text{prog},i}} \]

where \(k_{\text{prog,pop}}\) is the typical population progression rate and \(\eta_{\text{prog},i}\) represents individual deviation from that typical value.

Covariates may also explain differences in disease trajectories:

\[ k_{\text{prog},i} = k_{\text{prog,pop}} \left(\frac{\text{Covariate}_i}{\text{Reference}}\right)^\theta e^{\eta_i} \]

This allows a model to distinguish differences in disease natural history from differences in drug exposure and treatment response.

15 · Variability

15. Sources of Variability

Observed longitudinal outcomes usually contain several sources of variation.

SourceExample
Between-subject variabilityPatients differ in baseline disease state or natural progression rate.
PK variabilityClearance or volume differs between patients.
PD variabilityPatients differ in sensitivity to drug exposure.
Progression variabilityPatients have different natural disease trajectories.
Residual variabilityObserved measurements fluctuate around the model prediction.
Measurement errorThe biomarker or endpoint is measured imperfectly.

A useful disease-progression model attempts to separate these sources rather than treating all longitudinal variation as drug response.

16 · Model selection

16. Choosing the Disease-Progression Structure

Model selection should begin with the scientific question and the behavior of the observed endpoint.

  1. Define the disease endpoint. Is it a biomarker, continuous clinical measurement, categorical state, or event?
  2. Characterize natural history. Determine whether the endpoint is approximately stable, linear, nonlinear, or turnover-driven.
  3. Examine treatment timing. Determine whether drug effects appear rapidly or with delay.
  4. Consider plausible mechanisms. Decide whether treatment changes the disease level, slope, production, loss, or transition between states.
  5. Use the simplest adequate structure. Additional compartments or mechanisms should be supported by the data and scientific question.
  6. Evaluate diagnostics. Examine residuals, individual trajectories, parameter estimates, and predictive performance.
Avoid overinterpretation: a mathematically identifiable progression model is not automatically a validated biological mechanism. Different structures can sometimes produce similar predictions over the observed time window.
17 · Simulation

17. Using Disease-Progression Models for Simulation

Once adequately developed, a disease-progression PK/PD model can be used to simulate hypothetical treatment scenarios.

  • Different dose levels and dosing intervals.
  • Different exposure profiles.
  • Alternative treatment durations.
  • Different assumptions about progression rates.
  • Delayed or persistent drug effects.
  • Differences between treated and untreated trajectories.
  • Potential long-term consequences of modifying disease progression.

For chronic diseases, simulation can be especially informative because a modest change in progression rate can accumulate into a substantial difference between treated and untreated trajectories over longer periods.

18 · Interpretation

18. Important Limitations and Interpretation Issues

Disease-progression models can be powerful, but several issues require careful interpretation.

  • Natural history may be incompletely observed. The study duration may be too short to characterize long-term progression.
  • Progression and treatment effects can be confounded. Limited control data may make the two components difficult to separate.
  • Multiple models may fit the same data. A good statistical fit does not establish a unique biological mechanism.
  • Long-term extrapolation is model-dependent. Small differences in progression assumptions can produce large differences over extended periods.
  • Missing data can be informative. Dropout may depend on disease severity or treatment response.
  • Endpoint measurement matters. Measurement noise can obscure relatively slow progression.
  • Parameters have model-specific meanings. A progression rate in one structural model may not have the same interpretation in another.
Key interpretation point: a treatment that slows progression may produce little short-term change while producing a substantial separation between treated and untreated trajectories over time.
19 · Practical workflow

19. A Practical Workflow for Disease-Progression PK/PD Modeling

  1. Define the clinical question. Determine whether the goal is to characterize response, progression, delay, or long-term treatment effect.
  2. Define the endpoint. Identify the biomarker or clinical outcome that represents disease state.
  3. Explore longitudinal data. Plot individual trajectories and treatment groups over time.
  4. Characterize natural progression. Use control or placebo data where available.
  5. Build the PK component. Estimate or obtain individual concentration-time profiles.
  6. Develop the drug-effect model. Consider direct, indirect, turnover, or delayed-effect structures.
  7. Combine drug effect with disease progression. Specify how exposure changes the disease trajectory.
  8. Add variability and covariates. Explain important differences between patients where supported by data.
  9. Evaluate model adequacy. Use diagnostic plots, predictive checks, parameter plausibility, and sensitivity analyses.
  10. Simulate relevant scenarios. Use the model to investigate dosing and longer-term disease trajectories.

20. Key Takeaways

  • Disease-progression PK/PD models describe both the natural course of disease and the effect of drug exposure on that course.
  • The key conceptual distinction is between disease progression and drug effect.
  • A treatment does not necessarily have to reverse a disease marker to be effective; it may instead slow the rate of deterioration.
  • Progression can be represented using linear, nonlinear, turnover, logistic, transit, or state-transition models.
  • PK supplies the time-varying concentration \(C(t)\), which becomes an input to the drug-effect component.
  • Turnover and indirect-response models can explain delayed changes in biomarkers or clinical endpoints.
  • Control or placebo data can be particularly valuable for identifying the natural-history component of the model.
  • Population models can describe between-subject variability in baseline disease state, progression rate, drug exposure, and treatment response.
  • Different disease-progression models can sometimes describe the same observed data, so biological interpretation should be made cautiously.
  • Long-term predictions are especially sensitive to assumptions about the progression process.
  • A well-developed progression model can be used to simulate treatment trajectories and investigate how changes in disease progression accumulate over time.
  • The most useful model is the one that adequately represents the disease dynamics relevant to the scientific question without adding unsupported complexity.
Next step

Where to Go Next

A natural progression is to study indirect-response and turnover PK/PD models, followed by transit-compartment models, disease-progression models with delayed drug effects, tumor-growth inhibition models, biomarker production and loss models, and state-transition approaches.

These models build on the central idea introduced here: drug concentration changes over time, but the biological system receiving that signal may have its own dynamics. Understanding those dynamics is essential for separating immediate pharmacologic effects from longer-term changes in disease trajectory.

References

References

ReferenceRelevance
Jusko WJ, Ko HC. Physiologic indirect response models characterize diverse types of pharmacodynamic effects. Clinical Pharmacology & Therapeutics. 1994;56(4):406–419. Foundational framework for indirect-response and turnover approaches to pharmacodynamic modeling.
Dayneka NL, Garg V, Jusko WJ. Comparison of four basic models of indirect pharmacodynamic responses. Journal of Pharmacokinetics and Biopharmaceutics. 1993;21:457–478. Describes alternative mathematical structures for delayed pharmacodynamic responses.
Holford NH. A time to event tutorial for pharmacometricians. CPT: Pharmacometrics & Systems Pharmacology. Provides context for linking pharmacometric models with disease and clinical-event processes.
Holford NH, Sheiner LB. Understanding the dose-effect relationship: clinical application of pharmacokinetic-pharmacodynamic models. Clinical Pharmacokinetics. 1981;6:429–453. Classic discussion of PK/PD modeling and the interpretation of dose-effect relationships.
← Back to Pharmacokinetics Tutorials