Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Mapping Dosing Regimens to Predicted Dose Response

Dose-Response and Dose Selection · Tutorial 532 of 598

Mapping Dosing Regimens to Predicted Dose Response

Translate a proposed dosing regimen into predicted exposure and response while accounting for route, dosing interval, and relevant sources of variability.

Advanced 12 min read

What You'll Learn

  • Distinguish a regimen-level dose-response prediction from an exposure-response relationship.
  • Calculate steady-state exposure for intravenous bolus regimens with different intervals.
  • Represent oral absorption and bioavailability when predicting concentrations.
  • Propagate interindividual clearance variability into exposure and response distributions.
  • Diagnose when a chosen exposure summary or model assumption does not adequately describe a regimen.

From a Proposed Regimen to a Predicted Response

A dose amount alone does not specify the concentration-time profile experienced by a patient. The dosing interval determines how much drug remains when the next dose is given, and the route determines how drug enters the systemic circulation. These differences can matter even when two regimens deliver the same nominal amount per day. To predict response from a proposed regimen, a PK/PD model first maps that regimen to exposure and then maps exposure to a response.

This mapping also clarifies the distinction between two kinds of relationship. A dose-response relationship describes the predicted response associated with a dose or regimen, after accounting for the PK processes connecting dose to exposure. An exposure-response relationship describes the predicted response conditional on a specified exposure measure, such as average concentration or peak concentration. The latter does not, by itself, say which dose or route will produce that exposure.

Definition: A regimen is specified by at least the dose amount, dosing interval, and route; it may also specify formulation, infusion duration, and timing. A regimen-specific response prediction combines a PK model for the exposure generated by that regimen with a PD model relating exposure to response. If PK parameters or patient characteristics vary, the prediction can be a distribution of exposures and responses rather than one value.

A useful conceptual representation is \(D \rightarrow C(t) \rightarrow E\), where \(D\) represents the regimen, \(C(t)\) is the resulting concentration-time profile, and \(E\) is the predicted response. The arrows encode model assumptions: the PK model must represent the relevant route and dosing schedule, and the PD model must use an exposure measure and response relationship appropriate to the drug and endpoint.

Building the Regimen-Specific Exposure Prediction

For linear PK at steady state, the average concentration over a dosing interval is determined by the systemic amount delivered per interval and clearance. For an intravenous dose \(D\) every \(\tau\) hours, the systemic amount is \(D\). For an oral dose with bioavailability \(F\), it is \(FD\). The average concentration is therefore:

$$ C_{\mathrm{avg},ss} = \frac{F D}{\tau\,CL} $$

Here \(CL\) is clearance, and \(F=1\) for an intravenous dose. Units must be consistent: if \(D\) is in mg, \(\tau\) in hours, and \(CL\) in L/h, then \(C_{\mathrm{avg},ss}\) is in mg/L. This average does not determine the full profile. For a given average, a short interval may produce smaller peak-to-trough fluctuations than a longer interval.

For a one-compartment model with instantaneous intravenous boluses, elimination rate constant \(k=CL/V\), volume \(V\), and repeated dose interval \(\tau\), the steady-state concentration immediately after a dose and immediately before the next dose are:

$$ C_{\max,ss}=\frac{D/V}{1-e^{-k\tau}}, \qquad C_{\min,ss}=C_{\max,ss}e^{-k\tau}. $$

These expressions assume linear elimination, constant PK parameters, complete intravenous delivery, and doses repeated long enough to reach steady state. They do not apply unchanged to an oral regimen, because oral absorption takes time and may be incomplete.

Worked Example: Same Daily Intravenous Dose, Different Intervals

Hypothetical data. A one-compartment drug has \(V=50\) L and \(CL=5\) L/h. Compare 100 mg intravenously every 24 hours with 50 mg every 12 hours at steady state. Assume instantaneous boluses and linear elimination. The elimination rate constant is \(k=CL/V=5/50=0.1000\ \mathrm{h}^{-1}\).

For 100 mg every 24 hours, \(D/V=100/50=2.000\) mg/L and \(e^{-k\tau}=e^{-2.4}=0.09072\). Thus:

$$ C_{\max,ss} = \frac{2.000}{1-0.09072} = 2.19954\ \mathrm{mg/L}, \qquad C_{\min,ss} = 2.19954(0.09072) = 0.1995\ \mathrm{mg/L}. $$

For 50 mg every 12 hours, \(D/V=50/50=1.000\) mg/L and \(e^{-k\tau}=e^{-1.2}=0.301194\). Therefore:

$$ C_{\max,ss} = \frac{1.000}{1-0.301194} = \frac{1.000}{0.698806} = 1.4310\ \mathrm{mg/L}, \qquad C_{\min,ss} = 1.4310(0.301194) = 0.4310\ \mathrm{mg/L}. $$

Both regimens deliver 100 mg per day. Their average concentration is the same under these linear, steady-state assumptions: \((100\ \mathrm{mg}/24\ \mathrm{h})/(5\ \mathrm{L/h})=0.8333\ \mathrm{mg/L}\). The 12-hour regimen has a lower peak and a higher trough, although its average concentration is unchanged. If response depends on peak, trough, or time above a relevant concentration, equal daily dose and equal average concentration need not imply equal predicted response.

For oral dosing, the model must represent absorption as well as elimination. In a one-compartment model with first-order absorption rate constant \(k_a\), linear elimination rate constant \(k\), bioavailability \(F\), and repeated doses every \(\tau\), the steady-state concentration \(t\) hours after a dose, for \(0\leq t<\tau\), is:

$$ C_{ss}(t) = \frac{F D k_a}{V(k_a-k)} \left[ \frac{e^{-kt}}{1-e^{-k\tau}} - \frac{e^{-k_a t}}{1-e^{-k_a\tau}} \right], \qquad k_a\ne k. $$

This expression assumes first-order absorption, no lag time, one-compartment disposition, constant bioavailability, linear elimination, and repeated dosing until steady state. A different formulation or route can require a different input model; it should not be treated as a simple change in dose unless its effect on absorption and bioavailability is adequately represented.

Worked Example: Predicting Oral Concentration at a Specified Time

Hypothetical data. Consider 100 mg orally every 24 hours, with \(F=0.80\), \(V=50\) L, \(k=0.1000\ \mathrm{h}^{-1}\), and \(k_a=1.00\ \mathrm{h}^{-1}\). Assume the stated one-compartment, first-order absorption model and steady state. Predict concentration 6 hours after a dose.

The coefficient is \(FDk_a/[V(k_a-k)]=0.80(100)(1.00)/[50(1.00-0.1000)]=80/45=1.7778\) mg/L. Substituting \(t=6\) hours and \(\tau=24\) hours:

$$ C_{ss}(6) = 1.7778 \left[ \frac{e^{-0.1000(6)}}{1-e^{-0.1000(24)}} - \frac{e^{-1.00(6)}}{1-e^{-1.00(24)}} \right] \approx 1.0686\ \mathrm{mg/L} \approx 1.069\ \mathrm{mg/L}. $$

The predicted value is a concentration at a specified time, not an average concentration or a response prediction. Bioavailability reduces systemic input relative to the same nominal intravenous dose, while the absorption rate shapes when concentrations rise. A PD model is still needed to translate this profile or a suitable summary of it into response.

Mapping Exposure to Response, Including Variability

A PD model specifies how response changes with exposure. For example, a simple Emax relationship for a response that increases with concentration is:

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

Here \(E_0\) is the baseline response, \(E_{\max}\) is the maximum increase above baseline, and \(EC_{50}\) is the concentration producing half of that maximum increase. The formula presumes that the chosen concentration measure is the relevant predictor and that this direct relationship is suitable for the endpoint. If response is delayed, depends on concentration history, or is driven by peak or trough rather than average concentration, applying the formula to \(C_{\mathrm{avg},ss}\) may be inappropriate.

When individuals differ in PK parameters, the regimen generates an exposure distribution. Passing a population-average exposure through a nonlinear PD model generally does not give the average predicted response. Instead, the PK and PD calculations should be performed for representative individuals or parameter draws and then summarized across the target population. This distinction matters even if the model itself is deterministic for a given individual's parameters.

Worked Example: From a Regimen to a Population Response Distribution

Hypothetical data. Suppose patients have either \(CL=4\) L/h or \(CL=8\) L/h, with equal probability. They receive 80 mg orally every 24 hours, with \(F=0.80\). Assume linear steady-state PK and that response depends on average concentration through an Emax model with \(E_0=0\), \(E_{\max}=100\) percentage points, and \(EC_{50}=0.500\) mg/L.

The average systemic input rate is \(FD/\tau=0.80(80\ \mathrm{mg})/(24\ \mathrm{h})=2.6667\) mg/h. For a patient with \(CL=4\) L/h: \(C_{\mathrm{avg},ss}=2.6667/4=0.6667\) mg/L, giving \(E=100(0.6667)/(0.500+0.6667)=57.14\%\). For a patient with \(CL=8\) L/h: \(C_{\mathrm{avg},ss}=2.6667/8=0.3333\) mg/L, giving \(E=100(0.3333)/(0.500+0.3333)=40.00\%\).

With equal proportions in the two groups, the mean predicted response is \((57.14\%+40.00\%)/2=48.57\%\). The mean clearance is 6 L/h, which gives an average concentration of \(2.6667/6=0.4444\) mg/L and a response of \(100(0.4444)/(0.500+0.4444)=47.06\%\). That response is not the population mean prediction: the nonlinear PD transformation makes the average of individual predictions differ from the prediction at the average clearance.

This example produces a simple two-point response distribution. With continuously distributed covariates or PK parameters, the same logic applies by calculating concentrations and responses over the modeled population distribution. The resulting distribution reflects the variability included in the model; it should not be mistaken for uncertainty about the parameters unless parameter uncertainty has also been propagated.

Justification, Diagnostics, and Alternatives

The average-concentration equation follows from mass balance at steady state. Under linear PK, the amount eliminated over one interval equals the systemically available dose \(FD\). If concentration is \(C(t)\), the amount eliminated over the interval is \(CL\int_0^\tau C(t)\,dt\). Dividing by the interval and equating average elimination rate to average systemic input rate gives:

$$ CL\frac{1}{\tau}\int_0^\tau C(t)\,dt = \frac{FD}{\tau} \quad\Longrightarrow\quad C_{\mathrm{avg},ss} = \frac{1}{\tau}\int_0^\tau C(t)\,dt = \frac{FD}{\tau CL}. $$

The justification depends on steady state and linear clearance: the amount eliminated is proportional to concentration, with the same \(CL\) throughout the interval. Nonlinear elimination, time-varying clearance, nonadherence, or a changing disease state can invalidate this simple mapping.

Before interpreting a regimen-specific prediction, check whether the model reproduces observed concentration-time patterns for the route, formulation, and interval being considered. Systematic residual patterns, for example underprediction near peaks and overprediction during absorption, may indicate a misspecified input or absorption model. Check that dose units, time units, bioavailability, and route are implemented consistently. Predictions far outside the observed dose or exposure range are extrapolations and may be especially sensitive to model assumptions.

A mechanistic PK/PD model is not the only way to estimate dose-response. An empirical dose-response model can relate assigned dose directly to response and may be useful when the objective is prediction within the regimens studied. However, it may not explain how a new interval or route changes exposure, and it can confound differences in exposure with differences in regimen. An exposure-response model can describe response conditional on exposure, but does not by itself identify the exposure produced by a proposed dose. Integrating PK and PD makes the mapping between regimen, exposure, and response explicit, at the cost of assumptions in both parts of the model.

A common mistake is to assume that equal total daily dose means equal response. In the intravenous example, the two regimens had equal average concentration only under linear steady-state PK; their peaks and troughs differed. If PD depends on those profile features, the predicted responses may differ. Another mistake is to calculate response at the population-average exposure and call it the average response. With nonlinear PD, the worked example shows why predictions should instead be made across the relevant population distribution before summarizing.

Key takeaway: A regimen-specific dose-response prediction requires a route- and schedule-appropriate PK model, a justified exposure measure, and a suitable PD model. State which sources of variability are included, and distinguish the resulting predicted population distribution from parameter uncertainty.

Check Your Understanding

Try these questions, then compare your reasoning with the solutions.

  1. For a drug with linear steady-state PK, what additional information is needed beyond a dose amount to predict average concentration after oral dosing?
  2. Two intravenous regimens have equal daily dose and linear PK. Why can their predicted responses still differ?
  3. For a nonlinear Emax relationship, why is response at mean exposure not necessarily the mean response?
  4. What does an exposure-response model alone fail to specify when evaluating a new route?

Solutions. (1) At minimum, the dosing interval, clearance, and oral bioavailability are needed for the average-concentration calculation; a full profile additionally requires an appropriate absorption and disposition model. (2) A different interval changes the concentration-time profile, including peak and trough, which may matter if the PD relationship depends on those features rather than only average concentration. (3) In general, applying a nonlinear function to an average is not equivalent to averaging the function applied to individual values: \(f(\bar{C})\ne\overline{f(C)}\). (4) It does not specify the exposure produced by the new route; a route-appropriate PK model is needed to map the regimen to exposure.

Clinical Biostats