Calculating Accumulation and Peak–Trough Fluctuation

Why Accumulation and Fluctuation Are Different

Repeated doses do not simply reproduce the concentration profile of a single dose. Before one dose has been fully eliminated, the next dose adds to the residual drug, so concentrations build toward a repeating steady-state profile. Two useful questions are how much the peak has accumulated relative to a single dose, and how far the concentration falls from peak to trough during each dosing interval. These are related, but they are not the same quantity.

This tutorial considers repeated instantaneous intravenous bolus doses in a one-compartment model with linear, first-order elimination. It derives an accumulation ratio and several explicit peak-to-trough measures. Naming the measure matters: a peak-to-trough fold ratio is not a percentage fluctuation, and neither should be confused with the accumulation ratio.

Model and assumptions: A dose \(D\) is given as an instantaneous IV bolus every \(\tau\) time units. The apparent volume of distribution \(V\) and clearance \(CL\) are constant, and elimination is linear and first order, with \(k=CL/V\). Doses and intervals are identical, concentrations add across doses, and the regimen has continued long enough to reach steady state. There is no absorption phase, time-varying clearance, or nonlinear elimination in this model.

Concentration immediately after a single bolus is \(D/V\). Between doses it declines exponentially. The steady-state profile for repeated bolus dosing is the starting point for the calculations below: each new dose adds \(D/V\), while the concentration remaining from earlier doses decays by factors determined by \(k\) and \(\tau\).

Accumulation Ratio From a Geometric Series

Let \(q=e^{-k\tau}\), the fraction of concentration remaining after one complete dosing interval. Because \(k>0\) and \(\tau>0\), \(0<q<1\). Immediately after the \(n\)th dose, the current dose contributes \(D/V\), the preceding dose contributes \((D/V)q\), and each earlier dose contributes another factor of \(q\).

C_{\max,n}=\frac{D}{V}\sum_{j=0}^{n-1}q^j =\frac{D}{V}\frac{1-q^n}{1-q}

As \(n\) increases, \(q^n\) approaches zero. The steady-state peak is therefore \(D/V\) multiplied by the limiting geometric sum. The accumulation ratio compares this steady-state peak with the immediate post-dose concentration after one isolated dose.

Definition: For this linear repeated-bolus model, the steady-state accumulation ratio is \(R=C_{\max,ss}/(D/V)\). It is the factor by which the steady-state peak exceeds the immediate peak after a single dose. Thus \(R=1/(1-e^{-k\tau})\). At the same elapsed time after a dose, the steady-state concentration is also \(R\) times the concentration from a single dose, provided both profiles follow the same linear elimination model.

C_{\max,ss}=\frac{D/V}{1-q} \qquad R=\frac{1}{1-q}=\frac{1}{1-e^{-k\tau}}

The finite-dose expression also shows how the profile approaches steady state. The fraction of the eventual steady-state peak attained immediately after dose \(n\) is \(C_{\max,n}/C_{\max,ss}=1-q^n\). This statement concerns the idealized model and does not mean that every clinical regimen will achieve a desired concentration by a particular dose.

Peak, Trough, and Fluctuation at Steady State

At steady state, the trough immediately before the next dose is the peak multiplied by the fraction remaining through one interval. This gives the peak-to-trough fold ratio directly. A separate measure expresses the peak-to-trough difference relative to the average concentration over the interval.

C_{\min,ss}=qC_{\max,ss}

The mean concentration over a steady-state interval is the area under the concentration-time curve during that interval divided by \(\tau\). Under linear elimination, the amount eliminated over a steady-state interval equals the dose administered during it. Since the elimination rate is \(CL\) times concentration, the interval area is \(D/CL\).

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

The term “peak-to-trough fluctuation” is used with more than one convention. To avoid ambiguity, this tutorial reports both the fold ratio \(C_{\max,ss}/C_{\min,ss}\) and normalized fluctuation, defined as the peak-to-trough difference divided by the average concentration. The latter is often expressed as a percentage.

Definitions: The peak-to-trough fold ratio is \(F=C_{\max,ss}/C_{\min,ss}\). Normalized peak-to-trough fluctuation is \((C_{\max,ss}-C_{\min,ss})/C_{\mathrm{avg},ss}\); as a percentage, multiply this value by 100. These definitions use different denominators and should not be interchanged.

For this model, the main relationships follow by substitution. They show that the fold ratio depends on the fraction remaining, whereas normalized fluctuation has an especially simple expression in terms of the dimensionless product \(k\tau\).

F=\frac{C_{\max,ss}}{C_{\min,ss}}=\frac{1}{q}=e^{k\tau}

\frac{C_{\max,ss}-C_{\min,ss}}{C_{\mathrm{avg},ss}} =\frac{C_{\max,ss}(1-q)}{D/(CL\,\tau)} =\frac{D/V}{D/(CL\,\tau)} =k\tau

The last equality uses \(CL=kV\). The corresponding percentage fluctuation is \(100k\tau\%\). It can exceed 100%: the numerator is a peak-to-trough difference, not a probability or a proportion constrained to be at most one.

Worked Calculations

Worked Example: Accumulation and Concentrations for a 12-Hour Interval

Hypothetical regimen. An adult receives 100 mg as an IV bolus every 12 hours. Suppose \(V=20.0\) L and the elimination half-life is 8.00 hours. Calculate the steady-state accumulation ratio, peak, trough, average concentration, and normalized fluctuation.

First calculate \(k\) from the half-life, then calculate the fraction remaining over the 12-hour interval:

k=\frac{\ln(2)}{8.00\ \mathrm{h}}=0.08664\ \mathrm{h}^{-1} \qquad k\tau=(0.08664\ \mathrm{h}^{-1})(12.0\ \mathrm{h})=1.0397

q=e^{-1.0397}=0.3536 \qquad R=\frac{1}{1-0.3536}=1.547

The single-dose immediate peak is \(D/V=100\ \mathrm{mg}/20.0\ \mathrm{L}=5.000\ \mathrm{mg/L}\). Therefore:

C_{\max,ss}=R(D/V)=(1.547)(5.000\ \mathrm{mg/L})=7.735\ \mathrm{mg/L}

C_{\min,ss}=qC_{\max,ss}=(0.3536)(7.735\ \mathrm{mg/L})=2.735\ \mathrm{mg/L}

Clearance is \(CL=kV=(0.08664\ \mathrm{h}^{-1})(20.0\ \mathrm{L})=1.733\ \mathrm{L/h}\). Thus the average concentration and normalized fluctuation are:

C_{\mathrm{avg},ss}=\frac{100\ \mathrm{mg}}{(1.733\ \mathrm{L/h})(12.0\ \mathrm{h})} =4.809\ \mathrm{mg/L}

\frac{7.735-2.735\ \mathrm{mg/L}}{4.809\ \mathrm{mg/L}} =1.040 \approx k\tau

The normalized fluctuation is approximately 104.0%. The steady-state peak is 1.547 times the single-dose peak, while the steady-state peak-to-trough fold ratio is \(7.735/2.735=2.828\). These are three distinct summaries.

Worked Example: Interpreting a Peak-to-Trough Ratio of 2

Hypothetical profile. A repeated-bolus steady-state profile has a peak-to-trough fold ratio of 2. Under the stated model, determine the fraction remaining at the end of the interval, the accumulation ratio, and normalized fluctuation.

From \(F=1/q\), the fraction remaining is \(q=1/F=1/2=0.5000\). The accumulation ratio is then:

R=\frac{1}{1-q}=\frac{1}{1-0.5000}=2.000

Because \(q=e^{-k\tau}\), take the negative natural logarithm to obtain the dimensionless interval-elimination product:

k\tau=-\ln(0.5000)=\ln(2)=0.6931

Normalized fluctuation is therefore 0.6931, or 69.31%. The fold ratio is 2, the remaining fraction is 0.5000, and the normalized fluctuation is 69.31%; they are not interchangeable numbers. Since \(k=\ln(2)/t_{1/2}\), this particular fold ratio occurs when \(\tau=t_{1/2}\).

Worked Example: Changing the Interval Relative to Half-Life

Hypothetical comparison. Consider the same one-compartment drug with a 12.0-hour half-life, but compare dosing intervals of 6.00 and 24.0 hours. Calculate \(q\), \(R\), the peak-to-trough fold ratio, and normalized fluctuation for each interval. Dose size does not affect these dimensionless measures in the linear model.

Here \(k=\ln(2)/12.0=0.05776\ \mathrm{h}^{-1}\). For a 6.00-hour interval:

k\tau=(0.05776)(6.00)=0.3466 \qquad q=e^{-0.3466}=0.7071

R=\frac{1}{1-0.7071}=3.414 \qquad F=\frac{1}{0.7071}=1.414 \qquad 100k\tau=34.66\%

For a 24.0-hour interval, \(k\tau=(\ln 2)(24.0/12.0)=\ln(4)=1.3863\), so \(q=e^{-1.3863}=0.2500\). Therefore:

R=\frac{1}{1-0.2500}=1.333 \qquad F=\frac{1}{0.2500}=4.000 \qquad 100k\tau=138.6\%

The shorter interval leaves more drug in the body at the next dose and has greater accumulation, but its within-interval fluctuation is smaller. The longer interval has less accumulation and a larger peak-to-trough difference relative to the average. This comparison concerns dimensionless measures; absolute concentrations also depend on dose, volume, and clearance.

How Half-Life and Interval Control the Measures

The elimination half-life satisfies \(k=\ln(2)/t_{1/2}\). Substitution gives a direct relationship between the interval measured in half-lives and all three dimensionless summaries:

q=2^{-\tau/t_{1/2}} \qquad R=\frac{1}{1-2^{-\tau/t_{1/2}}} \qquad F=2^{\tau/t_{1/2}} \qquad \text{normalized fluctuation}=(\ln 2)\frac{\tau}{t_{1/2}}

For fixed half-life, increasing \(\tau\) decreases \(q\): less concentration remains until the next dose. It consequently decreases \(R\), increases \(F\), and increases normalized fluctuation. For fixed interval, increasing the half-life has the opposite effects because elimination is slower. The relevant comparison is the interval relative to half-life, not the interval alone.

If the dose per administration is held constant while the interval changes, the mean concentration \(D/(CL\tau)\) also changes. If instead a fixed total daily dose is divided among more frequent administrations, \(D/\tau\) is held constant and the mean concentration remains constant under this model, while the dimensionless fluctuation still changes with \(\tau/t_{1/2}\). These comparisons rely on linear pharmacokinetics and should not be assumed to hold unchanged when elimination is nonlinear.

Derivation Check, Diagnostics, and Limitations

The central identities can be checked without memorizing them. The peak is the geometric sum, the trough is the peak times \(q\), and the difference between peak and trough is exactly the concentration increment from one bolus:

C_{\max,ss}-C_{\min,ss} =C_{\max,ss}(1-q) =\frac{D/V}{1-q}(1-q) =\frac{D}{V}

This check explains why, at a fixed dose and volume, the absolute peak-to-trough difference is \(D/V\) in this idealized model even as the interval changes. The normalized fluctuation changes because the average concentration changes with the interval. As a second consistency check, \(F=e^{k\tau}\) must exceed 1 for positive \(k\) and \(\tau\), while \(R=1/(1-q)\) must also exceed 1.

When applying these calculations to data or model output, check that the dosing times are regular, the profile has reached a repeating pattern, and peak and trough are taken at consistent times relative to each dose. Confirm units: \(k\) must be an inverse time and \(k\tau\) must be dimensionless. A computed \(q\) outside \((0,1)\), a fold ratio at or below 1 for positive elimination, or a negative accumulation ratio signals an input, unit, or calculation error.

Departures from the expected profile may have biological as well as computational explanations. Absorption can delay the peak, multiple compartments can create distribution phases, missed or delayed doses can disrupt periodic steady state, and time-varying or nonlinear clearance can invalidate the geometric accumulation calculation. In those settings, use a model that represents the relevant process and calculate peaks and troughs from that model rather than forcing the one-compartment bolus formulas onto the data.

Common pitfall: A peak-to-trough ratio of 2 does not mean normalized fluctuation is 50%. It implies \(q=1/2\), so \(k\tau=-\ln(0.5)=\ln(2)=0.6931\), and normalized fluctuation is 69.31%. The 0.5000 value is the fraction remaining at the trough relative to the peak; the fluctuation uses the average concentration as its denominator.

The main alternative summaries answer different questions. The accumulation ratio compares a steady-state concentration with its single-dose counterpart. The fold ratio describes the peak relative to the trough. Normalized fluctuation scales the peak-to-trough difference by the average concentration, while the absolute difference retains concentration units. Reporting the definition and units prevents a numerically correct result from being misinterpreted.

Key takeaway: In a linear one-compartment repeated-bolus model, accumulation is governed by \(1/(1-e^{-k\tau})\), the peak-to-trough fold ratio by \(e^{k\tau}\), and normalized peak-to-trough fluctuation by \(k\tau\). The same interval-to-half-life relationship links all three, but each measures a different feature of the profile.

Check Your Understanding

Answer each question, then compare with the brief solutions below.

For \(k\tau=\ln(3)\), calculate \(q\), \(R\), the peak-to-trough fold ratio, and normalized fluctuation as a percentage.
A regimen has \(q=0.2500\). Calculate \(k\tau\), \(R\), and normalized fluctuation as a percentage.
For a fixed drug half-life, what happens to accumulation ratio and normalized fluctuation when the dosing interval is halved?
Why is the accumulation ratio not the same as the peak-to-trough fold ratio?
Name two assumptions whose failure could make these repeated-bolus calculations unsuitable for observed concentrations.

Solutions. (1) \(q=e^{-\ln(3)}=1/3=0.3333\); \(R=1/(1-1/3)=1.500\); \(F=1/q=3.000\); normalized fluctuation is \(100\ln(3)=109.9\%\). (2) \(k\tau=-\ln(0.2500)=\ln(4)=1.386\); \(R=1/(1-0.2500)=1.333\); normalized fluctuation is \(100(1.386)=138.6\%\). (3) The product \(k\tau\) is halved. Since \(R=1/(1-e^{-k\tau})\), the accumulation ratio increases; normalized fluctuation \(100k\tau\) decreases by half. (4) Accumulation compares the steady-state peak with a single-dose peak; the fold ratio compares the steady-state peak with the steady-state trough. (5) Examples include non-instantaneous absorption, multicompartment distribution, irregular dosing, failure to reach steady state, and nonlinear or time-varying elimination.