Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Transcellular and Paracellular Permeation

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

Transcellular and Paracellular Permeation

Compare the two epithelial permeation routes, calculate their contributions to flux, and recognize when a route-specific model is needed.

Advanced 11 min read

What You'll Learn

  • Distinguish transcellular permeation through epithelial cells from paracellular permeation between them.
  • Define route-specific permeability and total flux under steady-state passive transport assumptions.
  • Derive the area-weighted effective permeability for parallel routes.
  • Calculate route contributions and identify which pathway limits permeation.
  • Recognize conditions that invalidate a single effective-permeability parameter.
  • Compare a lumped permeability model with explicit route-specific representations.

Why Separate Epithelial Routes?

Intestinal drug movement across the epithelium can occur through cells or between them. The distinction matters because these routes do not have the same physical barriers. Transcellular permeation crosses the apical membrane, the cell interior, and the basolateral membrane. Paracellular permeation passes through the aqueous spaces and junctional pathway between adjacent cells. Lipophilic compounds may permeate cell membranes readily, whereas small hydrophilic compounds may rely more heavily on the paracellular route, subject to its size and charge selectivity.

A measured or modeled effective permeability can summarize both routes when their fluxes respond linearly to the same concentration difference. But the summary can conceal which route carries the drug, and it can fail when route contributions change with location, concentration, time, or experimental conditions. This tutorial formulates the route contributions and sets out when a single parameter is a defensible simplification.

Definition: Transcellular permeation is movement across the epithelial cell layer through the cells. Paracellular permeation is movement through the pathway between adjacent cells. A route-specific permeability \(P_r\) relates the steady passive flux through route \(r\) to the concentration difference driving that route. Permeability has units of length per time, such as cm/s.

Flux Through Parallel Routes

Consider a local epithelial region at steady state. Let \(C_L\) be the free, permeating drug concentration at the luminal boundary and \(C_B\) the corresponding concentration at the basolateral boundary. For passive linear transport, the route-specific flux per unit area of that route is proportional to the concentration difference \(C_L-C_B\). Flux is positive from lumen to blood when \(C_L>C_B\).

The two pathways are treated as parallel routes: drug can cross through either route between the same sides of the epithelium. If \(f_t\) and \(f_p\) are the fractions of the chosen reference area assigned to transcellular and paracellular routes, respectively, their sum is one in this two-route representation. The total flux per unit reference area is the sum of the area-weighted route fluxes.

$$ J_{\mathrm{total}} = f_t P_t(C_L-C_B) + f_p P_p(C_L-C_B) $$

Here \(J_{\mathrm{total}}\) has units of amount per area per time when concentration is expressed as amount per volume. If both routes have the same boundary concentrations and remain linear, the equation can be written using a single effective permeability.

$$ P_{\mathrm{eff}} = f_tP_t+f_pP_p, \qquad J_{\mathrm{total}} = P_{\mathrm{eff}}(C_L-C_B) $$

The fractions and route permeabilities must refer to a consistent area convention. In this formulation, \(P_t\) and \(P_p\) are permeabilities per unit area of their respective routes, while \(P_{\mathrm{eff}}\) is normalized to the chosen reference area. The \(f\) terms account for how much of that reference area each route occupies. In practice, route-specific surface areas may be difficult to measure; an empirical \(P_{\mathrm{eff}}\) can therefore summarize geometry and transport together.

Worked Example: Combining Two Passive Pathways

Suppose a hypothetical epithelial model assigns 98% of its reference area to the transcellular route and 2% to the paracellular route. Let \(P_t=1.20\times10^{-4}\) cm/s and \(P_p=3.00\times10^{-5}\) cm/s. Assume steady passive transport, the same \(C_L-C_B\) for both routes, and no route-specific accumulation.

The area-weighted effective permeability is \(P_{\mathrm{eff}}=f_tP_t+f_pP_p\). Substituting the values gives

$$ P_{\mathrm{eff}} = (0.98)(1.20\times10^{-4}) + (0.02)(3.00\times10^{-5}) = 1.176\times10^{-4} + 6.00\times10^{-7} = 1.182\times10^{-4}\ \mathrm{cm/s}. $$

The transcellular contribution is \(1.176\times10^{-4}\) cm/s and the paracellular contribution is \(6.00\times10^{-7}\) cm/s. The transcellular share is \(1.176\times10^{-4}/1.182\times10^{-4}=0.995\), or about 99.5%. The numerical result is an effective permeability for the stated reference area and assumptions, not a claim that all intestinal tissue has these route fractions.

What Is Inside the Transcellular Permeability?

The transcellular pathway includes sequential barriers. In a simple steady-state passive model, the same flux passes through the apical membrane, the cell interior, and the basolateral membrane. If each barrier has a linear permeability defined using the concentration difference across that barrier, the total concentration drop is the sum of the individual drops. This makes the permeation resistances add in series.

$$ \frac{1}{P_t} = \frac{1}{P_{\mathrm{apical}}} + \frac{1}{P_{\mathrm{cell}}} + \frac{1}{P_{\mathrm{basolateral}}} $$

This expression assumes that there is no net accumulation within the cell at steady state, that the barriers can each be represented by a linear permeability, and that the route is passive. It is a useful conceptual decomposition, not a requirement that every PBPK model resolve three separate transcellular parameters. A fitted transcellular permeability may instead lump these resistances.

The equation also explains why a high permeability of one barrier does not guarantee high permeability of the complete route. A relatively low-permeability membrane or slow intracellular step can dominate the total resistance. In contrast, the paracellular route is commonly represented as a separate pathway with its own effective permeability, reflecting such features as accessible junctional area and pathway restrictions.

Worked Example: A Transcellular Bottleneck

For a hypothetical compound, take \(P_{\mathrm{apical}}=2.00\times10^{-4}\) cm/s, \(P_{\mathrm{cell}}=1.00\times10^{-4}\) cm/s, and \(P_{\mathrm{basolateral}}=4.00\times10^{-4}\) cm/s. Assume the barriers are in steady passive series as described above.

The total transcellular resistance is

$$ \frac{1}{P_t} = \frac{1}{2.00\times10^{-4}} + \frac{1}{1.00\times10^{-4}} + \frac{1}{4.00\times10^{-4}} = 5000+10000+2500 = 17500\ \mathrm{s/cm}. $$

Therefore \(P_t=1/17500=5.71\times10^{-5}\) cm/s. The cell-interior term contributes the largest resistance, \(10000\) s/cm, so it is the principal bottleneck in this hypothetical calculation. Increasing the apical permeability alone would not remove that dominant resistance.

This serial calculation describes barriers within one route. It is different from combining transcellular and paracellular pathways, which are parallel alternatives and whose area-weighted fluxes add.

When the Route Contributions Differ

The algebraic sum remains useful when route permeabilities differ greatly. The relative contribution of a route is its flux divided by total flux, not simply its area fraction. For a common positive concentration difference, the transcellular fraction is

$$ \phi_t = \frac{f_tP_t}{f_tP_t+f_pP_p}, \qquad \phi_p = \frac{f_pP_p}{f_tP_t+f_pP_p}. $$

These fractions sum to one under the parallel, linear model. A route can occupy only a small area yet dominate if its route-specific permeability is sufficiently large. Conversely, an extensive route can carry little flux if its permeability is low.

Worked Example: A Small-Area Route Can Dominate

Consider a hypothetical small hydrophilic compound in an epithelial model with \(f_t=0.95\), \(f_p=0.05\), \(P_t=2.00\times10^{-7}\) cm/s, and \(P_p=2.00\times10^{-5}\) cm/s. Assume both routes experience the same concentration difference and behave linearly.

The weighted contributions are

$$ f_tP_t=(0.95)(2.00\times10^{-7})=1.90\times10^{-7}\ \mathrm{cm/s}, $$ $$ f_pP_p=(0.05)(2.00\times10^{-5})=1.00\times10^{-6}\ \mathrm{cm/s}. $$

Thus \(P_{\mathrm{eff}}=1.19\times10^{-6}\) cm/s. The paracellular contribution is \(1.00\times10^{-6}/1.19\times10^{-6}=0.840\), or 84.0% of total flux, despite representing only 5% of the reference area. The route with the smaller area dominates because its route-specific permeability is much higher in this example.

When One Effective Permeability Is Inadequate

A single \(P_{\mathrm{eff}}\) is a useful reduction when it preserves the flux needed by the model over the conditions of interest. It is inadequate if the reduction cannot reproduce relevant changes in flux or if the route distinction is itself needed for interpretation. Important warning signs include:

  • Different concentration dependencies. The linear sum assumes each route flux is proportional to the same concentration difference. If either route has a nonlinear relationship, a constant weighted permeability may fit only a limited concentration range.
  • Different driving concentrations. If the routes experience different local free concentrations or concentration differences, one common \((C_L-C_B)\) cannot be factored out of the flux sum. This can occur when the model resolves local boundary layers or pathway-specific conditions.
  • Time-dependent intracellular behavior. The simple serial transcellular expression is a steady-state result. If drug temporarily accumulates in cells, the transcellular flux can depend on time and intracellular state rather than only on the current lumen-to-blood concentration difference.
  • Changing route contributions. A fitted effective permeability may not transfer if route properties or accessible areas differ across intestinal regions, experimental systems, or conditions. Regional differences in overall permeability do not by themselves identify which route caused them.
  • Route-specific questions. A total flux may be predicted adequately while the model still cannot answer whether the flux is transcellular or paracellular. Route-specific data or assumptions are needed to make that distinction.

An effective permeability measured in an assay is therefore not automatically a molecular constant. It may summarize the compound, the epithelial preparation, the available pathways, and the assay conditions. The earlier material on effective permeability and intestinal absorption, and on translating permeability assays to human \(P_{\mathrm{eff}}\), provides the broader context for using such inputs in PBPK models.

Deriving the Parallel-Pathway Result

The single-parameter reduction follows directly from flux conservation across the two parallel routes. It does not require the routes to have equal permeability; it does require them to share the same driving concentration difference and to be linear under the conditions being modeled.

1
Write each route flux.
For route \(r\), assume \(J_r=P_r(C_L-C_B)\), with constant \(P_r\) and the same boundary concentrations for both routes.
2
Weight by route area.
The flux per unit reference area contributed by route \(r\) is \(f_rJ_r\), because that route occupies fraction \(f_r\) of the reference area.
3
Add the parallel fluxes.
Conservation gives \(J_{\mathrm{total}}=f_tP_t(C_L-C_B)+f_pP_p(C_L-C_B)\).
4
Factor the common driving difference.
\(J_{\mathrm{total}}=(f_tP_t+f_pP_p)(C_L-C_B)\). Defining the bracketed term as \(P_{\mathrm{eff}}\) yields the area-weighted result.

If the two routes instead have distinct driving differences, the correct expression remains \(J_{\mathrm{total}}=f_tP_t\Delta C_t+f_pP_p\Delta C_p\). Unless those differences are equal or have a known fixed relationship, this cannot be represented by the same constant \(P_{\mathrm{eff}}\) multiplying one common concentration difference. This is the precise mathematical reason a single parameter can fail.

Model Choices, Diagnostics, and Limitations

RepresentationWhat it capturesMain limitation
Single effective permeabilityTotal flux with a compact parameter when conditions support a common linear driving differenceDoes not identify route contributions and may not transfer across conditions
Area-weighted parallel routesSeparate transcellular and paracellular contributions under common-gradient passive assumptionsRequires route-specific permeabilities or their defensible estimates
Resolved serial transcellular barriersCan represent which transcellular barrier is limitingNeeds more parameters and steady-state assumptions; may be weakly identifiable from total flux alone

Useful diagnostics include checking whether a fitted permeability remains stable across concentration ranges, assay conditions, and relevant intestinal regions. Systematic residual patterns or a permeability estimate that shifts substantially across these settings may indicate that the lumped model is absorbing changing route contributions or other omitted processes. Such patterns are evidence to investigate, not proof that a particular route is responsible: total-flux data alone generally cannot uniquely separate \(f_tP_t\) from \(f_pP_p\).

A common mistake is to conclude that a pathway contributes little because it occupies little area. That reasoning confuses area with flux. The correction is to compare the weighted conductances \(f_tP_t\) and \(f_pP_p\); the pathway with the larger weighted conductance contributes more under the common-gradient model. Another mistake is to add serial permeabilities as if they were parallel. Serial barriers add as resistances, whereas parallel route fluxes add.

Key takeaway: Under steady, linear passive transport with a common concentration difference, transcellular and paracellular fluxes add, giving \(P_{\mathrm{eff}}=f_tP_t+f_pP_p\). Use that effective parameter only for conditions where the assumptions hold; preserve route-specific terms when their contributions or driving forces differ.

Check Your Understanding

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

  1. For \(f_t=0.80\), \(P_t=4.00\times10^{-5}\) cm/s, \(f_p=0.20\), and \(P_p=1.00\times10^{-5}\) cm/s, calculate \(P_{\mathrm{eff}}\) under the common-gradient model.
  2. In a parallel-pathway model, what determines whether the transcellular or paracellular route contributes more flux?
  3. Why do serial apical, intracellular, and basolateral resistances combine differently from parallel transcellular and paracellular routes?
  4. Give one condition under which a constant effective permeability cannot represent both route fluxes using a single common concentration difference.

Solutions. (1) \(P_{\mathrm{eff}}=(0.80)(4.00\times10^{-5})+(0.20)(1.00\times10^{-5})=3.20\times10^{-5}+2.00\times10^{-6}=3.40\times10^{-5}\) cm/s. (2) Compare the weighted conductances \(f_tP_t\) and \(f_pP_p\); the larger one carries more flux when both routes share the same positive concentration difference. (3) The same flux passes successively through serial barriers, so their concentration drops, or resistances \(1/P\), add. Parallel routes provide alternative paths, so their area-weighted fluxes add. (4) If the routes have different local concentration differences that vary independently, the sum \(f_tP_t\Delta C_t+f_pP_p\Delta C_p\) cannot generally be reduced to one constant permeability multiplying one common difference.