Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Oral Dose Entry and Gastric-Emptying Boundary Conditions

Compartmental Absorption and Transit Models · Tutorial 215 of 799

Oral Dose Entry and Gastric-Emptying Boundary Conditions

Learn to choose between an initial dose condition and a gastric-emptying input, then verify that the modeled dose is conserved.

Intermediate 10 min read

What You'll Learn

  • Distinguish an initial amount in the first intestinal segment from an input delivered over time.
  • Write a gastric-emptying input from an explicit stomach amount or a specified emptying profile.
  • Check that an input-rate function has the correct units and delivers the intended dose.
  • Account for direct luminal loading and gastric delivery without counting any dose twice.
  • Diagnose mass-balance errors caused by inconsistent initial conditions or boundary inputs.

Where an Oral Dose Enters a Segmental Model

A segmental intestinal model needs a precise description of how much drug is present at the start and how much arrives afterward. These are distinct modeling choices. A dose already present in the first modeled intestinal segment belongs in that segment’s initial condition. Drug that arrives later from the stomach belongs in an input-rate term at the segment’s inlet. Treating the same dose as both an initial amount and a continuing input counts it twice; treating it as neither omits it.

Here, amounts refer to the modeled drug mass, for example milligrams of active drug, on one consistent mass basis. Let \(A_i(t)\) denote the luminal drug amount in intestinal segment \(i\), and let \(u_1(t)\) denote the rate at which drug enters the first modeled segment. Thus \(A_i\) has units of amount and \(u_1\) has units of amount per time. Any losses or transfers within the intestinal model must be represented separately in its segment balances.

Definition: An initial luminal dose is represented by a nonzero initial amount, such as \(A_1(0)=D\), with no additional input for that same dose. A time-dependent dose input is represented by an inlet rate \(u_1(t)\), whose time integral gives the amount delivered. A dose may be partitioned among these representations only when the partition is explicit and the pieces sum to the administered modeled dose.

The starting condition and input are independent pieces of the model specification. A general first-segment balance can be written as

$$ \frac{dA_1}{dt} = u_1(t) + \text{other inflows} - \text{transfers out} - \text{other specified losses}. $$

The exact transfer and loss terms depend on the segmental mass balances. The important boundary-condition point is that \(u_1(t)\) represents drug crossing into the modeled intestinal system; it is not itself the amount currently in the segment. If the oral dose is directly loaded into segment 1 at time zero, set \(A_1(0)=D\) and set the dose-related part of \(u_1(t)\) to zero. If the dose is in a modeled stomach compartment, set its initial amount there and let gastric emptying generate \(u_1(t)\).

Gastric Emptying as an Inlet Boundary

There are two common bookkeeping choices. A model may represent the stomach explicitly, with a gastric amount \(G(t)\), or may prescribe the amount entering the intestine as a time-dependent boundary input. These choices can describe the same emptying process, but they should not both be applied to the same material in a way that duplicates delivery.

For an explicit stomach with first-order emptying, assume the dose is initially in the stomach, emptying is the only gastric loss, and the emptying rate constant \(k_g\) is constant and nonnegative. Then the gastric amount and intestinal input satisfy

$$ \frac{dG}{dt}=-k_gG,\qquad u_1(t)=k_gG(t),\qquad G(0)=D,\qquad A_1(0)=0. $$

Here \(k_g\) has units of inverse time. The solution is \(G(t)=D e^{-k_gt}\), so \(u_1(t)=k_gD e^{-k_gt}\). This is a particular emptying assumption, not a requirement for every oral-dose model. If emptying is delayed or follows a different time course, the input function must reflect that alternative assumption.

A prescribed gastric-emptying profile can instead be expressed through \(E(t)\), the cumulative fraction of the gastric dose emptied by time \(t\). For a dose \(D\) initially in the stomach, assume \(E(0)=0\), \(E(t)\) is nondecreasing, and its eventual value is 1 if the entire dose empties. If \(E\) is absolutely continuous on \([0,T]\), its ordinary derivative defines the inlet rate almost everywhere, and

$$ u_1(t)=D\frac{dE(t)}{dt}, \qquad \int_0^T u_1(t)\,dt=D E(T). $$

This form makes the amount delivered by any time \(T\) explicit. A profile that has a lag may have \(E(t)=0\) before the lag; a profile that leaves some material permanently in the stomach has a limiting value below 1, provided that the retained amount is accounted for in the model. If the profile has jumps or other non-absolutely-continuous increases, represent those portions as amount transfers or use an input measure that accounts for the full increase in \(E\), rather than silently approximating them by an ordinary finite rate.

Why the Boundary Condition Conserves Dose

For an explicit stomach and first intestinal segment, suppose the only transfer between them is gastric emptying. Write the first-segment balance schematically as \(dA_1/dt=k_gG-L_1(t)\), where \(L_1(t)\) collects the segment’s net losses, such as onward transit and absorption. Adding this equation to the gastric balance gives

$$ \frac{d}{dt}(G+A_1) = -k_gG+k_gG-L_1(t) = -L_1(t). $$

The gastric outflow disappears from the sum because it is an internal transfer: it leaves the stomach and enters segment 1 in equal amounts. With no intestinal losses, \(G(t)+A_1(t)\) remains equal to its initial value \(D\). With losses, the total amount in the stomach and intestine falls only by those losses. This cancellation is the basic mass-conservation check for an emptying boundary.

For a prescribed inlet rate and no explicit stomach state, integration of the first-segment balance shows that the dose delivered through the boundary by time \(T\) is \(\int_0^T u_1(t)\,dt\). The model must then account for any administered dose not yet delivered, for example through an explicitly tracked gastric remainder or through the stated cumulative emptying profile. Prescribing an input does not make undelivered dose disappear from the dose accounting.

1
Choose the dose location at time zero.
Put an amount in the stomach, in the first intestinal segment, or in explicitly identified fractions at different locations.
2
Specify the emptying assumption.
Use an explicit gastric state and its outflow, or prescribe a cumulative emptying profile and derive the inlet rate.
3
Check units and total delivery.
Confirm that the input has units of amount per time and that its integral matches the dose fraction assigned to that input.
4
Audit the dose ledger.
At each time, distinguish material already delivered from material still in the stomach; do not include either fraction twice.

Worked Example: A Dose Loaded Directly Into the First Segment

Suppose a hypothetical model begins with a 120 mg oral dose already dissolved in the first modeled intestinal segment. No gastric state is included, and there is no later dose-related input. Specify the dose conditions and check the initial dose accounting.

Because the entire modeled dose is present in segment 1 at time zero, set \(A_1(0)=120\) mg and \(u_1(t)=0\) mg/h for the dose input. The dose-related initial amount in all other modeled compartments is zero under this setup. The amount initially represented in the gut model is therefore \(120\) mg.

This condition does not assert that the dose was physiologically absorbed immediately. It states only that the model begins after the dose has reached the first represented segment. Subsequent transit and absorption are determined by the segmental balance terms. Adding a gastric-emptying input totaling another 120 mg would incorrectly make the model represent 240 mg.

Worked Example: First-Order Emptying With a Directly Loaded Fraction

Consider a hypothetical 200 mg dose partitioned by assumption into 60 mg loaded directly into segment 1 at time zero and 140 mg initially in a stomach compartment. Let gastric emptying be first-order with \(k_g=0.50\) h\(^{-1}\). Find the amount delivered from the stomach by 1.0 h, the total dose represented in the intestine by then before intestinal losses, and the gastric remainder.

The initial conditions are \(A_1(0)=60\) mg and \(G(0)=140\) mg. For the gastric fraction, \(G(t)=140e^{-0.50t}\) mg and \(u_1(t)=0.50(140)e^{-0.50t}\) mg/h. Integrating the input from 0 to 1.0 h gives

$$ \int_0^{1.0} u_1(t)\,dt = 140(1-e^{-0.50}) = 140(1-0.6065306597\ldots) = 55.0857076\ldots\ \text{mg}. $$

Thus the total dose delivered into the intestine by 1.0 h, before accounting for intestinal losses, is \(60+55.0857076\ldots=115.0857076\ldots\) mg, or 115.1 mg to four significant figures. The gastric remainder is \(140e^{-0.50}=84.9142924\ldots\) mg, or 84.91 mg. The unrounded dose ledger closes:

$$ \bigl[60+140(1-e^{-0.50})\bigr]+140e^{-0.50} =60+140=200\ \text{mg}. $$

The displayed rounded approximations, 115.09 mg delivered and 84.91 mg remaining, sum to 200.00 mg. The unrounded expressions sum exactly to 200 mg. This check concerns the dose entering the intestine and remaining in the stomach; any intestinal absorption or onward transit would change the amount currently in segment 1.

Worked Example: A Delayed, Prescribed Gastric Input

A hypothetical 100 mg dose is assigned to the stomach, with no emptying during the first 0.5 h. After that lag, assume first-order emptying with \(k_g=0.80\) h\(^{-1}\). Find the gastric input rate for \(t\geq0.5\) h, and calculate the amount delivered and the amount remaining in the stomach at 1.5 h.

The gastric amount remains 100 mg during the lag. For \(t\geq0.5\) h, it follows \(G(t)=100e^{-0.80(t-0.5)}\) mg. The input is zero before the lag and is \(u_1(t)=0.80(100)e^{-0.80(t-0.5)}\) mg/h afterward. At 1.5 h, one hour of post-lag emptying has elapsed:

$$ \text{delivered} = \int_{0.5}^{1.5}u_1(t)\,dt = 100(1-e^{-0.80}) = 55.0671\ldots\ \text{mg}. $$

The gastric remainder is \(100e^{-0.80}=44.9329\ldots\) mg, so delivered plus remaining is \(55.0671\ldots+44.9329\ldots=100.0000\ldots\) mg. The lag changes the timing, not the dose total. In a numerical implementation, the input must be zero through the lag and must use time since the end of the lag in the exponential; using time since zero would incorrectly empty the dose during the lag.

Diagnostics, Alternatives, and Limitations

A useful diagnostic is a dose ledger at selected times: initial amount, cumulative inlet delivery, amount remaining in the stomach, and any explicitly modeled losses or transfers. If no loss has yet occurred, the dose fractions in those locations should sum to the assigned dose. With losses, adding cumulative loss to the amounts still represented should recover the assigned dose, subject only to numerical integration error.

RepresentationWhat is specifiedKey check
Initial luminal amountA dose amount in the initial condition, such as \(A_1(0)=D\)Do not add another input for that same dose
Explicit stomach compartmentInitial gastric amount and an emptying outflow into segment 1Stomach outflow must equal the intestinal inlet transfer
Prescribed boundary inputAn inlet rate or cumulative emptying profileThe rate integral must equal the assigned delivered fraction

Common warning signs include an intestinal amount exceeding the administered dose before any external source or mass-generating process is specified, a nonzero input with units of amount rather than amount per time, or a cumulative input that never reaches its intended total. Another frequent error is to set \(A_1(0)=D\) and also use \(u_1(t)=k_gD e^{-k_gt}\). The latter input represents an additional dose; it is not a second way of expressing the same initial condition.

An explicit stomach compartment is useful when the remaining gastric amount matters or when its emptying law must interact with other modeled processes. A prescribed input is simpler when the gastric process is represented only through its delivery profile. Direct initial loading is appropriate when the model deliberately begins with the dose already in the first intestinal segment. These are alternative representations, not claims that the physiological process is identical in every setting. A fixed input cannot respond dynamically to changes in gastric amount unless that dependence is built into the input specification.

The first-order example assumes a constant emptying rate and an exponentially declining gastric remainder. It should not be used merely because it is mathematically convenient if the intended gastric delivery profile has a lag or another shape. Likewise, a prescribed profile is only as informative as its assumptions and supporting data. The introductory material on gastric emptying and intestinal transit-time distributions provides physiological context for choosing a profile; here the modeling requirement is to translate the chosen profile into a consistent boundary condition.

Key takeaway: Put each dose fraction in exactly one place at the start of the model or deliver it through exactly one inlet input. For gastric emptying, the amount leaving the stomach must enter the first intestinal segment in equal measure, and the cumulative input must agree with the assigned dose.

Check Your Understanding

For each item, identify the dose representation and check its mass accounting. Solutions are included below.

  1. A model starts with \(A_1(0)=75\) mg and has no dose-related input. How much dose is initially represented, and should gastric input be added for the same 75 mg?
  2. A 90 mg gastric dose empties according to \(E(t)=1-e^{-0.4t}\), with time in hours. What input rate follows, and how much has entered the intestine by 2 h?
  3. A model has 30 mg initially in segment 1 and an additional prescribed input whose integral over all time is 50 mg. What total dose is assigned across these two sources?
  4. For an explicit stomach with \(G(0)=80\) mg and no gastric loss other than emptying, what must be true of the gastric outflow and segment-1 inlet for the stomach-to-intestine transfer to conserve mass?

Solutions. 1. The initial dose represented is 75 mg; do not add gastric input for that same dose. 2. The rate is \(u_1(t)=90(0.4)e^{-0.4t}=36e^{-0.4t}\) mg/h. By 2 h, the delivered amount is \(90(1-e^{-0.8})\approx49.5604\) mg (49.6 mg to 3 significant figures), with \(90e^{-0.8}\approx40.4396\) mg remaining (40.4 mg to 3 significant figures); the unrounded amounts sum to 90 mg. 3. The assigned dose is \(30+50=80\) mg, assuming these are distinct fractions and no other dose source is present. 4. The outflow from the stomach must equal the inlet to segment 1 at each time; then the two transfer terms cancel when their mass balances are added.