Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Meal Timing and Gastric State Transitions

Food Effects and Gastric Emptying · Tutorial 287 of 799

Meal Timing and Gastric State Transitions

Represent meal timing with gastric state intervals and calculate how state transitions change the amount and rate of a dose leaving the stomach.

Advanced 10 min read

What You'll Learn

  • Define fed and fasted gastric-state intervals relative to meal and dose times.
  • Derive dose remaining and emptying rate under piecewise first-order emptying.
  • Calculate profiles when a meal occurs before or after dosing.
  • Handle multiple fed-state intervals without treating the whole period as continuously fed.
  • Check state assignments, mass balance, and assumptions behind abrupt transitions.

Why meal timing must be represented over time

A dose taken with food and a dose taken before a meal do not necessarily experience the same gastric conditions throughout their residence in the stomach. A meal eaten after dosing can change the emptying rate of drug still present in the stomach. Likewise, a dose taken after a meal may begin emptying in a fed state, then continue emptying after the stomach returns to a fasted state. The relevant question is therefore not only whether food was present at dosing, but which gastric state the dose experiences over time.

The earlier material on fed-state physiology and emptying models provides possible state-specific emptying parameters. Here, meal timing determines when those parameters apply. The examples use a first-order emptying hazard to make the transitions explicit; the same timeline logic can be used with other emptying models, provided their state dependence is defined appropriately.

Definition: Let \(q(t)\) denote the gastric state at time \(t\), with \(q(t)\) taking the value fasted or fed. A state transition is a specified time at which \(q(t)\) changes. In a piecewise model, meal timing defines intervals over which the state is treated as constant. For example, a meal at time \(t_m\) might initiate a fed interval at \(t_m+\delta_{\mathrm{on}}\), which ends at \(t_m+\delta_{\mathrm{off}}\). The onset and end times are model assumptions or inputs, not universal constants.

Set dosing time to \(t=0\), and let \(G(t)\) be the amount of the dose remaining in the stomach, in milligrams. In the first-order model, the state-specific emptying hazard is \(k_F\) in the fasted state and \(k_{\mathrm{fed}}\) in the fed state, both in inverse hours. The hazard is the instantaneous emptying rate divided by the amount still in the stomach. The amount leaving the stomach per unit time is \(k_{q(t)}G(t)\).

Model assumptions: The dose begins in the stomach at \(t=0\); the gastric amount is removed only by emptying; the state-specific hazard is constant within each state interval; and a state transition changes the hazard immediately without creating or removing dose mass. This first-order model is memoryless: the hazard depends on the current state, not on how long the dose has already been in the stomach.

Calculating a piecewise emptying profile

Within any interval where the state is constant, the gastric amount satisfies a first-order differential equation. If the state changes at time \(\tau\), calculate the amount at \(\tau\) using the state before the transition, then continue from that same amount using the new state-specific hazard. Thus, the remaining amount is continuous at the transition even though the emptying rate may change abruptly.

$$ \frac{dG(t)}{dt}=-k_{q(t)}G(t), \qquad G(t)=G(0)\exp\left(-\int_0^t k_{q(u)}\,du\right) $$

When the state is constant on successive intervals, the integral is a sum of hazard-times-duration terms. If interval \(j\) has hazard \(k_j\) and duration \(\Delta t_j\), then \(G(t)=G(0)\exp(-\sum_j k_j\Delta t_j)\). The cumulative amount emptied by time \(t\) is \(G(0)-G(t)\), and the instantaneous emptying rate within an interval is \(k_jG(t)\).

A short derivation

For a constant hazard \(k\), separation of variables gives \(dG/G=-k\,dt\). Integrating from the start \(t_a\) of an interval to time \(t\) gives \(\ln G(t)-\ln G(t_a)=-k(t-t_a)\), so \(G(t)=G(t_a)e^{-k(t-t_a)}\). At a transition, the amount at the end of the first interval becomes the initial amount for the next. Repeating this calculation across intervals multiplies the exponential terms, which is equivalent to adding their hazard-times-duration contributions in the exponent. This establishes the piecewise formula without resetting the dose at a meal.

Worked Example: A Meal Eaten After Dosing

Suppose a hypothetical 100 mg dose is taken at \(t=0\). It empties with \(k_F=0.80\ \mathrm{h}^{-1}\) while fasted. A meal at \(t=1\) hour initiates a fed interval lasting until \(t=5\) hours, during which \(k_{\mathrm{fed}}=0.15\ \mathrm{h}^{-1}\). Assume the state returns to fasted at \(t=5\) hours. Calculate the amount remaining at 1, 5, and 6 hours.

From \(t=0\) to \(t=1\), the dose is in the fasted state:

$$ G(1)=100e^{-(0.80)(1)}=44.93\ \mathrm{mg}. $$

From \(t=1\) to \(t=5\), it is in the fed state for 4 hours. The transition does not reset \(G\):

$$ G(5)=44.93e^{-(0.15)(4)} =44.93e^{-0.60} =24.66\ \mathrm{mg}. $$

At 5 hours, 75.34 mg of the original dose has emptied: \(100-24.66=75.34\ \mathrm{mg}\). After the return to fasted state, one further hour gives

$$ G(6)=24.66e^{-(0.80)(1)} =11.08\ \mathrm{mg}. $$

The amount remains continuous at the transitions, but the rate changes. Just before the meal, the rate is \((0.80)(44.93)=35.94\ \mathrm{mg/h}\); just after the fed-state transition, it is \((0.15)(44.93)=6.740\ \mathrm{mg/h}\). At 5 hours, the rate changes from \((0.15)(24.66)=3.699\ \mathrm{mg/h}\) to \((0.80)(24.66)=19.73\ \mathrm{mg/h}\). These are rates of gastric emptying, not concentrations in plasma.

Meals eaten before dosing

For a meal before dosing, the fed interval may already be underway at \(t=0\). In that case, the dose starts with the fed-state hazard, and a later transition to fasted state changes the hazard for whatever dose remains. The meal time itself is not the dose time: both must be placed on the same timeline.

Worked Example: Dosing During a Fed Interval

Consider a hypothetical 100 mg dose given at \(t=0\), one hour after a meal. Assume the fed interval runs from \(t=-1\) to \(t=3\) hours, with \(k_{\mathrm{fed}}=0.15\ \mathrm{h}^{-1}\). After \(t=3\), the state is fasted, with \(k_F=0.80\ \mathrm{h}^{-1}\). Find the amount remaining at 3 and 4 hours after dosing.

The dose is in the fed state for the first 3 hours after dosing:

$$ G(3)=100e^{-(0.15)(3)} =100e^{-0.45} =63.76\ \mathrm{mg}. $$

It then empties for one hour in the fasted state:

$$ G(4)=63.76e^{-(0.80)(1)} =28.65\ \mathrm{mg}. $$

The model predicts a transition in the emptying rate at \(t=3\): just before it, the rate is \((0.15)(63.76)=9.564\ \mathrm{mg/h}\); just after it, the rate is \((0.80)(63.76)=51.01\ \mathrm{mg/h}\). This example assumes that the stated fed interval is an appropriate representation of the meal's effect; it does not imply that every meal produces the same interval or hazards.

Multiple meals and separated fed intervals

If the modeled state returns to fasted between meals, each fed interval must be represented separately. Extending one fed interval from the first meal to the last would incorrectly apply the fed-state hazard during the intervening fasted period. The cumulative-hazard form makes the calculation straightforward: identify each interval, multiply its hazard by its duration, and add those contributions.

Worked Example: Two Meals During Gastric Residence

A hypothetical 100 mg dose is given at \(t=0\). Let \(k_F=0.60\ \mathrm{h}^{-1}\) and \(k_{\mathrm{fed}}=0.20\ \mathrm{h}^{-1}\). Suppose fed intervals occur from \(t=1\) to \(t=2\) and from \(t=4\) to \(t=5\), with fasted state at all other times. Calculate the amount remaining at 5 hours.

The dose experiences one hour fasted, one hour fed, two hours fasted, and one hour fed. Therefore:

$$ G(5)=100\exp[-(0.60)(1)-(0.20)(1)-(0.60)(2)-(0.20)(1)] =100e^{-2.20} =11.08\ \mathrm{mg}. $$

Thus \(100-11.08=88.92\ \mathrm{mg}\) has emptied by 5 hours. If the fed state were incorrectly treated as continuous from \(t=1\) to \(t=5\), the exponent would instead be \(-(0.60)(1)-(0.20)(4)=-1.40\), giving \(24.66\ \mathrm{mg}\) remaining. That different result comes from assigning the wrong state to the interval from 2 to 4 hours.

Choosing and checking the state representation

The piecewise model is useful when state-specific hazards and transition times are available or can be specified. Its simplicity also makes its assumptions visible. A static fed-or-fasted model assigns one state to the entire dose trajectory, often based only on the condition at dosing. A time-varying piecewise model instead updates the hazard when the defined state changes. A model with gradual transitions can represent a progressive change in gastric conditions rather than an abrupt switch, but it requires a time-varying hazard or additional state description. These are alternative representations, not interchangeable calculations: the appropriate choice depends on the physiological information and the intended level of detail.

A practical check is to write a timeline before calculating anything. Mark dosing time, meal times, the assumed onset of each fed interval, and the assumed return to fasted state. Then verify that every time after dosing belongs to exactly one state. For a first-order model with nonnegative hazards, \(G(t)\) must be nonincreasing and must never exceed \(G(0)\). The amount emptied plus the amount remaining must equal the initial gastric amount:

$$ \bigl[G(0)-G(t)\bigr]+G(t)=G(0). $$

At each transition, check that the calculated gastric amount is unchanged while the rate updates to the new hazard times that amount. A jump in the amount itself usually indicates that the interval calculation was restarted incorrectly or that a transition was treated as a mass transfer. If a numerical calculation produces increasing gastric mass under positive emptying hazards, check the sign in the differential equation and the time units. For example, hazards in \(\mathrm{h}^{-1}\) must be multiplied by durations in hours.

A further limitation is the memoryless assumption. If emptying depends on how long the dose has been in the stomach, a state-specific constant hazard may not capture the process. Lagged or transit-time descriptions may require tracking residence time or another model-specific quantity through the state change. A meal-state timeline alone does not specify those details. Nor does a gastric emptying profile by itself establish how much drug reaches the systemic circulation; that requires the relevant downstream absorption and disposition model.

Common pitfall: “The dose was taken fasted, so its entire emptying profile is fasted.” This is wrong whenever the modeled gastric state changes while dose remains in the stomach. The correction is to apply the fasted hazard only up to the transition, carry the remaining amount forward continuously, and then use the new state-specific hazard. The converse error—treating a dose as fed for its entire residence because a meal preceded dosing—is corrected in the same way.
Key takeaway: Meal timing affects gastric emptying through the sequence and duration of gastric states experienced by the dose. Place meals and dosing on one timeline, apply the appropriate hazard on each interval, and preserve the remaining dose across every transition.

Check Your Understanding

Use the piecewise first-order model and show how the state intervals determine each calculation.

  1. A 50 mg dose is given at \(t=0\). It is fasted with \(k_F=0.50\ \mathrm{h}^{-1}\) until a meal at \(t=2\), then fed with \(k_{\mathrm{fed}}=0.10\ \mathrm{h}^{-1}\) until \(t=4\). What amount remains at 4 hours?
  2. In question 1, is the amount remaining allowed to jump at \(t=2\)? Explain.
  3. A meal has a modeled fed interval from \(t=-2\) to \(t=1\). A dose is given at \(t=0\). Which state-specific hazard applies from dosing until \(t=1\)?
  4. For a positive, constant state-specific hazard, name two checks that should hold for a calculated gastric amount profile.

Solutions. (1) The dose spends 2 hours fasted and 2 hours fed, so \(G(4)=50\exp[-(0.50)(2)-(0.10)(2)]=50e^{-1.20}=15.06\ \mathrm{mg}\). (2) No. The amount at \(t=2\) is the starting amount for the fed interval; only the hazard and emptying rate change. (3) The fed-state hazard applies, because the dose is given within the specified fed interval. (4) \(G(t)\) is nonincreasing and nonnegative, and the amount emptied plus the amount remaining equals the initial amount. These checks follow from positive emptying hazards and the stated mass-balance assumption.