Comparing Tissues Requires More Than Comparing Blood Flow
Perfusion-limited tissues can respond very differently to the same blood concentration history. A tissue’s equilibrium tissue-to-blood partition coefficient describes how much drug it holds relative to blood at equilibrium. Its blood flow and volume, together with that partitioning, determine how quickly it approaches that equilibrium. Comparing tissues therefore requires keeping two questions separate: where does each tissue tend to equilibrate, and how fast does it get there?
The calculations here use hypothetical tissue parameters and a common blood input. They illustrate the model’s internal comparisons, not measured values for particular drugs or organs. The starting assumption is the perfusion-limited tissue model: tissue is well mixed, exchange with blood is fast enough that the venous concentration leaving the tissue is \(C_t/K_p\), and blood flow and tissue volume remain constant during the calculation. Local metabolism and other tissue losses are omitted.
With \(Q_i\) measured in L/h and \(V_i\) in L, the rate \(Q_i/V_i\) has units of inverse hours. The time constant has units of hours, while \(K_{p,i}\) and the equilibrium concentration ratio are dimensionless. A large partition coefficient raises the equilibrium tissue concentration for a given blood concentration, but it also increases the time constant when volume and flow are held fixed. A large flow tends to shorten the time constant, but flow alone does not determine it.
Why the Time Constant Has This Form
Rearranging the tissue equation shows that the rate at which a tissue approaches its equilibrium concentration is governed by \(Q_i/(V_iK_{p,i})\). If the blood concentration is held fixed at \(C_b\), define the deviation from equilibrium as \(D_i=C_{t,i}-K_{p,i}C_b\). Because \(C_b\) is constant, differentiating \(D_i\) gives
The solution is an exponential decline in the deviation: \(D_i(t)=D_i(0)e^{-t/\tau_i}\), where \(\tau_i=V_iK_{p,i}/Q_i\). At \(t=(\ln 2)\tau_i\), the remaining deviation is \(D_i(0)e^{-\ln 2}=D_i(0)/2\); this justifies the half-time formula. The half-time refers to halving the distance from equilibrium, not necessarily halving the tissue concentration itself.
Worked Example: Partitioning and Half-Times for Three Tissues
Suppose three hypothetical tissues have the following parameters. Muscle has \(V=28.0\) L, \(Q=50.4\) L/h, and \(K_p=2.00\). Adipose tissue has \(V=18.0\) L, \(Q=10.8\) L/h, and \(K_p=10.0\). A high-flow tissue has \(V=1.50\) L, \(Q=90.0\) L/h, and \(K_p=1.50\). Calculate each time constant and half-time.
For muscle, \(\tau=VK_p/Q=(28.0\ \mathrm{L})(2.00)/(50.4\ \mathrm{L/h})=1.11111\ \mathrm{h}\). Its half-time is \((\ln 2)(1.11111\ \mathrm{h})=0.7702\ \mathrm{h}\).
For adipose tissue, \(\tau=(18.0\ \mathrm{L})(10.0)/(10.8\ \mathrm{L/h})=16.6667\ \mathrm{h}\). Its half-time is \((\ln 2)(16.6667\ \mathrm{h})=11.55\ \mathrm{h}\).
For the high-flow tissue, \(\tau=(1.50\ \mathrm{L})(1.50)/(90.0\ \mathrm{L/h})=0.02500\ \mathrm{h}\). Its half-time is \((\ln 2)(0.02500\ \mathrm{h})=0.01733\ \mathrm{h}\), or \(1.040\) min.
At equilibrium, their tissue-to-blood concentration ratios are \(2.00\), \(10.0\), and \(1.50\), respectively. The high-flow tissue equilibrates fastest for these specified values, while adipose tissue has the longest half-time and the largest equilibrium ratio. These results follow from the full combination of \(V\), \(Q\), and \(K_p\); they should not be attributed to flow alone.
Responses to a Shared Declining Blood Concentration
A direct comparison can use one blood concentration history for every tissue. Suppose blood starts at \(C_0\) and declines exponentially as \(C_b(t)=C_0e^{-kt}\), where \(k\) is a positive rate constant. If tissue concentration is initially zero, the tissue equation becomes a first-order linear differential equation with a time-varying input. Multiplying by its integrating factor and integrating from zero to \(t\) gives
This expression is a response to a prescribed input, not an equilibrium calculation. The terms \(e^{-kt}\) and \(e^{-t/\tau_i}\) describe, respectively, the declining blood input and the tissue’s own equilibration rate. If the two rates coincide, the displayed denominator is zero but the response itself has a finite limit: \(C_{t,i}(t)=(Q_i/V_i)C_0t e^{-kt}\) when \(k=1/\tau_i\). The formula assumes the same blood input reaches each tissue and that the model parameters remain constant.
Worked Example: Tissue Concentrations at Two Hours
Use the three tissues above and the shared hypothetical blood input \(C_b(t)=(2.00\ \mathrm{mg/L})e^{-(0.500\ \mathrm{h}^{-1})t}\). Calculate each tissue concentration at \(t=2.00\) h. Retain unrounded values in intermediate calculations.
For muscle, \(\tau=1.11111\) h and \(k\tau=0.555556\). Substitution in the response formula gives
For adipose tissue, \(\tau=16.6667\) h and \(k\tau=8.33333\). The denominator is negative, as is the difference between the two exponential terms; their product is positive:
For the high-flow tissue, \(\tau=0.02500\) h and \(k\tau=0.01250\):
At two hours the blood concentration is \(2.00e^{-1}=0.7358\ \mathrm{mg/L}\). The corresponding equilibrium tissue concentrations at that instant would be \(1.472\ \mathrm{mg/L}\) for muscle, \(7.358\ \mathrm{mg/L}\) for adipose tissue, and \(1.104\ \mathrm{mg/L}\) for the high-flow tissue. None of those equilibrium values is a substitute for the calculated transient concentrations. In particular, adipose tissue has a high equilibrium ratio but has not approached the value implied by the current blood concentration during this declining input.
Reading the Comparison Correctly
The values in the worked comparison illustrate several important distinctions. First, \(K_p\) is the ratio at equilibrium under the model; it is not a rule that \(C_t(t)=K_pC_b(t)\) at every instant. During a changing input, tissue concentration reflects the preceding blood concentrations as well as the current one. Second, the time constant depends on both capacity, represented here by \(V K_p\), and delivery by flow \(Q\). Thus a high-flow tissue is not automatically the fastest tissue under every possible set of parameters.
One alternative is an equilibrium approximation: set \(C_t(t)\) equal to \(K_pC_b(t)\). This can be useful when the blood input changes slowly relative to the tissue time constant and the tissue has had enough time to approach equilibrium. It can mislead when input changes on a comparable or shorter timescale, or when tissue starts far from equilibrium. The dynamic perfusion-limited equation is the appropriate comparison when such lag matters. A full PBPK simulation may also couple the blood input to tissue return flows and other organs; the prescribed shared input used here is for isolating and comparing tissue responses, not for asserting that the three tissues jointly determine blood concentration.
Worked Example: A Constant Blood Concentration Step
Suppose instead that blood concentration is maintained at \(2.00\ \mathrm{mg/L}\) from time zero, and that a tissue starts at zero concentration. The solution is \(C_t(t)=K_pC_b(1-e^{-t/\tau})\). Calculate muscle and adipose concentrations at two hours, and compare them with equilibrium.
For muscle, \(C_t=(2.00)(2.00\ \mathrm{mg/L})(1-e^{-2.00/1.11111})=4.00(1-e^{-1.8})=3.339\ \mathrm{mg/L}\). Its equilibrium concentration is \(K_pC_b=(2.00)(2.00)=4.00\ \mathrm{mg/L}\), so it has reached about \(3.339/4.00=0.8347\), or \(83.47\%\), of equilibrium.
For adipose tissue, \(C_t=(10.0)(2.00\ \mathrm{mg/L})(1-e^{-2.00/16.6667})=20.0(1-e^{-0.12})=2.262\ \mathrm{mg/L}\). Its equilibrium concentration is \(20.0\ \mathrm{mg/L}\), so it has reached about \(2.262/20.0=0.1131\), or \(11.31\%\), of equilibrium. The same maintained blood concentration therefore produces different progress toward equilibrium at the same elapsed time. This step-input calculation and the declining-input example answer different questions; their tissue concentrations should not be interchanged.
Checks, Failure Modes, and Limitations
Several checks help catch errors before interpreting a tissue comparison. The time constant \(VK_p/Q\) must have units of time. For a positive step input and zero initial tissue concentration, the concentration should rise toward \(K_pC_b\) without exceeding it. For an exponentially declining input, a calculated concentration should be nonnegative; when \(k\tau>1\), both the denominator and the difference of exponentials in the response formula are negative, so their ratio remains positive. A negative result in that case usually signals a reversed subtraction or an inconsistent sign.
The model’s comparison is limited by its assumptions. It represents each tissue with one well-mixed concentration and assumes rapid enough exchange for the perfusion-limited relation to hold. If observed tissue behavior cannot be represented by that relation, these calculations should not be rescued by changing \(K_p\) or flow without evidence. Also, a transient ratio \(C_t/C_b\) may change over time and should not be reported as though it were the equilibrium partition coefficient.
- Do not rank tissues by flow alone. Compare \(V K_p/Q\) to assess their time constants.
- Do not confuse equilibrium concentration with a dynamic response. Use \(K_pC_b(t)\) only as the equilibrium target at that blood concentration.
- Do not round intermediate exponentials too early. Carry the specified parameters through the calculation and round the final concentration.
- Check the input assumption. A prescribed common blood curve isolates tissue response; it is not automatically the blood curve of a coupled whole-body model.
Check Your Understanding
Answer each question using the equations and hypothetical comparisons in this tutorial. Short solutions follow the questions.
- If \(V=12.0\) L, \(K_p=3.00\), and \(Q=18.0\) L/h, what are the time constant and half-time?
- For a blood concentration held at \(1.50\ \mathrm{mg/L}\), what is the equilibrium tissue concentration when \(K_p=4.00\)?
- Why does a large \(K_p\) not, by itself, guarantee the fastest approach to equilibrium?
- For an exponential blood input, what does it indicate if \(k\tau=1\) makes the closed-form expression appear undefined?
- At a particular time, can a transient \(C_t/C_b\) be assumed equal to \(K_p\)? Explain.
Solutions. (1) \(\tau=(12.0)(3.00)/18.0=2.00\) h and \(t_{1/2}=(\ln 2)(2.00)=1.386\) h. (2) \(C_t=K_pC_b=(4.00)(1.50)=6.00\ \mathrm{mg/L}\). (3) The time constant is \(V K_p/Q\), so volume and flow also matter; a larger \(K_p\) can lengthen the time constant if the other values are fixed. (4) The displayed form has a zero denominator, but the limit is finite: \(C_t(t)=(Q/V)C_0t e^{-kt}\). (5) No. \(K_p\) is the equilibrium ratio; during a changing input the tissue can lag behind blood, so the transient ratio may differ.