Deriving Regional Residence Times From Marker Measurements

Why Measure Residence Time With a Marker?

A regional residence time summarizes how long material takes to pass between two specified boundaries, such as the entrance to and exit from a gastrointestinal region. In PBPK modeling, it provides a physiological quantity that can be estimated from marker measurements rather than assumed from an arbitrary compartment rate. The estimate is meaningful only when the marker’s measured passage represents the material whose movement the model is intended to describe.

The earlier material on regional GI compartments and anatomical boundaries establishes why the boundaries must be named: changing the entry or exit location changes the region and therefore can change the measured residence time. Here the focus is the calculation from marker observations at those boundaries. A marker may be tracked as an amount present within the region, or through its rates of entry and exit. These measurements lead to equivalent calculations when mass is conserved and observations cover the complete passage.

Definition: The mean regional residence time is the amount-weighted average time that marker material spends between a specified entry boundary and exit boundary. For a finite marker passage, it can be calculated as the time integral of marker mass inside the region divided by the total marker amount entering, provided the initial and final inventories are zero and marker mass is conserved.

Inventory Area Gives Mean Residence Time

Let $M(t)$ be the marker mass inside the defined region at time $t$, in milligrams, and let $D$ be the total marker mass entering during the measurement, also in milligrams. The inventory-area estimate is

\bar{\tau}=\frac{\int_0^T M(t)\,dt}{D}.

The integral has units of milligram-hours, so division by milligrams gives hours. Intuitively, each marker particle contributes its time in the region to the accumulated inventory area: a particle remaining for longer contributes to $M(t)$ over more time. Dividing the total accumulated marker-time by the number of marker units gives their mean residence time.

For measured values at discrete times, the integral is usually approximated numerically. With observations $(t_i,M_i)$, the trapezoidal rule estimates the area between adjacent times by the average endpoint inventory multiplied by the time interval:

A \approx \sum_{i=0}^{k-1}\frac{M_i+M_{i+1}}{2}(t_{i+1}-t_i), \qquad \bar{\tau}\approx \frac{A}{D}.

This is an approximation to the area under the inventory curve, not a rule that the true curve is linear. It is most credible when sampling is frequent enough to capture important rises, peaks, and tails. The measured inventory must represent marker mass in the region. If only marker concentration is measured, concentration alone is not generally an inventory: it must be combined with an appropriate volume or other measurement that converts concentration to amount.

Worked Example: Residence Time From an Inventory Curve

In a hypothetical marker experiment, 100 mg enters a defined region at time zero. Marker mass inside the region is estimated at 100 mg at 0 h, 75 mg at 1 h, 50 mg at 2 h, 25 mg at 3 h, and 0 mg at 4 h. Use trapezoidal integration.

The four interval areas are $87.5$, $62.5$, $37.5$, and $12.5$ mg·h. For example, the first is $(100+75)(1)/2=87.5$ mg·h. Therefore,

A=87.5+62.5+37.5+12.5=200\ \text{mg}\cdot\text{h}, \qquad \bar{\tau}\approx\frac{200\ \text{mg}\cdot\text{h}}{100\ \text{mg}}=2.00\ \text{h}.

Under the stated measurements and trapezoidal approximation, the mean marker residence time is 2.00 h. The calculation assumes the observed inventory represents the complete passage, including the point at which no marker remains in the region.

Using Entry and Exit Measurements

Sometimes marker mass inside the region is not measured directly. Instead, entry and exit marker rates are observed. Let $q_{\mathrm{in}}(t)$ and $q_{\mathrm{out}}(t)$ denote marker mass flow rates across the entry and exit boundaries, in mg/h. If the initial inventory is zero and there is no marker generation, degradation, absorption, or unmeasured loss within the region, the inventory is the cumulative entry minus cumulative exit:

M(t)=\int_0^t\left[q_{\mathrm{in}}(s)-q_{\mathrm{out}}(s)\right]ds.

For a complete passage, the total amounts entering and exiting are equal, and the inventory is zero at the start and after the passage. The integral of inventory can then be written as the difference between the first time moments of the exit and entry rates. Dividing by the total marker amount gives the difference between mean exit time and mean entry time. This difference is the mean time spent in the region.

Result: If total entry and exit marker amounts are both $D$, and there is no initial or final marker inventory, then $$ \bar{\tau}=\frac{\int_0^T M(t)\,dt}{D} =\frac{\int_0^T t\,q_{\mathrm{out}}(t)\,dt}{D} -\frac{\int_0^T t\,q_{\mathrm{in}}(t)\,dt}{D}. $$ The time moments must refer to the same time origin and cover the complete entry and exit profiles.

A brief derivation makes the boundary conditions explicit. Write $M(t)$ as cumulative input minus cumulative output, then integrate over $0$ to $T$. Reversing the order of integration gives the time-weighted difference below. If total input and output equal $D$, the terms involving $T$ cancel.

\begin{aligned} \int_0^T M(t)\,dt &=\int_0^T (T-s)\left[q_{\mathrm{in}}(s)-q_{\mathrm{out}}(s)\right]ds\\ &=T(D_{\mathrm{in}}-D_{\mathrm{out}}) -\int_0^T s\,q_{\mathrm{in}}(s)\,ds +\int_0^T s\,q_{\mathrm{out}}(s)\,ds\\ &=\int_0^T s\,q_{\mathrm{out}}(s)\,ds -\int_0^T s\,q_{\mathrm{in}}(s)\,ds, \qquad D_{\mathrm{in}}=D_{\mathrm{out}}=D. \end{aligned}

Thus, a marker introduced over a spread of times must not be treated as though all marker entered at the start. Subtracting the entry mean time corrects for the input profile’s duration. For an effectively instantaneous input at time zero, its mean entry time is zero, and the mean exit time alone gives the regional mean residence time.

Worked Example: Recovering Residence Time From Exit Measurements

In a hypothetical experiment, an 80 mg marker pulse enters at time zero. Exit measurements indicate that 10 mg exits between 0 and 1 h, 25 mg between 1 and 2 h, 30 mg between 2 and 3 h, and 15 mg between 3 and 4 h. Assume the exit mass is approximately uniform within each interval, so each interval’s representative exit time is its midpoint.

The recovered amount is $10+25+30+15=80$ mg, matching the entry amount. The estimated exit time moment is $10(0.5)+25(1.5)+30(2.5)+15(3.5)=5+37.5+75+52.5=170$ mg·h. The input time moment is $80(0)=0$ mg·h. Hence,

\bar{\tau}\approx \frac{170-0}{80}=2.125\ \text{h}.

The estimated mean regional residence time is 2.125 h. This calculation relies on the midpoint approximation within each interval; it is not an exact moment calculation if the marker exits unevenly within those intervals.

Incomplete Recovery and Sparse Sampling

Incomplete recovery is not fixed by simply dividing by the amount observed at the exit. That calculation estimates the mean exit time among the marker that was recovered. It estimates the mean for all marker entering only if the unrecovered fraction can be justified as having the same timing distribution as the recovered fraction, or if its passage is otherwise accounted for. If marker is missing preferentially at particular times, the recovered-marker mean can be biased in either direction.

Worked Example: Interpreting Incomplete Recovery

Suppose 100 mg enters at time zero, but only 80 mg is detected exiting. The recovered amounts are 20 mg at a representative time of 1 h, 30 mg at 2 h, and 30 mg at 3 h. The recovered-marker mean exit time is

\bar{\tau}_{\mathrm{recovered}} =\frac{20(1)+30(2)+30(3)}{80} =\frac{170}{80} =2.125\ \text{h}.

This is a 2.125 h mean for the recovered 80 mg. The missing 20 mg has an unknown time contribution, so the result alone does not determine the mean for all 100 mg. For illustration, if that missing 20 mg were known to exit at 6 h, the complete mean would be $[170+20(6)]/100=290/100=2.90$ h. If the missing marker never exits because it is lost or transformed, a complete-marker exit residence time cannot be recovered from these observations without additional assumptions or measurements.

Sparse sampling can also change the area estimate, especially when the inventory curve bends sharply or has an unobserved peak or tail. It does not necessarily change an estimate in a predictable direction. Compare the following sparse observations with those in the first example: inventory is 100 mg at 0 h, 50 mg at 2 h, and 0 mg at 4 h. The sparse trapezoidal areas are $(100+50)(2)/2=150$ mg·h and $(50+0)(2)/2=50$ mg·h. Thus $A_{\mathrm{sparse}}=150+50=200$ mg·h and $\bar{\tau}_{\mathrm{sparse}}\approx200/100=2.00$ h. Here the sparse estimate agrees with the more frequent-sampling estimate. That agreement follows from these particular values; other curves can produce different estimates.

Worked Example: Checking a Sparse-Sampling Estimate

For the 100 mg input above, calculate the sparse estimate directly rather than assuming it must differ from the denser estimate.

From 0 to 2 h, the trapezoidal area is $(100+50)(2)/2=150$ mg·h. From 2 to 4 h, it is $(50+0)(2)/2=50$ mg·h. Therefore $A_{\mathrm{sparse}}=150+50=200$ mg·h, and

\bar{\tau}_{\mathrm{sparse}}\approx\frac{200\ \text{mg}\cdot\text{h}}{100\ \text{mg}}=2.00\ \text{h}.

The frequent and sparse trapezoidal estimates happen to agree at 2.00 h for these measurements. Agreement in this example is not evidence that sparse sampling is generally harmless: intermediate observations can reveal curvature that the two-interval calculation cannot detect.

Assumptions, Diagnostics, and Alternatives

The estimates above require a clearly defined region, a marker that moves with the material of interest, and measurements that correspond to marker mass rather than an unconverted concentration. For inventory integration, the initial and final inventories should be zero or explicitly accounted for. For the entry–exit moment method, entry and exit profiles must be measured over the full passage, with equal total amounts under the conservative-marker assumption. If the marker is absorbed, degraded, retained beyond the observation period, or otherwise lost, the mass-balance derivation does not apply without additional terms or assumptions.

Useful checks include comparing total entry and exit recovery, inspecting whether the calculated inventory stays nonnegative, and checking whether measurements show a plausible return to baseline. A negative inventory inferred from cumulative entry minus exit suggests measurement error, timing mismatch, or inconsistent flow estimates. A truncated exit tail can shift the exit time moment downward if late marker is omitted. Concentration sampled at a boundary should be converted to mass flow using the corresponding flow rate when that is the required measured quantity; concentration alone does not specify how much marker crossed per unit time.

The inventory-area method is direct when regional marker mass is measured repeatedly. The entry–exit moment difference is an alternative when boundary flux profiles are available, and it can correct for a non-instantaneous input. An exit-only mean is simpler but applies directly only when the input timing is effectively a pulse at the chosen time origin, or when the entry-time mean is subtracted separately. These are alternative ways to summarize the same complete passage, not interchangeable calculations when their assumptions differ.

A residence-time estimate is a summary of the measured marker passage, not proof that every drug molecule follows that passage. Marker mismatch, incomplete recovery, uncertain regional boundaries, and numerical integration all limit how directly the estimate can be used as a model input. Record the measurement basis, boundaries, recovery, sampling schedule, and any correction or approximation alongside a derived value.

Common Mistakes

Dividing by recovered marker and calling the result the full-population mean. This gives the mean for recovered marker unless the missing fraction’s timing is addressed.
Using concentration as though it were mass flow. A concentration profile alone does not tell how much marker crosses a boundary without the relevant flow information.
Ignoring the entry-time profile. With a spread-out input, the exit mean includes time spent before entry; subtract the mean entry time to estimate time between boundaries.
Assuming sparse sampling always underestimates residence time. The effect depends on the unobserved curve shape. Sparse and frequent estimates can agree, as the worked calculation demonstrates.
Stopping observation before the tail has passed. Unobserved late marker can make the measured exit profile incomplete and distort its mean.

Check Your Understanding

Try each question, then check the worked solutions below.

A 60 mg marker input produces an inventory curve with trapezoidal area 150 mg·h. What is the estimated mean residence time?
Why is boundary concentration alone insufficient to calculate the marker amount crossing per unit time?
A marker input mean time is 0.5 h and its complete exit mean time is 3.0 h. Under the mass-balance assumptions, what is the regional mean residence time?
For 100 mg entering, sparse inventory observations are 100 mg at 0 h, 50 mg at 2 h, and 0 mg at 4 h. Calculate the trapezoidal area and residence time.

Solutions. (1) $\bar{\tau}\approx150/60=2.50$ h. (2) Concentration is marker amount per volume; the amount crossing per unit time also depends on the volume flow across the boundary. (3) $3.0-0.5=2.50$ h, the difference between the mean exit and entry times. (4) The areas are $(100+50)(2)/2=150$ mg·h and $(50+0)(2)/2=50$ mg·h, totaling 200 mg·h; $200/100=2.00$ h.