Accounting for Elimination in the Global Mass Balance

Why Elimination Changes the Conservation Statement

In a closed perfusion model, drug can move among blood and tissues without changing the total amount inside the modeled system. Metabolism or excretion changes that accounting: parent drug leaves the system being tracked. Conservation still applies, but the balance must include the amount lost. The key distinction is between movement within the system and loss across its boundary.

For example, transfer from arterial blood into a tissue is not elimination if both blood and tissue are included in the model. By contrast, metabolism that converts parent drug to an untracked product is a loss from a parent-drug balance. Excretion of parent drug from the modeled body is also a loss. Summing compartment equations makes this distinction explicit: internal transfers cancel, while external inputs and losses remain.

Definition: Let \(A_i(t)\) be the amount of parent drug in modeled compartment \(i\) at time \(t\), and let \(A_{\mathrm{tot}}(t)=\sum_i A_i(t)\) be the total amount of parent drug in all included compartments. Let \(R_{\mathrm{in}}(t)\) be the total rate of parent-drug input across the system boundary, and let \(R_{\mathrm{elim},j}(t)\) be the rate of parent-drug loss through elimination pathway \(j\). The cumulative input and cumulative elimination by time \(t\) are the time integrals of their respective rates.

Amounts can be measured in mg and rates in mg/h, for example. The integration variable \(u\) below represents time, and the initial time is \(0\). Provided the rates are integrable over the interval, the global balance is

A_{\mathrm{tot}}(t) = A_{\mathrm{tot}}(0) + \int_0^t R_{\mathrm{in}}(u)\,du - \sum_j\int_0^t R_{\mathrm{elim},j}(u)\,du.

Equivalently, defining \(E_j(t)=\int_0^t R_{\mathrm{elim},j}(u)\,du\) and \(E_{\mathrm{tot}}(t)=\sum_j E_j(t)\), the identity is

A_{\mathrm{tot}}(0)+\int_0^t R_{\mathrm{in}}(u)\,du = A_{\mathrm{tot}}(t)+E_{\mathrm{tot}}(t).

This says that the initial amount plus all drug subsequently entering the system equals the parent drug still present plus all parent drug eliminated. It is an amount identity, not a claim that the amount remaining is constant. If the system has no input after time zero, the total amount of parent drug can only decrease when elimination is nonnegative.

Deriving the Global Balance

Assume each compartment has an amount balance: its rate of change equals rates entering it minus rates leaving it. Transfers between included compartments may be numerous, but each such transfer appears once as an outflow and once as an inflow. Elimination terms, in contrast, cross the boundary of the parent-drug system and appear only as losses.

Write compartment balances.
For each included compartment, include its external inputs, internal transfers, and any elimination from that compartment.

Sum over compartments.
Each internal transfer has equal magnitude and opposite signs in the two balances, so its contribution to the sum is zero.

Integrate over time.
The summed rate of change is the external input rate minus the total elimination rate. Integration from \(0\) to \(t\) gives the global amount identity.

For a more explicit justification, write \(J_{ik}(t)\) for the transfer rate from included compartment \(i\) to included compartment \(k\). This transfer contributes \(-J_{ik}(t)\) to the balance for \(i\) and \(+J_{ik}(t)\) to the balance for \(k\). Summing all compartment equations therefore cancels every internal transfer pair. If all external input rates sum to \(R_{\mathrm{in}}(t)\), and all elimination rates sum to \(\sum_jR_{\mathrm{elim},j}(t)\), the result is

\frac{dA_{\mathrm{tot}}(t)}{dt} = R_{\mathrm{in}}(t)-\sum_jR_{\mathrm{elim},j}(t).

Integrating both sides from \(0\) to \(t\) gives

A_{\mathrm{tot}}(t)-A_{\mathrm{tot}}(0) = \int_0^t R_{\mathrm{in}}(u)\,du - \sum_j\int_0^tR_{\mathrm{elim},j}(u)\,du.

Rearranging yields the stated identity. This derivation assumes that the amounts and rates refer to the same tracked drug species and system boundary, that internal transfers are represented consistently in both connected compartments, and that no unlisted source or sink crosses that boundary. A dose administered into an included compartment is an external input; transfer from one included compartment to another is not.

Worked Example: One Elimination Pathway During an Infusion

Hypothetical data. A modeled system contains 10.0 mg of parent drug initially. A constant infusion supplies 2.00 mg/h for 3.00 h, and the total parent-drug elimination rate is 1.00 mg/h throughout this interval. Find the cumulative input, cumulative elimination, and amount remaining at 3.00 h.

The cumulative input is the infusion rate multiplied by the duration:

\int_0^{3.00\ \mathrm{h}}2.00\ \mathrm{mg/h}\,du = 2.00\times3.00 = 6.00\ \mathrm{mg}.

Similarly, the cumulative elimination is

E_{\mathrm{tot}}(3.00\ \mathrm{h}) = \int_0^{3.00\ \mathrm{h}}1.00\ \mathrm{mg/h}\,du = 3.00\ \mathrm{mg}.

Substitute into the global balance:

A_{\mathrm{tot}}(3.00\ \mathrm{h}) = 10.0\ \mathrm{mg}+6.00\ \mathrm{mg}-3.00\ \mathrm{mg} = 13.0\ \mathrm{mg}.

The result is consistent with the stated rates: input exceeds elimination by 1.00 mg/h, so the total amount increases by 3.00 mg over 3.00 h. This calculation checks the global accounting; by itself, it does not determine how the 13.0 mg is distributed among compartments.

Elimination Rates and Cumulative Amounts

The global identity does not require a particular elimination model. It requires only that the elimination rate be defined on the same amount basis and integrated over time. For instance, a clearance-based parent-drug elimination rate can be written \(CL\,C\), where \(CL\) is clearance in volume per time and \(C\) is the relevant concentration in amount per volume. Their product has units of amount per time. A first-order pathway can instead be written \(kA\), where \(k\) has units of inverse time and \(A\) is the amount on which that pathway acts.

These expressions describe possible rates, not additional terms to add automatically. If a specified clearance already represents total systemic elimination, adding separate metabolic and excretory rates on top of it would count some loss twice. Each distinct pathway should be included once, and only if it removes parent drug from the chosen accounting system.

Worked Example: Blood Elimination Alongside Tissue–Blood Exchange

Hypothetical data. A model includes one tissue and its blood pool. Their initial amounts are 4.00 mg and 6.00 mg, respectively. During 2.00 h, an external infusion supplies 3.00 mg/h. Parent drug is eliminated from blood at a constant rate of 1.00 mg/h. Tissue–blood exchange occurs, but both compartments are included in the model. Calculate the final total amount.

The initial total is

A_{\mathrm{tot}}(0)=4.00\ \mathrm{mg}+6.00\ \mathrm{mg}=10.0\ \mathrm{mg}.

The external input over 2.00 h is \(3.00\times2.00=6.00\) mg. The eliminated amount is \(1.00\times2.00=2.00\) mg. Internal tissue–blood exchange contributes neither to external input nor to elimination because the transferred drug remains inside the modeled system. Thus,

A_{\mathrm{tot}}(2.00\ \mathrm{h}) = 10.0\ \mathrm{mg}+6.00\ \mathrm{mg}-2.00\ \mathrm{mg} = 14.0\ \mathrm{mg}.

The two compartment amounts at 2.00 h are not determined by this global calculation; only their sum is. Their individual values would require the compartment-specific model and its rates.

Worked Example: Separating Metabolic and Excretory Loss

Hypothetical data. A parent-drug system starts with 20.0 mg and has no further input. Metabolic loss and excretory loss are parallel first-order pathways with rates \(0.100A_{\mathrm{tot}}(t)\) and \(0.0500A_{\mathrm{tot}}(t)\), respectively, where the rate constants are in h\(^{-1}\). Find the amounts eliminated by each pathway and the parent drug remaining after 4.00 h.

The total loss rate is \(0.150A_{\mathrm{tot}}(t)\), so the parent amount obeys \(dA_{\mathrm{tot}}/dt=-0.150A_{\mathrm{tot}}\). Its solution over 4.00 h is

A_{\mathrm{tot}}(4.00\ \mathrm{h}) = 20.0e^{-(0.150)(4.00)} = 20.0e^{-0.600} = 10.98\ \mathrm{mg}.

For parallel first-order pathways acting on the same parent amount, the metabolic fraction of the total elimination rate is \(0.100/0.150=2/3\), and the excretory fraction is \(0.0500/0.150=1/3\). The total eliminated amount is

E_{\mathrm{tot}}(4.00\ \mathrm{h}) = 20.0-10.98 = 9.02\ \mathrm{mg}.

Therefore,

E_{\mathrm{met}}(4.00\ \mathrm{h}) = \frac{2}{3}(9.02) = 6.02\ \mathrm{mg}, \qquad E_{\mathrm{excr}}(4.00\ \mathrm{h}) = \frac{1}{3}(9.02) = 3.01\ \mathrm{mg}.

The pathway amounts sum to 9.03 mg using the displayed rounded values; before rounding, they sum to the total eliminated amount. The mass balance is \(20.0\ \mathrm{mg}=10.98\ \mathrm{mg}+9.02\ \mathrm{mg}\), with no input.

Choosing the Accounting Boundary

The phrase “eliminated amount” depends on what the model tracks. If the balance tracks parent drug only, conversion of parent drug into an untracked metabolite is a metabolic loss. If a metabolite is explicitly included as a state, the conversion is instead an internal transfer between tracked species; a separate loss term is needed only when that metabolite leaves the expanded system, such as through excretion. The accounting boundary must be stated before assigning a process to input, internal transfer, or elimination.

A useful alternative to integrating elimination after the fact is to add cumulative elimination as a state. For each pathway, define \(dE_j/dt=R_{\mathrm{elim},j}(t)\) and \(E_j(0)=0\). Then the global balance can be checked at every time as

A_{\mathrm{tot}}(t)+\sum_jE_j(t) = A_{\mathrm{tot}}(0)+\int_0^tR_{\mathrm{in}}(u)\,du.

This cumulative-state formulation is algebraically equivalent to integrating the rate separately. It can make pathway-specific accounting convenient, while a direct residual balance may be simpler when only the total loss is needed. Neither formulation repairs an incorrectly specified compartment equation or an inconsistent system boundary.

Diagnostics, Pitfalls, and Limitations

A global balance is a high-level accounting check, not a substitute for the full compartment model. If calculated amounts imply that cumulative elimination is negative despite all modeled elimination rates being nonnegative, or that drug remaining exceeds initial amount plus cumulative input, investigate the definitions and calculations. Also check units: an amount rate must be integrated over time to produce an amount. Mixing, for example, a rate in mg/h with time recorded in minutes requires a unit conversion.

Counting internal flow as elimination. A tissue-to-blood transfer is not a loss when both compartments are included. It cancels in the summed balance.
Omitting an external sink. If parent drug leaves through a pathway not represented in the elimination sum, the balance will not close.
Double counting total and component clearance. Do not add a total elimination term to pathway-specific terms if the total already includes those pathways.
Changing the tracked species without changing the boundary. Conversion to an untracked metabolite is a parent-drug loss; conversion between two tracked species is an internal transfer in the expanded balance.
Expecting the global balance to identify each compartment amount. It constrains the sum of amounts, not their distribution. Compartment-specific equations are still needed for that.

The identity also assumes that all relevant sources and sinks are represented. If an input is missing, a pathway is unmeasured, or the amounts combine incompatible species, apparent failure of the balance may reflect the accounting definition rather than the arithmetic. Conversely, a balance that closes does not establish that the assumed elimination mechanism or distribution among tissues is physiologically correct; different incorrect models can satisfy the same total identity.

Key takeaway: For a parent-drug system, initial amount plus cumulative external input equals parent drug remaining plus cumulative parent-drug elimination. Internal exchanges cancel only when both sides of each exchange are included in the same system.

Check Your Understanding

Answer each question, then compare with the solutions below.

A system starts with 8.0 mg, receives 5.0 mg, and eliminates 2.0 mg over an interval. What amount remains?
Why does transfer from an included tissue to an included blood pool not appear as a global elimination term?
A pathway has rate \(CL\,C\), with \(CL=0.40\) L/h and \(C=3.0\) mg/L maintained for 2.0 h. What amount is eliminated through this pathway?
If a parent-drug model tracks a metabolite as a state, how should conversion of parent drug to that metabolite be represented in the expanded balance?

Solutions. (1) The balance gives \(8.0+5.0-2.0=11.0\) mg remaining. (2) It is an internal transfer: the same rate is subtracted from the tissue and added to blood, so the terms cancel when the compartment equations are summed. (3) The rate is \(0.40\times3.0=1.2\) mg/h, so the eliminated amount is \(1.2\times2.0=2.4\) mg. (4) It is an internal transfer between tracked species, not a loss from the expanded system; elimination is counted when drug leaves that expanded system.