Tutorials › Model-Informed Drug Development (MIDD) › Joint Optimization of Dose and Dosing Interval

Dose Optimization and Titration · Tutorial 258 of 596

Joint Optimization of Dose and Dosing Interval

Calculate steady-state exposure across a dosing interval, then choose dose and interval together to meet clinical targets under an explicitly stated regimen objective.

Advanced 12 min read

What You'll Learn

  • State the assumptions behind a repeated-dose, one-compartment bolus model.
  • Derive steady-state concentration, peak, trough, and average exposure.
  • Convert peak and trough targets into feasible dose bounds at a fixed interval.
  • Determine when a dose interval makes the target constraints infeasible.
  • Optimize dose and interval for a stated regimen objective and explain the model's limitations.

Why Dose and Interval Must Be Chosen Together

A maintenance dose and its dosing interval jointly shape the exposure profile between doses. Raising the dose tends to raise concentrations throughout the interval; lengthening the interval gives more time for concentrations to fall before the next dose. Consequently, evaluating dose amount alone—or evaluating average exposure without checking the interval profile—can miss whether clinical targets are met at the times that matter.

The introductory material on decision objectives and constraints for regimen optimization establishes why an objective and its constraints must be specified. Here, the objective will be to use the longest feasible interval, thereby reducing dosing frequency, while keeping the steady-state trough above an efficacy target and the peak below a safety limit. A different objective could produce a different preferred regimen; there is no unique optimum without a stated objective.

Steady-State Profile Under a Simple PK Model

Begin with a linear, one-compartment model with instantaneous intravenous bolus doses. Let \(D\) be the dose in mg, \(\tau\) the dosing interval in hours, \(V\) the apparent volume of distribution in L, and \(k\) the first-order elimination rate constant in h\(^{-1}\). Let \(F\) be the fraction of the dose reaching the systemic circulation; for intravenous dosing, \(F=1\). Assume the same dose is given every \(\tau\) hours, pharmacokinetics are time invariant, elimination is first order, and the patient has reached steady state. These assumptions exclude absorption delay, nonlinear clearance, changing physiology, and missed doses.

Immediately after a dose at steady state, the concentration includes the residual concentration from all previous doses. Each previous dose contributes a fraction \(e^{-k\tau}\) of its preceding contribution. The resulting geometric series gives the profile for \(0\leq t<\tau\), where \(t=0\) is immediately after a dose.

$$ C_{\mathrm{ss}}(t) = \frac{FD}{V}e^{-kt}\sum_{j=0}^{\infty}e^{-jk\tau} = \frac{FD}{V}\frac{e^{-kt}}{1-e^{-k\tau}} $$

The series converges because \(k>0\) and \(\tau>0\), so \(0<e^{-k\tau}<1\). Concentration declines between instantaneous doses. Thus the maximum occurs at \(t=0\), just after dosing, and the minimum occurs just before the next dose, at \(t=\tau\).

Definition: In this model, the steady-state peak and trough concentrations are \(C_{\max,\mathrm{ss}}=\dfrac{FD/V}{1-e^{-k\tau}}\) and \(C_{\min,\mathrm{ss}}=\dfrac{FD/V}{e^{k\tau}-1}\), respectively. The steady-state average concentration over an interval is \(C_{\mathrm{avg,ss}}=\dfrac{FD}{CL\,\tau}\), where \(CL=kV\).

The trough expression follows by evaluating the profile at the end of the interval: multiply the peak by \(e^{-k\tau}\). The average follows because the area under the concentration curve over one steady-state interval is \(FD/CL\); dividing this interval AUC by \(\tau\) gives the average concentration. Peak, trough, and average answer different questions. Meeting an average-exposure target does not establish that a trough or peak constraint is met.

Worked Example: Checking One Dose–Interval Pair

Suppose a hypothetical adult receives \(D=100\) mg intravenously every \(\tau=12\) h. Let \(V=2.00\) L and \(k=0.100\) h\(^{-1}\). The clinical targets are a trough of at least 20.0 mg/L and a peak no greater than 80.0 mg/L. Check the regimen under the stated model.

The assumptions are those of the one-compartment, linear, instantaneous bolus model; the calculation also presumes dosing has continued long enough to attain steady state. First, \(e^{-k\tau}=e^{-1.2}\approx0.3012\), so \(1-e^{-k\tau}\approx0.6988\). The peak is

$$ C_{\max,\mathrm{ss}} = \frac{100\ \mathrm{mg}/2.00\ \mathrm{L}}{1-e^{-0.100(12)}} = \frac{50.0\ \mathrm{mg/L}}{0.6988} \approx71.55\ \mathrm{mg/L}. $$

For the trough, \(e^{1.2}\approx3.3201\), so \(e^{1.2}-1\approx2.3201\):

$$ C_{\min,\mathrm{ss}} = \frac{100\ \mathrm{mg}/2.00\ \mathrm{L}}{e^{0.100(12)}-1} = \frac{50.0\ \mathrm{mg/L}}{2.3201} \approx21.55\ \mathrm{mg/L}. $$

The calculated peak is below 80.0 mg/L, and the trough is above 20.0 mg/L, so this pair meets both specified constraints. As a consistency check, \(CL=kV=0.200\) L/h and \(C_{\mathrm{avg,ss}}=100/(0.200\times12)\approx41.67\) mg/L. The average lies between the peak and trough, as expected for a declining concentration profile.

Turning Target Constraints Into Dose Bounds

For fixed \(\tau\), both peak and trough are proportional to dose. Suppose the required trough is at least \(C_{\mathrm{trough,target}}\), while the permitted peak is at most \(C_{\mathrm{peak,max}}\). Rearranging the trough and peak formulas gives a lower and an upper dose bound. A feasible dose exists only if the lower bound does not exceed the upper bound.

$$ D_{\min}(\tau) = C_{\mathrm{trough,target}}\frac{V}{F}(e^{k\tau}-1), \qquad D_{\max}(\tau) = C_{\mathrm{peak,max}}\frac{V}{F}(1-e^{-k\tau}). $$

These bounds are meaningful when \(F>0\), \(V>0\), \(k>0\), and \(\tau>0\). A dose can satisfy both targets at that interval if and only if \(D_{\min}(\tau)\leq D_{\max}(\tau)\). Equality leaves exactly one dose satisfying both constraints; strict inequality leaves a range of feasible doses. This test avoids selecting a dose first and discovering afterward that the interval makes the targets incompatible.

Worked Example: Testing Feasibility at Two Intervals

Use the hypothetical parameters above, with \(V/F=2.00\) L, \(k=0.100\) h\(^{-1}\), \(C_{\mathrm{trough,target}}=20.0\) mg/L, and \(C_{\mathrm{peak,max}}=80.0\) mg/L. Calculate the feasible dose range at \(\tau=12\) h and check whether any dose works at \(\tau=16\) h.

At 12 h, the lower bound is

$$ D_{\min}(12) = (20.0\ \mathrm{mg/L})(2.00\ \mathrm{L})(e^{1.2}-1) = 40.0(2.3201) \approx92.80\ \mathrm{mg}. $$

The upper bound is

$$ D_{\max}(12) = (80.0\ \mathrm{mg/L})(2.00\ \mathrm{L})(1-e^{-1.2}) = 160(0.6988) \approx111.8\ \mathrm{mg}. $$

Therefore, any dose from about 92.80 mg through 111.8 mg satisfies both bounds at 12 h under this model. At 16 h, \(e^{1.6}\approx4.9530\), \(e^{1.6}-1\approx3.9530\), and \(e^{-1.6}\approx0.2019\). The bounds are

$$ D_{\min}(16) = (20.0)(2.00)(3.9530) \approx158.1\ \mathrm{mg}, \qquad D_{\max}(16) = (80.0)(2.00)(1-0.2019) \approx127.7\ \mathrm{mg}. $$

The lower bound exceeds the upper bound, so no dose can meet both targets at 16 h. Increasing the dose enough to reach the trough target would violate the peak limit; reducing it enough to respect the peak limit would miss the trough target. The units in both bound calculations reduce to mg.

Joint Optimization for the Longest Feasible Interval

The dose bounds can be used to optimize the interval instead of checking only a few candidate schedules. For the stated objective—maximize \(\tau\) subject to both exposure constraints—feasibility requires

$$ C_{\mathrm{trough,target}}\frac{V}{F}(e^{k\tau}-1) \leq C_{\mathrm{peak,max}}\frac{V}{F}(1-e^{-k\tau}). $$

Since \(V/F>0\), cancel it. Set \(x=e^{k\tau}\), which is greater than 1 for positive \(k\) and \(\tau\). Then \(1-e^{-k\tau}=(x-1)/x\), so the inequality becomes \(C_{\mathrm{trough,target}}(x-1)\leq C_{\mathrm{peak,max}}(x-1)/x\). Multiplying by \(x>0\) and dividing by \(x-1>0\) gives \(C_{\mathrm{trough,target}}x\leq C_{\mathrm{peak,max}}\). Thus, when the peak limit exceeds the trough target, the largest feasible interval under this model is

$$ \tau_{\max} = \frac{1}{k}\ln\left( \frac{C_{\mathrm{peak,max}}}{C_{\mathrm{trough,target}}} \right). $$

At this limiting interval, the lower and upper dose bounds coincide. For any shorter positive interval the ratio test permits both constraints to be met, subject to any additional dose restrictions; for any longer interval it does not. The result depends on the assumed objective and model. It is not a universal recommendation to choose the longest possible interval.

Worked Example: Finding the Longest Feasible Interval and Its Dose

For the hypothetical targets of 20.0 mg/L at trough and 80.0 mg/L at peak, with \(k=0.100\) h\(^{-1}\), find the longest interval and the dose that meets both targets at that boundary.

The peak-to-trough target ratio is \(80.0/20.0=4.00\), so

$$ \tau_{\max} = \frac{\ln(4.00)}{0.100\ \mathrm{h}^{-1}} = \frac{1.3863}{0.100\ \mathrm{h}^{-1}} \approx13.86\ \mathrm{h}. $$

At this interval, \(e^{k\tau}=4.00\). The dose needed to put the trough exactly at 20.0 mg/L is

$$ D = (20.0\ \mathrm{mg/L})(2.00\ \mathrm{L})(4.00-1) = 120\ \mathrm{mg}. $$

Checking the peak, \(e^{-k\tau}=1/4\), and \(C_{\max,\mathrm{ss}}=(120/2.00)/(1-1/4)=60.0/0.750=80.0\) mg/L. The boundary regimen therefore meets both targets exactly. This is the continuous-model optimum for minimizing dosing frequency under the specified constraints. It does not account for available tablet strengths, permitted dosing times, patient variability, or uncertainty in the PK parameters.

Diagnostics, Alternatives, and Limitations

A useful diagnostic is to check the profile's internal ordering: with positive dose and elimination, the peak must exceed the trough, and the interval average must lie between them. If a computed trough exceeds its peak, inspect whether the interval was placed in the exponent correctly. Also check units: \(k\tau\) must be dimensionless, \(D/V\) has concentration units, and the dose bounds must have dose units. The limiting interval formula requires \(C_{\mathrm{peak,max}}>C_{\mathrm{trough,target}}\); if the two targets are equal or the peak limit is lower, no positive interval satisfies both in this declining-profile model.

A common mistake is to calculate a dose from the trough target and assume the peak will consequently be acceptable. The trough calculation gives only a lower dose bound. The peak target gives an upper bound, and both must be checked at the same interval. Another error is to use \(e^{k\tau}\) where the peak formula requires \(1-e^{-k\tau}\), or to omit the subtraction of 1 in \(e^{k\tau}-1\). Keeping the derived lower and upper bounds visible makes these errors easier to catch.

An average-exposure-only approach is simpler when the clinical target is genuinely defined by interval AUC or average concentration, but it cannot guarantee a minimum trough or maximum peak. Conversely, this analytic calculation is convenient for a one-compartment linear bolus model; a more detailed PK model can calculate steady-state profiles numerically when absorption, multiple compartments, nonlinear elimination, or complex dosing makes the closed-form equations inappropriate. The model should be checked against observed data and the assumptions should be defensible before its calculated regimen is used to support a clinical decision.

Finally, the calculation is deterministic: the same parameters imply the same profile. A regimen meeting the targets for a typical parameter set need not meet them for every patient. Incorporating population variability or uncertainty requires additional analyses; this calculation alone does not establish a population success rate. Its proper role is to make the dose–interval consequences and the feasibility of specified targets explicit under a stated model.

Check Your Understanding

Use the one-compartment bolus model and show the essential calculation for each question. Answers follow the questions.

  1. If \(k=0.0800\) h\(^{-1}\), what is the fraction of a dose's concentration contribution remaining after 10.0 h?
  2. Why is a dose below \(D_{\min}(\tau)\) insufficient when the trough target is a minimum concentration?
  3. With \(V/F=3.00\) L, \(k=0.100\) h\(^{-1}\), \(\tau=8.00\) h, and a trough target of 15.0 mg/L, calculate the minimum dose.
  4. For peak and trough targets of 60.0 and 15.0 mg/L, respectively, with \(k=0.120\) h\(^{-1}\), calculate the longest feasible interval under the model.

Solutions. (1) The remaining fraction is \(e^{-k\tau}=e^{-0.0800(10.0)}=e^{-0.800}\approx0.4493\), or 44.93%. (2) A dose below \(D_{\min}(\tau)=C_{\mathrm{trough,target}}(V/F)(e^{k\tau}-1)\) produces a steady-state trough below the required minimum, by rearrangement of the trough equation. (3) \(D_{\min}=(15.0\ \mathrm{mg/L})(3.00\ \mathrm{L})(e^{0.800}-1)=45.0(2.2255-1)\approx55.15\) mg. (4) \(\tau_{\max}=\ln(60.0/15.0)/(0.120\ \mathrm{h}^{-1})=\ln(4)/0.120\approx11.55\) h. At this interval the two target constraints meet at their boundary; a longer interval cannot satisfy both under the stated assumptions.