Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › Free Versus Micelle-Bound Drug in Dissolution Models

Dissolution and Solubility Modeling · Tutorial 261 of 799

Free Versus Micelle-Bound Drug in Dissolution Models

Use free and micelle-bound concentrations to interpret dissolution, apparent solubility, and precipitation without confusing total dissolved drug with the pool that supplies the interfacial driving force.

Advanced 11 min read

What You'll Learn

  • Define free, micelle-bound, and total dissolved drug concentrations.
  • Calculate free concentration from linear or saturable micellar binding.
  • Relate free concentration to dissolution driving force and precipitation supersaturation.
  • Quantify errors from substituting total concentration for free concentration.
  • Recognize when equilibrium partitioning is inadequate and separate-pool dynamics are needed.

Why the Dissolved Pool Must Be Partitioned

A dissolution sample can contain drug molecules freely dispersed in the aqueous phase and drug associated with bile-salt micelles. Both contribute to a measurement of total dissolved drug, but they need not have the same availability at a solid surface or the same thermodynamic activity relevant to precipitation. Treating the measured total as though every molecule were free can therefore distort both dissolution and precipitation predictions.

The key distinction is between how much drug is dissolved in total and how much is present in the free pool. Micellar binding can increase apparent solubility by storing drug in the bound pool. If binding and release are sufficiently rapid, that pool can replenish free drug as it is removed. It does not follow that bound drug directly contributes to the interfacial concentration gradient or should be counted as free when evaluating supersaturation.

Definition: Let \(C_f\) be the free dissolved-drug concentration, \(C_b\) the micelle-bound dissolved-drug concentration, and \(C_T\) the total dissolved-drug concentration, all expressed on the same drug-equivalent concentration basis. Then \(C_T=C_f+C_b\). The free fraction is \(f_f=C_f/C_T\), and the bound fraction is \(f_b=C_b/C_T\), so \(f_f+f_b=1\) when \(C_T>0\).

These concentrations refer to dissolved drug; undissolved solid is not included in \(C_T\). This tutorial first treats the bulk medium as well mixed and assumes that micelle-bound drug does not directly cross the aqueous boundary layer to the solid surface. That is a modeling assumption, not a universal property of all micellar systems. The introductory material on dynamic micellar solubilization and drug exchange motivates why exchange timescale matters: an equilibrium partition is appropriate only when exchange is fast relative to the process being represented.

Equilibrium Partitioning and Its Assumptions

For a dilute system in which the micellar binding capacity is not appreciably saturated, a useful approximation is a linear partition relation \(C_b=\beta C_f\), where \(\beta\) is a nonnegative, dimensionless binding ratio for the specified medium. It summarizes the effects of micelle abundance and drug–micelle association under those conditions; it is not necessarily constant when medium composition changes.

$$ C_T=(1+\beta)C_f,\qquad C_f=\frac{C_T}{1+\beta},\qquad f_f=\frac{1}{1+\beta} $$

The linear relation assumes that the relevant micellar environment is effectively constant, that the bound pool responds rapidly enough to remain in equilibrium with the free pool, and that one effective bound pool is an adequate description. It should not be extrapolated beyond the range in which the relation holds. If binding sites or micellar capacity are finite, a saturable relation is more appropriate. One common equilibrium form is

$$ C_b=\frac{B_{\max}C_f}{K_d+C_f}, \qquad C_T=C_f+\frac{B_{\max}C_f}{K_d+C_f} $$

Here \(B_{\max}\) is the maximum concentration of drug that can be represented in the bound pool, and \(K_d\) is the free concentration at which the bound concentration is half of \(B_{\max}\). Both have concentration units. The relation is monotone in \(C_f\), so a nonnegative total concentration corresponds to one nonnegative free concentration. The capacity and affinity parameters must apply to the medium and conditions being modeled.

Worked Example: Recovering the Free Pool From a Linear Partition

Worked Example: Recovering the Free Pool From a Linear Partition

Suppose a hypothetical dissolution medium contains \(C_T=10.0\) mg/L total dissolved drug and has a linear binding ratio \(\beta=4.00\). Assume rapid equilibrium and that the linear relation is valid at this concentration. Calculate the free and bound concentrations and the free fraction.

From \(C_T=(1+\beta)C_f\),

$$ C_f=\frac{10.0\ \mathrm{mg/L}}{1+4.00} =2.00\ \mathrm{mg/L} $$

The bound concentration is the remainder of the total:

$$ C_b=C_T-C_f=10.0-2.00=8.00\ \mathrm{mg/L} $$

Thus \(f_f=2.00/10.0=0.200\), or 20.0%, and \(f_b=0.800\), or 80.0%. The total measurement is 10.0 mg/L, but under the stated partition model the free concentration is only 2.00 mg/L.

Which Concentration Drives Dissolution?

For a simple mass-transfer representation, let \(S_f\) be the free-drug concentration at equilibrium with the solid at its surface, and let \(G\) be an effective mass-transfer conductance. If the bulk free concentration is \(C_f\), the net dissolution rate can be written

$$ r_{\mathrm{diss}}=G(S_f-C_f) $$

With \(G\) in L/min and concentrations in mg/L, \(r_{\mathrm{diss}}\) is in mg/min. A positive value represents net dissolution; a negative value represents net transfer toward the solid under this simplified reversible-flux description. The equation assumes the same free-drug basis for the surface equilibrium concentration and bulk concentration, and it omits complications such as changing particle area or a spatially varying boundary layer.

The corresponding driving difference is \(S_f-C_f\), not \(S_f-C_T\), under the assumption that bound drug does not directly contribute to transport across the boundary layer. Bound drug can still matter indirectly: if exchange is fast, release from micelles can replenish free drug and sustain dissolution. That replenishment is captured through the relationship between \(C_T\) and \(C_f\), rather than by inserting \(C_T\) as the free concentration.

Worked Example: How Total Concentration Can Reverse a Dissolution Prediction

Consider hypothetical values \(S_f=3.00\) mg/L, \(C_T=10.0\) mg/L, and \(\beta=4.00\). Assume the linear equilibrium partition is valid, exchange is fast, and bound drug does not directly cross the boundary layer. Let \(G=0.100\) L/min. Determine the predicted net rate using the free pool, then compare it with the total-as-free calculation.

The first worked calculation gives \(C_f=2.00\) mg/L. Therefore,

$$ r_{\mathrm{diss}}=0.100\ \mathrm{L/min} \left(3.00-2.00\right)\ \mathrm{mg/L} =0.100\ \mathrm{mg/min} $$

This model predicts net dissolution. If total concentration were incorrectly substituted for free concentration, the calculation would instead be

$$ r_{\mathrm{wrong}}=0.100\ \mathrm{L/min} \left(3.00-10.0\right)\ \mathrm{mg/L} =-0.700\ \mathrm{mg/min} $$

The substitution changes the predicted direction of net transfer. That conclusion is not a small numerical correction: it is an artifact of using a pool that the assumed interfacial flux law does not treat as freely available.

Free Concentration and Precipitation

Supersaturation relevant to precipitation should likewise be evaluated using the concentration or activity of the precipitating free species, not automatically the total dissolved concentration. For the simple concentration-based model used here, let \(S_{\mathrm{eq}}\) be the equilibrium free-drug solubility for the specified solid form and medium. Define the free supersaturation ratio as \(C_f/S_{\mathrm{eq}}\). A ratio greater than 1 indicates that the free concentration exceeds this equilibrium reference; it does not, by itself, determine when precipitation begins.

For illustration, a precipitation-rate expression might use

$$ r_{\mathrm{prec}} =k_{\mathrm{prec}} \left[ \max\left(\frac{C_f}{S_{\mathrm{eq}}}-1,0\right) \right]^n $$

Here \(k_{\mathrm{prec}}\) has units of concentration per time, \(n\) is a specified exponent, and \(r_{\mathrm{prec}}\) is a concentration removal rate. This is an illustrative rate law, not a universal precipitation mechanism. Nucleation barriers, induction periods, and growth kinetics can alter the time course even when the free concentration is supersaturated. The important point for this tutorial is that replacing \(C_f\) with \(C_T\) can create a false supersaturation signal or exaggerate a real one.

Worked Example: Avoiding a False Precipitation Signal

In a hypothetical medium, let \(C_f=0.800\) mg/L, \(C_b=3.20\) mg/L, and \(S_{\mathrm{eq}}=1.00\) mg/L. Thus \(C_T=4.00\) mg/L. Under the concentration-based definition, calculate free and total supersaturation ratios and evaluate the illustrative rate law with \(n=2\) and \(k_{\mathrm{prec}}=0.0200\) mg/L/min.

The free supersaturation ratio is

$$ \frac{C_f}{S_{\mathrm{eq}}} =\frac{0.800}{1.00}=0.800 $$

Because this is below 1, the positive-part term is zero and the specified model predicts \(r_{\mathrm{prec}}=0\) mg/L/min. By contrast, the total-based ratio would be \(4.00/1.00=4.00\). Incorrectly using that value gives

$$ r_{\mathrm{wrong}} =0.0200\left[\max(4.00-1,0)\right]^2 =0.0200(3.00)^2 =0.180\ \mathrm{mg/L/min} $$

The difference is caused by the concentration substitution, not by a change in the stated precipitation model. Under the assumptions and rate law used here, the free pool is below the equilibrium reference even though total dissolved drug is above it.

Why Binding Can Increase Apparent Solubility

Micellar binding can raise the total concentration that coexists with a solid without raising the equilibrium free concentration. At equilibrium with the solid, set \(C_f=S_f\) in the partition relation. For linear binding, the apparent total solubility is \(S_{\mathrm{app}}=(1+\beta)S_f\). For saturable binding, it is \(S_{\mathrm{app}}=S_f+B_{\max}S_f/(K_d+S_f)\). These are total concentrations at the specified free equilibrium concentration; neither expression means the free solubility itself has increased by the same factor.

Worked Example: Solving a Saturable Partition

Suppose a hypothetical medium is described by \(B_{\max}=12.0\) mg/L and \(K_d=2.00\) mg/L. The measured total dissolved concentration is \(C_T=8.00\) mg/L. Assuming the saturable equilibrium relation applies, find \(C_f\), \(C_b\), and the free fraction.

Substitute the parameters into the total-concentration equation and let \(x=C_f\):

$$ 8.00=x+\frac{12.0x}{2.00+x} $$

Multiply both sides by \(2.00+x\):

$$ 8.00(2.00+x)=x(2.00+x)+12.0x $$

Expanding and rearranging gives \(x^2+6.00x-16.0=0\). The roots are \(x=2.00\) mg/L and \(x=-8.00\) mg/L. A concentration cannot be negative, so \(C_f=2.00\) mg/L. Then

$$ C_b=\frac{12.0(2.00)}{2.00+2.00} =6.00\ \mathrm{mg/L}, \qquad f_f=\frac{2.00}{8.00}=0.250 $$

The pools sum to \(2.00+6.00=8.00\) mg/L, as required. At a free equilibrium concentration of 1.00 mg/L, this same medium would have \(C_b=12.0(1.00)/(2.00+1.00)=4.00\) mg/L and apparent total solubility \(1.00+4.00=5.00\) mg/L. The total solubility is higher than the free equilibrium concentration because the bound pool stores additional dissolved drug.

Deriving the Effect of Binding on Total-Concentration Changes

When partitioning is at equilibrium, changes in total concentration do not translate one-for-one into changes in free concentration. For the saturable model, differentiate the total concentration with respect to \(C_f\):

$$ \frac{dC_T}{dC_f} =1+\frac{B_{\max}K_d}{(K_d+C_f)^2} $$

To see the result, write \(C_T=C_f+B_{\max}C_f/(K_d+C_f)\). The derivative of the bound term follows from the quotient rule:

$$ \frac{d}{dC_f}\left(\frac{B_{\max}C_f}{K_d+C_f}\right) = B_{\max}\frac{(K_d+C_f)-C_f}{(K_d+C_f)^2} = \frac{B_{\max}K_d}{(K_d+C_f)^2} $$

The derivative is at least 1 for nonnegative parameters. Thus a change in free concentration can require a larger change in total concentration because the bound pool changes at the same time. This relationship is a local equilibrium result; it does not by itself specify the time required for binding and unbinding. When exchange is slow relative to dissolution or precipitation, separate dynamic equations for \(C_f\) and \(C_b\) are needed rather than an instantaneous partition formula.

Diagnostics, Alternatives, and Limitations

A practical diagnostic is to check whether the model’s reported free and bound concentrations sum to the reported total, and whether each pool is nonnegative. Also check that the partition parameters and free-solubility reference apply to the same medium, temperature, pH, and drug species. A mismatch can appear as an implausible free fraction or as a predicted total concentration that cannot be reconciled with the assumed capacity.

Compare predictions based on \(C_f\) with those from a total-only calculation. A large difference in dissolution direction, supersaturation ratio, or precipitation rate is a warning that the pool distinction materially affects the prediction. It is not a reason to select whichever prediction better matches expectations; it is a reason to verify the partition assumption and the concentration basis of the measured data.

Modeling approachWhat it representsMain caution
Explicit free and bound pools with equilibrium partitioningTracks \(C_f\) and \(C_b\) while using a stated partition relationRequires equilibrium to be fast enough and the partition relation to apply
Total-only apparent-solubility representationSummarizes the total dissolved amount under specified conditionsCannot, by itself, identify the free concentration driving an interface process
Separate dynamic free and bound poolsAllows exchange and the two pool sizes to change over timeRequires exchange-rate information; more parameters may be difficult to identify

An explicit free–bound equilibrium model is a useful middle ground when partitioning is defensible but the detailed exchange timescale is not the target. A total-only representation may be adequate for reporting measured apparent solubility, but it is insufficient for a process whose driving force depends on free drug unless a validated conversion is supplied. A kinetic two-pool model is preferable when exchange is not rapid, but it should not be added without data or a justified timescale.

Finally, concentration is a practical proxy for chemical activity. Non-ideal solutions, changes in ionization, multiple drug species, or changes in micellar composition may require a more careful thermodynamic treatment. The simple formulas here also do not establish that micelle-bound drug can never reach a surface or contribute to transport: they state the assumption used in the illustrated flux and precipitation calculations. If evidence supports direct micellar transport or a different interfacial mechanism, that mechanism needs to be represented explicitly rather than hidden by substituting total for free concentration.

Common Mistake

  • Wrong reasoning: “The assay reports 10 mg/L dissolved drug, so the interfacial model should use 10 mg/L as the bulk concentration.”
  • Correction: The assay’s total may include micelle-bound drug. First determine \(C_f\) from an appropriate partition model or measurement; use that free concentration wherever the assumed driving force or supersaturation definition is based on free drug.
  • Another pitfall: “If apparent solubility rises, free solubility must rise by the same amount.” Apparent solubility can rise because the bound pool increases while the equilibrium free concentration remains fixed.
  • Check before interpreting: Confirm that total equals free plus bound, that the model’s equilibrium assumption is plausible, and that the free solubility and partition parameters refer to compatible conditions.
Key takeaway: Total dissolved concentration measures the sum of free and micelle-bound drug. Under a model in which free drug supplies the interfacial driving force and sets the precipitation supersaturation, use \(C_f\) for those calculations and represent the bound pool through an explicit, justified partition or exchange model.

Check Your Understanding

Answer each question using the assumptions stated in it. Solutions are included below the questions.

  1. Under linear binding with \(\beta=3.00\), what are \(C_f\) and \(C_b\) if \(C_T=12.0\) mg/L?
  2. If \(S_f=2.00\) mg/L and the values in Question 1 apply, what is the sign and magnitude of \(S_f-C_f\)?
  3. For \(C_f=1.50\) mg/L, \(C_b=4.50\) mg/L, and \(S_{\mathrm{eq}}=2.00\) mg/L, calculate the free supersaturation ratio. Would substituting total concentration indicate the same status relative to equilibrium?
  4. Give one condition under which the equilibrium partition relation should not be used to collapse the free and bound pools into a single instantaneous relationship.

Solutions. (1) \(C_f=12.0/(1+3.00)=3.00\) mg/L, and \(C_b=12.0-3.00=9.00\) mg/L. (2) \(S_f-C_f=2.00-3.00=-1.00\) mg/L, indicating a negative dissolution driving difference under the stated flux convention. (3) \(C_T=1.50+4.50=6.00\) mg/L, so the free ratio is \(1.50/2.00=0.750\), below equilibrium. The total-based ratio is \(6.00/2.00=3.00\), which would incorrectly indicate a concentration above equilibrium under a free-concentration criterion. (4) If micelle–drug exchange is not fast relative to dissolution or precipitation, the instantaneous equilibrium relation may be inadequate; separate dynamic pools are then needed.