Tutorials › Nonlinear Mixed-Effects Models (NLME) › Parallel Absorption Pathways and Fractional Dose Allocation

Absorption Models · Tutorial 76 of 797

Parallel Absorption Pathways and Fractional Dose Allocation

Learn to represent parallel absorption pathways with valid dose fractions and calculate the resulting systemic input without confusing dose allocation with bioavailability.

Intermediate 10 min read

What You'll Learn

  • Define parallel absorption pathways using dose fractions and pathway-specific absorption rates.
  • Check that pathway fractions form a valid allocation of the administered dose.
  • Calculate each pathway’s systemic input rate and the total input rate.
  • Distinguish allocated dose from the amount that ultimately reaches systemic circulation.
  • Identify common allocation, unit, and parameter-interpretation errors.

Why Represent More Than One Absorption Pathway?

An oral dose does not always behave as though all of it enters the body through one route with one absorption rate. A model may instead assign portions of the administered dose to distinct pathways. For example, one portion might be absorbed relatively quickly and another more slowly. The pathways are parallel because each portion proceeds through its own absorption process; one pathway does not have to finish before another begins.

The central accounting question is simple but essential: how much of the administered dose is assigned to each pathway, and do those allocations add up correctly? The rate parameters describe how quickly the allocated amounts leave their respective absorption depots. They do not determine how much dose was allocated in the first place.

Definition: For an administered dose \(D\), a model with \(J\) parallel pathways assigns fraction \(f_j\) of the dose to pathway \(j\), where \(f_j\geq 0\) and \(\sum_{j=1}^{J} f_j=1\). The amount assigned to pathway \(j\) is \(f_jD\). In a first-order pathway, \(k_{a,j}\) is that pathway’s absorption rate constant, with units of inverse time.

Allocation, Absorption Rate, and Systemic Input

For a single dose at time zero, let \(A_j(t)\) be the amount remaining in the absorption depot for pathway \(j\) at time \(t\). Under first-order absorption, the depot amount decreases in proportion to the amount still present. The systemic input rate is the amount entering systemic circulation per unit time.

$$ A_j(0)=f_jD,\qquad \frac{dA_j(t)}{dt}=-k_{a,j}A_j(t),\qquad R_{\mathrm{sys}}(t)=\sum_{j=1}^{J}F_j k_{a,j}A_j(t). $$

Here \(F_j\) is the fraction of the amount leaving pathway \(j\)'s depot that reaches systemic circulation. If the model assumes the same bioavailability \(F\) for every pathway, use \(F_j=F\). If complete systemic availability is assumed, use \(F_j=1\). These are assumptions about systemic availability, not about dose allocation: the allocation fractions still obey \(\sum_j f_j=1\).

The distinction matters because an allocated amount need not equal the amount that ultimately reaches systemic circulation. The dose assigned to pathway \(j\) is \(f_jD\), while the eventual systemic amount from that pathway is \(F_jf_jD\), under this model. If no pathway-specific loss is assumed, set \(F_j=1\); do not add an availability factor without stating that assumption.

A Short Mass-Balance Justification

The first-order depot equation has the solution \(A_j(t)=f_jD e^{-k_{a,j}t}\). Multiplying the amount leaving the depot per unit time by \(F_j\) gives the systemic input from that pathway. Integrating that input over time verifies the total systemic amount and makes the role of the fraction explicit.

$$ \begin{aligned} R_j(t)&=F_j k_{a,j}f_jD e^{-k_{a,j}t},\\ \int_0^\infty R_j(t)\,dt &=F_jf_jD\int_0^\infty k_{a,j}e^{-k_{a,j}t}\,dt\\ &=F_jf_jD\left[-e^{-k_{a,j}t}\right]_{0}^{\infty}\\ &=F_jf_jD. \end{aligned} $$

This calculation assumes \(k_{a,j}>0\), so that the depot empties over an indefinitely long interval, and a constant \(F_j\). Summing over pathways gives total systemic input \(D\sum_j f_jF_j\). If availability is common across pathways, this reduces to \(FD\), because the allocation fractions sum to one. The result checks the amount entering the systemic circulation; it is not a concentration profile.

Worked Example: Two Pathways With a Common Availability

Suppose a hypothetical \(500\)-mg dose is split between two first-order pathways. The fractions are \(f_1=0.70\) and \(f_2=0.30\), the absorption rates are \(k_{a,1}=0.40\ \mathrm{h}^{-1}\) and \(k_{a,2}=0.08\ \mathrm{h}^{-1}\), and both pathways have availability \(F=0.80\). Check the allocation, calculate the initial systemic input rate, and calculate the eventual systemic amount.

Allocation check. The fractions are nonnegative and \(0.70+0.30=1.00\), so they allocate the full dose. The pathway depot amounts at time zero are \(0.70(500)=350\) mg and \(0.30(500)=150\) mg.

Initial input rate. At time zero, each depot’s first-order outflow is its rate constant multiplied by its depot amount. Applying the common availability gives

$$ \begin{aligned} R_{\mathrm{sys}}(0) &=F\left(k_{a,1}f_1D+k_{a,2}f_2D\right)\\ &=0.80\left[(0.40\ \mathrm{h}^{-1})(350\ \mathrm{mg}) +(0.08\ \mathrm{h}^{-1})(150\ \mathrm{mg})\right]\\ &=0.80(140+12)\ \mathrm{mg/h}\\ &=121.6\ \mathrm{mg/h}. \end{aligned} $$

Eventual systemic amount. With common availability, the total is \(FD=(0.80)(500\ \mathrm{mg})=400\) mg. The two allocated amounts sum to \(500\) mg, but the model sends \(400\) mg into systemic circulation because \(20\%\) is unavailable.

Enforcing a Valid Allocation

Fractions are constrained quantities, not unconstrained coefficients. For a complete allocation, each must be between zero and one and their sum must be one. A set of numbers that violates these conditions does not describe a division of the administered dose. It should not be repaired by silently changing the interpretation or rescaling the values after calculations have begun.

Worked Example: Correcting an Invalid Proposed Split

A hypothetical \(400\)-mg dose is proposed to have fractions \(0.65\) and \(0.45\) in two pathways. Determine whether this is a valid allocation, then give one explicitly corrected allocation and its pathway amounts.

Both proposed fractions are nonnegative, but \(0.65+0.45=1.10\), not \(1\). They assign \(110\%\) of the administered dose: \(0.65(400)=260\) mg and \(0.45(400)=180\) mg, totaling \(440\) mg. That is not a valid complete allocation.

One possible correction, if the intended first fraction remains \(0.65\), is to assign the remainder to pathway 2: \(f_2=1-0.65=0.35\). The corrected amounts are \(0.65(400)=260\) mg and \(0.35(400)=140\) mg, which sum to \(400\) mg. This is a revised modeling decision, not an arithmetic consequence that makes the original pair valid.

If part of the dose is intentionally not assigned to the modeled pathways, represent that remainder explicitly rather than hiding it. For example, an unmodeled fraction \(f_{\mathrm{other}}\) can be included so that the fractions across all specified portions sum to one. The meaning of that remainder must be clear; it should not be confused with bioavailability, which describes how much of an allocated amount reaches systemic circulation.

Pathway-Specific Availability and Total Input

When availability differs by pathway, the total systemic amount is a weighted sum, not generally a single common availability multiplied by the dose. Each pathway contributes according to both its dose fraction and its availability. The instantaneous rate also depends on the pathway’s absorption rate constant. Keeping these roles separate prevents a frequent bookkeeping error: treating \(f_j\), \(F_j\), and \(k_{a,j}\) as interchangeable kinds of fractions or rates.

Worked Example: Three Pathways With Different Availability

Consider a hypothetical \(600\)-mg dose allocated across three pathways with fractions \(0.25\), \(0.50\), and \(0.25\). Let the respective absorption rates be \(0.50\), \(0.10\), and \(0.025\ \mathrm{h}^{-1}\), and let the pathway availabilities be \(0.80\), \(0.90\), and \(0.60\). Calculate the initial systemic input rate and eventual systemic amount.

Allocation check and depot amounts. The fractions are nonnegative and \(0.25+0.50+0.25=1.00\). The amounts allocated are \(150\) mg, \(300\) mg, and \(150\) mg.

Initial systemic input. For each pathway, multiply its depot amount by its rate constant and availability, then add the pathway-specific rates:

$$ \begin{aligned} R_{\mathrm{sys}}(0) &=(0.80)(0.50\ \mathrm{h}^{-1})(150\ \mathrm{mg})\\ &\quad +(0.90)(0.10\ \mathrm{h}^{-1})(300\ \mathrm{mg})\\ &\quad +(0.60)(0.025\ \mathrm{h}^{-1})(150\ \mathrm{mg})\\ &=60+27+2.25\\ &=89.25\ \mathrm{mg/h}. \end{aligned} $$

Eventual systemic amount. Weight each availability by its allocated amount: \((0.80)(150)+(0.90)(300)+(0.60)(150)=120+270+90=480\) mg. Thus \(480\) mg reaches systemic circulation under these assumptions, even though the three pathways together were assigned the full \(600\)-mg dose.

Choosing and Checking a Parallel-Pathway Representation

A parallel model is useful when the scientific model calls for portions of one dose to enter through distinct absorption processes. It differs from a single first-order pathway, which assigns the entire dose to one absorption depot and one rate constant. It also differs from sequential absorption, where material passes through one modeled process before another. A transit representation instead describes a distributed delay through linked stages; it is not equivalent to splitting the dose into independent pathway allocations.

These alternatives are not interchangeable merely because each can produce input spread over time. The choice should reflect the intended structure and available information. This tutorial specifies the dose split and input accounting, not how to derive concentration profiles or determine which pathway structure best fits data.

Before using a parallel-pathway specification, check the following:

  • Fractions: verify \(f_j\geq0\) for every pathway and \(\sum_j f_j=1\) when all of the dose is allocated.
  • Amounts: confirm that \(f_jD\) has units of dose and that the allocated amounts sum to \(D\).
  • Rates: ensure each \(k_{a,j}\) has units of inverse time; multiplying it by an amount gives an amount per time.
  • Availability: state whether \(F_j\) is common, pathway-specific, or set to one. Do not apply availability twice.
  • Total input: integrate or otherwise check the pathway contributions against \(D\sum_j f_jF_j\) under the stated assumptions.

A model can be mathematically valid yet poorly supported by observations. For example, an allocation fraction near zero makes its pathway contribute very little, so its rate may be difficult to distinguish from other model features. A fitted fraction outside \([0,1]\), or fractions that do not sum to one under a complete-allocation specification, indicates a parameterization or implementation problem rather than a plausible dose split. These checks do not establish that multiple pathways are needed; they only verify that the specified allocation and input accounting are coherent.

Common pitfall: A learner may see fractions \(0.65\) and \(0.45\), call them a \(65\%/45\%\) split, and proceed because each number is individually between zero and one. The correction is to check the whole allocation: \(0.65+0.45=1.10\), so these values allocate \(110\%\) of the dose. Revise the allocation explicitly or include a clearly defined additional portion; do not silently renormalize it.

The main limitation of this representation is that its arithmetic does not establish that the pathways correspond to distinct biological mechanisms. Different combinations of fractions, rates, and availability can sometimes imply similar overall input behavior. The model also relies on the specified dose allocation and availability assumptions; changing those assumptions changes the total systemic amount. Keep those limitations separate from the basic mass-balance check: a valid allocation is necessary for the intended model, but it is not evidence by itself that the model is appropriate.

Check Your Understanding

Answer each question using the stated assumptions. Brief solutions follow the questions.

  1. A \(250\)-mg dose is allocated with fractions \(0.40\) and \(0.60\). What amount is assigned to each pathway, and is the allocation valid?
  2. For two pathways, \(f_1=0.30\), \(f_2=0.70\), \(F_1=0.50\), and \(F_2=1.00\). What total systemic amount is expected from a \(200\)-mg dose?
  3. Why is \(k_{a,j}\) not a fraction of the dose? State its units and how it enters the initial pathway input rate.
  4. A proposed complete allocation has fractions \(0.20\), \(0.30\), and \(0.40\). What is wrong, and what remainder would be needed to make the allocation complete?

Solutions. (1) The amounts are \(0.40(250)=100\) mg and \(0.60(250)=150\) mg; the fractions are nonnegative and sum to one, so the allocation is valid. (2) The systemic amount is \(200[(0.30)(0.50)+(0.70)(1.00)]=200(0.85)=170\) mg. (3) \(k_{a,j}\) is a rate constant with units of inverse time, not a dose share. The initial systemic contribution is \(F_jk_{a,j}f_jD\), which has units of amount per time. (4) The fractions sum to \(0.90\), so \(10\%\) is unassigned. For a complete allocation, include an explicitly defined remainder fraction of \(0.10\), or revise the listed fractions so their sum is one.