Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Effective Permeability and Intestinal Absorption

Intestinal Permeability and First-Pass Gut Metabolism · Tutorial 268 of 799

Effective Permeability and Intestinal Absorption

Learn how effective permeability converts a luminal concentration gradient and intestinal surface area into enterocyte uptake in a PBPK model.

Advanced 10 min read

What You'll Learn

  • Define effective intestinal permeability and state its units and operational meaning.
  • Distinguish an in vivo effective parameter from apparent permeability measured in an assay.
  • Formulate luminal-to-enterocyte flux with a concentration gradient and intestinal surface area.
  • Derive the first-order absorption rate that follows under sink and well-mixed assumptions.
  • Check mass balance, units, and model behavior when diagnosing permeability implementations.
  • Recognize when a single effective-permeability term is an inadequate representation.

From Dissolved Drug to Enterocyte Uptake

Dissolution provides drug in the intestinal lumen, but dissolved drug is not necessarily absorbed at the same rate that it becomes available. The rate of transfer across the luminal boundary depends on both the concentration driving force and how readily the intestinal interface permits transport. In a PBPK model, effective intestinal permeability represents that transfer in a form that can be combined with the segment’s luminal concentration and available surface area.

The term “effective” matters. It is a model-scale parameter for uptake from the lumen under a specified representation of the gut, not a universal molecular constant. It should not be confused with the apparent permeability reported by a cell or membrane assay. The latter describes transport in that assay’s particular setup; it is not automatically the parameter required by an intestinal PBPK compartment.

Definition: Effective intestinal permeability, \(P_{\mathrm{eff}}\), is an operational permeability coefficient used to describe net drug transfer from the intestinal lumen toward the enterocyte compartment per unit area and concentration difference. Its usual units are length per time, such as cm/s. It summarizes the net uptake behavior represented at the model boundary; it is not, by itself, a measure of dissolution, systemic bioavailability, or gut metabolic escape.

For the uptake boundary, let \(C_{\mathrm{L,free}}\) be the free dissolved drug concentration in the lumen. Let \(C_{\mathrm{E}}^*\) be the enterocyte-side concentration expressed on the same concentration basis as the luminal driving force. This equivalent concentration allows the two sides of the boundary to be compared without assuming that their raw compartment concentrations have identical reference scales. With \(A_{\mathrm{int}}\) denoting the effective intestinal surface area, the net flux density and total transfer rate are

$$ J_{\mathrm{L}\rightarrow\mathrm{E}} = P_{\mathrm{eff}}\left(C_{\mathrm{L,free}}-C_{\mathrm{E}}^*\right), \qquad R_{\mathrm{L}\rightarrow\mathrm{E}} = A_{\mathrm{int}}J_{\mathrm{L}\rightarrow\mathrm{E}}. $$

Here \(J_{\mathrm{L}\rightarrow\mathrm{E}}\) has units of amount per area per time, and \(R_{\mathrm{L}\rightarrow\mathrm{E}}\) has units of amount per time. A positive concentration difference gives net uptake; a negative difference gives net transfer in the opposite direction. Under a sink approximation, enterocyte-side concentration is negligible relative to luminal concentration, so \(C_{\mathrm{E}}^*\approx 0\) and \(R_{\mathrm{L}\rightarrow\mathrm{E}}\approx P_{\mathrm{eff}}A_{\mathrm{int}}C_{\mathrm{L,free}}\). The sink approximation is a modeling assumption, not part of the definition of permeability.

How Permeability Enters a PBPK Mass Balance

Suppose a luminal segment is well mixed, has fixed volume \(V_{\mathrm{L}}\), and contains free dissolved amount \(A_{\mathrm{L,free}}\), with \(C_{\mathrm{L,free}}=A_{\mathrm{L,free}}/V_{\mathrm{L}}\). If the enterocyte-side concentration is negligible and no other process changes the free luminal amount during the interval considered, the luminal mass balance is

$$ \frac{dA_{\mathrm{L,free}}}{dt} = -R_{\mathrm{L}\rightarrow\mathrm{E}} = -\frac{P_{\mathrm{eff}}A_{\mathrm{int}}}{V_{\mathrm{L}}}A_{\mathrm{L,free}}. $$

This equation shows why permeability is often said to control the absorption rate: with area and volume fixed, it sets the first-order rate constant for removal from the free luminal pool. It does not imply that total luminal drug always follows a single exponential. Dissolution, precipitation, transit, binding, and other represented processes can replenish or remove drug and alter the time course. In a coupled PBPK model, those processes must appear as their own mass-balance terms.

Model assumption: The first-order form \(R_{\mathrm{L}\rightarrow\mathrm{E}}=P_{\mathrm{eff}}A_{\mathrm{int}}C_{\mathrm{L,free}}\) requires a sink on the enterocyte side, a well-mixed luminal concentration, and fixed \(P_{\mathrm{eff}}\) and \(A_{\mathrm{int}}\) over the interval. It describes transfer out of the lumen, not the subsequent fate of drug in the enterocyte.

Worked Example: Calculating Luminal-to-Enterocyte Flux

Consider hypothetical values for a well-mixed intestinal segment: \(P_{\mathrm{eff}}=2.0\times10^{-4}\) cm/s, \(A_{\mathrm{int}}=600\) cm\(^2\), \(C_{\mathrm{L,free}}=0.015\) mg/cm\(^3\), and \(C_{\mathrm{E}}^*=0.005\) mg/cm\(^3\). Assume the effective permeability is constant over the calculation interval and that the stated concentrations are on a common driving-force basis.

The net concentration difference is \(0.015-0.005=0.010\) mg/cm\(^3\). Thus the flux density is

$$ J_{\mathrm{L}\rightarrow\mathrm{E}} = (2.0\times10^{-4}\ \mathrm{cm/s})(0.010\ \mathrm{mg/cm^3}) = 2.0\times10^{-6}\ \mathrm{mg/(cm^2\,s)}. $$

Multiplying by surface area gives the transfer rate:

$$ R_{\mathrm{L}\rightarrow\mathrm{E}} = (600\ \mathrm{cm^2})(2.0\times10^{-6}\ \mathrm{mg/(cm^2\,s)}) = 1.2\times10^{-3}\ \mathrm{mg/s}. $$

The positive rate indicates net transfer from lumen toward enterocytes. In an hour, this rate would correspond to \(0.0012\times3600=4.32\) mg only if the rate remained constant for the entire hour. In a dynamic model, concentrations and therefore the rate can change as drug is absorbed or replenished.

Deriving the Sink-Condition Absorption Time Course

A useful check on the compartment equation is to solve it under its stated assumptions. Let \(k_{\mathrm{abs}}=P_{\mathrm{eff}}A_{\mathrm{int}}/V_{\mathrm{L}}\). This quantity has units of inverse time, since permeability has units of length per time and area divided by volume has units of inverse length. The mass balance becomes \(dA_{\mathrm{L,free}}/dt=-k_{\mathrm{abs}}A_{\mathrm{L,free}}\).

1
Separate variables.
For positive luminal amount, divide by \(A_{\mathrm{L,free}}\) and multiply by \(dt\): \(dA_{\mathrm{L,free}}/A_{\mathrm{L,free}}=-k_{\mathrm{abs}}dt\).
2
Integrate from the initial amount.
Integrate from time 0 and amount \(A_{\mathrm{L,free}}(0)\) to time \(t\) and amount \(A_{\mathrm{L,free}}(t)\). This gives \(\ln[A_{\mathrm{L,free}}(t)/A_{\mathrm{L,free}}(0)]=-k_{\mathrm{abs}}t\).
3
Exponentiate and interpret.
Therefore \(A_{\mathrm{L,free}}(t)=A_{\mathrm{L,free}}(0)e^{-k_{\mathrm{abs}}t}\). The fraction remaining is \(e^{-k_{\mathrm{abs}}t}\), and the fraction removed by uptake is \(1-e^{-k_{\mathrm{abs}}t}\), provided no other process changes this pool.

Worked Example: Checking the Segmental Absorption Rate

Use hypothetical parameters \(P_{\mathrm{eff}}=1.5\times10^{-4}\) cm/s, \(A_{\mathrm{int}}=800\) cm\(^2\), and \(V_{\mathrm{L}}=200\) cm\(^3\). Assume a well-mixed lumen, sink conditions, and fixed parameters for 30 minutes.

The first-order rate constant is

$$ k_{\mathrm{abs}} = \frac{(1.5\times10^{-4}\ \mathrm{cm/s})(800\ \mathrm{cm^2})}{200\ \mathrm{cm^3}} = 6.0\times10^{-4}\ \mathrm{s^{-1}}. $$

Thirty minutes is \(1{,}800\) s, so the modeled fraction remaining is \(e^{-(6.0\times10^{-4})(1800)}=e^{-1.08}\approx0.340\). The fraction removed from the free luminal pool by uptake is therefore \(1-0.340=0.660\), or about 66.0%. This is not automatically the fraction of an administered dose absorbed: the calculation excludes any drug that is not in the free luminal pool and assumes no competing processes or replenishment.

Effective Permeability and Apparent Assay Permeability

Apparent permeability, commonly denoted \(P_{\mathrm{app}}\), is calculated from transport measured in a specified in vitro assay. For a typical donor-to-receiver assay, its operational calculation is

$$ P_{\mathrm{app}} = \frac{1}{A_{\mathrm{assay}}C_{\mathrm{donor}}} \frac{dQ_{\mathrm{receiver}}}{dt}, $$

where \(A_{\mathrm{assay}}\) is the assay membrane area, \(C_{\mathrm{donor}}\) is the donor concentration used for normalization, and \(dQ_{\mathrm{receiver}}/dt\) is the measured rate of appearance in the receiver. It also has units of length per time. Matching units does not make \(P_{\mathrm{app}}\) interchangeable with \(P_{\mathrm{eff}}\): the assay and PBPK parameter refer to different experimental or model systems, with potentially different barriers, surface-area definitions, concentration conditions, and sources of loss or retention.

Worked Example: Calculating an Assay Value Without Equating It to Peff

In a hypothetical assay, the measured receiver appearance rate is \(0.012\) \(\mu\)g/s, the membrane area is \(1.5\) cm\(^2\), and the donor concentration is \(0.40\) \(\mu\)g/cm\(^3\). Suppose the assay definition above applies and the rate is stable over the measurement interval.

Substitution gives

$$ P_{\mathrm{app}} = \frac{0.012\ \mathrm{\mu g/s}} {(1.5\ \mathrm{cm^2})(0.40\ \mathrm{\mu g/cm^3})} = \frac{0.012}{0.60}\ \mathrm{cm/s} = 0.020\ \mathrm{cm/s}. $$

This calculation yields the apparent permeability for that assay setup. It does not establish that the intestinal PBPK model should use \(P_{\mathrm{eff}}=0.020\) cm/s. Translating assay results into a human effective parameter requires additional assumptions and evidence; the numerical value cannot be transferred merely because both quantities have permeability units.

Diagnostics, Alternatives, and Limitations

A permeability implementation should pass simple physical and numerical checks before it is used to interpret absorption. First, confirm units: \(P_{\mathrm{eff}}A_{\mathrm{int}}C\) must have units of amount per time. Second, verify signs and mass conservation: a positive luminal-to-enterocyte transfer rate must reduce luminal amount by exactly the amount added to the receiving compartment at that interface. Third, inspect limiting behavior. With zero free luminal concentration and a sink assumption, uptake must be zero; with zero permeability or zero area, transfer must also be zero. If the model predicts uptake when all three conditions eliminate the driving force, its flux expression or bookkeeping is inconsistent.

A common warning sign is an apparent absorption rate that changes when concentration units are changed, even though the physical state is unchanged. This often indicates that permeability, concentration, or area units have been mixed. Another is an implausibly large uptake rate caused by applying an effective area twice—for example, using an area-scaled transfer coefficient and then multiplying by area again. Finally, if luminal amount falls faster than enterocyte amount rises, or vice versa, check whether another explicitly modeled process accounts for the difference; otherwise the interface violates mass conservation.

The main alternative to a permeability-driven flux is a direct first-order absorption term, such as \(R=k_{\mathrm{abs}}A_{\mathrm{L,free}}\). Under the well-mixed sink assumptions, these forms are mathematically equivalent when \(k_{\mathrm{abs}}=P_{\mathrm{eff}}A_{\mathrm{int}}/V_{\mathrm{L}}\). The permeability form makes the roles of area and concentration explicit; a direct rate constant can be simpler when only the lumped removal rate is needed. Neither form is automatically more mechanistic: the choice depends on what the model needs to represent and which parameters can be supported.

An effective coefficient is a lumped description. It cannot, on its own, identify which physical barrier controls uptake, describe distinct transport routes, or determine the drug’s later enterocyte disposition. A fixed value may also be inadequate if the represented permeability, area, or concentration driving force changes substantially with location or conditions. Those limitations should be addressed through a richer model only when the data and intended predictions justify the added structure; an effective parameter should not be interpreted as a direct measurement of every underlying process.

Common pitfall: “The assay reports \(P_{\mathrm{app}}\) in cm/s, so I can use that number directly as \(P_{\mathrm{eff}}\).” The units are compatible, but the quantities are defined for different systems. Use \(P_{\mathrm{app}}\) as an assay result; use a justified effective parameter in the PBPK flux equation. A numerical conversion requires its own rationale and evidence.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A segment has \(P_{\mathrm{eff}}=1.0\times10^{-4}\) cm/s, \(A_{\mathrm{int}}=500\) cm\(^2\), and concentration difference \(0.008\) mg/cm\(^3\). What is the net uptake rate?
  2. Under sink conditions, what happens to the predicted uptake rate if free luminal concentration is doubled while the other inputs are fixed?
  3. Why does an assay \(P_{\mathrm{app}}\) value not automatically provide the PBPK model’s \(P_{\mathrm{eff}}\), even when both are reported in cm/s?
  4. For a well-mixed, fixed-volume luminal segment under sink conditions, state the expression for the fraction removed by uptake after time \(t\).

Solutions. (1) \(R=P_{\mathrm{eff}}A_{\mathrm{int}}\Delta C=(1.0\times10^{-4})(500)(0.008)=4.0\times10^{-4}\) mg/s. The units reduce to mg/s. (2) The sink-condition rate is proportional to \(C_{\mathrm{L,free}}\), so it doubles. (3) The assay and PBPK parameter describe different systems and may use different barriers, areas, and concentration conditions; equal units do not establish equal values. (4) With \(k_{\mathrm{abs}}=P_{\mathrm{eff}}A_{\mathrm{int}}/V_{\mathrm{L}}\), the fraction removed is \(1-e^{-k_{\mathrm{abs}}t}\), assuming uptake is the only process changing the free luminal pool.