Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Deriving Membrane Flux From Permeability

Permeability-Limited and Multi-Compartment Tissues · Tutorial 73 of 799

Deriving Membrane Flux From Permeability

Derive membrane flux from a concentration gradient, calculate transfer rates with consistent units, and distinguish permeability from the permeability-surface-area product.

Intermediate 10 min read

What You'll Learn

  • Derive passive membrane flux from Fick’s first law for a thin, homogeneous membrane.
  • Identify when flux is proportional to the concentration difference across a membrane.
  • Distinguish permeability, with units of length per time, from permeability-surface-area clearance, with units of volume per time.
  • Calculate flux per area and total transfer rate using consistent concentration and area units.
  • Recognize when partitioning, unstirred layers, or active transport make the simple relation inadequate.

From a Membrane Concentration Gradient to a Transfer Rate

When a drug crosses a cell membrane, the concentration difference between the two sides is not by itself a transfer rate. The rate also depends on how readily the drug crosses the barrier and how much membrane is available. A useful starting point is therefore to separate three quantities: flux through a unit of membrane area, the total membrane area, and the concentration difference that drives net transfer.

The relation derived here describes passive transfer under specified assumptions. It is a constitutive relation for membrane crossing, not a complete tissue model: it does not, on its own, describe how concentrations change over time in the surrounding compartments.

Definition: Let \(C_1\) and \(C_2\) be the free drug concentrations on the two sides of a membrane, expressed on the same concentration basis. Define \(J_{1\to2}\) as the net amount crossing per unit membrane area per unit time, positive from side 1 to side 2. Under a linear passive-transfer model, \(J_{1\to2}=P(C_1-C_2)\), where \(P\) is the membrane permeability. For a uniform membrane area \(A\), the net transfer rate is \(A J_{1\to2}=PA(C_1-C_2)\).

The sign convention carries information: if \(C_1>C_2\), net flux is from side 1 to side 2; if \(C_1<C_2\), it is from side 2 to side 1. If the two concentrations are equal under the model’s assumptions, the net flux is zero. This does not mean that individual molecules stop moving; it means opposing molecular movements balance.

Deriving the Linear Permeability Relation

Consider a flat membrane of thickness \(L\). Let \(x\) measure distance through the membrane, from side 1 at \(x=0\) to side 2 at \(x=L\). Fick’s first law states that diffusive flux within the membrane is proportional to the negative concentration gradient:

$$ J=-D_m\frac{dC_m}{dx}, $$

where \(C_m\) is drug concentration in the membrane and \(D_m\) is the drug’s diffusion coefficient there. The minus sign indicates movement down the concentration gradient. Suppose the membrane is homogeneous, \(D_m\) is constant, and the membrane concentration profile is at steady state. The gradient is then constant across the membrane, so:

$$ J=\frac{D_m(C_{m,1}-C_{m,2})}{L}. $$

Here \(C_{m,1}\) and \(C_{m,2}\) are the membrane concentrations at its two interfaces. If drug partitions between the aqueous solution and membrane with the same partition factor \(K_{m/w}\) at both interfaces, then \(C_{m,1}=K_{m/w}C_1\) and \(C_{m,2}=K_{m/w}C_2\). Substitution gives:

$$ J=\frac{D_mK_{m/w}}{L}(C_1-C_2). $$

This has the form \(J=P(C_1-C_2)\), with effective permeability \(P=D_mK_{m/w}/L\). This expression shows why permeability is not simply a diffusion coefficient: it also reflects membrane partitioning and thickness. Permeability may be measured or estimated directly; the expression in terms of \(D_m\), \(K_{m/w}\), and \(L\) is a simplified physical derivation.

Assumptions for the simple relation: Transfer is passive and non-saturable; the membrane is homogeneous and of constant thickness; its effective diffusion coefficient and partitioning are constant; the membrane concentration profile is approximately steady; and there is no important membrane metabolism, active transport, or additional transport resistance. \(C_1\) and \(C_2\) must be free concentrations on compatible reference bases. If these assumptions do not hold, a single constant \(P\) multiplying the concentration difference may not represent transfer adequately.

If the two sides have different partitioning into the membrane, the interface concentrations need not be \(K_{m/w}\) times the same factor on both sides. With interface factors \(K_1\) and \(K_2\), the derivation instead gives \(J=(D_m/L)(K_1C_1-K_2C_2)\). In that case, equal bulk concentrations do not necessarily imply zero net flux. More generally, the thermodynamic driving force is a difference in chemical potential; a concentration difference is a suitable linear driving force only when the concentration bases and partition conditions make it so.

Permeability, Flux, and Permeability-Surface-Area Clearance

Permeability \(P\) has units of length per time. For example, if \(J\) is measured in amount per area per time and concentration in amount per volume, then \(P=J/(C_1-C_2)\) has units of volume per area per time, which reduces to length per time. It describes transfer per unit area for a given concentration difference.

Multiplying permeability by area gives \(PA\), often called the permeability-surface-area product or permeability-surface-area clearance. Denote it here by \(CL_{PS}=PA\). Its units are volume per time, and the total net transfer rate is:

$$ \text{net transfer rate}=CL_{PS}(C_1-C_2)=PA(C_1-C_2). $$

This is “clearance-like” because a volume-per-time coefficient multiplies a concentration to produce an amount-per-time rate. It is not automatically systemic drug clearance or elimination: it quantifies exchange across the specified membrane under the stated conditions. Nor should \(P\) and \(CL_{PS}\) be interchanged. The former is an area-normalized membrane property; the latter combines that property with the area available for transfer.

QuantityMeaningTypical units
\(P\)Permeability; transfer capacity per unit area and concentration differencecm/s
\(J\)Net flux through each unit of areaamount/(cm²·s)
\(A\)Membrane area available for transfercm²
\(CL_{PS}=PA\)Permeability-surface-area productcm³/s, equivalent to mL/s
\(PA(C_1-C_2)\)Total net transfer rate across the specified areaamount/s

Worked Calculations

Worked Example: Deriving Permeability and Flux From Membrane Properties

Suppose a hypothetical drug has membrane diffusion coefficient \(D_m=2.00\times10^{-6}\) cm²/s, membrane-to-water partition factor \(K_{m/w}=0.500\), and membrane thickness \(L=1.00\times10^{-4}\) cm. A membrane area of 20.0 cm² separates solutions with free concentrations \(C_1=6.00\) μmol/L and \(C_2=2.00\) μmol/L. Assume the membrane is homogeneous, the profile is steady, and the two solutions use the same concentration basis.

First calculate the effective permeability:

$$ P=\frac{D_mK_{m/w}}{L} =\frac{(2.00\times10^{-6}\ \text{cm}^2/\text{s})(0.500)} {1.00\times10^{-4}\ \text{cm}} =0.0100\ \text{cm/s}. $$

The concentration difference is \(4.00\) μmol/L. Since \(1\) L equals \(1000\) cm³, this is \(0.00400\) μmol/cm³. Thus the flux per area is:

$$ J=P(C_1-C_2) =(0.0100\ \text{cm/s})(0.00400\ \text{μmol/cm}^3) =4.00\times10^{-5}\ \text{μmol/(cm}^2\text{s)}. $$

The total transfer rate is the flux times area:

$$ AJ=(20.0\ \text{cm}^2)(4.00\times10^{-5}\ \text{μmol/(cm}^2\text{s)}) =8.00\times10^{-4}\ \text{μmol/s}. $$

Equivalently, \(CL_{PS}=PA=(0.0100\ \text{cm/s})(20.0\ \text{cm}^2)=0.200\) cm³/s, or 0.200 mL/s. Multiplying by \(0.00400\) μmol/cm³ gives the same \(8.00\times10^{-4}\) μmol/s. The positive sign indicates net transfer from side 1 to side 2.

Worked Example: Converting a Permeability-Surface-Area Product

For a hypothetical barrier, let \(P=0.00200\) cm/s, \(A=100\) cm², \(C_1=25.0\) μmol/L, and \(C_2=10.0\) μmol/L. Assume passive linear transfer and compatible free-concentration measurements.

The permeability-surface-area product is:

$$ CL_{PS}=PA=(0.00200\ \text{cm/s})(100\ \text{cm}^2) =0.200\ \text{cm}^3/\text{s} =0.200\ \text{mL/s}. $$

Convert the driving concentration difference before multiplying. It is \(15.0\) μmol/L, or \(0.0150\) μmol/cm³. Therefore:

$$ \text{net transfer rate} =(0.200\ \text{cm}^3/\text{s})(0.0150\ \text{μmol/cm}^3) =0.00300\ \text{μmol/s} =0.180\ \text{μmol/min}. $$

The result is a rate of net transfer across this barrier, not a permeability: it includes the specified 100 cm² area. Reporting \(0.200\) cm³/s as \(P\) would be a dimensional and conceptual error.

Worked Example: Interpreting a Reversed Concentration Gradient

Consider a hypothetical membrane with \(P=0.00500\) cm/s and \(A=40.0\) cm². Let \(C_1=3.00\) μmol/L and \(C_2=9.00\) μmol/L. Use positive flux to mean transfer from side 1 to side 2, and assume the same linear passive model.

The signed concentration difference is \(C_1-C_2=-6.00\) μmol/L, equal to \(-0.00600\) μmol/cm³. Hence:

$$ J_{1\to2} =(0.00500\ \text{cm/s})(-0.00600\ \text{μmol/cm}^3) =-3.00\times10^{-5}\ \text{μmol/(cm}^2\text{s)}. $$

The total rate is \((40.0\ \text{cm}^2)J_{1\to2}=-0.00120\) μmol/s. The negative sign does not mean an impossible negative amount; it indicates net transfer in the direction opposite to the chosen positive direction, from side 2 toward side 1. The magnitude of the net rate is \(0.00120\) μmol/s.

When the Relation Is Useful—and When It Can Mislead

The simple relation is useful when a membrane barrier is the rate-limiting resistance and the relevant concentrations are known on both sides. It provides a direct way to relate a measured or estimated permeability to net transfer and clarifies which parameters change flux: increasing \(P\), available area, or the driving concentration difference increases the rate in direct proportion, provided the model remains valid.

The earlier material on perfusion-limited tissue models treats a different approximation: blood delivery is sufficiently rapid that tissue exchange is not limited by passage across a membrane. A permeability-based description is appropriate when membrane crossing itself may limit equilibration. These are alternative limiting descriptions, not interchangeable labels for the same parameter. A real tissue may involve both blood delivery and membrane crossing; deciding which is limiting requires evidence about its concentration profiles and transport conditions.

Several observations should prompt scrutiny of the constant-\(P\) approximation. A measured rate that is not proportional to the concentration difference may indicate saturable transport, concentration-dependent binding, active transport, or changing membrane properties. An apparent permeability that varies with experimental mixing or stirring can signal an unstirred aqueous layer adding resistance. If membrane thickness or composition varies substantially across the surface, a single uniform \(P\) may conceal spatial heterogeneity. These effects do not prove that the membrane relation is useless, but they require a more suitable effective model or a better-defined experimental system.

A further diagnostic is the predicted equilibrium condition. Under the simple same-basis model, the net flux must be zero when \(C_1=C_2\), and its direction must reverse when the concentration difference reverses. If the system instead has unequal equilibrium concentrations because of partitioning, using \(C_1-C_2\) alone gives the wrong driving force; concentrations must be adjusted to a common thermodynamic reference or the partition terms must be represented explicitly.

Common Pitfall

Common mistake: “The measured permeability is \(0.20\) mL/s, so \(P=0.20\) mL/s.” This treats \(PA\), which has units of volume per time, as if it were \(P\). Correction: permeability has units of length per time, such as cm/s. If \(0.20\) mL/s is the measured product and the relevant area is 20 cm², then \(P=(0.20\ \text{cm}^3/\text{s})/(20\ \text{cm}^2)=0.010\) cm/s. Always identify whether a reported parameter includes membrane area before interpreting or using it.

A related practical error is multiplying \(P A\) by a concentration expressed in μmol/L while leaving \(PA\) in cm³/s. The units will not cancel correctly unless the concentration is converted to μmol/cm³ or the volume units are made consistent. A quick dimensional check should yield amount per time for a transfer rate.

Check Your Understanding

Try each question first, then compare with the concise solutions below.

  1. A membrane has \(P=0.00400\) cm/s and \(C_1-C_2=5.00\) μmol/L. What is the flux per area in μmol/(cm²·s)?
  2. If \(P=0.00300\) cm/s and \(A=50.0\) cm², calculate \(CL_{PS}\) in mL/min.
  3. For a fixed membrane and area, what happens to the signed net flux if the concentration difference changes from \(C_1-C_2=+4\) to \(-4\), with the same units?
  4. State one condition under which \(J=P(C_1-C_2)\) is not an adequate description, and explain why.

Solutions. (1) Convert \(5.00\) μmol/L to \(0.00500\) μmol/cm³. Then \(J=(0.00400)(0.00500)=2.00\times10^{-5}\) μmol/(cm²·s). (2) \(CL_{PS}=PA=(0.00300)(50.0)=0.150\) cm³/s; multiplying by 60 gives \(9.00\) mL/min. (3) The flux reverses direction with the same magnitude, because it is proportional to the signed concentration difference. (4) For example, saturable transport violates the constant-permeability assumption: the rate may no longer increase in direct proportion to the concentration difference.