Drug Dosing Calculator: One- and Two-Compartment Pharmacokinetics

Simulate repeated oral doses or IV boluses with one- and two-compartment pharmacokinetics. Compare accumulation, steady-state peaks, troughs and clearance scenarios.

#accumulation #IV bolus #one-compartment model #oral absorption #peak-trough #pharmacokinetics #two-compartment model

Caution

Do not use this simulator to choose or change a medicine dose. Use a clinician-validated dosing method and measured patient data for treatment decisions.

This drug dosing regimen calculator explores how repeated equal doses produce accumulation, peaks and troughs in linear pharmacokinetic models. It is intended for engineering students, biomedical modelers and readers studying compartment models. It provides a reproducible way to compare a one-compartment pharmacokinetic model with a two-compartment model, and to see why a concentration after several doses can still differ from the eventual steady-state profile.

Choose the disposition model and either instantaneous IV bolus administration or oral administration with first-order absorption. Enter the dose amount, dosing interval and number of doses. The first dose occurs at time zero; subsequent doses occur at equal intervals. The displayed regimen ends one interval after the last administered dose. Steady-state results describe continued repetition of the same regimen indefinitely, even when the finite simulation contains only one dose.

For one compartment, enter volume of distribution and systemic clearance. For two compartments, enter central volume, peripheral volume and intercompartmental clearance. The model allows elimination only from the central compartment. The oral depot is separate from these disposition compartments: choosing two compartments with oral dosing therefore introduces three drug-amount states. Oral dosing also requires an absorption rate and a bioavailability fraction. If parameters are already apparent values such as CL/F and V/F, use F = 1 consistently rather than correcting for bioavailability again.

Read the steady-state peak, predose trough, mean concentration and area under the curve per interval together. Peak accumulation compares the steady-state peak with the maximum reached during the first dosing interval; those peaks need not occur at the same time. The interval comparison chart separates the first interval, the final simulated interval and steady state. The terminal half-life includes slow absorption when it controls the observed tail. Use the play button to follow the concentration through the finite regimen. Pharmacokinetic units remain mg, L, h and mg/L in both unit-system views.

The clearance multiplier offers a simple sensitivity study. For example, a value of 0.5 represents a persistent halving of clearance from the start of dosing while all other parameters remain fixed. It does not identify interacting drugs, estimate inhibition strength or represent an interaction developing during treatment. The comparison-band defaults are illustrative reference concentrations, not established therapeutic or toxicity thresholds.

The tool uses analytic exponential responses and superposition rather than time-stepped Euler integration. It does not fit patient measurements, optimize prescriptions, model missed or unequal doses, loading doses, infusions, delayed absorption, nonlinear elimination, active metabolites, protein binding, organ-specific physiology or pharmacodynamic effects. It also excludes population variability and changing clearance. The default values describe an illustrative system, not a named medicine. Clinical treatment requires a drug-specific validated model, appropriate patient information and professional interpretation; changing a medicine dose from these curves can cause harm.

How it is calculated

The model follows linear mass balances. For the central and peripheral amounts, X˙c=I−(k10+k12)Xc+k21Xp\dot X_c=I-(k_{10}+k_{12})X_c+k_{21}X_p and X˙p=k12Xc−k21Xp\dot X_p=k_{12}X_c-k_{21}X_p, with C=Xc/VcC=X_c/V_c, k10=CL/Vck_{10}=CL/V_c, k12=Q/Vck_{12}=Q/V_c and k21=Q/Vpk_{21}=Q/V_p. In one-compartment mode, the exchange terms vanish. Oral input has a depot with X˙g=−kaXg\dot X_g=-k_aX_g and I=kaXgI=k_aX_g; each oral dose adds FDFD to this available-dose depot.

Bourne, Two Compartment Body Model (2014), PDF pp. 3–4, gives the eigenvalue relations α+β=k10+k12+k21\alpha+\beta=k_{10}+k_{12}+k_{21} and αβ=k10k21\alpha\beta=k_{10}k_{21}. The IV response is D/Vc [wae−αt+wbe−βt]D/V_c\,[w_ae^{-\alpha t}+w_be^{-\beta t}], where wa=(α−k21)/(α−β)w_a=(\alpha-k_{21})/(\alpha-\beta) and wb=1−waw_b=1-w_a. One-compartment mode uses one exponential with rate CL/VCL/V. Oral responses convolve each mode with first-order absorption. Equal rates use the continuous limit kate−katk_at e^{-k_at}; a small-argument exponential series avoids cancellation.

Responses are summed for each administered dose. At steady state, geometric sums replace the finite sums (Bourne, Multiple IV Bolus Dose Administration, 2014, PDF pp. 5–8; Introducing Mathematical Biology, repeated oral doses). A logarithmic time scan brackets each peak before bounded minimization refines it. Trough means the concentration immediately before the next dose. The first-interval maximum is searched only between zero and the next scheduled dose.

Conservation gives AUCss,τ=FD/CLAUC_{ss,\tau}=FD/CL and Cavg,ss=FD/(CLτ)C_{avg,ss}=FD/(CL\tau), with F=1F=1 for IV boluses. The clearance multiplier applies for the entire regimen. The observed terminal half-life is ln⁡(2)/min⁡(ka,β)\ln(2)/\min(k_a,\beta) for oral dosing and ln⁡(2)/β\ln(2)/\beta for IV dosing. Chart intervals include both sides of each IV jump.

Assumptions

Equal positive doses begin at time zero, with no pre-existing drug and no loading dose. Parameters remain constant; kinetics are linear and time invariant. Each disposition compartment is well mixed. Distribution is reversible through a common intercompartmental clearance; elimination occurs only centrally. Bioavailability and first-order absorption are constant, with no lag.

Defaults are illustrative: central volume 40 L, clearance 3 L/h and absorption rate 1.2 h⁻¹ follow Table 4 of WiredWhite post 44027; the 200 mg amount, 12 h interval, 60 L peripheral volume and 8 L/h distribution clearance are author-selected scenarios. No default is a drug-specific recommendation. The comparison lines of 1 and 8 mg/L are illustrative values used in that article.

The schematic is conceptual; its sizes do not represent anatomical volumes. Its small concentration curve is normalized to the plotted scale. Playback changes only the cursor, concentration readout and diagram.

Where it applies

For educational and model-development comparisons of fixed-parameter, one- or two-compartment linear PK with equal oral doses or instantaneous IV boluses. No prescription optimization or clinical safety classification is performed. A comparison band does not establish therapeutic efficacy or toxicity. Final-interval values are not automatically steady-state values.

The terminal disposition half-life differs from the central elimination half-life in a two-compartment model. Slow oral absorption can determine the observable terminal tail (flip-flop kinetics). Apparent parameters (CL/F, V/F and Q/F) must be used consistently with F = 1.

The manual clearance scenario is inspired by the sensitivity questions in WiredWhite post 48377; its multi-drug inhibition, binding and pharmacodynamic models are outside this tool. There is no patient variability, uncertainty interval, parameter estimation, lag, infusion, nonlinear elimination, changing clearance or pharmacodynamic prediction. Use a validated drug-specific method for clinical interpretation.

Source consistency: WiredWhite post 44027, Results and Discussion, describes “peaks exceeding the maximum safe concentration (Cmax)”. Its Tables 1 and 4 specify six 5 mg doses, V = 40 L and a 1–8 mg/L comparison band. Even retaining all six doses centrally would give only 30/40 = 0.75 mg/L, so that narrative cannot follow the listed regimen. This tool uses the stated PK parameters as illustrations and does not treat those plots or concentration claims as validation.

References

  • [WiredWhite, Computational Optimization of Drug Dosing Regimens Through Pharmacokinetic-Pharmacodynamic Modeling Using Matlab, post 44027](https://wiredwhite.com/drug-dosing-optimization-matlab/). Table 4 supplies the illustrative one-compartment PK parameters; this tool uses analytic responses rather than the article's Euler integrator.
  • [WiredWhite, Optimizing Multi-Drug Therapy: Drug Interaction Simulation Tools in Python, post 48377](https://wiredwhite.com/drug-interaction-simulation-python/). Background for clearance sensitivity; no drug-specific interaction factors are taken from this article.
  • [David W. A. Bourne, PHAR 7633 Chapter 14, Multiple IV Bolus Dose Administration, 2014](https://www.boomer.org/c/p4/c14/c14.pdf). PDF pp. 3–8: accumulation, finite geometric sums and steady-state peak–trough examples.
  • [David W. A. Bourne, PHAR 7633 Chapter 19, Multi-Compartment Pharmacokinetic Models, 2014](https://boomer.org/c/p4/c19/c19.pdf). PDF pp. 3–4: two-compartment mass balance and worked 500 mg bolus example.
  • [Introducing Mathematical Biology, Single and repeated oral doses](https://cuny.manifoldapp.org/read/introducing-mathematical-biology/section/1e98ad1a-6c1a-4f2e-a6a3-e76b53e99cf7). First-order oral superposition and worked repeated-dose trough and mean concentrations.

Inputs

InputSymbolUnitDefault
Disposition model
The oral absorption depot is additional to the selected disposition compartments.
modelTwo compartments
Dose routerouteOral, first-order absorption
Amount per dose
Illustrative amount, without a named drug or patient recommendation.
D_dosemg200
Dosing intervaldtauh12
Number of doses
Equal doses at 0, interval, 2 × interval, and so on. Plot ends one interval after the last dose.
N_dose6
Central volume (V for one compartment)VcL40
Baseline systemic clearanceCLL/h3
Peripheral volume
An illustrative second-compartment volume; not a fitted patient value.
VpL60
Intercompartmental clearance
Bidirectional distribution clearance Q; k₁₂ = Q/Vc and k₂₁ = Q/Vp.
Q_distL/h8
Oral bioavailability fraction
Use F = 1 for an apparent-parameter model with V/F and CL/F; do not apply F twice.
F_bio1
First-order absorption rate
No absorption lag. Equal absorption and disposition rates are supported.
kah^-11.2
Clearance multiplier
1 leaves clearance unchanged; 0.5 halves it. This is a user-defined sensitivity scenario, not an interaction prediction.
f_cl1
Lower comparison concentration
Illustrative reference line only; not a validated therapeutic limit.
C_lowermg/L1
Upper comparison concentration
Enter drug-specific evidence separately. These lines do not establish safety or efficacy.
C_uppermg/L8
Time cursor ▶
Moves the marker on the schematic; concentrations and chart curves are computed independently of playback.
t_playh26.4

Results

With the default inputs above, the calc gives:

ResultValue
Steady-state peak7.3102 mg/L
Steady-state trough4.224 mg/L
Steady-state average5.5556 mg/L
Peak accumulation factor2.2267
Peak time after dose at steady state1.557 h
Steady-state AUC per interval66.6667 mg·h/L
Final-interval peak6.4773 mg/L
Final-interval predose concentration3.59 mg/L
Terminal half-life (including absorption)26.4898 h
Exposure / baseline-clearance exposure1

Checked against sources

Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).

IV doses every one half-life

Source: Bourne, PHAR 7633 Chapter 14 (2014), PDF pp. 3–4, 100 mg every 6 h, V = 25 L and half-life = 6 h. Published example gives first peak 4 mg/L, third interval 7 to 3.5 mg/L, and steady state 8 to 4 mg/L.

Inputs: Disposition model = One compartment; Dose route = Instantaneous IV bolus; Amount per dose = 100 mg; Central volume (V for one compartment) = 25 L; Baseline systemic clearance = 2.888 L/h; Dosing interval = 6 h; Number of doses = 3; Time cursor ▶ = 6 h; the others as above.

ResultSourceThis calc
Peak in the first dosing interval4 mg/L4 mg/L
Steady-state peak8 mg/L8 mg/L
Steady-state trough4 mg/L4 mg/L
Peak accumulation factor22
Final-interval peak7 mg/L7 mg/L
Final-interval predose concentration3.5 mg/L3.5 mg/L

Published IV plateau: 4 h half-life, 6 h interval

Source: Bourne, PHAR 7633 Chapter 14 (2014), PDF pp. 7–8: 100 mg every 6 hours, V = 10 L and 4 h half-life; published rounded peak 15.5 and trough 5.48 mg/L.

Inputs: Disposition model = One compartment; Dose route = Instantaneous IV bolus; Amount per dose = 100 mg; Central volume (V for one compartment) = 10 L; Baseline systemic clearance = 1.733 L/h; Dosing interval = 6 h; Time cursor ▶ = 6 h; the others as above.

ResultSourceThis calc
Steady-state peak15.5 ± 0.04 mg/L15.469 mg/L
Steady-state trough5.48 ± 0.015 mg/L5.469 mg/L

Two-compartment 500 mg bolus at 1.5 h

Source: Bourne, PHAR 7633 Chapter 19 (2014), PDF p. 4. Published Vc = 90.5 kg × 0.39 L/kg = 35.295 L; k10 = 0.192, k12 = 1.86, k21 = 1.68 h⁻¹. Inputs CL = k10 Vc, Q = k12 Vc, Vp = Q/k21 express that same published case in clearance form. Expected concentration 5.59 mg/L is the published answer. Its last symbolic line writes alpha and beta as concentration coefficients, while its numerical substitution correctly uses A and B; this tool uses A and B.

Inputs: Dose route = Instantaneous IV bolus; Amount per dose = 500 mg; Central volume (V for one compartment) = 35.295 L; Baseline systemic clearance = 6.777 L/h; Intercompartmental clearance = 65.649 L/h; Peripheral volume = 39.077 L; Number of doses = 1; Time cursor ▶ = 1.5 h; the others as above.

ResultSourceThis calc
Fast disposition rate3.643 h⁻¹3.643 h⁻¹
Terminal disposition rate0.08853 h⁻¹0.08853 h⁻¹
Fast available-dose impulse coefficient7.824 ± 0.01 mg/L7.824 mg/L
Slow available-dose impulse coefficient6.343 ± 0.01 mg/L6.342 mg/L
Concentration at the play cursor5.59 ± 0.015 mg/L5.586 mg/L

Published oral trough and average concentration

Source: Introducing Mathematical Biology, Single and repeated oral doses, Example problems: k = 0.25 h⁻¹, a = 0.4 h⁻¹, V = 25 L, 400 mg every 6 h. Published trough is 8.0 mg/L and mean is 10.67 mg/L. The text explicitly identifies its 13.3 mg/L peak estimate as an overestimate; it is not used as an exact peak reference.

Inputs: Disposition model = One compartment; Amount per dose = 400 mg; Central volume (V for one compartment) = 25 L; Baseline systemic clearance = 6.25 L/h; First-order absorption rate = 0.4 h⁻¹; Dosing interval = 6 h; Time cursor ▶ = 6 h; the others as above.

ResultSourceThis calc
Steady-state trough8 ± 0.05 mg/L7.998 mg/L
Steady-state average10.67 mg/L10.667 mg/L

Regression: oral peak just before the final probe

Source: Independent reviewer regression supplied by the user: for D = 200 mg, V = 40 L, ka = CL/V = 1 h⁻¹, the exact single-dose equal-rate limit is C(t) = 5 t exp(−t) mg/L. Its derivative is 5(1−t) exp(−t), so the peak is at t = 1 h with C = 5/e = 1.8393972058572117 mg/L, inside the 1.09 h interval. The reviewer gives the corrected accumulation factor as 2.61667. These are analytic regression values, not a new published worked example; the general oral model is described in Introducing Mathematical Biology, https://cuny.manifoldapp.org/read/introducing-mathematical-biology/section/1e98ad1a-6c1a-4f2e-a6a3-e76b53e99cf7.

Inputs: Disposition model = One compartment; Baseline systemic clearance = 40 L/h; First-order absorption rate = 1 h⁻¹; Dosing interval = 1.09 h; Number of doses = 1; Time cursor ▶ = 1 h; the others as above.

ResultSourceThis calc
Time of peak in the first interval1 h1 h
Time of peak in the final interval1 h1 h
Peak in the first dosing interval1.839 mg/L1.839 mg/L
Final-interval peak1.839 mg/L1.839 mg/L
Peak accumulation factor2.6172.617

Related Resources

Responses

Your email address will not be published. Required fields are marked *

L ading...