Two Different Questions: Where and How Fast?
A tissue model may specify that a drug has a preferred equilibrium concentration in one compartment relative to another. That relationship answers where the concentrations tend in the long run; it does not say how quickly the drug crosses the membrane. A finite membrane permeability answers the second question. Keeping these roles separate is essential when interpreting tissue concentration profiles: the same eventual partitioning can accompany very different transient behavior.
Consider an extracellular compartment with volume \(V_e\) and concentration \(C_e\), and an intracellular compartment with volume \(V_i\) and concentration \(C_i\). The compartments exchange drug across a membrane. For the derivation and examples below, assume fixed positive volumes, a closed system with no input or elimination, and a constant, positive equilibrium partition ratio \(K\). Assume also that exchange is linear in the concentration difference expressed relative to that ratio. These assumptions isolate membrane exchange; they are not a claim that every tissue has constant volumes or no other processes.
For this two-compartment model, define the net flux from extracellular to intracellular space as \(J=PS(C_e-C_i/K)\). A positive \(J\) means drug is moving into the intracellular compartment. The factor \(1/K\) makes zero flux correspond to the specified equilibrium relationship \(C_i=K C_e\). With this sign convention, a concentration difference from equilibrium can drive flux in either direction.
The model assumes that the concentration-based flux law applies throughout the interval being considered, and that \(PS\) and \(K\) do not change during that interval. It also treats each compartment as well mixed, so a single concentration represents each compartment. If binding, active transport, changing membrane properties, or spatial gradients are important, this simple model may not adequately describe the observed profiles.
Mass Conservation and the Equilibrium Concentrations
The amount of drug in the two compartments is \(A=V_e C_e+V_i C_i\). In this closed system, the exchange flux leaves one compartment at exactly the rate it enters the other:
Thus, \(A\) remains constant. At equilibrium, \(J=0\), so \(C_i=K C_e\). Combining this relationship with the conserved amount gives \(A=(V_e+K V_i)C_e\). The equilibrium concentrations are therefore:
These expressions show what sets the long-time concentrations under the stated closed-system assumptions: the total amount, the compartment volumes, and the equilibrium ratio. The exchange capacity \(PS\) does not appear. Changing \(PS\) alone cannot change these equilibrium concentrations; it changes how long the model takes to approach them.
Worked Example: Finding Equilibrium From Total Amount
Suppose a hypothetical closed tissue model has \(V_e=0.50\) L, \(V_i=1.50\) L, total amount \(A=4.00\) mg, and equilibrium partition ratio \(K=2.00\). Assume the volumes and partition ratio are constant and there is no input or elimination. Find the equilibrium concentrations.
The denominator is \(V_e+K V_i=0.50+(2.00)(1.50)=3.50\) L. Therefore:
The amounts at these rounded concentrations are approximately \((0.50)(1.143)=0.572\) mg and \((1.50)(2.286)=3.429\) mg, totaling 4.001 mg because of rounding. Using the unrounded concentrations, the amounts sum to exactly 4.00 mg. The higher intracellular concentration does not imply that more drug is in the intracellular compartment solely because of concentration: amount depends on both concentration and volume.
Deriving the Approach to Equilibrium
To describe the transient, define \(q=C_e-C_i/K\), the concentration difference that drives the modeled flux. At equilibrium \(q=0\). Differentiating \(q\) and substituting the two compartment equations gives a single differential equation:
For constant positive \(PS\), \(V_e\), \(V_i\), and \(K\), this equation has the exponential solution below. The coefficient \(\lambda\) has units of inverse time.
Because \(A\) is conserved, the deviations of \(C_e\) and \(C_i\) from their equilibrium concentrations must offset one another in amount. Combining that constraint with \(q=C_e-C_i/K\) yields the individual concentrations:
For example, if \(q(0)>0\), then \(C_e\) begins above its equilibrium value and \(C_i\) below its equilibrium value. As \(q(t)\) decays toward zero, the extracellular concentration falls and the intracellular concentration rises. If \(q(0)<0\), the directions reverse. This is a relaxation toward the same equilibrium partition relationship, not a changing value of \(K\).
The exponential form follows by separating variables when \(q\ne 0\): \(dq/q=-\lambda\,dt\). Integrating from 0 to \(t\) gives \(\ln|q(t)|-\ln|q(0)|=-\lambda t\), and exponentiating gives \(q(t)=q(0)e^{-\lambda t}\). If \(q(0)=0\), the system is already at equilibrium, \(J=0\), and the same solution gives \(q(t)=0\) throughout.
Fast Exchange: Equilibrium Is Approached Quickly
When \(PS\) is large relative to the compartment volumes and observation period, \(\lambda\) is large and the exponential term decays rapidly. After a short transient, concentrations are close to their equilibrium values. An instantaneous-equilibrium approximation sets \(C_i=K C_e\) at all modeled times; it can be useful when exchange is sufficiently fast for the question being asked. It is not the same as saying that \(K\) itself is large.
Worked Example: A Fast Membrane Exchange
Use a hypothetical model with \(V_e=0.25\) L, \(V_i=0.75\) L, \(K=2.00\), \(A=2.00\) mg, and \(PS=10.0\) L/h. At time zero, let \(C_e=8.00\) mg/L and \(C_i=0\) mg/L; these initial concentrations contain \((0.25)(8.00)+(0.75)(0)=2.00\) mg. Assume no other processes act over the 0.100 h interval.
The equilibrium concentrations are \(C_e^*=2.00/[0.25+(2.00)(0.75)]=1.143\) mg/L and \(C_i^*=2C_e^*=2.286\) mg/L. Initially \(q(0)=8.00-0/2=8.00\) mg/L. The model’s decay coefficient and concentration difference at 0.100 h are:
The concentration coefficients are \(K V_i/(V_e+K V_i)=1.50/1.75\) and \(K V_e/(V_e+K V_i)=0.500/1.75\). Applying them gives:
The two concentrations are already close to the equilibrium values, but are not exactly at equilibrium: the remaining \(q\) is about 0.07523 mg/L. The total amount remains 2.00 mg, up to rounding. Treating the compartments as instantly equilibrated would be an approximation at this observation time, not a different equilibrium partition ratio.
Slow Exchange: The Long-Time Ratio May Be Far Away
If \(PS\) is very small, little exchange occurs over a fixed observation interval. The equilibrium ratio still defines the eventual state for every positive \(PS\) in this closed model, but it may be a poor description of concentrations measured before that state is approached. At exactly \(PS=0\), there is no exchange: each compartment retains its initial concentration, and the system does not move toward the partition relationship unless it already satisfies it.
Worked Example: A Slow Membrane Exchange
Use the same hypothetical volumes, \(K\), total amount, and initial concentrations as in the fast-exchange example, but set \(PS=0.00100\) L/h. Keep the no-input, no-elimination assumption and evaluate the concentrations after 0.100 h. The equilibrium values remain 1.143 mg/L extracellularly and 2.286 mg/L intracellularly because \(PS\) does not enter their formulas.
Now the decay coefficient is \(\lambda=(0.00100)(4.00+0.666667)=0.004667\ \mathrm{h^{-1}}\). Thus:
Using the same concentration coefficients as before:
After this short interval, the extracellular concentration remains close to its initial value and the intracellular concentration remains close to zero, despite the equilibrium ratio \(K=2.00\). The small positive intracellular concentration is consistent with a small but nonzero exchange flux. It would be incorrect to conclude from this early concentration profile that the equilibrium ratio is close to zero.
Interpreting the Limits and Checking the Model
The limiting cases clarify the distinction between equilibrium and rate. For any positive finite \(PS\), the model gives \(q(t)\to0\) as \(t\to\infty\), so the concentrations tend toward \(C_i=K C_e\). As \(PS\to\infty\), the decay is arbitrarily rapid; in the idealized limit, an initial mismatch disappears immediately after time zero. As \(PS\to0\), the rate of change tends to zero. At \(PS=0\), an initially unequal partition remains unequal indefinitely. Thus, a zero-permeability model and a very-low-permeability model can look similar over a short interval but differ in their long-time behavior.
A useful diagnostic is to check whether a proposed trajectory respects the assumptions and mathematical constraints. In the closed model, total amount \(V_e C_e+V_i C_i\) must remain constant. The sign of \(J\) must agree with \(q\): if \(C_e>C_i/K\), flux is inward, so \(C_e\) decreases and \(C_i\) increases. For nonnegative initial concentrations and positive parameters, a computed trajectory should not produce negative concentrations. A violation can indicate a sign error, inconsistent volumes, a unit mismatch, or use of the formula outside its assumptions.
Units provide another check. Since \(J\) is amount per time and \(C_e-C_i/K\) is amount per volume, \(PS\) must have units of volume per time. Then \(PS(1/V_e+1/(K V_i))\) has units of inverse time, as required for the exponent \(-\lambda t\) to be dimensionless. The partition ratio \(K\) is dimensionless when both concentrations use the same units.
| Model choice | What it represents | Important limitation |
|---|---|---|
| Finite-rate exchange | A concentration difference drives a permeability-dependent flux over time. | Requires an exchange law and parameters such as \(PS\); this simple form assumes constant values and well-mixed compartments. |
| Instantaneous equilibrium | Imposes \(C_i=K C_e\) at each modeled time. | Cannot represent a measurable delay in crossing; it is an approximation when exchange is fast enough for the intended analysis. |
| No exchange | Sets \(PS=0\), so neither compartment changes through membrane flux. | Cannot approach the equilibrium ratio from an initially unequal state. |
Perfusion-limited tissue models use a different simplification: they treat tissue and blood as equilibrating sufficiently rapidly for the relevant model purpose. A permeability-limited model retains a finite membrane flux when that assumption is not suitable. Neither approach makes the partition ratio a rate parameter. The choice between them depends on whether finite exchange delay matters for the question and data at hand.
Common Mistake
A frequent error is to treat \(K\) as a measure of how quickly drug enters the cell. For example, reasoning that a larger \(K\) must mean faster intracellular accumulation confuses a concentration ratio at zero net flux with a kinetic property. In this model, changing \(K\) changes the equilibrium condition and also appears in the mathematical expression for \(\lambda\), but the rate still depends on \(PS\), volumes, and the specified exchange law. A large partition ratio does not by itself establish a large permeability or rapid equilibration.
A related error is to impose \(C_i=K C_e\) while also claiming to model a substantial transient intracellular lag. The algebraic equality removes that lag by assumption. To represent delayed equilibration, retain the differential equations and finite flux; to use instantaneous equilibrium, state that the approximation is intended and check that the unresolved transient is unimportant for the application.
Check Your Understanding
Answer each question, then compare with the solutions below.
- In a closed two-compartment model, what determines the equilibrium concentrations, and which parameter controls the flux rate?
- If \(K=3\), \(V_e=1.0\) L, \(V_i=1.0\) L, and \(A=8.0\) mg, calculate \(C_e^*\) and \(C_i^*\).
- For positive finite \(PS\), what happens to \(q=C_e-C_i/K\) as time increases? What changes if \(PS=0\) and \(q(0)\ne0\)?
- A model predicts that \(C_e\) and \(C_i\) change through exchange, but \(V_eC_e+V_iC_i\) increases despite no input. Give a likely interpretation of this result.
Solutions. (1) The equilibrium concentrations depend on the conserved amount \(A\), the volumes, and \(K\); the finite flux rate depends on \(PS\) as well as the volumes and \(K\). (2) The denominator is \(1.0+(3)(1.0)=4.0\) L, so \(C_e^*=8.0/4.0=2.0\) mg/L and \(C_i^*=3(2.0)=6.0\) mg/L. These concentrations imply amounts of 2.0 mg and 6.0 mg, summing to 8.0 mg. (3) Since \(q(t)=q(0)e^{-\lambda t}\) with \(\lambda>0\), \(q\) tends to zero asymptotically. If \(PS=0\), then \(\lambda=0\) and \(q(t)=q(0)\): the initial mismatch persists. (4) Exchange alone must conserve the total amount. An increasing total amount suggests a sign or implementation error, an unmodeled input, or an incorrect accounting of compartment volumes or amounts.