Coupling Dissolution Input to Luminal Transit

Why Dissolution Must Be Coupled to Luminal Transit

A luminal absorption model follows dissolved drug as it moves through gut segments, is absorbed, or exits a segment by transit. For a formulation that dissolves over time, the dissolved mass does not necessarily enter the first segment all at once. Instead, dissolution supplies dissolved drug to the segment where it occurs, while transit and absorption remove drug already dissolved there. The amount available for absorption at any time therefore depends on all three processes.

This tutorial represents dissolution as a specified, time-varying input to one or more segments. The dissolution input is prescribed rather than predicted from a modeled solid-drug state. That distinction keeps the model focused on how a source of dissolved drug interacts with luminal transit and absorption.

Definition: For segment \(i\), \(L_i(t)\) is the dissolved drug mass in the lumen at time \(t\), in mg. \(S_i(t)\) is the dissolution input rate into that segment, in mg/h. \(k_{a,i}\) is its first-order absorption rate constant, and \(k_{tr,i}\) is its first-order transit rate constant; both have units of h\(^{-1}\). The segmental absorption and transit fluxes are \(k_{a,i}L_i(t)\) and \(k_{tr,i}L_i(t)\), respectively, in mg/h.

A prescribed source must be physically compatible with the available formulation dose. If the total dose available to dissolve is \(D\) mg, then the cumulative dissolution input cannot exceed \(D\):

S_i(t) \geq 0, \qquad \sum_{i=1}^{m}\int_0^\infty S_i(t)\,dt \leq D.

The inequality allows some dose to remain undissolved or otherwise unavailable to this model. If the formulation is assumed to dissolve completely, the total integral equals \(D\). The input functions also specify where and when dissolved drug appears. They do not by themselves explain the physical mechanism that determines that pattern.

Coupled Segmental Mass Balances

Consider \(m\) sequential gut segments. Drug entering segment \(i\) from the preceding segment arrives at rate \(k_{tr,i-1}L_{i-1}(t)\). The first segment can also receive dissolved drug from a boundary input, such as drug dissolved before entering that segment. Denote this additional input by \(G_1(t)\); set \(G_1(t)=0\) if there is no such input. For later segments, \(G_i(t)=0\). The mass balances are:

\frac{dL_1}{dt} = S_1(t)+G_1(t) - \left(k_{a,1}+k_{tr,1}\right)L_1(t),

\frac{dL_i}{dt} = S_i(t)+k_{tr,i-1}L_{i-1}(t) - \left(k_{a,i}+k_{tr,i}\right)L_i(t), \qquad i=2,\ldots,m.

The outflow from the final segment, \(k_{tr,m}L_m(t)\), leaves this sequential model. The equations assume well-mixed dissolved drug within each segment; first-order absorption and transit; nonnegative, specified dissolution inputs; and no degradation or return of drug to the lumen. Rate constants are treated as constant in the equations above. If they vary with time, the same flux accounting applies, but the resulting equations generally need to be solved with the time-varying rates.

These equations show the competition directly. Increasing \(S_i(t)\) adds dissolved mass, whereas increasing either \(k_{a,i}\) or \(k_{tr,i}\) increases removal from segment \(i\). Transit removal supplies the next segment; absorption removal does not. Consequently, dissolution in an upstream segment can contribute to absorption both there and downstream.

Deriving the Single-Segment Response

The integrating-factor solution makes the effect of a time-varying dissolution input explicit. For one segment with no incoming boundary flux, write \(K=k_a+k_{tr}\). The equation is \(dL/dt=S(t)-KL(t)\). Multiplying by \(e^{Kt}\) and applying the product rule gives:

\frac{d}{dt}\left(e^{Kt}L(t)\right)=e^{Kt}S(t).

Integrating from time 0 to time \(t\), and then solving for \(L(t)\), yields:

L(t)=L(0)e^{-Kt}+\int_0^t S(\tau)e^{-K(t-\tau)}\,d\tau.

Each amount dissolved at time \(\tau\) contributes to luminal mass at time \(t\) after being depleted by both absorption and transit. The instantaneous fluxes are \(k_aL(t)\) and \(k_{tr}L(t)\), so their cumulative amounts through time \(t\) are \(k_a\int_0^t L(\tau)d\tau\) and \(k_{tr}\int_0^t L(\tau)d\tau\). This solution relies on a linear, well-mixed segment with constant rates; it does not require the dissolution rate itself to be constant.

Worked Example: Constant Dissolution Over a Finite Interval

Suppose a hypothetical 20 mg dose dissolves into one segment at \(S(t)=10\) mg/h from 0 to 2 h, and \(S(t)=0\) thereafter. Initially \(L(0)=0\). Let \(k_a=0.6\) h\(^{-1}\) and \(k_{tr}=0.4\) h\(^{-1}\). The dissolution input is dose-consistent because \(10\text{ mg/h}\times2\text{ h}=20\text{ mg}\). Calculate the dissolved mass and cumulative absorption and transit through 2 h.

Here \(K=0.6+0.4=1.0\) h\(^{-1}\). For \(0\leq t\leq2\), the single-segment solution is:

L(t)=\frac{10}{1.0}\left(1-e^{-1.0t}\right) =10\left(1-e^{-t}\right)\text{ mg}.

At 2 h, \(L(2)=10(1-e^{-2})=8.647\) mg. Cumulative absorption is the integral of the absorption flux:

A(2)=0.6\int_0^2 10(1-e^{-t})\,dt =6\left[2-(1-e^{-2})\right] =6.812\text{ mg}.

Cumulative transit is:

T(2)=0.4\int_0^2 10(1-e^{-t})\,dt =4\left[2-(1-e^{-2})\right] =4.541\text{ mg}.

The mass balance is \(20.000=8.647+6.812+4.541\) mg, up to rounding: input equals mass still in the segment plus absorption plus transit. The input stops at 2 h, but the 8.647 mg remaining continues to be absorbed and transit afterward.

How Sequential Transit Changes Eventual Absorption

With constant rates, a dissolved amount entering a segment faces competing first-order absorption and transit. In segment \(i\), the probability of absorption before transit is \(q_{a,i}=k_{a,i}/(k_{a,i}+k_{tr,i})\), and the probability of transit first is \(q_{tr,i}=k_{tr,i}/(k_{a,i}+k_{tr,i})\). These proportions follow from integrating the corresponding fluxes while an amount decays at total rate \(k_{a,i}+k_{tr,i}\).

Let \(F_i\) be the eventual fraction absorbed from a dissolved amount entering segment \(i\), assuming it encounters the modeled segments and rates until it is absorbed or exits the final segment. It can be absorbed locally, or transit onward and be absorbed later. Thus \(F_i=q_{a,i}+q_{tr,i}F_{i+1}\), with \(F_{m+1}=0\) because drug leaving the last segment is not absorbed in the model. Repeated substitution gives:

F_i = 1-\prod_{j=i}^{m}q_{tr,j} = 1-\prod_{j=i}^{m} \frac{k_{tr,j}}{k_{a,j}+k_{tr,j}}.

This result assumes constant rates, sequential transfer, and no other loss or return pathway. It gives an eventual fraction, not the time course of absorption. In this setting, a source \(S_i(t)\) contributes eventual absorbed mass \(F_i\int_0^\infty S_i(t)dt\), provided rates remain constant as the drug moves through the model. Timing still matters for the luminal mass and absorption profile at finite times.

Worked Example: Dissolution in an Upstream Segment

Consider two hypothetical segments. In segment 1, \(k_{a,1}=0.5\) h\(^{-1}\) and \(k_{tr,1}=0.5\) h\(^{-1}\). In segment 2, \(k_{a,2}=1.0\) h\(^{-1}\) and \(k_{tr,2}=0.5\) h\(^{-1}\). An 8 mg dose dissolves into segment 1 at 8 mg/h from 0 to 1 h. There is no initial dissolved mass.

During the first hour, \(L_1(t)=8(1-e^{-t})\) mg. Therefore \(L_1(1)=8(1-e^{-1})=5.057\) mg. Absorption and transit fluxes in segment 1 are equal because their rate constants are equal. Each cumulative amount through 1 h is:

0.5\int_0^1 8(1-e^{-t})\,dt =4e^{-1} =1.472\text{ mg}.

The first-hour mass check is \(8.000=5.057+1.472+1.472\) mg. Some mass remains in segment 1 at 1 h, so the first-hour fluxes alone do not describe its eventual fate.

For drug entering segment 1, \(q_{tr,1}=0.5/(0.5+0.5)=0.5\). In segment 2, \(q_{tr,2}=0.5/(1.0+0.5)=1/3\). Hence \(F_1=1-(0.5)(1/3)=5/6\). The full 8 mg input therefore yields \(8(5/6)=6.667\) mg eventual absorption, and \(8(1/6)=1.333\) mg eventual exit from the final segment. This result includes drug that remains in segment 1 after dissolution stops and later transits to segment 2.

Worked Example: Exponentially Declining Dissolution in a Distal Segment

Suppose a hypothetical 6 mg amount dissolves directly into segment 2 according to \(S_2(t)=6e^{-t}\) mg/h for \(t\geq0\). Its total input is \(\int_0^\infty6e^{-t}dt=6\) mg. Let \(k_{a,2}=0.8\) h\(^{-1}\), \(k_{tr,2}=0.2\) h\(^{-1}\), and \(L_2(0)=0\). Calculate the mass and fluxes at 1 h, and the eventual absorbed amount.

The total removal rate is 1.0 h\(^{-1}\). Substituting the time-varying source in the single-segment solution gives:

L_2(t)=\int_0^t6e^{-\tau}e^{-(t-\tau)}\,d\tau =6te^{-t}\text{ mg}.

Thus \(L_2(1)=6e^{-1}=2.207\) mg. The amount dissolved by 1 h is \(6(1-e^{-1})=3.793\) mg. Cumulative absorption is:

A(1)=0.8\int_0^1 6te^{-t}\,dt =4.8\left(1-\frac{2}{e}\right) =1.268\text{ mg}.

Cumulative transit is \(T(1)=0.2\int_0^1 6te^{-t}dt=0.3171\) mg. The first-hour balance is \(3.793=2.207+1.268+0.317\) mg. Eventually, the absorbed fraction is \(F_2=0.8/(0.8+0.2)=0.8\), so eventual absorption is \(6(0.8)=4.800\) mg. An instantaneous 6 mg input would have the same eventual absorbed amount under these constant rates, but a different early absorption profile.

Diagnostics, Alternatives, and Limitations

A useful diagnostic is to track input, luminal mass, absorption, and transit together. Summing the segmental balances cancels internal transit transfers: mass leaving one segment and entering the next is not a loss from the whole modeled gut. The cumulative input must equal the sum of all segmental luminal masses, cumulative absorption, and cumulative exit from the final segment, with any initial luminal mass included. A discrepancy indicates an implementation or accounting error.

Also inspect whether each \(S_i(t)\) is nonnegative, has units of mass per time, and integrates to no more than the dose available for dissolution. A source that continues after its prescribed dose has been exhausted creates drug. A source assigned to a segment with no plausible formulation delivery may also misrepresent the intended model. A negative luminal mass or absorption flux is inconsistent with these mass balances and usually signals a numerical or coding problem.

A simple alternative is an instantaneous-dissolution approximation: place the full dissolved dose in the model at time zero. This can be reasonable when dissolution is rapid relative to transit and absorption, but it cannot represent delayed input or ongoing dissolution. A prescribed time-varying source allows the timing and location of dissolution to be specified while keeping the model linear. Its limitation is that the dissolution profile is supplied rather than predicted; it cannot respond within this formulation to changes in luminal conditions or to depletion of a modeled solid pool.

Finally, the eventual-fraction calculation is not a substitute for a time-course calculation. Two sources with the same total mass and location can have the same eventual absorbed fraction under constant rates while producing different absorption timing. When the rates vary over time, the constant-rate segment fractions and product formula need not apply; solve the coupled balances with the specified time-varying rates instead.

Common pitfall: Treating dissolution as an extra amount to add to a segment without also including it in the mass-balance input, or counting transit out of one segment as both a whole-gut loss and an input to the next segment. Dissolution adds mass to the dissolved pool; internal transit transfers it; absorption and final-segment exit remove it from the modeled luminal system.

Check Your Understanding

Answer each question, then compare with the solutions below.

A segment has \(k_a=0.3\) h\(^{-1}\) and \(k_{tr}=0.7\) h\(^{-1}\). What fraction of a dissolved amount entering it is absorbed before transit?
A 12 mg dose dissolves into one segment at 3 mg/h for 4 h. Is the input dose-consistent? What is the segmental mass-balance equation if its initial dissolved mass is zero?
In a two-segment model, the transit-first fractions are \(0.4\) and \(0.25\). What fraction entering segment 1 is eventually absorbed, assuming the stated constant-rate sequential model?
Why can two dissolution profiles with the same total input have equal eventual absorption but different absorption by 1 h?

Solutions. (1) The total removal rate is \(0.3+0.7=1.0\) h\(^{-1}\), so the fraction absorbed first is \(0.3/1.0=0.300\). (2) The input is \(3\text{ mg/h}\times4\text{ h}=12\) mg, so it is dose-consistent. The balance is \(dL/dt=3-(k_a+k_{tr})L\) during the 4 h input interval, with \(L(0)=0\); afterward the source is zero. (3) The eventual absorbed fraction is \(1-(0.4)(0.25)=0.900\). (4) Dissolution timing changes when dissolved mass is available to be absorbed, so the integrated absorption flux by 1 h can differ. Under constant rates and the same source location, the eventual fraction depends on the total input and segmental rates, not on the timing of that input.