Week 5: Volume of Distribution — Mechanisms, Determinants, and Mathematical Foundations
Lecture Notes
1 Overview
This lecture provides a comprehensive summary of the volume of distribution, consolidating material from the preceding sessions. The discussion covers drug transport across membranes, the rate and extent of tissue distribution, protein binding, and the mathematical frameworks that underpin compartmental and physiologically-based pharmacokinetic (PBPK) modeling. Volume of distribution and clearance are the two fundamental pharmacokinetic parameters; without a firm grasp of both, it is not possible to meaningfully interpret or model pharmacokinetic data.
2 1. Drug Transport Across Membranes
2.1 1.1 Types of Membrane Transport
Biological membranes are multi-layered structures. Each layer comprises cells with distinct properties, and a drug must traverse these layers to reach systemic circulation. Transport occurs via two primary mechanisms:
- Paracellular transport: The drug passes through the gaps between cells, moving from the lumen into the blood.
- Transcellular transport: The drug passes directly through the cells themselves. This can occur by:
- Passive diffusion — driven by a concentration gradient.
- Facilitated (active) transport — mediated by carrier proteins and may require cellular energy.
2.2 1.2 Fick’s First Law of Diffusion
Passive diffusion is the movement of a substance down its concentration gradient. The net rate of transport across a membrane is described by:
\[\text{Rate} = P \times SA \times (C_1 - C_2)\]
where: - \(P\) = permeability (units: cm/min, a velocity representing the ease of membrane penetration) - \(SA\) = surface area (units: cm²) - \(C_1 - C_2\) = concentration gradient
The product \(P \times SA\) has units of cm³/min, which is volume per unit time — equivalent to clearance. This is a foundational insight: the permeability-surface area product is clearance in its most fundamental form.
Extending this: since \(P \times SA\) is volume/time and the concentration gradient has units of amount/volume, the overall rate has units of amount/time — consistent with a mass flow rate. Units must always reconcile; this is an essential discipline in pharmacokinetic modeling.
Fick’s law can also be expressed as:
\[J = -\frac{D \cdot K_{\text{ow}} \cdot \Delta C}{h}\]
where \(J\) is flux (the rate of movement per unit area), \(D\) is the diffusion coefficient, \(K_{\text{ow}}\) is the octanol–water partition coefficient, \(\Delta C\) is the concentration gradient, and \(h\) is the membrane thickness.
In the context of absorption from the gut, the rate of drug absorption into a compartment is commonly written as:
\[\frac{dG}{dt} = -k_a \cdot G\]
where \(k_a\) is the first-order absorption rate constant. This parameter encapsulates the entire complexity of Fick’s law under the simplifying assumption of first-order kinetics. Advanced gut absorption models (e.g., GastroPlus/ACAT, Simcyp/ADAM) represent the full mechanistic detail that \(k_a\) compresses into a single parameter.
2.3 1.3 Drug Properties Determining Permeability
Three drug properties govern membrane permeability:
1. Molecular size Smaller molecules traverse membranes more readily. As a rule of thumb, for two drugs of equal lipophilicity, doubling the molecular weight reduces permeability approximately 300-fold. This is why small molecules remain highly tractable from a drug discovery standpoint.
2. Lipophilicity Lipophilicity is measured as the octanol–water partition coefficient (\(K_{\text{ow}}\)): the drug is added to a mixture of octanol and water, the vessel is shaken, and the resulting concentration ratio is measured. Greater lipophilicity increases permeability. Molecular size and lipophilicity are interdependent — a small, polar molecule may have limited membrane penetration despite its size.
3. Degree of ionization The ionization state of a drug is governed by the Henderson–Hasselbalch equation. For a weak acid:
\[\text{pH} = \text{p}K_a + \log_{10}\frac{[A^-]}{[HA]}\]
For a weak base, the ratio is inverted (ionized in the denominator). The \(\text{p}K_a\) is a drug-specific property — it is the pH at which the concentrations of ionized and unionized forms are equal (i.e., the ratio equals 1, and \(\log_{10}(1) = 0\), so \(\text{pH} = \text{p}K_a\)).
As the pH of the environment changes along the gastrointestinal tract (from acidic in the stomach to neutral/alkaline in the intestine), the fraction of drug in each ionization state shifts accordingly. The Henderson–Hasselbalch equation allows formulation scientists to design drug release profiles that leverage these pH transitions.
3 2. Carrier-Mediated Transport
3.1 2.1 Equilibrating vs. Concentrating Transporters
Beyond passive diffusion, drugs may be transported by carrier proteins. These fall into two classes:
Equilibrating (passive facilitated) transporters - Driven by concentration gradient - Move drug bidirectionally across membranes - No energy (ATP) requirement - Net movement continues until equilibrium is reached on both sides of the membrane
Concentrating (active) transporters - ATP-dependent - Can move drug against a concentration gradient - Directionally specific (either uptake into or efflux out of cells) - Examples: P-glycoprotein (P-gp, an efflux transporter); organic anion-transporting polypeptides (OATPs, uptake transporters)
These two transporter classes often operate in tandem. An equilibrating transporter attempts to establish equilibrium, while a concentrating transporter actively counteracts this by driving drug out of (or into) a compartment. Simultaneously, systemic clearance is continuously removing drug. The interplay of these processes gives rise to the mono-, bi-, or poly-exponential concentration–time profiles observed in practice.
4 3. Rate of Tissue Distribution
4.1 3.1 Determinants of Tissue Distribution
Three factors govern the rate at which a drug distributes into tissues:
- Tissue perfusion (blood flow to the organ)
- Binding to plasma proteins and tissue components
- Permeability of tissue membranes to the drug
4.2 3.2 Perfusion Rate-Limited Distribution
Perfusion rate-limited distribution applies when membrane permeability is not a barrier — membranes are leaky, and the drug is a small, highly lipophilic molecule. Under these conditions, the rate-limiting step is the delivery of drug to the tissue via blood flow, not its penetration across the membrane.
Tissue perfusion is defined as:
\[Q_T = \frac{Q}{V_T}\]
where \(Q\) is blood flow (mL/min) and \(V_T\) is tissue volume (mL). Perfusion therefore has units of mL/min per mL of tissue, or min⁻¹.
For reference, hepatic blood flow is approximately 1.5 L/min; lung perfusion is among the highest (~10 mL/min/mL of tissue), while adipose tissue perfusion is low (~0.03 mL/min/mL).
Mathematical treatment
Consider an organ with blood flow \(Q\), arterial drug concentration \(C_a\), and venous drug concentration \(C_v\). Under the assumption that tissue membranes present no barrier, extracellular and intracellular spaces are treated as a single well-mixed compartment.
The amount of drug in the tissue is:
\[A_T = V_T \cdot C_T = V_T \cdot K_P \cdot C_v\]
where \(K_P\) is the equilibrium tissue-to-venous concentration ratio (the tissue partition coefficient):
\[K_P = \frac{C_T}{C_v}\]
\(K_P\) is defined using the venous concentration because venous blood is the only measurable surrogate for tissue drug levels in practice.
The fractional rate of exit from the tissue, \(k_T\), is:
\[k_T = \frac{\text{Rate of exit}}{A_T} = \frac{Q \cdot C_v}{V_T \cdot K_P \cdot C_v} = \frac{Q}{V_T \cdot K_P}\]
This can be rewritten as:
\[k_T = \frac{Q_T}{K_P} = \frac{\text{perfusion}}{K_P}\]
Analogously, for whole-body elimination:
\[k_{el} = \frac{CL}{V_d}\]
Thus \(k_T\) is the tissue-level equivalent of the whole-body elimination rate constant — it is a first-order rate constant with units of time⁻¹.
Distribution half-life
\[t_{1/2,\text{dist}} = \frac{0.693}{k_T} = \frac{0.693 \cdot K_P \cdot V_T}{Q}\]
This is the time for tissue concentration to reach 50% of its equilibrium value. Two key relationships follow:
- Higher blood flow \(Q\) → shorter distribution half-life → drug equilibrates more rapidly with the tissue.
- Lower \(K_P\) (low tissue affinity for the drug) → shorter distribution half-life → drug does not accumulate in the tissue.
4.3 3.3 Permeability Rate-Limited Distribution
When the flux across tissue membranes is much less than the rate of drug delivery by blood flow, distribution becomes permeability rate-limited. This is predominantly relevant for drugs distributed into the CNS, where the blood–brain barrier imposes significant permeability restrictions. Lipophilicity and the Henderson–Hasselbalch equation are the primary governing factors in these situations.
5 4. Extent of Distribution and Apparent Volume of Distribution
5.1 4.1 Definition of Apparent Volume of Distribution
At equilibrium, the apparent volume of distribution \(V_d\) is defined as:
\[V_d = \frac{A}{C}\]
where \(A\) is the total amount of drug in the body and \(C\) is the plasma (or blood) drug concentration. \(V_d\) is a proportionality constant relating the amount of drug in the body to the measured plasma concentration — it does not correspond to any real anatomical volume. Observed values range from approximately 3 L to 40,000 L in humans, which underscores its apparent nature.
The fraction of total body drug that resides within plasma is:
\[f_{\text{plasma}} = \frac{V_P}{V_d}\]
Since plasma volume is approximately 3 L, and \(V_d\) can be orders of magnitude larger, only a small fraction of total body drug is typically present in plasma at any given time.
5.2 4.2 Derivation of Total Apparent Volume
Partitioning the drug between plasma and a single tissue compartment:
| Compartment | Amount | Volume |
|---|---|---|
| Plasma | \(V_P \cdot C\) | \(V_P\) |
| Tissue | \(V_T \cdot K_P \cdot C\) | \(V_T\) |
Total amount in body: \(A = V_P \cdot C + V_T \cdot K_P \cdot C\)
Dividing both sides by \(C\):
\[V_d = V_P + V_T \cdot K_P\]
For the whole body (summing over all tissue compartments \(i\)):
\[V_d = V_P + \sum_i V_{T,i} \cdot K_{P,i}\]
This equation makes the “apparent” nature of \(V_d\) explicit: \(K_P\) scales the tissue volume contribution. In obese individuals, for example, adipose tissue volume is increased and adipose \(K_P\) values for lipophilic drugs can be substantial (e.g., 5×), causing \(V_d\) to increase markedly. Consequently, a larger dose may be required to maintain a target plasma concentration, because more drug is sequestered in fat.
6 5. Plasma Protein Binding
6.1 5.1 Key Binding Proteins
- Acidic drugs bind primarily to albumin
- Basic drugs bind primarily to alpha-1 acid glycoprotein (AAG)
- Some drugs bind to multiple plasma proteins
6.2 5.2 Fraction Unbound in Plasma (\(f_u\))
Three relevant concentration metrics:
- \(C_{BL}\): concentration in whole blood
- \(C\): total concentration in plasma (bound + unbound)
- \(C_u\): concentration of unbound drug in plasma water
The degree of plasma protein binding is defined as \(C_{\text{bound}} / C\), but because bound drug is difficult to measure directly, the fraction unbound is the more practically useful quantity:
\[f_u = \frac{C_u}{C}\]
It follows that \(f_u = 1 - f_{\text{bound}}\).
Only unbound drug is available to interact with target tissues and receptors; only unbound drug can cross membranes and exert pharmacological effect.
6.3 5.3 Association Constant and Binding Equilibrium
For a drug \(D\) binding to a protein \(P\) at a single binding site:
\[D + P \rightleftharpoons DP\]
The association constant \(K_a\) is:
\[K_a = \frac{[DP]}{[D][P]}\]
where \([DP]\) is the concentration of drug–protein complex (bound drug), \([D]\) is the unbound drug concentration, and \([P]\) is the concentration of unoccupied protein.
Since the fraction of unoccupied protein relates to total protein by \(P/P_T = f_u^P\) (the fraction of unoccupied binding sites), and given that \(C_u = f_u \cdot C\), the fundamental equation for fraction unbound is:
\[f_u = \frac{1}{1 + K_a \cdot f_u^P \cdot P_T}\]
This equation has far-reaching implications:
- If total protein \(P_T\) increases (e.g., during infection, when AAG is elevated as an acute-phase reactant), the denominator increases and \(f_u\) decreases — less free drug is available.
- At high drug concentrations, binding sites can become saturated, causing \(f_u\) to increase non-linearly with dose.
- These non-linearities in protein binding are a common source of unexpected pharmacokinetic behavior.
7 6. Tissue Binding and the Volume of Distribution
7.1 6.1 Role of Tissue Binding
Since tissue drug concentrations cannot be measured directly in clinical settings, tissue binding is inferred from plasma measurements. Starting from the total amount of drug in the body:
\[V_d \cdot C = V_P \cdot C + V_{T,W} \cdot C_{T,W}\]
Defining the tissue fraction unbound \(f_{uT} = C_{uT}/C_{T,W}\) and using \(f_u = C_u/C\), this leads to:
\[V_d = V_P + V_{T,W} \cdot \frac{f_u}{f_{uT}}\]
This is a key equation: the ratio \(f_u / f_{uT}\) drives the tissue contribution to volume of distribution. If a drug binds more avidly to tissue proteins than to plasma proteins (\(f_{uT} < f_u\)), the apparent volume of distribution will be large.
7.2 6.2 Active Transport and Volume of Distribution
When active transporters are present, at equilibrium:
\[P_{\text{up}} \cdot SA \cdot C_u = P_{\text{ef}} \cdot SA \cdot C_{uT}\]
where \(P_{\text{up}}\) is the uptake permeability and \(P_{\text{ef}}\) is the efflux permeability. This gives:
\[\frac{C_{uT}}{C_u} = \frac{P_{\text{up}}}{P_{\text{ef}}}\]
Incorporating this into the volume equation:
\[V_d = V_P + \sum_i V_{T,i} \cdot \frac{f_u}{f_{uT,i}} \cdot \frac{P_{\text{up},i}}{P_{\text{ef},i}}\]
This general equation encompasses both passive and active distribution. For passive-only distribution, \(P_{\text{up}} / P_{\text{ef}} = 1\), and the equation reduces to the earlier form. This expression is the foundation of PBPK modeling, where software systems use literature-curated values for \(V_{T,i}\), \(K_{P,i}\), and transporter parameters, derived from drug properties such as lipophilicity, molecular weight, \(\text{p}K_a\), and molecular size.
8 7. Summary: Connecting the Concepts
The four lectures on volume of distribution addressed two complementary questions:
- Rate of distribution — by what processes does drug enter tissues, and how quickly?
- Extent of distribution — how much drug distributes into tissues, and what determines this?
The rate is governed by perfusion (blood flow), permeability (size, lipophilicity, ionization), and transporter activity. The extent is captured by \(V_d\), which reflects plasma protein binding, tissue binding, and transporter-mediated accumulation.
The mathematical framework developed here — Fick’s first law, the Henderson–Hasselbalch equation, perfusion rate equations, the protein binding equilibrium, and the volume partitioning equation — recurs throughout pharmacokinetics. These same equations reappear in the context of clearance, and later in population PK and PBPK modeling. Mastery of these fundamentals is essential before advancing to those topics.
8.1 Key Equations Reference
| Concept | Equation |
|---|---|
| Net diffusion rate | \(\text{Rate} = P \cdot SA \cdot (C_1 - C_2)\) |
| Henderson–Hasselbalch (acid) | \(\text{pH} = \text{p}K_a + \log_{10}\frac{[A^-]}{[HA]}\) |
| Tissue amount | \(A_T = V_T \cdot K_P \cdot C_v\) |
| Fractional tissue exit rate | \(k_T = Q / (V_T \cdot K_P)\) |
| Distribution half-life | \(t_{1/2,\text{dist}} = 0.693 \cdot K_P \cdot V_T / Q\) |
| Apparent volume (one tissue) | \(V_d = V_P + V_T \cdot K_P\) |
| Fraction unbound | \(f_u = 1 / (1 + K_a \cdot f_u^P \cdot P_T)\) |
| Volume with tissue binding | \(V_d = V_P + V_{T,W} \cdot (f_u / f_{uT})\) |
| Volume with active transport | \(V_d = V_P + \sum_i V_{T,i} \cdot (f_u / f_{uT,i}) \cdot (P_{\text{up},i} / P_{\text{ef},i})\) |