Why Dosing Schedule Matters
Two regimens can produce similar drug exposure over a day yet different responses over time. Continuous exposure may sustain an adaptation process that reduces drug effect, while drug-free intervals allow that adaptation to recede. The size of this benefit depends on how quickly adaptation forms and recovers, as well as on the exposure pattern itself. A schedule comparison is therefore a dynamic prediction, not a rule that intermittent dosing is always preferable.
This tutorial uses a compact state-variable model to simulate that comparison. The adaptation state represents a fractional loss of response; it is not itself a measured receptor concentration or a direct clinical endpoint. Earlier material on state-variable models of pharmacodynamic tolerance and exposure-driven adaptation dynamics motivates this model structure. Here the emphasis is on calculating its consequences across dosing schedules.
The exposure input may be defined from concentration, for example \(u(t)=C(t)/(EC_{50}+C(t))\), where \(C(t)\) is concentration and \(EC_{50}\) is the concentration giving half-maximal input. In the worked schedule calculations below, exposure is idealized as a square wave: \(u=1\) while dosing and \(u=0\) during a holiday. This isolates adaptation dynamics from concentration changes. It assumes that concentration rises and falls rapidly relative to the schedule intervals, and that the active-dose input is the same from dose to dose.
The state remains in its stated range if it starts there. At \(A=0\), the derivative is \(k_{\mathrm{form}}u\), which is nonnegative; at \(A=1\), it is \(-k_{\mathrm{rec}}\), which is negative. Thus the model does not predict adaptation below zero or above one. The response multiplier \(M\) is correspondingly between zero and one.
Solutions Within Dosing and Holiday Intervals
When \(u\) is constant, the adaptation state approaches a schedule-specific equilibrium. Define \(\lambda=k_{\mathrm{form}}u+k_{\mathrm{rec}}\). Setting the derivative to zero gives the equilibrium \(A_{\mathrm{eq}}=k_{\mathrm{form}}u/\lambda\). Subtracting this equilibrium from the differential equation shows how quickly the state approaches it:
This is the exact update over an interval of length \(\Delta\) with constant input. During a holiday, \(u=0\), so \(A_{\mathrm{eq}}=0\) and \(\lambda=k_{\mathrm{rec}}\). Therefore, after \(h\) drug-free time units, the state is multiplied by \(e^{-k_{\mathrm{rec}}h}\). A holiday recovers only the fraction allowed by its duration and the recovery rate; it does not automatically reset adaptation to zero.
Worked Example: Adaptation During Continuous Exposure
Suppose \(k_{\mathrm{form}}=0.12\ \mathrm{h}^{-1}\), \(k_{\mathrm{rec}}=0.08\ \mathrm{h}^{-1}\), and \(u=1\) continuously. Assume \(A(0)=0\). Calculate the adaptation state and response multiplier after 8 hours, then find their long-run values.
The equilibrium and approach rate are $$ A_{\mathrm{eq}}=\frac{0.12}{0.12+0.08}=0.600, \qquad \lambda=0.200\ \mathrm{h}^{-1}. $$ Using the constant-input solution at 8 hours, $$ A(8)=0.600+(0-0.600)e^{-(0.200)(8)} =0.600\bigl(1-e^{-1.6}\bigr) \approx0.4789. $$ The response multiplier is \(M(8)=1-0.4789=0.5211\). Thus, under these model assumptions, the response at the same exposure is about 52.11% of the response without adaptation after 8 hours.
As time tends to infinity, \(A(t)\) approaches \(0.600\), and \(M(t)\) approaches \(0.400\). The state approaches but does not exceed this equilibrium under continuous exposure. This calculation is conditional on the specified rate constants and the idealized constant exposure input.
Intermittent Dosing and Periodic Adaptation
With repeated on/off dosing, adaptation generally does not return to zero before each dose. Instead, after enough cycles it approaches a periodic pattern: the state at the start of a dosing interval repeats from one cycle to the next. This periodic state is useful for comparing the response in an established regimen with the response during the first dose, when the initial adaptation may be lower.
For an on interval of length \(t_{\mathrm{on}}\) at constant input \(u_{\mathrm{on}}\), write \(A_{\mathrm{on,eq}}=k_{\mathrm{form}}u_{\mathrm{on}}/(k_{\mathrm{form}}u_{\mathrm{on}}+k_{\mathrm{rec}})\) and \(\lambda_{\mathrm{on}}=k_{\mathrm{form}}u_{\mathrm{on}}+k_{\mathrm{rec}}\). If \(A_0\) is the state at the start of dosing, the end-of-dose state is
For an off interval of length \(t_{\mathrm{off}}\), the next cycle begins at \(A_{\mathrm{next}}=A_1e^{-k_{\mathrm{rec}}t_{\mathrm{off}}}\). At the periodic state, \(A_{\mathrm{next}}=A_0\). Solving this recurrence gives the start-of-dose adaptation:
The denominator is positive because \(0<r<1\) and \(0<q<1\), so \(rq<1\). This periodic solution follows directly by substituting \(A_1=A_{\mathrm{on,eq}}(1-q)+qA_0\) into \(A_0=rA_1\) and rearranging. It describes the recurring state only after the regimen has approached periodic behavior; it is not generally the state on the first dose.
Worked Example: A Periodic 4-Hour-On, 8-Hour-Off Schedule
Use \(k_{\mathrm{form}}=0.12\ \mathrm{h}^{-1}\), \(k_{\mathrm{rec}}=0.08\ \mathrm{h}^{-1}\), and \(u=1\) during each 4-hour dosing interval. Set \(u=0\) during each 8-hour holiday. Find the periodic adaptation at the start and end of dosing, and the average response multiplier during dosing.
The on-interval equilibrium is \(A_{\mathrm{on,eq}}=0.600\), and \(\lambda_{\mathrm{on}}=0.200\ \mathrm{h}^{-1}\). Thus $$ q=e^{-(0.200)(4)}=e^{-0.8}\approx0.4493, \qquad r=e^{-(0.08)(8)}=e^{-0.64}\approx0.5273. $$ The periodic start-of-dose state is $$ A_0=\frac{(0.5273)(0.600)(1-0.4493)}{1-(0.5273)(0.4493)} \approx0.2283. $$ At the end of the 4-hour dosing interval, $$ A_1=0.600+(0.2283-0.600)(0.4493)\approx0.4330. $$ The state then decays during the 8-hour holiday: \(A_1r\approx(0.4330)(0.5273)=0.2283\), which returns to the calculated start-of-dose state, as required for a periodic cycle.
During dosing, \(A(t)=0.600+(A_0-0.600)e^{-0.200t}\). The average adaptation over the 4-hour on interval is $$ \bar{A}_{\mathrm{on}} =\frac{1}{4}\int_0^4 A(t)\,dt =0.600+(0.2283-0.600)\frac{1-e^{-0.8}}{(0.200)(4)} \approx0.3442. $$ Therefore, the average response multiplier during dosing is \(1-\bar{A}_{\mathrm{on}}\approx0.6558\). This is the average response relative to the unadapted response at the same on-dose exposure; it is not the average response over the full 12-hour cycle.
Drug Holidays and Recovery on Resumption
A holiday can be evaluated by propagating the adaptation state through the off interval, then calculating the response multiplier when exposure resumes. This makes the distinction between recovery and complete reversal explicit. With the fractional-loss model, the immediate response multiplier on resumption is \(1-A_{\mathrm{holiday}}\), provided the exposure input is again the same as the reference exposure.
Worked Example: A 12-Hour Holiday
Suppose continuous exposure has brought the state to \(A(0)=0.600\), and the drug is then stopped for 12 hours. Let \(k_{\mathrm{rec}}=0.08\ \mathrm{h}^{-1}\). Find the state at the end of the holiday and the response multiplier when the same exposure resumes.
During the holiday, \(u=0\), so $$ A(12)=0.600e^{-(0.08)(12)} =0.600e^{-0.96} \approx0.2297. $$ When exposure resumes, the initial response multiplier is $$ M=1-A(12)=1-0.2297\approx0.7703. $$ Thus this model predicts recovery from an adapted multiplier of \(1-0.600=0.400\) to approximately \(0.7703\) at the beginning of resumed exposure. It does not predict an immediate return to the unadapted multiplier of 1.000, because adaptation remains at about 0.2297.
Comparing Schedules and Checking the Simulation
Schedule comparisons require a clearly stated comparison target. The continuous-exposure example has a long-run response multiplier of 0.400, while the periodic intermittent example has an average on-dose multiplier of about 0.6558. These numbers describe response preservation during active exposure. They do not establish that intermittent dosing produces more total therapeutic response, because the intermittent schedule has fewer hours of exposure. A meaningful comparison may require matching cumulative exposure, dose intensity, or another clinically relevant constraint.
| Schedule feature | Model calculation | What it summarizes |
|---|---|---|
| Continuous exposure | Approach \(A_{\mathrm{eq}}\) at the constant-input rate \(\lambda\) | Long-run adaptation and response under sustained input |
| Repeated on/off dosing | Apply the on and off updates in sequence until the state repeats | Adaptation range and response during an established dosing cycle |
| Drug holiday | Multiply the starting state by \(e^{-k_{\mathrm{rec}}h}\) | Residual adaptation and response at resumption |
For a numerical simulation with changing exposure, calculate or interpolate \(u(t)\) from the PK profile, then integrate the adaptation differential equation alongside the relevant exposure trajectory. The exact constant-input updates above are preferable within intervals where \(u\) is genuinely constant; a numerical solver may be needed when concentration varies continuously. As a basic check, a simulation should respect \(0\leq A(t)\leq1\), approach the appropriate equilibrium under constant input, and decay during a zero-input holiday. The simulated periodic start and end states should agree across successive cycles within the chosen numerical tolerance.
A common mistake is to treat every missed dose or holiday as a complete reset, then calculate the next response as if \(A=0\). The correction is to propagate the state through the actual off-time using \(A(t+h)=A(t)e^{-k_{\mathrm{rec}}h}\); the recovery rate and holiday duration jointly determine the residual state. Another mistake is to compare response multipliers during active dosing and call the larger one a superior overall regimen. A larger multiplier can indicate response preservation during exposure, but it says nothing by itself about total benefit, safety, adherence, or the consequences of lower exposure.
The square-wave input is a useful teaching and screening approximation, but its limitations matter. Real concentration may persist after the last dose, so a nominal holiday may still have \(u(t)>0\); in that case adaptation does not follow pure exponential recovery at rate \(k_{\mathrm{rec}}\). The response multiplier \(1-A\) also assumes adaptation acts as a proportional loss. Other PD structures may represent adaptation as a potency shift, an efficacy change, or a separate opposing response. These alternatives can yield different schedule predictions even when the adaptation trajectory is similar. Schedule simulation supports a conditional model-based comparison; it cannot, by itself, establish which mechanism is operating or prove that a regimen is clinically preferable.
Check Your Understanding
Answer each question, then check the solution immediately below it.
- A drug-free interval begins at \(A=0.500\). If \(k_{\mathrm{rec}}=0.10\ \mathrm{h}^{-1}\), what is \(A\) after 5 hours?
- For \(k_{\mathrm{form}}=0.12\ \mathrm{h}^{-1}\), \(k_{\mathrm{rec}}=0.08\ \mathrm{h}^{-1}\), and constant \(u=0.5\), find the equilibrium adaptation state and the response multiplier at equilibrium.
- Why is the first intermittent dose not necessarily representative of the average response during an established periodic regimen?
- A schedule simulation predicts \(A<0\) during a long holiday. Is this consistent with the model, and what should be checked?
Solutions. (1) During a holiday, \(A(5)=0.500e^{-(0.10)(5)}=0.500e^{-0.5}\approx0.3033\). (2) The equilibrium is \(A_{\mathrm{eq}}=(0.12)(0.5)/[(0.12)(0.5)+0.08]=0.060/0.140\approx0.4286\). The equilibrium response multiplier is \(1-0.4286=0.5714\). (3) A first dose may begin with little or no adaptation, whereas repeated cycles approach a nonzero periodic start-of-dose state. The within-dose adaptation trajectory, and hence its average response, can therefore differ. (4) No. The stated model preserves \(A\geq0\); a negative simulated state suggests a numerical or implementation problem, such as an unstable integration step or incorrect update equation. With a zero-input holiday, the exact update is \(A(t+h)=A(t)e^{-k_{\mathrm{rec}}h}\), which remains nonnegative.