Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Limiting Cases for Fast Transit and Fast Absorption

Compartmental Absorption and Transit Models · Tutorial 230 of 799

Limiting Cases for Fast Transit and Fast Absorption

Learn how the ratio of absorption to transit rates determines where a dose is absorbed, and how limiting-case calculations provide checks on segmental simulations.

Intermediate 11 min read

What You'll Learn

  • Derive the absorbed and transited fractions from competing first-order rates.
  • Predict dose allocation when transit is much faster than absorption.
  • Predict dose allocation when absorption is much faster than transit.
  • Extend single-segment limits to a sequential chain of gut segments.
  • Check simulation outputs against limiting-case calculations and mass balance.
  • Recognize when a rate comparison alone is insufficient to establish a limit.

Why Limiting Cases Are Useful

A segmental absorption model assigns drug two possible fates while it is in a gut segment: it may be absorbed, or it may transit to the next segment. The relative rates of those competing processes determine the fraction taking each route. Limiting cases make that competition easier to interpret. If transit is much faster than absorption, little drug is absorbed in each segment. If absorption is much faster than transit, most drug is absorbed before it can move downstream.

These predictions are useful both for understanding model behavior and for checking numerical simulations. A simulation that predicts substantial absorption in a segment where transit is overwhelmingly faster, for example, deserves careful investigation. The comparisons here assume dissolved drug is available to undergo absorption; they do not model dissolution as a separate rate-limiting process.

Definition: In a segment with constant first-order absorption rate \(k_a\) and first-order transit rate \(k_{tr}\), the fast-transit limit is the regime \(k_{tr} \gg k_a\), and the fast-absorption limit is the regime \(k_a \gg k_{tr}\). The symbols \(k_a\) and \(k_{tr}\) have units of inverse time. These comparisons concern the rates’ relative sizes; their sum also determines how quickly the competing events occur.

One Segment: Deriving the Competing Fractions

Let \(A(t)\) be the dissolved drug amount in one segment at time \(t\), with initial amount \(A_0\). Assume the segment is well mixed for this model, both exit processes are first order with constant rates, and there are no other gains or losses from the segment. Absorption removes drug from the lumen, while transit transfers it to the next segment. The amount remaining therefore obeys

$$ \frac{dA}{dt}=-(k_a+k_{tr})A,\qquad A(0)=A_0. $$

Solving this first-order equation gives \(A(t)=A_0e^{-(k_a+k_{tr})t}\). The absorption and transit fluxes are \(k_aA(t)\) and \(k_{tr}A(t)\), respectively. Integrating each flux from entry until the amount has left the segment gives the corresponding amounts:

$$ A_{\mathrm{abs}}=\int_0^\infty k_a A_0e^{-(k_a+k_{tr})t}\,dt =A_0\frac{k_a}{k_a+k_{tr}}, \qquad A_{\mathrm{tr}}=\int_0^\infty k_{tr} A_0e^{-(k_a+k_{tr})t}\,dt =A_0\frac{k_{tr}}{k_a+k_{tr}}. $$

The two fractions sum to one. This is a direct consequence of the assumptions that these are the only two exits and that both rates remain constant while drug is in the segment. The characteristic time for disappearance from the segment is \(1/(k_a+k_{tr})\); the allocation between absorption and transit depends on the ratio of the rates, not on that time alone.

Key result: The fraction absorbed in one segment is \(k_a/(k_a+k_{tr})\), and the fraction transited is \(k_{tr}/(k_a+k_{tr})\). In the fast-transit limit, the absorbed fraction is approximately \(k_a/k_{tr}\). In the fast-absorption limit, the transited fraction is approximately \(k_{tr}/k_a\).

Fast Transit Across a Sequential Chain

For a chain of segments, use \(k_{a,j}\) and \(k_{tr,j}\) for the absorption and transit rates in segment \(j\). If the dose entering segment 1 is \(D\), and there is no bypass or return flow, the fraction of that dose reaching segment \(j\) is the product of the transit fractions in all preceding segments. Thus the amount absorbed in segment \(j\) is

$$ A_{\mathrm{abs},j} = D\left[\prod_{i=1}^{j-1} \frac{k_{tr,i}}{k_{a,i}+k_{tr,i}}\right] \frac{k_{a,j}}{k_{a,j}+k_{tr,j}}. $$

An empty product for \(j=1\) is defined as 1. The amount remaining after the final segment, if transit out of that segment is included, is \(D\prod_{i=1}^{N} k_{tr,i}/(k_{a,i}+k_{tr,i})\), where \(N\) is the number of segments. The absorbed amounts plus that final transited amount sum to \(D\).

If every segment has \(k_{a,j}/k_{tr,j}=\epsilon\), the total fraction absorbed over \(N\) segments is

$$ F_{\mathrm{abs}}=1-\left(\frac{1}{1+\epsilon}\right)^N. $$

For fixed \(N\) and small \(\epsilon\), this is approximately \(N\epsilon\) when \(N\epsilon\) is also small. In particular, with a fixed number of segments, the total absorbed fraction tends to zero as \(k_a/k_{tr}\) tends to zero. A very large number of segments can change the practical outcome: small absorption opportunities can accumulate across the chain. The single-segment limit should therefore not be mistaken for a claim that fast transit always means zero total absorption in any possible model.

Worked Example: Small Absorption in Three Fast-Transit Segments

Consider a hypothetical 100 mg dissolved dose entering three segments. In each segment, let \(k_a=0.020\ \mathrm{h}^{-1}\) and \(k_{tr}=1.50\ \mathrm{h}^{-1}\). Transit is 75 times faster than absorption, so this is a fast-transit regime. Assume constant rates and no other loss or gain.

In one segment, the absorbed fraction is \(0.020/(1.50+0.020)=1/76\), and the transit fraction is \(1.50/1.52=75/76\). Therefore the first segment absorbs \(100/76=1.31579\) mg. The second segment receives \(100(75/76)\) mg and absorbs \(100(75/76)(1/76)=1.29848\) mg. The third absorbs \(100(75/76)^2(1/76)=1.28139\) mg.

The total absorbed amount is computed from the unrounded fractions: \(100[1-(75/76)^3]=3.89566\) mg, or 3.896% of the entering dose. The amount transited beyond segment 3 is \(100(75/76)^3=96.1043\) mg. The two amounts sum to 100 mg. The small but nonzero total absorption illustrates why the fast-transit approximation means “small per segment,” not “exactly none.”

Fast Absorption and the Location of Uptake

When \(k_a\gg k_{tr}\), the single-segment absorbed fraction \(k_a/(k_a+k_{tr})\) approaches one. For a finite rate ratio, it is not exactly one: some drug still transits. In a sequential chain, the small fraction that escapes absorption in an early segment can reach later segments, but the amount reaching each successive segment generally decreases sharply when absorption dominates in each segment.

For equal rates across \(N\) segments, write \(\delta=k_{tr}/k_a\). The fraction that transits one segment is \(\delta/(1+\delta)\), so the fraction remaining after transit through all \(N\) segments is \([\delta/(1+\delta)]^N\). The total fraction absorbed is one minus that residual fraction. As \(\delta\) tends to zero, nearly all the dose is absorbed, and nearly all uptake occurs in the first segment. This location prediction depends on the sequential-entry assumption: a dose delivered directly to a later segment does not first pass through earlier segments.

Worked Example: Uptake Is Nearly Complete Before the Third Segment

For a hypothetical 100 mg dose entering three sequential segments, take \(k_a=9.0\ \mathrm{h}^{-1}\) and \(k_{tr}=1.0\ \mathrm{h}^{-1}\) in every segment. The rates are constant, absorption is the only competing loss, and all material that transits enters the next segment.

Each segment absorbs \(9/(9+1)=0.9\) of the amount entering it and transits \(0.1\). Segment 1 absorbs \(100(0.9)=90.0\) mg, leaving 10.0 mg to enter segment 2. Segment 2 absorbs \(10.0(0.9)=9.00\) mg, leaving 1.00 mg to enter segment 3. Segment 3 absorbs \(1.00(0.9)=0.900\) mg, leaving \(0.100\) mg to transit beyond the chain.

The total absorbed amount is \(90.0+9.00+0.900=99.9\) mg, and the final transited amount is 0.100 mg. The mass balance closes at 100 mg. Although this is a fast-absorption regime, the 0.100 mg that exits the third segment shows that “nearly complete” should not be reported as exact completion.

Mixed Regimes and an Analytical Simulation Check

A chain need not have the same rate ordering everywhere. Segment geometry, effective permeability, and other modeled conditions can yield fast transit in one segment and fast absorption in another. The general product formula allows those segments to be handled in sequence without replacing each finite rate by an exact zero or infinity.

Worked Example: Fast Transit Followed by Fast Absorption

Suppose a hypothetical 50 mg dissolved dose enters segment 1, where \(k_{a,1}=0.10\ \mathrm{h}^{-1}\) and \(k_{tr,1}=1.00\ \mathrm{h}^{-1}\). It then enters segment 2, where \(k_{a,2}=4.00\ \mathrm{h}^{-1}\) and \(k_{tr,2}=0.20\ \mathrm{h}^{-1}\). Assume no other exits, constant rates, and complete transfer of the segment-1 transit amount into segment 2.

In segment 1, the absorbed fraction is \(0.10/1.10=1/11\). The absorbed amount is \(50/11=4.54545\) mg, and \(50(10/11)=45.4545\) mg enters segment 2.

In segment 2, the absorbed fraction is \(4.00/4.20=20/21\). The absorbed amount there is \((500/11)(20/21)=10000/231=43.2900\) mg. The amount transited beyond segment 2 is \((500/11)(1/21)=500/231=2.16450\) mg. Total absorption is \(4.54545+43.2900=47.8355\) mg, using the unrounded amounts, and the remaining 2.16450 mg brings the total to 50.0000 mg. The first segment has a fast-transit ordering, whereas the second has a fast-absorption ordering; neither finite rate ratio is treated as an exact limiting value.

Analytical predictions are also useful when a numerical solver integrates the amount over time. For example, with \(A_0=10.0\) mg, \(k_a=100\ \mathrm{h}^{-1}\), and \(k_{tr}=1.00\ \mathrm{h}^{-1}\), the eventual absorbed amount is \(10(100/101)=9.90099\) mg and the transited amount is \(10(1/101)=0.0990099\) mg. A simulation intended to run until the segment is effectively empty should approach those values. If it does not, check that the rates, units, initial amount, and flux definitions match the equations. Also check whether the simulation has run long enough: the remaining amount at time \(t\) is \(10e^{-101t}\) mg, so a short run may not yet have allocated all the dose.

Assumptions, Diagnostics, and Common Pitfalls

The calculations above are exact for constant, competing first-order rates in each segment, with absorption and transit as the only exits from the available luminal amount. The sequential-chain formula additionally assumes that all material transited from one segment enters the next, without a delay, bypass, or return path. If a model includes time-varying rates, additional losses, incomplete transfer, or separate solid and dissolved states, these formulas may not describe its full system. They can still provide a check for a simplified component, but the assumptions must be matched before comparing results.

A limiting comparison is most informative when stated as a ratio. For instance, \(k_{tr}/k_a=10\) gives an absorbed fraction of \(1/11\), not zero. Whether that is sufficiently small for a particular purpose depends on the needed accuracy and on how many segments the dose traverses. Likewise, \(k_a/k_{tr}=10\) gives a per-segment absorption fraction of \(10/11\), not exactly one. Limits are approximations to finite-rate behavior.

The sum \(k_a+k_{tr}\) controls the disappearance timescale. If both rates are multiplied by the same factor, their fractions are unchanged, even though both events occur faster. Thus a claim about fast absorption relative to transit requires comparing \(k_a\) with \(k_{tr}\), not merely observing that \(k_a\) is numerically large in isolation. Ensure rates are expressed in compatible time units before taking their ratio.

ApproachWhat it providesImportant limitation
Limiting-case calculationA quick prediction of dose allocation when one rate greatly exceeds the otherApproximations may be poor when the rates are comparable or when many small opportunities accumulate
Finite-rate analytical calculationExact competing fractions and sequential allocations under the stated first-order assumptionsDoes not represent omitted processes or time-varying rates
Numerical simulationTime courses and more complex model structuresCan be affected by implementation errors, unit mismatches, insufficient run time, or numerical error

A common mistake is to set the slower process to zero as soon as the other rate is called “fast.” For example, reasoning that \(k_{tr}\gg k_a\) means no absorption ignores the finite fraction \(k_a/(k_a+k_{tr})\), which can accumulate over multiple segments. The correction is to calculate the exact fraction for the specified rates, then decide whether its limiting approximation is adequate. Similarly, fast absorption does not mean that all drug is absorbed immediately at segment entry: absorption is a first-order process over time, and the competing transit fraction remains nonzero unless the rate ratio is taken to its limiting value.

Key diagnostic: Under the stated assumptions, the amounts absorbed across all segments plus the amount transited beyond the last segment must equal the entering dose. A failure of this balance, or a segment allocation inconsistent with \(k_a/(k_a+k_{tr})\), signals that the calculation, implementation, or assumed model structure needs review.

Check Your Understanding

Answer each question using the competing-rate assumptions described above. Brief solutions follow the questions.

  1. In one segment, \(k_a=0.050\ \mathrm{h}^{-1}\) and \(k_{tr}=0.95\ \mathrm{h}^{-1}\). What fraction is absorbed there?
  2. For three identical segments with \(k_a/k_{tr}=0.01\), write the exact total absorbed fraction and explain why it is not exactly 0.03.
  3. In one segment, both rates are multiplied by four. What happens to the absorbed fraction and to the disappearance timescale?
  4. A simulation reports 100% absorption in a segment with finite positive \(k_{tr}\). Is this the exact competing-rate prediction?

Solutions. (1) The absorbed fraction is \(0.050/(0.050+0.95)=0.0500\), or 5.00%. (2) The exact fraction is \(1-[1/(1.01)]^3\), approximately 0.0294. It is close to 0.03, but lower because each segment’s exact transit fraction is \(1/1.01\), and later segments receive only what escaped earlier. (3) The absorbed fraction is unchanged because the ratio of the rates is unchanged; the disappearance rate \(k_a+k_{tr}\) is four times larger, so the characteristic time \(1/(k_a+k_{tr})\) is one-quarter as long. (4) No. With finite positive transit and absorption rates, \(k_a/(k_a+k_{tr})<1\). The result may reflect rounding, a different model structure, or an implementation or reporting problem.