Tutorials › Model-Informed Drug Development (MIDD) › Selecting Titration Step Sizes and Adjustment Intervals

Dose Optimization and Titration · Tutorial 264 of 596

Selecting Titration Step Sizes and Adjustment Intervals

Use pharmacokinetic and pharmacodynamic time courses to choose defensible dose increments and the time to wait before reassessing a titration step.

Advanced 12 min read

What You'll Learn

  • Calculate the fraction of a pharmacokinetic dose change expressed after a specified number of elimination half-lives.
  • Relate a permitted exposure change to a dose increment using a local dose-exposure relationship.
  • Derive the response time course for linked first-order PK and PD equilibration.
  • Choose a provisional adjustment interval while accounting for both accumulation and equilibration.
  • Identify when simplified time-course calculations need to be replaced by regimen-specific simulation.

Why Step Size and Timing Must Be Chosen Together

A titration rule specifies how dose changes in response to an observation. Its practical performance also depends on how large each change is and how long the patient remains on that dose before the next adjustment. A large increase can move exposure rapidly toward a useful range, but it can also create an excessive transient exposure or response. An assessment made too soon may reflect only part of the previous dose change; another increase at that point can compound the transient.

The calculations below provide a model-based way to screen candidate increments and intervals. They are not universal clinical rules. The choice depends on the exposure and response metric used by the titration rule, the acceptable transient, the dosing regimen, and the model’s assumptions. In particular, repeated oral doses produce peak–trough fluctuations that a single smooth time constant does not fully describe.

Time Courses After a Dose Change

Definition: A titration step changes the prescribed dose from \(D_{\mathrm{old}}\) to \(D_{\mathrm{new}}\). The adjustment interval is the time from that change until the next planned dose decision. The residual fraction at time \(t\) is the fraction of the difference between the old and new steady-state values that remains unexpressed at that time. An adjustment interval is adequate only relative to a specified exposure or response criterion; there is no single settling time that applies to every drug and endpoint.

For a one-compartment model with first-order elimination, the elimination rate constant is \(k=\mathrm{CL}/V\), and the half-life is \(t_{1/2}=\ln(2)/k\). With repeated doses at a fixed interval \(\tau\), the concentration at a fixed point in the dosing cycle approaches its steady-state value geometrically. After \(m\) intervals, the fraction of the new steady-state change expressed is \(1-\exp(-km\tau)\), provided the doses and disposition are linear and unchanged apart from the dose amount. The remaining fraction is \(\exp(-km\tau)\).

This same expression describes the change from one steady-state dose to another at a fixed phase of the dosing cycle under linear pharmacokinetics: by superposition, the difference from the new steady state decays at the elimination rate. It does not mean that concentration is constant within an interval. Peak and trough concentrations can remain clinically important even when the interval-averaged or same-phase change is mostly expressed.

A common pharmacodynamic approximation is a first-order equilibration model. If the exposure has stabilized and the corresponding equilibrium response is \(R_{\mathrm{new}}\), then the response \(R(t)\) approaches it at rate \(k_{\mathrm{PD}}\). Its residual fraction is \(\exp(-k_{\mathrm{PD}}t)\), with response half-life \(\ln(2)/k_{\mathrm{PD}}\). If exposure itself is still changing, the PD input is not an immediate step to its new equilibrium value; the coupled time course must be considered.

Worked Example: Estimating Pharmacokinetic Settling Time

Suppose a hypothetical drug has an elimination half-life of 24 hours and is given once daily. The titration rule uses a same-phase exposure measure, and the planning criterion is to express at least 90% of the steady-state change before reassessment. Assume linear one-compartment elimination and unchanged dosing interval.

The residual fraction after \(m\) daily intervals is \( \exp(-km\tau)=2^{-m\tau/t_{1/2}}=2^{-m} \), because \(\tau=t_{1/2}=24\) hours. The criterion is \(1-2^{-m}\geq 0.90\), or \(2^{-m}\leq 0.10\).

$$ m \geq \frac{\ln(0.10)}{\ln(0.5)}=3.322 $$

A planned assessment after a whole number of daily intervals therefore requires four intervals. After three days, the expressed fraction is \(1-2^{-3}=0.875\), or 87.5%, which does not meet the criterion. After four days it is \(1-2^{-4}=0.9375\), or 93.75%. Under these assumptions, four days is the first daily assessment time meeting the chosen 90% PK criterion. This calculation alone does not establish that the PD response has settled.

Choosing a Dose Increment

A dose step should be evaluated against the change it is predicted to cause in a decision-relevant quantity: for example, steady-state exposure, a response measure, or a modeled safety endpoint. If exposure \(E\) is approximately proportional to dose in the local range, define the local slope \(g=\Delta E/\Delta D\). A prespecified maximum permitted exposure change \(\Delta E_{\max}\) then gives the screening constraint \(|\Delta D|\leq\Delta E_{\max}/|g|\). This is a model-based bound, not a claim that one fixed percentage increase is safe for all drugs or patients.

When dose proportionality is not credible, calculate the predicted change directly as \(E(D_{\mathrm{new}})-E(D_{\mathrm{old}})\) using the exposure model. A nonlinear dose–exposure or exposure–response relationship can make equal dose increments produce unequal changes in exposure or response. The step may also need to respect available strengths and rounding constraints, while ensuring the rounded dose still meets the modeled bound. The dose-optimization objectives and constraints established for the regimen remain relevant: a convenient increment is not justified if its predicted transient or steady-state consequences are unacceptable.

Worked Example: Applying an Exposure-Based Step Limit

In a hypothetical linear range, a model predicts steady-state exposure \(E=0.08D\), where dose \(D\) is in milligrams and exposure is in arbitrary units. A patient is taking 50 mg, and the selected target exposure is 8.0 units. For this illustration, the prespecified maximum exposure increase at one titration step is 1.2 units. Available doses are in 5 mg increments.

At 50 mg, predicted exposure is \(0.08(50)=4.0\) units. The exposure slope is \(g=0.08\) units per mg, so the maximum unrounded dose increment is \( \Delta D_{\max}=1.2/0.08=15\) mg. The largest available increment meeting that bound is 15 mg. The next dose is therefore 65 mg, with predicted exposure \(0.08(65)=5.2\) units; the increase is \(5.2-4.0=1.2\) units.

If the same bound is applied at subsequent steps, the sequence could be 50, 65, 80, 95, and 100 mg. Predicted exposures are 4.0, 5.2, 6.4, 7.6, and 8.0 units, respectively. The final step is 5 mg, corresponding to a 0.4-unit increase. This arithmetic demonstrates consistency with the stated exposure cap; it does not establish that the cap is clinically appropriate. That requires justification from the exposure–response and safety objectives for the drug.

Combining PK Accumulation With PD Equilibration

When both exposure and response change over time, a useful simplified model represents the normalized exposure change \(C(t)\) as a first-order approach to its new level, and the normalized response change \(R(t)\) as a first-order equilibration toward \(C(t)\): \(dC/dt=k_{\mathrm{PK}}(1-C)\) and \(dR/dt=k_{\mathrm{PD}}(C-R)\), with \(C(0)=R(0)=0\). Here, 0 and 1 denote the old and new steady-state levels on normalized scales. This model assumes a linear PD link over the change considered, constant rate constants, and no extra delay or effect-site process beyond the stated PD equilibration. It is a smooth approximation to the exposure trajectory; a regimen with substantial peak–trough variation should be assessed with its actual dosing model.

The PK equation gives \(C(t)=1-\exp(-k_{\mathrm{PK}}t)\). To obtain the response, define the response deviation from its new equilibrium as \(Y(t)=R(t)-1\). Since 1 is constant, \(dY/dt=dR/dt=k_{\mathrm{PD}}(C-R)\). Substituting \(C(t)\) and solving the resulting linear equation gives the residual response fraction below. This derivation makes explicit that the response is driven by the changing exposure, not by an instantaneous jump to the new PD equilibrium.

$$ 1-R(t)= \frac{k_{\mathrm{PD}}\exp(-k_{\mathrm{PK}}t) -k_{\mathrm{PK}}\exp(-k_{\mathrm{PD}}t)} {k_{\mathrm{PD}}-k_{\mathrm{PK}}}, \qquad k_{\mathrm{PD}}\ne k_{\mathrm{PK}} $$

When the two rate constants are equal to \(k\), the limiting residual is \((1+kt)\exp(-kt)\). These expressions can be used to check a stated residual-response criterion under the simplified model. They also show why simply taking the larger of the separate PK and PD settling times need not guarantee that the coupled response has settled: the PD response continues to lag while the exposure is moving.

Worked Example: Finding an Interval for a Coupled PK–PD Change

Consider a hypothetical drug with a PK half-life of 18 hours and a PD equilibration half-life of 36 hours. Use the coupled first-order model above, and require the residual response deviation to be no more than 10%. For the numerical calculation, use the smooth exposure approximation; the final clinical interval would also need checking against the actual dosing schedule and concentration metric.

The rate constants are \(k_{\mathrm{PK}}=\ln(2)/18=0.03851\) per hour and \(k_{\mathrm{PD}}=\ln(2)/36=0.01925\) per hour. At 144 hours, the exponentials are \(\exp(-k_{\mathrm{PK}}t)=1/256=0.003906\) and \(\exp(-k_{\mathrm{PD}}t)=1/16=0.06250\). Substitution gives \(1-R(144)=0.1211\), or 12.11%. Six days does not meet the 10% criterion.

At 168 hours, the corresponding exponentials are approximately 0.001550 and 0.03937. Therefore, \(1-R(168)\approx[(0.01925)(0.001550)-(0.03851)(0.03937)]/(0.01925-0.03851)=0.0772\), or 7.72%. Seven days meets the criterion in this model. The interval is therefore at least seven days when assessments are planned in whole days. This is a model-based result for the stated assumptions, not a general recommendation for a drug with these half-lives.

Practical Checks, Alternatives, and Limitations

Before adopting an increment and interval, check that the calculation concerns the same metric used by the titration rule. A trough-based rule cannot be justified solely by a calculation for average exposure if peak–trough differences are material. Check the model’s dose proportionality in the proposed range, whether the half-life or clearance changes with dose or time, and whether the response model represents the endpoint and its delay. If relevant, simulate the exact dosing times, dose history, and observation times, and inspect both transient and steady-state predictions.

A useful diagnostic is to compare predicted and observed changes after a step at the intended assessment time. Persistent drift may indicate that the interval is too short, the assumed time constants are wrong, or the response model omits a delay. A response that changes much more or less than predicted can also indicate nonlinearity, dose-dependent disposition, adherence or administration differences, or measurement error. These patterns do not identify a cause by themselves; they signal that the simplified calculation should not be treated as sufficient evidence for the next step.

The principal alternative to half-life-based screening is a regimen-specific PK/PD simulation that predicts the exact metric at candidate assessment times. It can account for dose timing, peaks and troughs, nonlinear kinetics, covariates, and a nonlinear response model when these are supported by the available model. A simpler half-life calculation is transparent and useful for an initial screen, but can mislead when the assumptions are poor. Conversely, simulation is only as reliable as its structural assumptions and parameter values; more detail does not remove uncertainty.

Ultimately, step size and interval are separate design choices that interact. A smaller step may reduce the consequence of a transient overshoot but can require more adjustments; a longer interval can allow more of the exposure and response change to appear but may delay reaching a useful regimen. The selected choices should be justified against the objective and constraints of the titration rule, and checked over plausible model behavior rather than inferred from half-life alone.

Common pitfall: “The drug has reached four half-lives, so the response must be settled.” Four half-lives leave about 6.25% of a simple first-order change, but that statement applies only to the modeled process. If the PD response equilibrates more slowly, or if the exposure has substantial peak–trough fluctuation, the response criterion may still fail. Calculate or simulate the response metric itself.

Check Your Understanding

Try each question, then check the solutions below.

  1. A hypothetical drug has a 12-hour half-life. Under linear first-order elimination, what fraction of a steady-state change is expressed after 36 hours?
  2. In a dose-proportional range, exposure is \(0.05D\) units when \(D\) is in milligrams. If the maximum permitted exposure increase per step is 0.75 units, what is the largest dose increment under this constraint?
  3. Why does a 90% PK settling criterion not by itself establish that the PD response has settled?
  4. In the coupled model, suppose \(k_{\mathrm{PK}}=k_{\mathrm{PD}}=k\). What is the response residual, and why is the equal-rate expression needed?

Solutions: (1) Thirty-six hours is three half-lives, so the expressed fraction is \(1-2^{-3}=0.875\), or 87.5%. (2) The local slope is 0.05 units per mg, so \(\Delta D_{\max}=0.75/0.05=15\) mg. (3) PK settling describes the exposure change; PD equilibration can be slower and is driven by the exposure trajectory. The response must be checked with an appropriate PD model or direct response calculation. (4) The residual is \((1+kt)\exp(-kt)\). The unequal-rate formula has denominator \(k_{\mathrm{PD}}-k_{\mathrm{PK}}\), which is zero when the rates are equal; taking its limit gives the stated expression.